SEARCH · Search NASA

Results for “electrostatic interactions”

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.

SourceResearch connectionFrom yearThrough year

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

At least 91 records · Page 5

All sources NASA NTRS DOE OSTI

DOE OSTI2017

Crosslinked anion exchange membranes with connected cations

The selective transport of ions has crucial importance in biological systems as well as modern-day energy devices, such as batteries and fuel cells, and water purification membranes. Control over ion movement can be exerted by ligation, ion channel dimensions, solvation, and electrostatic interactions. Polyelectrolyte hydrogels can provide aligned pathways for counter ion transport but lack mechanical integrity, while polyelectrolyte membranes typically suffer from the absence of an ion transport channel network. To develop polymer membranes for improved ion transport, we present the design of a novel material that combines the advantages of aligned pathways found in polyelectrolyte hydrogel and mechanical robustness in conventional membranes. The ionic species were organized via controlled copolymerization of a quaternizable monomer. Additionally, dimensional stability was then incorporated through a cast/crosslinking method to lock in the network of connected cationic groups. This strategy resulted in dramatically enhanced ion transport, as characterized by ionic conductivities (>80 mS/cm for Cl-, and approximate to 200 mS/cm for OH-).

Zhang, Wenxu↗

DOE OSTI2018

In-Situ 27Al NMR Spectroscopy of Aluminate in Sodium Hydroxide Solutions Above and Below Saturation with Respect to Gibbsite

Aluminum hydroxide (Al(OH)3, gibbsite) dissolution and precipitation processes in alkaline environments play a commanding role in aluminum refining and nuclear waste processing, yet mechanistic aspects underlying sluggish kinetics during crystallization have remained obscured due to a lack of in situ probes capable of isolating incipient ion pairs. At a molecular level Al is cycling between tetrahedral (Td) coordination in solution to octahedral (Oh) in the solid. We explored dissolution of Al(OH)3 that prepared variably saturated aluminate (Al(OH)4 -)-containing solutions under alkaline conditions (pH > 13) with in-situ 27Al magic angle spinning (MAS)-nuclear magnetic resonance (NMR) spectroscopy, and interrogated the results with ab initio molecular dynamics (AIMD) simulations complemented with chemical shift calculations based on high level theory. The collective results highlight the overall stability of the solvation structure for Td Al in the Al(OH)4 - oxyanion as a function of both temperature and Al concentration. The observed chemical shift did not change significantly even when the Al concentration in solution became supersaturated upon cooling and limited precipitation of the octahedral Al(OH)3 phase occurred. However, subtle changes in Al(OH)4 - speciation correlated with the dissolution/precipitation reaction were found. AIMDinformed chemical shift calculations indicate that measurable perturbations should begin when the Al(OH)4 -?Na+ distance is less than 6 Å, increasing dramatically at shorter distances, coinciding with appreciable changes to the electrostatic interaction and reorganization of the Al(OH)4 - solvation shell. The integrated findings thus suggest that, at conditions incipient to and concurrent with gibbsite crystallization, nominally expected contact ion pairs are insignificant and instead medium-range (4 – 6 Å) solvent-separated Al(OH)4 -?Na+ pairs predominate. Moreover, the fact that these medium-range interactions bear directly on resulting gibbsite characteristics was demonstrated by detailed microscopic and X-ray diffraction analysis, and by progressive changes in the FWHM of the Oh resonance, as measured by in-situ NMR. Sluggish gibbsite crystallization may arise from the activation energy associated with disrupting this robust medium-range ion pair interaction.

Graham, Trenton R.↗

DOE OSTI2017

The origin of selectivity in the conversion of glucose to fructose and mannose in Sn-BEA and Na-exchanged Sn-BEA zeolites

We investigate the isomerization and epimerization of glucose to fructose and mannose in Sn-BEA andNa-exchanged Sn-BEA using density-functional theory calculations on periodic BEA crystals. We comparereaction pathways both in the absence and presence of water molecules in the vicinity of the active siteand find that water effectively determines the selectivity in Na-Sn-BEA. We identify two competingepimerization pathways, one involving direct 1,2-carbon shift and the other involving 1,2-hydride shiftvia fructose. In Sn-BEA, the isomerization to fructose is the kinetically dominant pathway, while mannoseis formed via the indirect 1,2-hydride shift epimerization pathway. In Na-Sn-BEA, the kinetically domi-nant pathway is epimerization to mannose via the direct 1,2-carbon shift pathway (Bilik mechanism)only in the presence of water solvent in the vicinity of the active site, whereas isomerization is preferredin the absence of water. We argue that polar water molecules that coordinate around the Na cation screenstrong electrostatic interactions between Na+and the glucose backbone that are responsible for thestrong inhibition of the 1,2-carbon shift mechanism in the absence of water. In Sn-BEA, the presenceof water does not influence the selectivity. Our calculations resolve for the first time the role of waterand Na cations in the catalytic activity of Sn-BEA, and rationalize the experimental data.

Li, Sha↗

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↗

DOE OSTI2018

Impact of long-range electrostatic and dispersive interactions on theoretical predictions of adsorption and catalysis in zeolites

In this article, we review the importance of long-range zeolite framework interactions in theoretical predictions for a variety of zeolite-catalyzed processes, and we show why such interactions must be determined accurately in order to reproduce experimentally measured adsorption and activation energies. We begin with an overview of the different strategies that have been used to account for long-range coulombic and dispersive interactions of zeolite framework atoms with species adsorbed at an active site. These methods include full periodic-DFT calculations and multi-layer hybrid techniques. Electrostatic interactions are observed to have a more significant impact than dispersive interactions on the geometries of ion-pair transition states and adsorbed species. Stabilization of the TS relative to reactant complexes is also dictated by electrostatic interactions. Dispersion effects are found to significantly stabilize both transition and reactant states for adsorbed species, especially those which have dimensions that provide good fits within the zeolite pore or cavity. We also show that the relevance of particular active site configurations can be missed, if the effects of long-range interactions are neglected. As a case in point, we demonstrate that a site previously considered inactive for ethane dehydrogenation, [GaH 2 ] + may in fact be more active than previously thought, when the impact of long-range interactions on the predicted activation energy is taken into account. In conclusion, the use of hybrid quantum mechanics/molecular mechanics approaches on extended, finite zeolite clusters has emerged as an accurate, highly cost-effective, and versatile alternative towards overcoming some of the present-day limitations of periodic calculations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2018

Probing Chromophore Energetics and Couplings for Singlet Fission in Solar Cell Applications: Final technical report

The overall goal of this project is to provide a molecular-level understanding of the mechanisms and requirements for efficient exciton multiplication through singlet fission processes. Such understanding would open the paths for improving alternative sources of energy including an improved efficiency of dye-sensitized solar cells. Using high-resolution anion photoelectron spectroscopy, electronic structures of chromophores are explored by probing directly the singlet and triplet states involved in singlet fission, determining precisely their energetics, vibrational structures and couplings. To carry out these experiments, we built a flexible-source instrument that allows for cryogenic cooling of precursor anions, providing improved spectral resolutions. Moreover, our velocity map imaging spectrometer includes novel design considerations that provide better focus with less stringent parameters. Using this instrument, we studied the three-ring polycyclic aromatic hydrocarbon molecules, anthracene and fluoranthene, the latter can be considered as a covalently-bonded benzyl and naphthyl system. Starting from the radical anion ground state, we accessed the singlet and triplet states at the same time, and the results revealed intriguing vibronic coupling signatures in both electronic states, which may possibly be related to singlet fission mechanisms. We further studied aromatic oxide molecules, phenoxy, 1-naphthoxy, and 2-naphthoxy, whose large molecular dipole moments supported highly delocalized dipole bound states. An intrinsic property of singlet fission is the delocalization of excitation to nearby chromophore molecules. In an effort to extend our studies to solvated chromophore clusters, we development new instrumental capabilities that allowed for formation of weakly bound clusters around any ion by utilizing liquid-nitrogen cooled octupole ion traps. In characterizing this new approach for gas phase cluster formation, we studied the structures of a variety of solvated clusters, from solvated ionic liquid cation, where intermolecular interactions are dominated by electrostatic interactions, to solvated peptides, where hydrogen-bonding with solvent molecules can directly alter the peptide structure. Overall, we found that the formation of weakly-bound clusters in a cryogenic ion trap can be a viable approach for accessing structures that better represent those found in condensed phase. In the process, we also implemented a straightforward double-resonance setup and established a protocol for isotopic substitution using cryogenic ion traps, both of which aided in decongesting crowded experimental spectra that often accompany these larger solvated systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Tunable exciton binding energy in 2D hybrid layered perovskites through donor–acceptor interactions within the organic layer

The strength of electrostatic interactions within semiconductors strongly affects their performance in optoelectronic devices. An important target is the tuning of a material’s exciton binding energy—the energy binding an electron–hole pair through the electrostatic Coulomb force—independent of its electronic band gap. Here, we report on the doping of a family of two-dimensional hybrid perovskites, in which inorganic lead halide sheets alternate with naphthalene-based organic layers, with tetrachloro-1,2-benzoquinone (TCBQ). For four out of seven n = 1 perovskites, the incorporation of the electron-accepting TCBQ dopant into the organic sublattice containing the electron-donating naphthalene species enabled the tuning of the materials’ 1s exciton binding energy. The naphthalene–TCBQ electron donor–acceptor interactions increased the electrostatic screening of the exciton, in turn lowering its binding energy relative to the undoped perovskite—by almost 50% in one system. Structural and optical characterization showed that the inorganic lattice is not significantly perturbed even though the layer-to-layer spacing increases upon molecular dopant incorporation.

36 MATERIALS SCIENCE↗

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↗

DOE OSTI2022

Driving force and pathway in polyelectrolyte complex coacervation

There is notable discrepancy between experiments and coarse-grained model studies regarding the thermodynamic driving force in polyelectrolyte complex coacervation: experiments find the free energy change to be dominated by entropy, while simulations using coarse-grained models with implicit solvent usually report a large, even dominant energetic contribution in systems with weak to intermediate electrostatic strength. Here, using coarse-grained, implicit-solvent molecular dynamics simulation combined with thermodynamic analysis, we study the potential of mean force (PMF) in the two key stages on the coacervation pathway for symmetric polyelectrolyte mixtures: polycation–polyanion complexation and polyion pair–pair condensation. We show that the temperature dependence in the dielectric constant of water gives rise to a substantial entropic contribution in the electrostatic interaction. By accounting for this electrostatic entropy, which is due to solvent reorganization, we find that under common conditions (monovalent ions, room temperature) for aqueous systems, both stages are strongly entropy-driven with negligible or even unfavorable energetic contributions, consistent with experimental results. Furthermore, for weak to intermediate electrostatic strengths, this electrostatic entropy, rather than the counterion-release entropy, is the primary entropy contribution. From the calculated PMF, we find that the supernatant phase consists predominantly of polyion pairs with vanishingly small concentration of bare polyelectrolytes, and we provide an estimate of the spinodal of the supernatant phase. Finally, we show that prior to contact, two neutral polyion pairs weakly attract each other by mutually induced polarization, providing the initial driving force for the fusion of the pairs.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Substituent and Heteroatom Effects on π–π Interactions: Evidence That Parallel-Displaced π-Stacking is Not Driven by Quadrupolar Electrostatics

Stacking interactions are a recurring motif in supramolecular chemistry and biochemistry, where a persistent theme is a preference for parallel-displaced aromatic rings rather than face-to-face π-stacking. This is usually explained in terms of quadrupole–quadrupole interactions between the arene moieties but that interpretation is inconsistent with accurate calculations, which reveal that the quadrupolar picture is qualitatively wrong. At typical π-stacking distances, quadrupolar electrostatics may differ in sign from an exact calculation based on charge densities of the interacting arenes. We apply symmetry-adapted perturbation theory to dimers composed of substituted benzene and various aromatic heterocycles, which display a wide range of electrostatic interactions, and we investigate the interplay of Pauli repulsion, dispersion, and electrostatics as it pertains to parallel-displaced π-stacking. Profiles of energy components along cofacial slip-stacking coordinates support a prominent role for the “van der Waals model” (dispersion in competition with Pauli repulsion), even for polar monomers where electrostatic interactions are significant. While electrostatic interactions are necessary to explain the optimal face-to-face π-stacking distance and to account for the relative orientation of one polar arene with respect to another, we find no evidence to support continued invocation of quadrupolar electrostatics as a basis for π-stacking. Our results suggest that a driving force for offset-stacking exists even in the absence of electrostatic interactions. Consequently, tuning electrostatics via functionalization does not guarantee that slip-stacking can be avoided. This has implications for rational design of soft materials and other supramolecular architectures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

Interfacial Adhesion Mechanism of Anionic Polyelectrolyte Brushes Induced by Oppositely Charged Macromolecular Counterions

Abstract Polyelectrolyte brushes are widely used as model systems for investigating electrostatic interactions at soft interfaces and to achieve exceptional lubrication properties. Previous studies have primarily focused on their behavior in the presence of multivalent counterions, which induce brush collapse, ionic crosslinking, and pronounced changes in interfacial structures. However, interactions between polyelectrolyte brushes and oppositely charged macromolecular counterions, such as polycations, remain poorly understood. Here, we prepare well-defined polystyrene sulfonate (PSS) brushes fabricated via surface-initiated grafting and employ surface forces apparatus measurements to investigate their interactions with oppositely charged polycations. In contrast to multivalent ions, polycations do not induce noticeable brush collapse, but instead generate significant adhesion between symmetric PSS brush layers. This adhesion increases with both contact time and applied load, eventually reaching a steady plateau, indicating that the interaction is governed not simply by electrostatic screening but by the gradual formation of polycation-mediated bridging under confinement. Furthermore, the introduction of monovalent Na+ ions disrupts the adhesive interaction even at very low concentrations, suggesting that the bridging function of the adsorbed polycation is highly sensitive to competitive ionic screening. In comparison, measurements with trivalent counterions reveal the expected brush collapse behavior but minimal dependence of adhesion on contact time, highlighting a clear mechanistic distinction from the polymeric counterion case. Collectively, these results demonstrate that molecular size, configurational restriction, and confinement-induced rearrangement, rather than charge valency alone, govern the interaction behavior of macromolecular counterions at brush interfaces. This work provides new insight into the molecular origins of adhesion and the regulation of interfacial interactions in charged polymer brush systems.

Park, Jinwoo [Argonne National Laboratory , , , ,;↗

DOE OSTI2025

Interactions between Per- and Polyfluoroalkyl Substances (PFAS) at the Water–Air Interface

Per- and polyfluoroalkyl substances (PFAS)--so-called “forever chemicals”--contaminate the drinking water of about 100 million people in the U.S. alone and are inefficiently removed by standard treatment techniques. A key property of these compounds that underlies their fate and transport and the efficacy of several promising remediation approaches is that they accumulate at the water−air interface. This phenomenon remains incompletely understood, particularly under conditions relevant to natural and treatment systems where water−air interfaces often carry significant loads of other organic contaminants or natural organic matter. To understand the impact of organic loading on PFAS adsorption, we carried out molecular dynamics simulations of PFAS at varying interfacial densities. We find that adsorbed PFAS form strong mutual interactions (attraction between perfluoroalkyl chains and electrostatic interactions among charged head groups) that give rise to ordered interfacial coatings. Furthermore, these interactions often involve near-cancellation of hydrophobic attraction and Coulomb repulsion. Our findings explain an apparent paradox whereby PFAS adsorption isotherms often suggest minimal mutual interactions while simultaneously displaying a high sensitivity to the composition and density of interfacial coatings. Consideration of the compounds present with PFAS at the interface has the potential to allow for more accurate predictions of fate and transport and the design of more efficient remediation approaches.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Unravelling the Glycan Code: Molecular Dynamics and Quantum Chemistry Reveal How O-Glycan Functional Groups Govern OgpA Selectivity in Mucin Degradation by Akkermansia muciniphila

Mucins, heavily O-glycosylated glycoproteins, are a key component of mucus, and certain gut microbiota, including Akkermansia muciniphila , can utilise mucin glycans as a carbon source. Akkermansia muciniphila produces the O-glycopeptidase enzyme OgpA, which cleaves peptide bonds at the N-terminus of serine (Ser) or threonine (Thr) residues carrying O-glycan substitutions, with selectivity influenced by the O-glycan functional groups. Using molecular dynamics (MD) simulations and quantum chemistry calculations, we explored how different O-glycan groups affect OgpA's selectivity. Our results show that peptides bind to the enzyme via hydrogen bonds, π–π interactions, van der Waals forces and electrostatic interactions, with key residues, including Tyr90, Val138, Gly176, Tyr210 and Glu91, playing important roles. The primary determinant of selectivity is the interaction between the peptide's functional group and the enzyme's binding cavity, while peptide–enzyme interface interactions are secondary. Quantum chemistry calculations reveal that OgpA prefers peptides with a lower electrophilic character. This study provides new insights into mucin degradation by gut microbiota enzymes, advancing our understanding of this critical biological process.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRSDate not supplied

Electrostatic Charging of the Lunar Surface

Lunar regolith dust particles accumulate charge and interact electrostatically with rover wheels, astronaut boots, and equipment. We have developed instrumentation for in situ measurements of the electrostatic charge developed by the interactions of lunar regolith dust with the space-rated materials on these devices. This instrument is also capable of measuring the distribution of electric fields on or near the lunar surface and the ion currents present near the lunar surface. We also report on our efforts to characterize the charging behavior of lunar dust in low gravity environments. This behavior is nonintuitive due to complex interactions between individual dust grains. We are developing an experiment to study this interaction in a microgravity vacuum environment. Better understanding of this interaction will allow for improved dust mitigation on the lunar surface.

electrostatics↗

NASA NTRSDate not supplied

Electrostatic charging of the lunar regolith

electrostatics↗

NASA NTRS1971

Long-range electrostatic and electromagnetic interactions between atoms

Long-range electrostatic and electromagnetic interactions between two neutral hydrogenic atoms in ground state

Cortez, J. L. M.↗

DOE OSTI2024

Elucidate the molecular basis of ampholytic chitosan as a high-performance cryoprotectant to myosin denaturation: The importance of saccharide charges

The uses of charged poly/oligo-saccharides enable the retarding of protein denaturation against various environmental stresses during food storage and manufacturing. However, at subzero temperatures, the molecular basis of such stabilization behaviors, i.e., cryoprotections, remain less explored. Here, in this study, we introduced an ampholytic saccharide, carboxymethyl chitooligosaccharide (CMCO) that effectively inhibited the freezing-induced myosin denaturation. The in-depth cryoprotective mechanism was systematically investigated by using molecular dynamic simulation and multispectral characterizations. Results showed that CMCO may interact with myosin through hydrogen bonding and electrostatic interactions, which caused the expelling of water at protein surfaces and the reduced conformational flexibility of myosin molecules. Due to this water replacement event, both secondary and ternary structures of myosin became freezing-resistant, leading to the inhibited protein aggregations and retained functionalities, such as solubility, Ca 2+ -ATPase activity, and gelling properties. Moreover, cryoprotective behaviors of CMCO were charge-dependent. CMCO with a higher degree of carboxymethyl substitution (DS: 1.2) was inclined to bind and stabilize myosin molecules better than the low-DS (DS: 0.8) one, even though both outperformed other cryoprotective saccharides. Therefore, this investigation not only introducing a high-performance myosin cryoprotectant, but also elaborated the cryoprotective mechanism of ampholytic saccharides.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2023

Molecular Dynamics Simulation Study of the Protonation State Dependence of Glutamic Acid Transport through a Cyclic Peptide Nanotube

The effect of the protonation state of glutamic acid on its translocation through cyclic peptide nanotubes (CPNs) was assessed by using molecular dynamics (MD) simulations. Anionic (GLU–), neutral zwitterionic (GLU0), and cationic (GLU+) forms of glutamic acid were selected as three different protonation states for an analysis of energetics and diffusivity for acid transport across a cyclic decapeptide nanotube. Based on the solubility-diffusion model, permeability coefficients for the three protonation states of the acid were calculated and compared with experimental results for CPN-mediated glutamate transport through CPNs. Potential of mean force (PMF) calculations reveal that, due to the cation-selective nature of the lumen of CPNs, GLU–, so-called glutamate, shows significantly high free energy barriers, while GLU+ displays deep energy wells and GLU0 has mild free energy barriers and wells inside the CPN. The considerable energy barriers for GLU– inside CPNs are mainly attributed to unfavorable interactions with DMPC bilayers and CPNs and are reduced by favorable interactions with channel water molecules through attractive electrostatic interactions and hydrogen bonding. Unlike the distinct PMF curves, position-dependent diffusion coefficient profiles exhibit comparable frictional behaviors regardless of the charge status of three protonation states due to similar confined environments imposed by the lumen of the CPN. The calculated permeability coefficients for the three protonation states clearly demonstrate that glutamic acid has a strong protonation state dependence for its transport through CPNs, as determined by the energetics rather than the diffusivity of the protonation state. In addition, the permeability coefficients also imply that GLU– is unlikely to pass through a CPN due to the high energy barriers inside the CPN, which is in disagreement with experimental measurements, where a considerable amount of glutamate permeating through the CPN was detected. To resolve the discrepancy between this work and the experimental observations, several possibilities are proposed, including a large concentration gradient of glutamate between the inside and outside of lipid vesicles and bilayers in the experiments, the glutamate activity difference between our MD simulations and experiments, an overestimation of energy barriers due to the artifacts imposed in MD simulations, and/or finally a transformation of the protonation state from GLU– to GLU0 to reduce the energy barriers. Altogether, our study demonstrates that the protonation state of glutamic acid has a strong effect on the transport of the acid and suggests a possible protonation state change for glutamate permeating through CPNs.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗