Source author record

Martin Z. Bazant

Martin Z. Bazant appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

physics.chem-ph physics.flu-dyn cond-mat.mtrl-sci cond-mat.soft cond-mat.stat-mech cond-mat.mes-hall Biological Physics math-ph math.MP physics.app-ph Applications cond-mat cond-mat.other math.AP math.PR physics.class-ph physics.comp-ph physics.hist-ph

Catalog footprint

What is connected

63works

18topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2022arXiv

Gelation, Clustering and Crowding in the Electrical Double Layer of Ionic Liquids

Understanding the bulk and interfacial properties of super-concentrated electrolytes, such as ionic liquids (ILs), has attracted significant attention lately for their promising applications in supercapacitors and batteries. Recently, McEldrew \textit{et al.} developed a theory for reversible ion associations in bulk ILs, which accounted for the formation of all possible Cayley tree clusters and a percolating ionic network (gel). Here we adopt and develop this approach to understand the associations of ILs in the electrical double layer at electrified interfaces. With increasing charge of the electrode, the theory predicts a transition from a regime dominated by a gelled or clustered state to a regime dominated by free, crowded ions. This transition from gelation to crowding is conceptually similar to the overscreening to crowding transition.

preprint2022arXiv

Interfacial resistive switching by multiphase polarization in ion-intercalation nanofilms

Nonvolatile resistive-switching (RS) memories promise to revolutionize hardware architectures with in-memory computing. Recently, ion-interclation materials have attracted increasing attention as potential RS materials for their ion-modulated electronic conductivity. In this Letter, we propose RS by multiphase polarization (MP) of ion-intercalated thin films between ion-blocking electrodes, in which interfacial phase separation triggered by an applied voltage switches the electron-transfer resistance. We develop an electrochemical phase-field model for simulations of coupled ion-electron transport and ion-modulated electron-transfer rates and use it to analyze the MP switching current and time, resistance ratio, and current-voltage response. The model is able to reproduce the complex cyclic voltammograms of lithium titanate (LTO) memristors, which cannot be explained by existing models based on bulk dielectric breakdown. The theory predicts the achievable switching speeds for multiphase ion-intercalation materials and could be used to guide the design of high-performance MP-based RS memories.

preprint2022arXiv

Polar Liquids at Charged Interfaces: A Dipolar Shell Theory

The structure of polar liquids and electrolytic solutions, such as water and aqueous electrolytes, at interfaces underlies numerous phenomena in physics, chemistry, biology, and engineering. In this work, we develop a continuum theory that captures the essential features of dielectric screening by polar liquids at charged interfaces, including oscillations in charge and mass, starting from the molecular properties of the solvent. The theory predicts an anisotropic dielectric tensor of interfacial polar liquids previously studied in molecular dynamics simulations. We explore the effect of the interfacial polar liquid properties on the capacitance of the electrode/electrolyte interface and on hydration forces between two plane-parallel polarized surfaces. In the linear response approximation, we obtain simple formulas for the characteristic decay lengths of molecular and ionic profiles at the interface.

preprint2021arXiv

Correlative image learning of chemo-mechanics in phase-transforming solids

Constitutive laws underlie most physical processes in nature. However, learning such equations in heterogeneous solids (e.g., due to phase separation) is challenging. One such relationship is between composition and eigenstrain, which governs the chemo-mechanical expansion in solids. In this work, we developed a generalizable, physically-constrained image-learning framework to algorithmically learn the chemo-mechanical constitutive law at the nanoscale from correlative four-dimensional scanning transmission electron microscopy and X-ray spectro-ptychography images. We demonstrated this approach on Li$_X$FePO$_4$, a technologically-relevant battery positive electrode material. We uncovered the functional form of composition-eigenstrain relation in this two-phase binary solid across the entire composition range (0 $\leq$ X $\leq$ 1), including inside the thermodynamically-unstable miscibility gap. The learned relation directly validates Vegard's law of linear response at the nanoscale. Our physics-constrained data-driven approach directly visualizes the residual strain field (by removing the compositional and coherency strain), which is otherwise impossible to quantify. Heterogeneities in the residual strain arise from misfit dislocations and were independently verified by X-ray diffraction line profile analysis. Our work provides the means to simultaneously quantify chemical expansion, coherency strain and dislocations in battery electrodes, which has implications on rate capabilities and lifetime. Broadly, this work also highlights the potential of integrating correlative microscopy and image learning for extracting material properties and physics.

preprint2021arXiv

Dip-coating of bidisperse particulate suspensions

Dip-coating consists in withdrawing a substrate from a bath to coat it with a thin liquid layer. This process is well-understood for homogeneous fluids, but heterogeneities such as particles dispersed in the liquid lead to more complex situations. Indeed, particles introduce a new length scale, their size, in addition to the thickness of the coating film. Recent studies have shown that at first order, the thickness of the coating film for monodisperse particles can be captured by an effective capillary number based on the viscosity of the suspension, providing that the film is thicker than the particle diameter. However, suspensions involved in most practical applications are polydisperse, characterized by a wide range of particle sizes, introducing additional length scales. In this study, we investigate the dip coating of suspensions having a bimodal size distribution of particles. We show that the effective viscosity approach is still valid in the regime where the coating film is thicker than the diameter of the largest particles, although bidisperse suspensions are less viscous than monodisperse suspensions of the same solid fraction. We also characterize the intermediate regime that consists of a heterogeneous coating layer and where the composition of the film is different from the composition of the bath. A model to predict the probability of entraining the particles in the liquid film depending on their sizes is proposed and captures our measurements. In this regime, corresponding to a specific range of withdrawal velocities, capillarity filters the large particles out of the film.

preprint2021arXiv

Growth morphology and symmetry selection of interfacial instabilities in anisotropic environments

The displacement of a fluid by another less viscous one in a quasi-two dimensional geometry typically leads to complex fingering patterns. In an isotropic system, dense-branching growth arises, which is characterized by repeated tip-splitting of evolving fingers. When anisotropy is present in the interfacial dynamics, the growth morphology changes to dendritic growth characterized by regular structures. We introduce anisotropy by engraving a six-fold symmetric lattice of channels on a Hele-Shaw cell. We show that the morphology transition in miscible fluids depends not only on the previously reported degree of anisotropy set by the lattice topography, but also on the viscosity ratio between the two fluids. Remarkably, the viscosity ratio and the degree of anisotropy also govern the global features of the dendritic patterns, inducing a systematic change from six-fold towards twelve-fold symmetric dendrites. Varying either control parameter provides a new method to tune the symmetry of complex patterns, which may also have relevance for analogous phenomena of gradient-driven interfacial dynamics, such as directional solidification or electrodeposition.

preprint2021arXiv

Ion Clusters and Networks in "Water-in-Salt Electrolytes"

Water-in-salt electrolytes (WiSEs) are a class of super-concentrated electrolytes that have shown much promise in replacing organic electrolytes in lithium-ion batteries. At the extremely high salt concentrations of WiSEs, ionic association is more complicated than the simple ion pair description. In fact, large branched clusters can be present in WiSEs, and past a critical salt concentration, an infinite percolating ionic network can form spontaneously. In this work, we simplify our recently developed thermodynamic model of reversible ionic aggregation and gelation, tailoring it specifically for WiSEs. Our simplified theory only has a handful of parameters, all of which can be readily determined from simulations. Our model is able to quantitatively reproduce the populations of ionic clusters of different sizes as a function of salt concentration, the critical salt concentration for ionic gelation, and the fraction of ions incorporated into the ionic gel, as observed from molecular simulations of three different lithium-based WiSEs. The extent of ionic association and gelation greatly affects the effective ionic strength of solution, the coordination environment of active cations that is known to govern the chemistry of the solid-electrolyte interface, and the thermodynamic activity of all species in the electrolyte.

preprint2020arXiv

Analysis, Design, and Generalization of Electrochemical Impedance Spectroscopy (EIS) Inversion Algorithms

We introduce a framework for analyzing and designing EIS inversion algorithms. Our framework stems from the observation of four features common to well-defined EIS inversion algorithms, namely (1) the representation of unknown distributions, (2) the minimization of a metric of error to estimate parameters arising from the chosen representation, subject to constraints on (3) the complexity control parameters, and (4) a means for choosing optimal control parameter values. These features must be present to overcome the ill-posed nature of EIS inversion problems. We review three established EIS inversion algorithms to illustrate the pervasiveness of these features, and show the utility of the framework by resolving ambiguities concerning three more algorithms. Our framework is then used to design the generalized EIS inversion (gEISi) algorithm, which uses Gaussian basis function representation, modality control parameter, and cross-validation for choosing the optimal control parameter value. The gEISi algorithm is applicable to the generalized EIS inversion problem, which allows for a wider range of underlying models. We also considered the construction of credible intervals for distributions arising from the algorithm. The algorithm is able to accurately reproduce distributions which have been difficult to obtain using existing algorithms. It is provided gratis on the repository https://github.com/suryaeff/gEISi.git.

preprint2020arXiv

Continuum theory of electrostatic correlations at charged surfaces

The standard model for diffuse charge phenomena in colloid science, electrokinetics and biology is the Poisson-Boltzmann mean-field theory, which breaks down for multivalent ions and large surface charge densities due to electrostatic correlations. In this paper, we formulate a predictive continuum theory of correlated electrolytes based on two extensions of the Bazant-Storey-Kornyshev (BSK) framework: (i) a physical boundary condition enforcing continuity of the Maxwell stress at a charged interface, which upholds the Contact Theorem for dilute primitive-model electrolytes, and (ii) scaling relationships for the correlation length, for a one-component plasma at a charged plane and around a cylinder, as well as a dilute z:1 electrolyte screening a planar surface. In these cases, the theory accurately reproduces Monte Carlo simulation results from weak to strong coupling, and extensions are possible for more complex models of electrolytes and ionic liquids.

preprint2020arXiv

Dielectric breakdown by electric-field induced phase separation

The control of the dielectric and conductive properties of device-level systems is important for increasing the efficiency of energy- and information-related technologies. In some cases, such as neuromorphic computing, it is desirable to increase the conductivity of an initially insulating medium by several orders of magnitude, resulting in effective dielectric breakdown. Here, we show that by tuning the value of the applied electric field in systems { with variable permittivity and electric conductivity}, e.g. ion intercalation materials, we can vary the device-level electrical conductivity by orders of magnitude. We attribute this behavior to the formation of filament-like conductive domains that percolate throughout the system, { which form only when the electric conductivity depends on the concentration}. We conclude by discussing the applicability of our results in neuromorphic computing devices and Li-ion batteries.

preprint2020arXiv

Equivalent circuit model for electrosorption with redox active materials

Electrosorption is a promising technique for brackish water deionization and waste water remediation. Faradaic materials with redox activity have recently been shown to enhance both the adsorption capacity and the selectivity of electrosorption processes. Development of the theory of electrosorption with redox active materials can provide a fundamental understanding of the electrosorption mechanism and a means to extract material properties from small-scale experiments for process optimization and scale-up. Here, we present an intuitive, physics-based equivalent circuit model to describe the electrosorption performance of redox active materials, which is able to accurately fit experimental cyclic voltammetry measurements. The model can serve as an efficient and easy-to-implement tool to evaluate properties of redox active materials and help to distinguish between the transport-limited and reaction-limited regimes in electrosorption processes. And the extracted intrinsic material properties can be further incorporated into process models under lower supporting electrolyte concentrations for realistic electrosorption applications.

preprint2020arXiv

Interfacial layering in the electric double layer of ionic liquids

Ions in ionic liquids and concentrated electrolytes reside in a crowded, strongly-interacting environment, leading to the formation of discrete layers of charges at interfaces. Here, we propose a continuum theory that captures the transition from overscreening-- alternating layers of excess charge at low surface potential, to overcrowding-- the formation of dense layers of charge of the same sign at high surface potential. The model outputs slowly-decaying oscillations in the charge density with a wavelength of single ion diameters, as shown by analysis of the gradient expansion. The gradient expansion suggests a new structure for partial differential equations describing the electrostatic potential at charged interfaces. We find quantitative agreement between the theory and molecular simulations in the differential capacitance and concentration profiles.

preprint2020arXiv

Theory of Faradaically Modulated Redox Active Electrodes for Electrochemically Mediated Selective Adsorption Processes

Electrochemically mediated selective adsorption is an emerging electrosorption technique that utilizes Faradaically enhanced redox active electrodes, which can adsorb ions not only electrostatically, but also electrochemically. The superb selectivity (>100) of this technique enables selective removal of toxic or high-value target ions under low energy consumption. Here, we develop a general theoretical framework to describe the competitive electrosorption phenomena involving multiple ions and surface-bound redox species. The model couples diffusion, convection and electromigration with competitive surface adsorption reaction kinetics, consistently derived from non-equilibrium thermodynamics. To optimize the selective removal of the target ions, design criteria were derived analytically from physically relevant dimensionless groups and time scales, where the propagation of the target anions concentration front is the limiting step. Detailed computational studies are reported for three case studies that cover a wide range of inlet concentration ratios between the competing ions. And in all three cases, target anions in the electrosorption cell forms a self-sharpening reaction-diffusion wave front. Based on the model, a three-step stop-flow operation scheme with a pure stripping solution of target anions is proposed that optimizes the ion adsorption performance and increases the purity of the regeneration stream to almost 100%, which is beneficial for downstream processing.

preprint2020arXiv

Theory of Ion Aggregation and Gelation in Super-Concentrated Electrolytes

In concentrated electrolytes with asymmetric or irregular ions, such as ionic liquids and solvent-in-salt electrolytes, ion association is more complicated than simple ion-pairing. Large branched aggregates can form at significant concentrations at even moderate salt concentrations. When the extent of ion association reaches a certain threshold, a percolating ionic gel networks can form spontaneously. Gelation is a phenomenon that is well known in polymer physics, but it is practically unstudied in concentrated electrolytes. However, despite this fact, the ion-pairing description is often applied to these systems for the sake of simplicity. In this work, drawing strongly from established theories in polymer physics, we develop a simple thermodynamic model of reversible ionic aggregation and gelation in concentrated electrolytes accounting for the competition between ion solvation and ion association. Our model predicts the populations of ionic clusters of different sizes as a function of salt concentration, it captures the onset of ionic gelation and also the post-gel partitioning of ions into the gel. We discuss the applicability of our model, as well as the implications of its predictions on thermodynamic, transport, and rheological properties.

preprint2020arXiv

Tuning the stability of Electrochemical Interfaces by Electron Transfer reactions

The morphology of interfaces is known to play fundamental role on the efficiency of energy-related applications, such light harvesting or ion intercalation. Altering the morphology on demand, however, is a very difficult task. Here, we show ways the morphology of interfaces can be tuned by driven electron transfer reactions. By using non-equilibrium thermodynamic stability theory, we uncover the operating conditions that alter the interfacial morphology. We apply the theory to ion intercalation and surface growth where electrochemical reactions are described using Butler-Volmer or coupled ion-electron transfer kinetics. The latter connects microscopic/quantum mechanical concepts with the morphology of electrochemical interfaces. Finally, we construct non-equilibrium phase diagrams in terms of the applied driving force (current/voltage) and discuss the importance of engineering the density of states of the electron donor in applications related to energy harvesting and storage, electrocatalysis and photocatalysis.

preprint2019arXiv

Deionization Shock Driven by Electroconvection in a Circular Channel

In a circular channel passing over-limiting current (faster than diffusion), transient vortices of bulk electroconvection are observed in salt-depleted region within the horizontal plane. The spatiotemporal evolution of the salt concentration is directly visualized, revealing the propagation of a deionization shock wave driven by bulk electroconvection up to millimeter scales. This novel mechanism leads to quantitatively similar dynamics as for deionization shocks in charged porous media, which are driven instead by surface conduction and electro-osmotic flow at micron to nanometer scales. The remarkable generality of deionization shocks under over-limiting current could be used to manipulate ion transport in complex geometries for desalination and water treatment.

preprint2019arXiv

Ionic activity in concentrated electrolytes: solvent structure effect revisited

We revisit the role of the local solvent structure on the activity coefficient of electrolytes with a general non-local dielectric function approach. We treat the concentrated electrolyte as a dielectric medium and suggest an interpolated formula for the dielectric response. The pure water limit is calibrated based on MD simulations and experimental data. Solving our model around a central ion, we find strong over-screening and oscillations in the potential, which are absent in the standard "primitive model" predictions. We obtain mathematically tractable closed-form expressions for the activity coefficients and show reasonable agreement with experimental data.

preprint2019arXiv

Physics of Electrostatic Projection Revealed by High-Speed Video Imaging

Processes based on electrostatic projection are used extensively in industry, e.g. for mineral separations, electrophotography or manufacturing of coated abrasives, such as sandpaper. Despite decades of engineering practice, there are still unanswered questions. In this paper, we present a comprehensive experimental study of projection process of more than 1500 individual spherical alumina particles with a nominal size of 500 $μ$m, captured by high-speed video imaging and digital image analysis. Based on flight trajectories of approximately 1100 projected particles, we determined the acquired charge and dynamics as a function of relative humidity (RH) and electric field intensity and compared the results with classical theories. For RH levels of 50\% and above, more than 85\% of disposed particles were projected, even when the electric field intensity was at its minimum level. This suggests that, beyond a critical value of electric field intensity, relative humidity plays a more critical role in the projection process. We also observed that the charging time is reduced dramatically for RH levels of 50\% and above, possibly due to the build-up of thin water films around the particles which can facilitate charge transfer. In contrast, projected particles at 30\% RH level exhibited an excessive amount of electric charge, between two to four times than that of saturation value, which might be attributed to triboelectric charging effects. Finally, the physics of electrostatic projection is compared and contrasted with those of induced-charge electrokinetic phenomena, which share similar field-square scaling, as the applied field acts on its own induced charge to cause particle motion.

preprint2019arXiv

The Heat of Nervous Conduction: A Thermodynamic Framework

Early recordings of nervous conduction revealed a notable thermal signature associated with the electrical signal. The observed production and subsequent absorption of heat arise from physicochemical processes that occur at the cell membrane level during the conduction of the action potential. In particular, the reversible release of electrical energy stored as a difference of potential across the cell membrane appears as a simple yet consistent explanation for the heat production, as proposed in the "Condenser Theory." However, the Condenser Theory has not been analyzed beyond the analogy between the cell membrane and a parallel-plate capacitor, i.e. a condenser, which cannot account for the magnitude of the heat signature. In this work, we use a detailed electrostatic model of the cell membrane to revisit the Condenser Theory. We derive expressions for free energy and entropy changes associated with the depolarization of the membrane by the action potential, which give a direct measure of the heat produced and absorbed by neurons. We show how the density of surface charges on both sides of the membrane impacts the energy changes. Finally, considering a typical action potential, we show that if the membrane holds a bias of surface charges, such that the internal side of the membrane is 0.05 C m$^{-2}$ more negative than the external side, the size of the heat predicted by the model reaches the range of experimental values. Based on our study, we identify the change in electrical energy of the membrane as the primary mechanism of heat production and absorption by neurons during nervous conduction.

preprint2016arXiv

Dendrite Suppression by Shock Electrodeposition in Charged Porous Media

It is shown that surface conduction can stabilize electrodeposition in random, charged porous media at high rates, above the diffusion-limited current. After linear sweep voltammetry and impedance spectroscopy, copper electrodeposits are visualized by scanning electron microscopy and energy dispersive spectroscopy in two different porous separators (cellulose nitrate, polyethylene), whose surfaces are modified by layer-by-layer deposition of positive or negative charged polyelectrolytes. Above the limiting current, surface conduction inhibits growth in the positive separators and produces irregular dendrites, while it enhances growth and suppresses dendrites behind a deionization shock in the negative separators, also leading to improved cycle life. The discovery of stable uniform growth in the random media differs from the non-uniform growth observed in parallel nanopores and cannot be explained by classical quasi-steady leaky membrane models, which always predict instability and dendritic growth. Instead, the experimental results suggest that transient electro-diffusion in random porous media imparts the stability of a deionization shock to the growing metal interface behind it. Shock electrodeposition could be exploited to enhance the cycle life and recharging rate of metal batteries or to accelerate the fabrication of metal matrix composite coatings.

preprint2016arXiv

Exact solutions and physical analogies for unidirectional flows

Unidirectional flow is the simplest phenomenon of fluid mechanics. Its mathematical description, the Dirichlet problem for Poisson's equation in two dimensions with constant forcing, arises in many physical contexts, such as the torsion of elastic beams, first solved by de Saint-Venant for complex shapes. Here, the literature is unified and extended by identifying seventeen physical analogies for unidirectional flow and describing their common mathematical structure. Besides classical analogies in fluid and solid mechanics, applications are discussed in stochastic processes (first passage in two dimensions), pattern formation (river growth by erosion), and electrokinetics (ion transport in nanochannels), which also involve Poisson's equation with non-constant forcing. Methods are given to construct approximate geometries that admit exact solutions, by adding harmonic functions to quadratic forms or by truncating eigenfunction expansions. Exact solutions for given geometries are also derived by conformal mapping. We prove that the remarkable geometrical interpretation of Poiseuille flow in an equilateral triangular pipe (the product of the distances from an interior point to the sides) is only shared by parallel plates and unbounded equilateral wedges (with the third side hidden behind the apex). We also prove Onsager reciprocity for linear electrokinetic phenomena in straight pores of arbitrary shape and surface charge, based on the mathematics of unidirectional flow.

preprint2016arXiv

Membraneless flow battery leveraging flow-through heterogeneous porous media for improved power density and reduced crossover

A key factor preventing the market penetration of renewable, intermittent energy sources, such as solar, wind and wave, is the lack of cost-effective energy storage options to counteract intermittency. Here, we propose and demonstrate a novel flow battery architecture that replaces traditional ion-exchange membranes with less expensive heterogeneous flow-through porous media. We present an experimentally-validated model which demonstrates that our architecture promises reduced crossover of reactive species compared to typical membraneless systems employing co-laminar flows in open channels. In addition, our prototype battery exhibits significantly improved power density (0.925 W/cm2) and maximum current density (3 A/cm2) compared to previous membraneless systems.

preprint2016arXiv

Performance and Degradation of A Lithium-Bromine Rechargeable Fuel Cell Using Highly Concentrated Catholytes

Lithium-air batteries have been considered as ultimate solutions for the power source of long-range electrified transportation, but state-of-the-art prototypes still suffer from short cycle life, low efficiency and poor power output. Here, a lithium-bromine rechargeable fuel cell using highly concentrated bromine catholytes is demonstrated with comparable specific energy, improved power density, and higher efficiency. The cell is similar in structure to a hybrid-electrolyte Li-air battery, where a lithium metal anode in nonaqueous electrolyte is separated from aqueous bromine catholytes by a lithium-ion conducting ceramic plate. The cell with a flat graphite electrode can discharge at a peak power density around 9mW cm-2 and in principle could provide a specific energy of 791.8 Wh kg-1, superior to most existing cathode materials and catholytes. It can also run in regenerative mode to recover the lithium metal anode and free bromine with 80-90% voltage efficiency, without any catalysts. Degradation of the solid electrolyte and the evaporation of bromine during deep charging are challenges that should be addressed in improved designs to fully exploit the high specific energy of liquid bromine. The proposed system offers a potential power source for long-range electric vehicles, beyond current Li-ion batteries yet close to envisioned Li-air batteries.

preprint2016arXiv

Resistive Switching in Aqueous Nanopores by Shock Electrodeposition

Solid-state programmable metallization cells have attracted considerable attention as memristive elements for Redox-based Resistive Random Access Memory (ReRAM) for low-power and low-voltage applications. In principle, liquid-state metallization cells could offer the same advantages for aqueous systems, such as biomedical lab-on-a-chip devices, but robust resistive switching has not yet been achieved in liquid electrolytes, where electrodeposition is notoriously unstable to the formation of fractal dendrites. Here, the recently discovered physics of shock electrodeposition are harnessed to stabilize aqueous copper growth in polycarbonate nanopores, whose surfaces are modified with charged polymers. Stable bipolar resistive switching is demonstrated for 500 cycles with <10s retention times, prior to any optimization of the geometry or materials.

preprint2015arXiv

Asymmetric collapse by dissolution or melting in a uniform flow

An advection--diffusion-limited dissolution model of an object being eroded by a two-dimensional potential flow is presented. By taking advantage of the conformal invariance of the model, a numerical method is introduced that tracks the evolution of the object boundary in terms of a time-dependent Laurent series. Simulations of a variety of dissolving objects are shown, which shrink and then collapse to a single point in finite time. The simulations reveal a surprising exact relationship whereby the collapse point is the root of a non-analytic function given in terms of the flow velocity and the Laurent series coefficients describing the initial shape. This result is subsequently derived using residue calculus. The structure of the non-analytic function is examined for three different test cases, and a practical approach to determine the collapse point using a generalized Newton--Raphson root-finding algorithm is outlined. These examples also illustrate the possibility that the model breaks down in finite time prior to complete collapse, due to a topological singularity, as the dissolving boundary overlaps itself rather than breaking up into multiple domains (analogous to droplet pinch-off in fluid mechanics). In summary, the model raises fundamental mathematical questions about broken symmetries in finite-time singularities of both continuous and stochastic dynamical systems.

preprint2014arXiv

Charge Transfer Kinetics at the Solid-Solid Interface in Porous Electrodes

Interfacial charge transfer is widely assumed to obey Butler-Volmer kinetics. For certain liquid-solid interfaces, Marcus-Hush-Chidsey theory is more accurate and predictive, but it has not been applied to porous electrodes. Here we report a simple method to extract the charge transfer rates in carbon-coated LiFePO4 porous electrodes from chronoamperometry experiments, obtaining curved Tafel plots that contradict the Butler-Volmer equation but fit the Marcus-Hush-Chidsey prediction over a range of temperatures. The fitted reorganization energy matches the Born solvation energy for electron transfer from carbon to the iron redox site. The kinetics are thus limited by electron transfer at the solid-solid (carbon-LixFePO4) interface, rather than by ion transfer at the liquid-solid interface, as previously assumed. The proposed experimental method generalizes Chidsey's method for phase-transforming particles and porous electrodes, and the results show the need to incorporate Marcus kinetics in modeling batteries and other electrochemical systems.

preprint2014arXiv

Experimental Verification of Overlimiting Current by Surface Conduction and Electro-osmotic Flow in Microchannels

Possible mechanisms of overlimiting current in unsupported electrolytes, exceeding diffusion limitation, have been intensely studied for their fundamental significance and applications to desalination, separations, sensing, and energy storage. In bulk membrane systems, the primary physical mechanism is electro-convection, driven by electro-osmotic instability on the membrane surface. It has recently been predicted that confinement by charged surfaces in microchannels or porous media favors two new mechanisms, electro-osmotic flow (EOF) and surface conduction (SC), driven by large electric fields in the depleted region acting on the electric double layers on the sidewalls. Here, we provide the first direct evidence for the transition from SC to EOF above a critical channel height, using in situ particle tracking and current-voltage measurements in a micro/nanofluidic device. The dependence of the over-limiting conductance on channel depth (d) is consistent with theoretical predictions, scaling as d^-1 for SC and d^4/5 for EOF with a transition around d=8um. This complete picture of surface-driven over-limiting current can guide engineering applications of ion concentration polarization phenomena in microfluidics and porous media.

preprint2014arXiv

Homogenization of the Poisson-Nernst-Planck Equations for Ion Transport in Charged Porous Media

Effective Poisson-Nernst-Planck (PNP) equations are derived for macroscopic ion transport in charged porous media under periodic fluid flow by an asymptotic multi-scale expansion with drift. The microscopic setting is a two-component periodic composite consisting of a dilute electrolyte continuum (described by standard PNP equations) and a continuous dielectric matrix, which is impermeable to the ions and carries a given surface charge. Four new features arise in the upscaled equations: (i) the effective ionic diffusivities and mobilities become tensors, related to the microstructure; (ii) the effective permittivity is also a tensor, depending on the electrolyte/matrix permittivity ratio and the ratio of the Debye screening length to the macroscopic length of the porous medium; (iii) the microscopic fluidic convection is replaced by a diffusion-dispersion correction in the effective diffusion tensor; and (iv) the surface charge per volume appears as a continuous "background charge density", as in classical membrane models. The coefficient tensors in the upscaled PNP equations can be calculated from periodic reference cell problems. For an insulating solid matrix, all gradients are corrected by the same tensor, and the Einstein relation holds at the macroscopic scale, which is not generally the case for a polarizable matrix, unless the permittivity and electric field are suitably defined. In the limit of thin double layers, Poisson's equation is replaced by macroscopic electroneutrality (balancing ionic and surface charges). The general form of the macroscopic PNP equations may also hold for concentrated solution theories, based on the local-density and mean-field approximations. These results have broad applicability to ion transport in porous electrodes, separators, membranes, ion-exchange resins, soils, porous rocks, and biological tissues.

preprint2014arXiv

Multicomponent Gas Diffusion in Porous Electrodes

Multicomponent gas transport is investigated with unprecedented precision by AC impedance analysis of porous YSZ anode-supported solid oxide fuel cells. A fuel gas mixture of H2-H2O-N2 is fed to the anode, and impedance data are measured across the range of hydrogen partial pressure (10-100%) for open circuit conditions at three temperatures (800C, 850C and 900C) and for 300mA applied current at 800C. For the first time, analytical formulae for the diffusion resistance (Rb) of three standard models of multicomponent gas transport (Fick, Stefan-Maxwell, and Dusty Gas) are derived and tested against the impedance data. The tortuosity is the only fitting parameter since all the diffusion coefficients are known. Only the Dusty Gas model leads to a remarkable data collapse for over twenty experimental conditions, using a constant tortuosity consistent with permeability measurements and the Bruggeman relation. These results establish the accuracy of the Dusty Gas model for multicomponent gas diffusion in porous media and confirm the efficacy of electrochemical impedance analysis to precisely determine transport mechanisms.

preprint2014arXiv

Phase Separation Dynamics in Isotropic Ion-Intercalation Particles

Lithium-ion batteries exhibit complex nonlinear dynamics, resulting from diffusion and phase transformations coupled to ion intercalation reactions. Using the recently developed Cahn-Hilliard reaction (CHR) theory, we investigate a simple mathematical model of ion intercalation in a spherical solid nanoparticle, which predicts transitions from solid-solution radial diffusion to two-phase shrinking-core dynamics. This general approach extends previous Li-ion battery models, which either neglect phase separation or postulate a spherical shrinking-core phase boundary, by predicting phase separation only under appropriate circumstances. The effect of the applied current is captured by generalized Butler-Volmer kinetics, formulated in terms of diffusional chemical potentials, and the model consistently links the evolving concentration profile to the battery voltage. We examine sources of charge/discharge asymmetry, such as asymmetric charge transfer and surface "wetting" by ions within the solid, which can lead to three distinct phase regions. In order to solve the fourth-order nonlinear CHR initial-boundary-value problem, a control-volume discretization is developed in spherical coordinates. The basic physics are illustrated by simulating many representative cases, including a simple model of the popular cathode material, lithium iron phosphate (neglecting crystal anisotropy and coherency strain). Analytical approximations are also derived for the voltage plateau as a function of the applied current.

preprint2014arXiv

Phase Transformation Dynamics in Porous Battery Electrodes

Porous electrodes composed of multiphase active materials are widely used in Li-ion batteries, but their dynamics are poorly understood. Two-phase models are largely empirical, and no models exist for three or more phases. Using a modified porous electrode theory based on non-equilibrium thermodynamics, we show that experimental phase behavior can be accurately predicted from free energy models, without artificially placing phase boundaries or fitting the open circuit voltage. First, we simulate lithium intercalation in porous iron phosphate, a popular two-phase cathode, and show that the zero-current voltage gap, sloping voltage plateau and under-estimated exchange currents all result from size-dependent nucleation and mosaic instability. Next, we simulate porous graphite, the standard anode with three stable phases, and reproduce experimentally observed fronts of color-changing phase transformations. These results provide a framework for physics-based design and control for electrochemical systems with complex thermodynamics.

preprint2014arXiv

Simple Formula for Marcus-Hush-Chidsey Kinetics

The Marcus-Hush-Chidsey (MHC) model is well known in electro-analytical chemistry as a successful microscopic theory of outer-sphere electron transfer at metal electrodes, but it is unfamiliar and rarely used in electrochemical engineering. One reason may be the difficulty of evaluating the MHC reaction rate, which is defined as an improper integral of the Marcus rate over the Fermi distribution of electron energies. Here, we report a simple analytical approximation of the MHC integral that interpolates between exact asymptotic limits for large overpotentials, as well as for large or small reorganization energies, and exhibits less than 5\% relative error for all reasonable parameter values. This result enables the MHC model to be considered as a practical alternative to the ubiquitous Butler-Volmer equation for improved understanding and engineering of electrochemical systems.

preprint2014arXiv

Water Purification by Shock Electrodialysis: Deionization, Filtration, Separation, and Disinfection

The development of energy and infrastructure efficient water purification systems are among the most critical engineering challenges facing our society. Water purification is often a multi-step process involving filtration, desalination, and disinfection of a feedstream. Shock electrodialysis (shock ED) is a newly developed technique for water desalination, leveraging the formation of ion concentration polarization (ICP) zones and deionization shock waves in microscale pores near to an ion selective element. While shock ED has been demonstrated as an effective water desalination tool, we here present evidence of other simultaneous functionalities. We show that, unlike electrodialysis, shock ED can thoroughly filter micron-scale particles and aggregates of nanoparticles present in the feedwater. We also demonstrate that shock ED can enable disinfection of feedwaters, as approximately $99\%$ of viable bacteria (here \textit{E. coli}) in the inflow were killed or removed by our prototype. Shock ED also separates positive from negative particles, contrary to claims that ICP acts as a virtual barrier for all charged particles. By combining these functionalities (filtration, separation and disinfection) with deionization, shock ED has the potential to enable more compact and efficient water purification systems.

preprint2013arXiv

Nonlinear Dynamics of Ion Concentration Polarization in Porous Media: The Leaky Membrane Model

The conductivity of highly charged membranes is nearly constant, due to counter-ions screening pore surfaces. Weakly charged porous media, or "leaky membranes", also contain a significant concentration of co-ions, whose depletion at high current leads to ion concentration polarization and conductivity shock waves. To describe these nonlinear phenomena the absence of electro-osmotic flow, a simple Leaky Membrane Model is formulated, based on macroscopic electroneutrality and Nernst-Planck ionic fluxes. The model is solved in cases of unsupported binary electrolytes: steady conduction from a reservoir to a cation-selective surface, transient response to a current step, steady conduction to a flow-through porous electrode, and steady conduction between cation-selective surfaces in cross flow. The last problem is motivated by separations in leaky membranes, such as shock electrodialysis. The article begins with a tribute to Neal Amundson, whose pioneering work on shock waves in chromatography involved similar mathematics.

preprint2013arXiv

Overlimiting Current and Shock Electrodialysis in Porous Media

Most electrochemical processes, such as electrodialysis, are limited by diffusion, but in porous media, surface conduction and electro-osmotic flow also contribute to ionic fluxes. In this paper, we report experimental evidence for surface-driven over-limiting current (faster than diffusion) and deionization shocks (propagating salt removal) in a porous medium. The apparatus consists of a silica glass frit (1 mm thick with 500 nm mean pore size) in an aqueous electrolyte (CuSO$_4$ or AgNO$_3$) passing ionic current from a reservoir to a cation-selective membrane (Nafion). The current-voltage relation of the whole system is consistent with a proposed theory based on the electro-osmotic flow mechanism over a broad range of reservoir salt concentrations (0.1 mM - 1.0 M), after accounting for (Cu) electrode polarization and pH-regulated silica charge. Above the limiting current, deionized water ($\approx 10 μ$ $M$) can be continuously extracted from the frit, which implies the existence of a stable shock propagating against the flow, bordering a depleted region that extends more than 0.5mm across the outlet. The results suggest the feasibility of "shock electrodialysis" as a new approach to water desalination and other electrochemical separations.

preprint2013arXiv

Particle-Level Modeling of the Charge-Discharge Behavior of Nanoparticulate Phase-Separating Li-Ion Battery Electrodes

In nanoparticulate phase-separating electrodes, phase separation inside the particles can be hindered during their charge/discharge cycles even when a thermodynamic driving force for phase separation exists. In such cases, particles may (de)lithiate discretely in a process referred to as mosaic instability. This instability could be the key to elucidating the complex charge/discharge dynamics in nanoparticulate phase-separating electrodes. In this paper, the dynamics of the mosaic instability is studied using Smoothed Boundary Method simulations at the particle level, where the concentration and electrostatic potential fields are spatially resolved around individual particles. Two sets of configurations consisting of spherical particles with an identical radius are employed to study the instability in detail. The effect of an activity-dependent exchange current density on the mosaic instability, which leads to asymmetric charge/discharge, is also studied. While we show that our model reproduces the results of a porous-electrode model for the simple setup studied here, it is a powerful framework with the capability to predict the detailed dynamics in three-dimensional complex electrodes and provides further insights into the complex dynamics that result from the coupling of electrochemistry, thermodynamics, and transport kinetics.

preprint2013arXiv

Theory of Chemical Kinetics and Charge Transfer based on Nonequilibrium Thermodynamics

Classical theories of chemical kinetics assume independent reactions in dilute solutions, whose rates are determined by mean concentrations. In condensed matter, strong interactions alter chemical activities and create inhomogeneities that can dramatically affect the reaction rate. The extreme case is that of a reaction coupled to a phase transformation, whose kinetics must depend on the order parameter -- and its gradients, at phase boundaries. This Account presents a general theory of chemical kinetics based on nonequilibrium thermodynamics. The reaction rate is a nonlinear function of the thermodynamic driving force (free energy of reaction) expressed in terms of variational chemical potentials. The Cahn-Hilliard and Allen-Cahn equations are unified and extended via a master equation for non-equilibrium chemical thermodynamics. For electrochemistry, both Marcus and Butler-Volmer kinetics are generalized for concentrated solutions and ionic solids. The theory is applied to intercalation dynamics in the phase separating Li-ion battery material Li$_x$FePO$_4$.

preprint2013arXiv

Theory of Nucleation in Phase-separating Nanoparticles

The basic physics of nucleation in solid \hl{single-crystal} nanoparticles is revealed by a phase-field theory that includes surface energy, chemical reactions and coherency strain. In contrast to binary fluids, which form arbitrary contact angles at surfaces, complete "wetting" by one phase is favored at binary solid surfaces. Nucleation occurs when surface wetting becomes unstable, as the chemical energy gain (scaling with area) overcomes the elastic energy penalty (scaling with volume). The nucleation barrier thus decreases with the area-to-volume ratio and vanishes below a critical size, and nanoparticles tend to transform in order of increasing size, leaving the smallest particles homogeneous (in the phase of lowest surface energy). The model is used to simulate phase separation in realistic nanoparticle geometries for \ce{Li_XFePO4}, a popular cathode material for Li-ion batteries, and collapses disparate experimental data for the nucleation barrier, with no adjustable parameters. Beyond energy storage, the theory generally shows how to tailor the elastic and surface properties of a solid nanostructure to achieve desired phase behavior.

preprint2012arXiv

Effects of electrostatic correlations on electrokinetic phenomena

Classical theory of the electric double layer is based on the fundamental assumption of a dilute solution of point ions. There are a number of situations such as high applied voltages, high concentration of electrolytes, systems with multivalent ions, or solvent-free ionic liquids where the classical theory is often applied but the fundamental assumptions cannot be justified. Perhaps the most basic assumption underlying continuum models in electrokinetics is the mean-field approximation, that the electric field acting on each discrete ion is self-consistently determined by the local mean charge density. This paper considers situations where the mean-field approximation breaks down and electrostatic correlations become important. A fourth-order modified Poisson equation is developed that accounts for electrostatic correlations and captures the essential features in a simple continuum framework. The theory is derived variationally as a gradient approximation for non-local electrostatics, in which the dielectric permittivity becomes a differential operator. The only new parameter is a characteristic length scale for correlated ion pairs. The model is able to capture subtle aspects of more detailed simulations based on Monte Carlo, molecular dynamics, or density functional theory and allows for the straightforward calculation of electrokinetic flows in correlated liquids, for the first time. Departures from classical Helmholtz-Smoluchowski theory are controlled by the dimensionless ratio of the correlation length to the Debye screening length. Charge-density oscillations tend to reduce electro-osmotic flow and streaming current, and over-screening of the surface charge can lead to flow reversal. These effects also help to explain the apparent charge-induced thickening of double layers in induced-charge electrokinetic phenomena.

preprint2012arXiv

Nonequilibrium Thermodynamics of Porous Electrodes

We reformulate and extend porous electrode theory for non-ideal active materials, including those capable of phase transformations. Using principles of non-equilibrium thermodynamics, we relate the cell voltage, ionic fluxes, and Faradaic charge-transfer kinetics to the variational electrochemical potentials of ions and electrons. The Butler-Volmer exchange current is consistently expressed in terms of the activities of the reduced, oxidized and transition states, and the activation overpotential is defined relative to the local Nernst potential. We also apply mathematical bounds on effective diffusivity to estimate porosity and tortuosity corrections. The theory is illustrated for a Li-ion battery with active solid particles described by a Cahn-Hilliard phase-field model. Depending on the applied current and porous electrode properties, the dynamics can be limited by electrolyte transport, solid diffusion and phase separation, or intercalation kinetics. In phase-separating porous electrodes, the model predicts narrow reaction fronts, mosaic instabilities and voltage fluctuations at low current, consistent with recent experiments, which could not be described by existing porous electrode models.

preprint2012arXiv

Phase Separation with Anisotropic Coherency Strain

We consider the effects of anisotropic coherency strain (due to lattice mismatch) on phase separation in intercalation materials, motivated by the high-rate Li-ion battery material \fepo. Using a phase-field model coupled to elastic stresses, we analyze spinodal decomposition (linear instability of the homogeneous state) as well as nonlinear evolution of the phase pattern at constant mean filling. We consider fully anisotropic coherency strain and focus on the novel case of simultaneous expansion and contraction along different crystal axis, as in the case of \fepo, which leads to tilted, striped phase boundaries in equilibrium.

preprint2012arXiv

Theory of SEI Formation in Rechargeable Batteries: Capacity Fade, Accelerated Aging and Lifetime Prediction

Cycle life is critically important in applications of rechargeable batteries, but lifetime prediction is mostly based on empirical trends, rather than mathematical models. In practical lithium-ion batteries, capacity fade occurs over thousands of cycles, limited by slow electrochemical processes, such as the formation of a solid-electrolyte interphase (SEI) in the negative electrode, which compete with reversible lithium intercalation. Focusing on SEI growth as the canonical degradation mechanism, we show that a simple single-particle model can accurately explain experimentally observed capacity fade in commercial cells with graphite anodes, and predict future fade based on limited accelerated aging data for short times and elevated temperatures. The theory is extended to porous electrodes, predicting that SEI growth is essentially homogeneous throughout the electrode, even at high rates. The lifetime distribution for a sample of batteries is found to be consistent with Gaussian statistics, as predicted by the single-particle model. We also extend the theory to rapidly degrading anodes, such as nanostructured silicon, which exhibit large expansion on ion intercalation. In such cases, large area changes during cycling promote SEI loss and faster SEI growth. Our simple models are able to accurately fit a variety of published experimental data for graphite and silicon anodes.

preprint2011arXiv

Coherency strain and the kinetics of phase separation in LiFePO4

A theoretical investigation of the effects of elastic coherency on the thermodynamics, kinetics, and morphology of intercalation in single LiFePO4 nanoparticles yields new insights into this important battery material. Anisotropic elastic stiffness and misfit strains lead to the unexpected prediction that low-energy phase boundaries occur along {101} planes, while conflicting reports of phase boundary orientations are resolved by a partial loss of coherency in the {100} direction. Elastic relaxation near surfaces leads to the formation of a striped morphology, whose characteristic length scale is predicted by the model and yields an estimate of the interfacial energy. The effects of coherency strain on solubility and galvanostatic discharge are studied with a reaction-limited phase-field model, which quantitatively captures the influence of misfit strain, particle size, and temperature on solubility seen in experiments. Coherency strain strongly suppresses phase separation during discharge, which enhances rate capability and extends cycle life. The effects of elevated temperature and the feasibility of nucleation are considered in the context of multi-particle cathodes.

preprint2011arXiv

Corrections to the Central Limit Theorem for Heavy-Tailed Probability Densities

Classical Edgeworth expansions provide asymptotic correction terms to the Central Limit Theorem (CLT) up to an order that depends on the number of moments available. In this paper, we provide subsequent correction terms beyond those given by a standard Edgeworth expansion in the general case of regularly varying distributions with diverging moments (beyond the second). The subsequent terms can be expressed in a simple closed form in terms of certain special functions (Dawson's integral and parabolic cylinder functions), and there are qualitative differences depending on whether the number of moments available is even, odd or not an integer, and whether the distributions are symmetric or not. If the increments have an even number of moments, then additional logarithmic corrections must also be incorporated in the expansion parameter. An interesting feature of our correction terms for the CLT is that they become dominant outside the central region and blend naturally with known large-deviation asymptotics when these are applied formally to the spatial scales of the CLT.

preprint2011arXiv

Desalination shocks in microstructures

Salt transport in bulk electrolytes is limited by diffusion and convection, but in microstructures with charged surfaces (e.g. microfluidic devices, porous media, soils, or biological tissues) surface conduction and electro-osmotic flow also contribute to ionic fluxes. For small applied voltages, these effects lead to well known linear electrokinetic phenomena. In this paper, we predict some surprising nonlinear dynamics that can result from the competition between bulk and interfacial transport at higher voltages. When counter-ions are selectively removed by a membrane or electrode, a "desalination shock" can propagate through the microstructure, leaving in its wake an ultrapure solution, nearly devoid of co-ions and colloidal impurities. We elucidate the basic physics of desalination shocks and develop a mathematical theory of their existence, structure, and stability, allowing for slow variations in surface charge or channel geometry. Via asymptotic approximations and similarity solutions, we show that desalination shocks accelerate and sharpen in narrowing channels, while they decelerate and weaken, and sometimes disappear, in widening channels. These phenomena may find applications in separations (desalination, decontamination, biological assays) and energy storage (batteries, supercapacitors) involving electrolytes in microstructures.

preprint2011arXiv

Double layer in ionic liquids: Overscreening vs. crowding

We develop a simple Landau-Ginzburg-type continuum theory of solvent-free ionic liquids and use it to predict the structure of the electrical double layer. The model captures overscreening from short-range correlations, dominant at small voltages, and steric constraints of finite ion sizes, which prevail at large voltages. Increasing the voltage gradually suppresses overscreening in favor of the crowding of counterions in a condensed inner layer near the electrode. The predicted ion profiles and capacitance-voltage relations are consistent with recent computer simulations and experiments on room-temperature ionic liquids, using a correlation length of order the ion size.

preprint2011arXiv

Overlimiting current in a microchannel

We revisit the classical problem of diffusion-limited ion transport to a membrane (or electrode) by considering the effects of charged side walls. Using simple mathematical models and numerical simulations, we identify three basic mechanisms for over-limiting current in a microchannel: (i) {\it surface conduction} carried by excess counterions, which dominates for very thin channels, (ii) convection by {\it electro-osmotic flow} on the side walls, which dominates for thicker channels and transitions to (iii) {\it electro-osmotic instability} on the membrane end in very thick channels. These intriguing electrokinetic phenomena may find applications in biological separations, water desalination, and electrochemical energy storage.

preprint2011arXiv

Suppression of Phase Separation in LiFePO4 Nanoparticles During Battery Discharge

Using a novel electrochemical phase-field model, we question the common belief that LixFePO4 nanoparticles separate into Li-rich and Li-poor phases during battery discharge. For small currents, spinodal decomposition or nucleation leads to moving phase boundaries. Above a critical current density (in the Tafel regime), the spinodal disappears, and particles fill homogeneously, which may explain the superior rate capability and long cycle life of nano-LiFePO4 cathodes.

preprint2011arXiv

Theory of Sorption Hysteresis in Nanoporous Solids: I. Snap-Through Instabilities

The sorption-desorption hysteresis observed in many nanoporous solids, at vapor pressures low enough for the the liquid (capillary) phase of the adsorbate to be absent, has long been vaguely attributed to changes in the nanopore structure, but no mathematically consistent explanation has been presented. The present work takes an analytical approach to account for discrete molecular forces in the nanopore fluid and proposes two related mechanisms that can explain the hysteresis at low vapor pressure without assuming any change in the nanopore structure. The first mechanism, presented in Part I, consists of a series of snap-through instabilities during the filling or emptying of non-uniform nanopores or nanoscale asperities. The instabilities are caused by non-uniqueness in the misfit disjoining pressures engendered by a difference between the nanopore width and an integer multiple of the thickness of a monomolecular adsorption layer. The second mechanism, presented in Part II, consists of molecular coalescence within a partially filled surface, nanopore or nanopore network. This general thermodynamic instability is driven by attractive intermolecular forces within the adsorbate and forms the basis to develop a unified theory of both mechanisms. The ultimate goals of the theory are to predict the fluid transport in nanoporous solids from microscopic first principles, and to determine the pore size distribution and internal surface area from sorption tests.

preprint2011arXiv

Theory of Sorption Hysteresis in Nanoporous Solids: II. Molecular condensation

Motivated by the puzzle of sorption hysteresis in Portland cement concrete or cement paste, we develop in Part II of this study a general theory of vapor sorption and desorption from nanoporous solids, which attributes hysteresis to hindered molecular condensation with attractive lateral interactions. The classical mean-field theory of van der Waals is applied to predict the dependence of hysteresis on temperature and pore size, using the regular solution model and gradient energy of Cahn and Hilliard. A simple "hierarchical wetting" model for thin nanopores is developed to describe the case of strong wetting by the first monolayer, followed by condensation of nanodroplets and nanobubbles in the bulk. The model predicts a larger hysteresis critical temperature and enhanced hysteresis for molecular condensation across nanopores at high vapor pressure than within monolayers at low vapor pressure. For heterogeneous pores, the theory predicts sorption/desorption sequences similar to those seen in molecular dynamics simulations, where the interfacial energy (or gradient penalty) at nanopore junctions acts as a free energy barrier for snap-through instabilities. The model helps to quantitatively understand recent experimental data for concrete or cement paste wetting and drying cycles and suggests new experiments at different temperatures and humidity sweep rates.

preprint2010arXiv

Diffusion-Limited Aggregation on Curved Surfaces

We develop a general theory of transport-limited aggregation phenomena occurring on curved surfaces, based on stochastic iterated conformal maps and conformal projections to the complex plane. To illustrate the theory, we use stereographic projections to simulate diffusion-limited-aggregation (DLA) on surfaces of constant Gaussian curvature, including the sphere ($K>0$) and pseudo-sphere ($K<0$), which approximate "bumps" and "saddles" in smooth surfaces, respectively. Although curvature affects the global morphology of the aggregates, the fractal dimension (in the curved metric) is remarkably insensitive to curvature, as long as the particle size is much smaller than the radius of curvature. We conjecture that all aggregates grown by conformally invariant transport on curved surfaces have the same fractal dimension as DLA in the plane. Our simulations suggest, however, that the multifractal dimensions increase from hyperbolic ($K<0$) to elliptic ($K>0$) geometry, which we attribute to curvature-dependent screening of tip branching.

preprint2010arXiv

Effective slip boundary conditions for arbitrary periodic surfaces: The surface mobility tensor

In a variety of applications, most notably microfluidic design, slip-based boundary conditions have been sought to characterize fluid flow over patterned surfaces. We focus on laminar shear flows over surfaces with periodic height fluctuations and/or fluctuating Navier scalar slip properties. We derive a general formula for the "effective slip", which describes equivalent fluid motion at the mean surface as depicted by the linear velocity profile that arises far from it. We show that the slip and the applied stress are related linearly through a tensorial mobility matrix, and the method of domain perturbation is then used to derive an approximate formula for the mobility law directly in terms of surface properties. The specific accuracy of the approximation is detailed, and the mobility relation is then utilized to address several questions, such as the determination of optimal surface shapes and the effect of random surface fluctuations on fluid slip.

preprint2009arXiv

Anisotropic electro-osmotic flow over super-hydrophobic surfaces

Patterned surfaces with large effective slip lengths, such as super-hydrophobic surfaces containing trapped gas bubbles, have the potential to greatly enhance electrokinetic phenomena. Existing theories assume either homogeneous flat surfaces or patterned surfaces with thin double layers (compared to the texture correlation length) and thus predict simple surface-averaged, isotropic flows (independent of orientation). By analyzing electro-osmotic flows over striped slip-stick surfaces with arbitrary double-layer thickness, we show that surface anisotropy generally leads to a tensorial electro-osmotic mobility and subtle, nonlinear averaging of surface properties. Interestingly, the electro-osmotic mobility tensor is not simply related to the hydrodynamic slip tensor, except in special cases. Our results imply that significantly enhanced electro-osmotic flows over super-hydrophobic surfaces are possible, but only with charged liquid-gas interfaces.

preprint2009arXiv

Effective zero-thickness model for a conductive membrane driven by an electric field

The behavior of a conductive membrane in a static (DC) electric field is investigated theoretically. An effective zero-thickness model is constructed based on a Robin-type boundary condition for the electric potential at the membrane, originally developed for electrochemical systems. Within such a framework, corrections to the elastic moduli of the membrane are obtained, which arise from charge accumulation in the Debye layers due to capacitive effects and electric currents through the membrane and can lead to an undulation instability of the membrane. The fluid flow surrounding the membrane is also calculated, which clarifies issues regarding these flows sharing many similarities with flows produced by induced charge electro-osmosis (ICEO). Non-equilibrium steady states of the membrane and of the fluid can be effectively described by this method. It is both simpler, due to the zero thickness approximation which is widely used in the literature on fluid membranes, and more general than previous approaches. The predictions of this model are compared to recent experiments on supported membranes in an electric field.

preprint2009arXiv

Size-dependent spinodal and miscibility gaps for intercalation in nano-particles

Using a recently-proposed mathematical model for intercalation dynamics in phase-separating materials [Singh, Ceder, Bazant, Electrochimica Acta 53, 7599 (2008)], we show that the spinodal and miscibility gaps generally shrink as the host particle size decreases to the nano-scale. Our work is motivated by recent experiments on the high-rate Li-ion battery material LiFePO4; this serves as the basis for our examples, but our analysis and conclusions apply to any intercalation material. We describe two general mechanisms for the suppression of phase separation in nano-particles: (i) a classical bulk effect, predicted by the Cahn-Hilliard equation, in which the diffuse phase boundary becomes confined by the particle geometry; and (ii) a novel surface effect, predicted by chemical-potential-dependent reaction kinetics, in which insertion/extraction reactions stabilize composition gradients near surfaces in equilibrium with the local environment. Composition-dependent surface energy and (especially) elastic strain can contribute to these effects but are not required to predict decreased spinodal and miscibility gaps at the nano-scale.

preprint2009arXiv

Strongly nonlinear dynamics of electrolytes in large ac voltages

We study the response of a model micro-electrochemical cell to a large ac voltage of frequency comparable to the inverse cell relaxation time. To bring out the basic physics, we consider the simplest possible model of a symmetric binary electrolyte confined between parallel-plate blocking electrodes, ignoring any transverse instability or fluid flow. We analyze the resulting one-dimensional problem by matched asymptotic expansions in the limit of thin double layers and extend previous work into the strongly nonlinear regime, which is characterized by two novel features - significant salt depletion in the electrolyte near the electrodes and, at very large voltage, the breakdown of the quasi-equilibrium structure of the double layers. The former leads to the prediction of "ac capacitive desalination", since there is a time-averaged transfer of salt from the bulk to the double layers, via oscillating diffusion layers. The latter is associated with transient diffusion limitation, which drives the formation and collapse of space-charge layers, even in the absence of any net Faradaic current through the cell. We also predict that steric effects of finite ion sizes (going beyond dilute solution theory) act to suppress the strongly nonlinear regime in the limit of concentrated electrolytes, ionic liquids and molten salts. Beyond the model problem, our reduced equations for thin double layers, based on uniformly valid matched asymptotic expansions, provide a useful mathematical framework to describe additional nonlinear responses to large ac voltages, such as Faradaic reactions, electro-osmotic instabilities, and induced-charge electrokinetic phenomena.

preprint2009arXiv

Topology and shape optimization of induced-charge electro-osmotic micropumps

For a dielectric solid surrounded by an electrolyte and positioned inside an externally biased parallel-plate capacitor, we study numerically how the resulting induced-charge electro-osmotic (ICEO) flow depends on the topology and shape of the dielectric solid. In particular, we extend existing conventional electrokinetic models with an artificial design field to describe the transition from the liquid electrolyte to the solid dielectric. Using this design field, we have succeeded in applying the method of topology optimization to find system geometries with non-trivial topologies that maximize the net induced electro-osmotic flow rate through the electrolytic capacitor in the direction parallel to the capacitor plates. Once found, the performance of the topology optimized geometries has been validated by transferring them to conventional electrokinetic models not relying on the artificial design field. Our results show the importance of the topology and shape of the dielectric solid in ICEO systems and point to new designs of ICEO micropumps with significantly improved performance.

preprint2008arXiv

Tensorial hydrodynamic slip

We describe a tensorial generalization of the Navier slip boundary condition and illustrate its use in solving for flows around anisotropic textured surfaces. Tensorial slip can be derived from molecular or microstructural theories or simply postulated as an constitutive relation, subject to certain general constraints on the interfacial mobility. The power of the tensor formalism is to capture complicated effects of surface anisotropy, while preserving a simple fluid domain. This is demonstrated by exact solutions for laminar shear flow and pressure-driven flow between parallel plates of arbitrary and different textures. From such solutions, the effects of rotating a texture follow from simple matrix algebra. Our results may be useful to extracting local slip tensors from global measurements, such as the permeability of a textured channel or the force required to move a patterned surface, in experiments or simulations.

preprint2006arXiv

A stochastic flow rule for granular materials

There have been many attempts to derive continuum models for dense granular flow, but a general theory is still lacking. Here, we start with Mohr-Coulomb plasticity for quasi-2D granular materials to calculate (average) stresses and slip planes, but we propose a "stochastic flow rule" (SFR) to replace the principle of coaxiality in classical plasticity. The SFR takes into account two crucial features of granular materials - discreteness and randomness - via diffusing "spots" of local fluidization, which act as carriers of plasticity. We postulate that spots perform random walks biased along slip-lines with a drift direction determined by the stress imbalance upon a local switch from static to dynamic friction. In the continuum limit (based on a Fokker-Planck equation for the spot concentration), this simple model is able to predict a variety of granular flow profiles in flat-bottom silos, annular Couette cells, flowing heaps, and plate-dragging experiments -- with essentially no fitting parameters -- although it is only expected to function where material is at incipient failure and slip-lines are inadmissible. For special cases of admissible slip-lines, such as plate dragging under a heavy load or flow down an inclined plane, we postulate a transition to rate-dependent Bagnold rheology, where flow occurs by sliding shear planes. With different yield criteria, the SFR provides a general framework for multiscale modeling of plasticity in amorphous materials, cycling between continuum limit-state stress calculations, meso-scale spot random walks, and microscopic particle relaxation.

preprint2004arXiv

Conformal mapping methods for interfacial dynamics

The article provides a pedagogical review aimed at graduate students in materials science, physics, and applied mathematics, focusing on recent developments in the subject. Following a brief summary of concepts from complex analysis, the article begins with an overview of continuous conformal-map dynamics. This includes problems of interfacial motion driven by harmonic fields (such as viscous fingering and void electromigration), bi-harmonic fields (such as viscous sintering and elastic pore evolution), and non-harmonic, conformally invariant fields (such as growth by advection-diffusion and electro-deposition). The second part of the article is devoted to iterated conformal maps for analogous problems in stochastic interfacial dynamics (such as diffusion-limited aggregation, dielectric breakdown, brittle fracture, and advection-diffusion-limited aggregation). The third part notes that all of these models can be extended to curved surfaces by an auxilliary conformal mapping from the complex plane, such as stereographic projection to a sphere. The article concludes with an outlook for further research.

preprint2003arXiv

Conformal mapping of some non-harmonic functions in transport theory

Conformal mapping has been applied mostly to harmonic functions, i.e. solutions of Laplace's equation. In this paper, it is noted that some other equations are also conformally invariant and thus equally well suited for conformal mapping in two dimensions. In physics, these include steady states of various nonlinear diffusion equations, the advection-diffusion equations for potential flows, and the Nernst-Planck equations for bulk electrochemical transport. Exact solutions for complicated geometries are obtained by conformal mapping to simple geometries in the usual way. Novel examples include nonlinear advection-diffusion layers around absorbing objects and concentration polarizations in electrochemical cells. Although some of these results could be obtained by other methods, such as Boussinesq's streamline coordinates, the present approach is based on a simple unifying principle of more general applicability. It reveals a basic geometrical equivalence of similarity solutions for a broad class of transport processes and paves the way for new applications of conformal mapping, e.g. to non-Laplacian fractal growth.

preprint1997arXiv

Environment Dependent Interatomic Potential for Bulk Silicon

We use recent theoretical advances to develop a new functional form for interatomic forces in bulk silicon. The theoretical results underlying the model include a novel analysis of elastic properties for the diamond and graphitic structures and inversions of ab initio cohesive energy curves. The interaction model includes two-body and three-body terms which depend on the local atomic environment through an effective coordination number. This formulation is able to capture successfully: (i) the energetics and elastic properties of the ground state diamond lattice; (ii) the covalent re-hybridization of undercoordinated atoms; (iii) and a smooth transition to metallic bonding for overcoordinated atoms. Because the essential features of chemical bonding in the bulk are built into the functional form, this model promises to be useful for describing interatomic forces in silicon bulk phases and defects. Although this functional form is remarkably realistic by usual standards, it contains a small number of fitting parameters and requires computational effort comparable to the most efficient existing models. In a companion paper, a complete parameterization of the model is given, and excellent performance for condensed phases and bulk defects is demonstrated.

preprint1997arXiv

Interatomic Potential for Silicon Defects and Disordered Phases

We develop an empirical potential for silicon which represents a considerable improvement over existing models in describing local bonding for bulk defects and disordered phases. The model consists of two- and three-body interactions with theoretically motivated functional forms that capture chemical and physical trends as explained in a companion paper. The numerical parameters in the functional form are obtained by fitting to a set of ab initio results from quantum mechanical calculations based on density functional theory in the local density approximation, which include various bulk phases and defect structures. We test the potential by applying it to the relaxation of point defects, core properties of partial dislocations and the structure of disordered phases, none of which are included in the fitting procedure. For dislocations, our model makes predictions in excellent agreement with ab initio and tight-binding calculations. It is the only potential known to describe both the 30$^\circ$- and 90$^\circ$-partial dislocations in the glide set {111}. The structural and thermodynamic properties of the liquid and amorphous phases are also in good agreement with experimental and ab initio results. Our potential is the first capable of simulating a quench directly from the liquid to the amorphous phase, and the resulting amorphous structure is more realistic than with existing empirical preparation methods. These advances in transferability come with no extra computational cost, since force evaluation with our model is faster than with the popular potential of Stillinger-Weber, thus allowing reliable atomistic simulations of very large atomic systems.

Martin Z. Bazant

What is connected

Connect this record

See the researcher in context

Building this map preview

63 published item(s)

Gelation, Clustering and Crowding in the Electrical Double Layer of Ionic Liquids

Interfacial resistive switching by multiphase polarization in ion-intercalation nanofilms

Polar Liquids at Charged Interfaces: A Dipolar Shell Theory

Correlative image learning of chemo-mechanics in phase-transforming solids

Dip-coating of bidisperse particulate suspensions

Growth morphology and symmetry selection of interfacial instabilities in anisotropic environments

Ion Clusters and Networks in "Water-in-Salt Electrolytes"

Analysis, Design, and Generalization of Electrochemical Impedance Spectroscopy (EIS) Inversion Algorithms

Continuum theory of electrostatic correlations at charged surfaces

Dielectric breakdown by electric-field induced phase separation

Equivalent circuit model for electrosorption with redox active materials

Interfacial layering in the electric double layer of ionic liquids

Theory of Faradaically Modulated Redox Active Electrodes for Electrochemically Mediated Selective Adsorption Processes

Theory of Ion Aggregation and Gelation in Super-Concentrated Electrolytes

Tuning the stability of Electrochemical Interfaces by Electron Transfer reactions

Deionization Shock Driven by Electroconvection in a Circular Channel

Ionic activity in concentrated electrolytes: solvent structure effect revisited

Physics of Electrostatic Projection Revealed by High-Speed Video Imaging

The Heat of Nervous Conduction: A Thermodynamic Framework

Dendrite Suppression by Shock Electrodeposition in Charged Porous Media

Exact solutions and physical analogies for unidirectional flows

Membraneless flow battery leveraging flow-through heterogeneous porous media for improved power density and reduced crossover

Performance and Degradation of A Lithium-Bromine Rechargeable Fuel Cell Using Highly Concentrated Catholytes

Resistive Switching in Aqueous Nanopores by Shock Electrodeposition

Asymmetric collapse by dissolution or melting in a uniform flow

Charge Transfer Kinetics at the Solid-Solid Interface in Porous Electrodes

Experimental Verification of Overlimiting Current by Surface Conduction and Electro-osmotic Flow in Microchannels

Homogenization of the Poisson-Nernst-Planck Equations for Ion Transport in Charged Porous Media

Multicomponent Gas Diffusion in Porous Electrodes

Phase Separation Dynamics in Isotropic Ion-Intercalation Particles

Phase Transformation Dynamics in Porous Battery Electrodes

Simple Formula for Marcus-Hush-Chidsey Kinetics

Water Purification by Shock Electrodialysis: Deionization, Filtration, Separation, and Disinfection

Nonlinear Dynamics of Ion Concentration Polarization in Porous Media: The Leaky Membrane Model

Overlimiting Current and Shock Electrodialysis in Porous Media

Particle-Level Modeling of the Charge-Discharge Behavior of Nanoparticulate Phase-Separating Li-Ion Battery Electrodes

Theory of Chemical Kinetics and Charge Transfer based on Nonequilibrium Thermodynamics

Theory of Nucleation in Phase-separating Nanoparticles

Effects of electrostatic correlations on electrokinetic phenomena

Nonequilibrium Thermodynamics of Porous Electrodes

Phase Separation with Anisotropic Coherency Strain

Theory of SEI Formation in Rechargeable Batteries: Capacity Fade, Accelerated Aging and Lifetime Prediction

Coherency strain and the kinetics of phase separation in LiFePO4

Corrections to the Central Limit Theorem for Heavy-Tailed Probability Densities

Desalination shocks in microstructures

Double layer in ionic liquids: Overscreening vs. crowding

Overlimiting current in a microchannel

Suppression of Phase Separation in LiFePO4 Nanoparticles During Battery Discharge

Theory of Sorption Hysteresis in Nanoporous Solids: I. Snap-Through Instabilities

Theory of Sorption Hysteresis in Nanoporous Solids: II. Molecular condensation

Diffusion-Limited Aggregation on Curved Surfaces

Effective slip boundary conditions for arbitrary periodic surfaces: The surface mobility tensor

Anisotropic electro-osmotic flow over super-hydrophobic surfaces

Effective zero-thickness model for a conductive membrane driven by an electric field

Size-dependent spinodal and miscibility gaps for intercalation in nano-particles

Strongly nonlinear dynamics of electrolytes in large ac voltages

Topology and shape optimization of induced-charge electro-osmotic micropumps

Tensorial hydrodynamic slip

A stochastic flow rule for granular materials

Conformal mapping methods for interfacial dynamics

Conformal mapping of some non-harmonic functions in transport theory

Environment Dependent Interatomic Potential for Bulk Silicon

Interatomic Potential for Silicon Defects and Disordered Phases