Search NASA⌕ Search

SEARCH · Search NASA

Results for “Poisson-Boltzmann equation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

7 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

A model for discrete fracture-clay rock interaction incorporating electrostatic effects on transport

Abstract A model based on the code CrunchClay is presented for a fracture-clay matrix system that takes electrostatic effects on transport into account. The electrostatic effects on transport include those associated with the development of a diffusion potential as captured by the Nernst-Planck equation, and the formation of a diffuse layer bordering negatively charged clay particles within which partial anion exclusion occurs. The model is based on a dual continuum formulation that accounts for diffuse layer and bulk water pore space, providing a more flexible framework than is found in the classical mean electrostatic potential models. The diffuse layer model is obtained by volume averaging ion concentrations in the Poisson-Boltzmann equation, but also includes the treatment of longitudinal transport within this continuum. The calculation of transport within the bulk and diffuse layer porosity is based on a new formulation for the Nernst-Planck equation that considers averaging of diffusion coefficients and accumulation factors at grid cell interfaces. Equations for function residuals and the associated Jacobian matrix are presented such that the system of nonlinear differential-algebraic equations can be solved with Newton’s method. As an example, we consider a 2D system with a single discrete fracture within which flow and advective transport occurs that is coupled to diffusion in the clay-rich matrix. The simulation results demonstrate the lack of retardation for anions (e.g., 36 Cl − ) of the contaminant plume within the fracture flow system because they are largely excluded from the charged clay rock, while the migration of cations (e.g., 90 Sr ++ ) is more strongly attenuated. The diffusive loss of divalent cations in particular from the fracture is accentuated by their accumulation in the diffuse layer within the clay-rich matrix.

58 GEOSCIENCES↗
NASA NTRS1998

Locations of Halide Ions in Tetragonal Lysozyme Crystals

Anions play an important role in the crystallization of lysozyme, and are known to bind to the crystalline protein. Previous studies employing X-ray crystallography had found one chloride ion binding site in the tetragonal crystal form of the protein and four nitrate ion binding sites in the monoclinic form. Studies using other approaches have reported more chloride ion binding sites, but their locations were not known. Knowing the precise location of these anions is also useful in determining the correct electrostatic fields surrounding the protein. In the first part of this study the anion positions in the tetragonal form were determined from the difference Fourier map obtained from the lysozyme crystals grown in bromide and chloride solutions under identical conditions. The anion locations were then obtained from standard crystallographic methods and five possible anion binding sites were found in this manner. The sole chloride ion binding site found in previous studies was confirmed. The remaining four sites were new ones for tetragonal lysozyme crystals. However, three of these new sites and the previously found one corresponded to the four unique binding sites found for nitrate ions in monoclinic crystals. This suggests that most of the anion binding sites in lysozyme remain unchanged, even when different anions and different crystal forms of lysozyme are employed. It is unlikely that there are many more anions in the tetragonal lysozyme crystal structure. Assuming osmotic equilibrium it can be shown that there are at most three more anions in the crystal channels. Some of the new anion binding sites found in this study were, as expected, in pockets containing basic residues. However, some of them were near neutral, but polar, residues. Thus, the study also showed the importance of uncharged, but polar groups, on the protein surface in determining its electrostatic field. This was important for the second part of this study where the electrostatic field surrounding the protein was accurately determined. This was achieved by solving the linearized version of the Poisson-Boltzmann equation for the protein in solution. The solution was computed employing the commercial code Delphi which uses a finite difference technique. This has recently become available as a module in the general protein visualization code Insight II. Partial charges were assigned to the polar groups of lysozyme for the calculations done here. The calculations showed the complexity of the electrostatic field surrounding the protein. Although most of the region near the protein surface had a positive field strength, the active site cleft was negatively charged and this was projected a considerable distance. This might explain the occurrence of "head-to-side" interactions in the formation of lysozyme aggregates in solution. Pockets of high positive field strength were also found in the vicinity of the anion locations obtained from the crystallographic part of this study, confirming the validity of these calculations. This study clearly shows not only the importance of determining the counterion locations in protein crystals and the electrostatic fields surrounding the protein, but also the advantage of performing them together.

Lim, Kap↗
NASA NTRS1996

Electrical double layers at the oil/water interface

This review presents the historical development and current status of the theory of the electrical double layer at a liquid/liquid interface. It gives rigorous thermodynamic definitions of all basic concepts related to liquid interfaces and to the electrical double layer. The difference between the surface of a solid electrode and the interface of two immiscible electrolyte solutions (ITIES) is analyzed in connection to their electrical properties. The most important classical relationships for the electrical double layer are presented and critically discussed. The generalized adsorption isotherm is derived. After a short review of the classical Gouy-Chapman and Verwey-Niessen models, more recent developments of the double layer theory are presented. These include effects of variable dielectric permittivity, nonlocal electrostatics, hydration forces, the modified Poisson-Boltzmann equation and the ion-dipole plasma. The relative merits of different theories are estimated by comparing them with computer simulation of the ITIES and electrical double layer. Special attention is given to the structure of ITIES and its variation due to adsorption of ions and amphiphilic molecules.

Non-NASA Center↗
DOE OSTI2024

Nonlinear Poisson–Boltzmann solutions for charged parallel plates: When opposite charges repel

I present an exact solution of the Poisson–Boltzmann equation for two parallel plates and discuss the solution properties. I discuss in more detail plates with opposite charges: In this case, there are two critical separations, L c,1 < L c,2 . For separations less than L c,1 , the force between plates is repulsive. It switches to attractive at L c,1 , but with the electric potential having the same sign on both plates. For L > L c,2 , the force remains attractive, and the potential at the plates has the same sign as the charge on each plate. I also describe charge regulation, determined by pK a , and provide formulas for both the critical distance where oppositely charged plates repel and their charging process. Finally, the implications of these results for the nanoparticle assembly, as driven by electrostatic interactions, are also discussed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Proton diffusion and hydrogen/deuterium exchange in amorphous solid water at temperatures from 114 to 134 K

The reaction coefficient for hydrogen/deuterium (H/D) exchange and the diffusion of hydrated excess protons within amorphous solid water (ASW) are characterized as a function of temperature. For these experiments, water films are deposited on a Pt(111) substrate at 108 K, and reactions with pre-adsorbed hydrogen atoms produce hydrated protons. Upon heating, protons diffuse within the water, and H/D exchange occurs when they encounter D2O probe molecules deposited in the films. The time-dependent concentration of D2O is monitored with infrared spectroscopy, and it indicates the protons diffusion from the substrate and establish an equilibrium distribution prior to significant H/D exchange for temperatures 114 K ≤T≤ 134 K. By controlling the distance between the D2O molecules and the substrate, we probe the distribution of protons within the film. It decays as x−2 for the examined range of x (12–52 nm) due to the electric field that develops between the diffusing protons and their image charges in the metal substrate. This agrees with the theoretical distance scaling for the equilibrated proton concentration in a dielectric near a metal boundary. From the proton concentration and the measured D2O decay rate, a lower bound for the proton diffusion coefficient ranging from 10−20 m2/s at 114 K to 10−18 m2/s at 134 K is estimated. The diffusion coefficient has an activation energy of 0.40 eV, which is comparable to energies reported for molecular translations and rotations of H2O, suggesting they may play a critical role in the proton diffusion mechanism within ASW.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Breakdown of electroneutrality in nanopores

Ion transport in extremely narrow nanochannels has gained increasing interest in recent years due to unique physical properties at the nanoscale and the technological advances that allow us to study them. It is tempting to approach this confined regime with the theoretical tools and knowledge developed for membranes and microfluidic devices, and naively apply continuum models, such as the Poisson-Nernst-Planck and Navier-Stokes equations. However, it turns out that some of the most basic principles we take for granted in larger systems, such as the complete screening of surface charge by counter-ions, can break down under extreme confinement. Additionally, we show that in a truly one-dimensional system of ions interacting with three-dimensional electrostatic interactions, the screening length is exponentially large, and can easily exceed the macroscopic length of a nanotube. Without screening, electroneutrality breaks down within the nanotube, with fundamental consequences for ion transport and electrokinetic phenomena. In this work, we build a general theoretical framework for electroneutrality breakdown in nanopores, focusing on the most interesting case of a one-dimensional nanotube, and show how it provides an elegant interpretation for the peculiar scaling observed in experimental measurements of ionic conductance in carbon nanotubes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗