Abstract:This work presents a few variational multiscale models for charge transport in complex physical, chemical and biological systems and engineering devices, such as fuel cells, solar cells, battery cells, nanofluidics, transistors and ion channels. An essential ingredient of the present models, introduced in an earlier paper (Bulletin of Mathematical Biology, 72, 1562-1622, 2010), is the use of differential geometry theory of surfaces as a natural means to geometrically separate the macroscopic domain from the microscopic domain, meanwhile, dynamically couple discrete and continuum descriptions. Our main strategy is to construct the total energy functional of a charge transport system to encompass the polar and nonpolar free energies of solvation, and chemical potential related energy. By using the Euler-Lagrange variation, coupled Laplace-Beltrami and Poisson-Nernst-Planck (LB-PNP) equations are derived. The solution of the LB-PNP equations leads to the minimization of the total free energy, and explicit profiles of electrostatic potential and densities of charge species. To further reduce the computational complexity, the Boltzmann distribution obtained from the Poisson-Boltzmann (PB) equation is utilized to represent the densities of certain charge species so as to avoid the computationally expensive solution of some Nernst-Planck (NP) equations. Consequently, the coupled Laplace-Beltrami and Poisson-Boltzmann-Nernst-Planck (LB-PBNP) equations are proposed for charge transport in heterogeneous systems. A major emphasis of the present formulation is the consistency between equilibrium LB-PB theory and non-equilibrium LB-PNP theory at equilibrium. Another major emphasis is the capability of the reduced LB-PBNP model to fully recover the prediction of the LB-PNP model at non-equilibrium settings. To account for the fluid impact on the charge transport, we derive coupled Laplace-Beltrami, Poisson-Nernst-Planck and Navier-Stokes equations from the variational principle for chemo-electro-fluid systems. A number of computational algorithms is developed to implement the proposed new variational multiscale models in an efficient manner. A set of ten protein molecules and a realistic ion channel, Gramicidin A, are employed to confirm the consistency and verify the capability. Extensive numerical experiment is designed to validate the proposed variational multiscale models. A good quantitative agreement between our model prediction and the experimental measurement of current-voltage curves is observed for the Gramicidin A channel transport. This paper also provides a brief review of the field.

Model of substance transport through the membrane of the last universal common ancestor (LUCA)

The Unfinished Reconstructed Nature of the Last Universal Common Ancestor

Evolution of Cellular Organization Along the First Branches of the Tree of Life

The Last Universal Common Ancestor of Ribosome-Encoding Organisms: Portrait of LUCA

Modeling of SGLT1 in Reconstituted Systems Reveals Apparent Ion-Dependencies of Glucose Uptake and Strengthens the Notion of Water-Permeable Apo States

Microscopic Characterization of Membrane Transporter Function by In Silico Modeling and Simulation

Ion transport through cell membrane channels

Functional Blood-Brain Barrier Model with Tight Connected Minitissue by Liquid Substrates Culture.

The First Universal Common Ancestor (FUCA) as the Earliest Ancestor of LUCA’s (Last UCA) Lineage

Quantitative description of ion transport via plasma membrane of yeast and small cells

The nature of the last universal common ancestor and its impact on the early Earth system

A unified evolutionary origin for the ubiquitous protein transporters SecY and YidC

A Model of Ion Transport Processes along and Across the Neuronal Membrane.

Looking for the Last Universal Common Ancestor (LUCA)

The metabolic network of the last bacterial common ancestor

A Simple Quantitative Model of Neuromodulation. Part I: Ion Flow Through Neural Ion Channels

The Origin(s) of Luca: Computer Simulation of a New Theory

Non-Markovian intracellular transport with sub-diffusion and run-length dependent detachment rate

A unifying model for membrane protein biogenesis

Membrane tubes with active pumping: water transport, vacuole formation and osmoregulation

Variational multiscale models for charge transport