Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “solution species”

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.

At least 721 records · Page 40
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Optimizing CO 2 -Loaded Aqueous Amine Solutions for Higher Electrocatalytic CO 2 Reduction Activity

The activity of aqueous-based carbon dioxide reduction (CO 2 R) reactions is often limited by the solubility of CO 2 . The addition of amines can increase the total dissolved carbon in water through the formation of bicarbonate and carbamate species, which has been used to a great effect to capture CO 2 from dilute streams. Here, in this study, we explore the effect of 12 primary and secondary amines of varying Brønsted basicity, steric profile, and hydrogen-bonding capabilities on the aqueous CO 2 R to CO activity of a molecular Ni(cyclam)Cl 2 catalyst with a Hg electrode. Addition of some of the amines results in greater activity and selectivity for CO production compared to equivalent aqueous solutions without added amines. Under optimal conditions (0.4 M 3-amino-propionitrile), there is an over sevenfold increase in partial current density and greater selectivity for CO compared to equivalent conditions with no amine. Interestingly, the increase in activity did not correlate to any single property across the 12 amines. To elucidate the effect of the amine additives on catalysis, we used vapor–liquid equilibrium modeling (VLE), 13 C NMR spectroscopy, and computational analysis to determine the carbon speciation of the solutions. These results indicate that for amines without ethylalcohol functionalities, CO 2 R activity correlates with carbamate concentration, which is in turn governed by amine basicity and steric effects. However, this correlation does not persist for amines with ethylalcohol functionalities, which can form more stable carbamates through intramolecular-hydrogen bonding. These studies demonstrate that amine additives can enhance aqueous CO 2 R activity and selectivity and describe amine properties that lead to these higher performance metrics.

Amines↗
NASA NTRS1998

Improved Modeling of Finite-Rate Turbulent Combustion Processes in Research Combustors

The objective of this thesis is to further develop and test a stochastic model of turbulent combustion in recirculating flows. There is a requirement to increase the accuracy of multi-dimensional combustion predictions. As turbulence affects reaction rates, this interaction must be more accurately evaluated. In this work a more physically correct way of handling the interaction of turbulence on combustion is further developed and tested. As turbulence involves randomness, stochastic modeling is used. Averaged values such as temperature and species concentration are found by integrating the probability density function (pdf) over the range of the scalar. The model in this work does not assume the pdf type, but solves for the evolution of the pdf using the Monte Carlo solution technique. The model is further developed by including a more robust reaction solver, by using accurate thermodynamics and by more accurate transport elements. The stochastic method is used with Semi-Implicit Method for Pressure-Linked Equations. The SIMPLE method is used to solve for velocity, pressure, turbulent kinetic energy and dissipation. The pdf solver solves for temperature and species concentration. Thus, the method is partially familiar to combustor engineers. The method is compared to benchmark experimental data and baseline calculations. The baseline method was tested on isothermal flows, evaporating sprays and combusting sprays. Pdf and baseline predictions were performed for three diffusion flames and one premixed flame. The pdf method predicted lower combustion rates than the baseline method in agreement with the data, except for the premixed flame. The baseline and stochastic predictions bounded the experimental data for the premixed flame. The use of a continuous mixing model or relax to mean mixing model had little effect on the prediction of average temperature. Two grids were used in a hydrogen diffusion flame simulation. Grid density did not effect the predictions except for peak temperature and tangential velocity. The hybrid pdf method did take longer and required more memory, but has a theoretical basis to extend to many reaction steps which cannot be said of current turbulent combustion models.

VanOverbeke, Thomas J.↗
DOE OSTI2024

Experimental and Computational Study of Pyrogenic Carbonaceous Matter Facilitated Hydrolysis of 2,4-Dinitroanisole (DNAN)

This study investigated the reaction pathway of 2,4-dinitroanisole (DNAN) on the pyrogenic carbonaceous matter (PCM) to assess the scope and mechanism of PCM-facilitated surface hydrolysis. DNAN degradation was observed at pH 11.5 and 25 °C with a model PCM, graphite, whereas no significant decay occurred without graphite. Experiments were performed at pH 11.5 due to the lack of DNAN decay at pH below 11.0, which was consistent with previous studies. Graphite exhibited a 1.78-fold enhancement toward DNAN decay at 65 °C and pH 11.5 relative to homogeneous solution by lowering the activation energy for DNAN hydrolysis by 54.3 ± 3.9%. This is supported by our results from the computational modeling using Car–Parrinello simulations by ab initio molecular dynamics/molecular mechanics (AIMD/MM) and DFT free energy simulations, which suggest that PCM effectively lowered the reaction barriers by approximately 8 kcal mol –1 compared to a homogeneous solution. Quaternary ammonium (QA)-modified activated carbon performed the best among several PCMs by reducing DNAN half-life from 185 to 2.5 days at pH 11.5 and 25 °C while maintaining its reactivity over 10 consecutive additions of DNAN. We propose that PCM can affect the thermodynamics and kinetics of hydrolysis reactions by confining the reaction species near PCM surfaces, thus making them less accessible to solvent molecules and creating an environment with a weaker dielectric constant that favors nucleophilic substitution reactions. Nitrite formation during DNAN decay confirmed a denitration pathway, whereas demethylation, the preferred pathway in homogeneous solution, produces 2,4-dinitrophenol (DNP). Denitration catalyzed by PCM is advantageous to demethylation because nitrite is less toxic than DNAN and DNP. These findings provide critical insights for reactive adsorbent design that has broad implications for catalyst design and pollutant abatement.

2,4-dinitroanisole (DNAN)↗
DOE OSTI2022

Theoretical Pathway toward Improved Reverse Osmosis Membrane Selectivity for Neutral Solutes: Inspiration from Gas Separations

We report highly rigid membrane materials with tailored structures have exhibited permeabilities and selectivities that exceed the polymer upper bound in gaseous and organic solvent separations for challenging mixtures containing species that are similar in size and shape. One potential question is whether such membrane materials can maintain meaningful guest diffusivities in situations where the microporous spaces are essentially full of guest molecules. Here, we use a simplified transition state theory approach to estimate the diffusivity of water and small organics within a microporous membrane. The transition state theory model is parameterized using experimental values from zeolites and carbon molecular sieve materials found in the literature. We demonstrate the differences in transport and Maxwell–Stefan diffusivities based on guest species loading with different isotherm behaviors. These calculations theorize a path forward for highly selective reverse osmosis membranes for aqueous phase separations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1994

Computation of iodine species concentrations in water

During an evaluation of the use of iodine as a water disinfectant and the development of methods for measuring various iodine species in water onboard Space Freedom, it became necessary to compute the concentration of the various species based on equilibrium principles alone. Of particular concern was the case when various amounts of iodine, iodide, strong acid, and strong base are added to water. Such solutions can be used to evaluate the performance of various monitoring methods being considered. The authors of this paper present an overview of aqueous iodine chemistry, a set of nonlinear equations which can be used to model the above case, and a computer program for solving this system of equations using the Newton-Raphson method. The program was validated by comparing results over a range of concentrations and pH values with those previously presented by Gottardi for a given pH. Use of this program indicated that there are multiple roots to many cases and selecting an appropriate initial guess is important. Comparison of program results with laboratory results for the case when only iodine is added to water indicates the program gives high pH values for the iodine concentrations normally used for water disinfection. Extending the model to include the effects of iodate formation results in the computer pH values being closer to those observed, but the model with iodate does not agree well for the case in which base is added in addition to iodine to raise the pH. Potential explanations include failure to obtain equilibrium conditions in the lab, inaccuracies in published values for the equilibrium constants, and inadequate model of iodine chemistry and/or the lack of adequate analytical methods for measuring the various iodine species in water.

Schultz, John R.↗
DOE OSTI2020

A Principle for Highly Active Metal Oxide Catalysts via NaCl-Based Solid Solution

Toward the preparation of industrial metal oxide catalysts, sacrificial organic templates, excessive solvents, complex impregnation, and drying steps are generally required. Here, we report a versatile rule for the simple synthesis of highly porous metal oxides with well-dispersed noble metal species. Porous metal oxides (Co 3 O 4 , Fe x O y , and Cr 2 O 3 ) are obtained with some surface areas (e.g., Cr 2 O3: 224 m 2 ·g -1 ) beyond the record value. Surprisingly, small noble metal nanoparticles (e.g., Pd: 3.1 and Pt: 3.2 nm) could be incorporated by this solid-state process simultaneously. Corresponding Rh-Co 3 O 4 , Pd-Fe x O y , and Pt-Cr 2 O 3 exhibit excellent performance: CH 4 combustion (T90 = ~360°C and thermal stability: >100 h at 680°C), hydrogenation of nitrobenzene and derivatives (turnover number [TON] = 2.49 × 104, 300 mmol per run), and reversed water gas shift (RWGS) reaction (44% CO 2 conversion with ~98% CO selectivity and thermal stability: >100 h at 500°C), respectively. Therefore, current principle via a NaCl-based solid solution could provide a solid-state, fast, and efficient route for processing metal oxide catalysts.

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

In-situ polarization modulation IRRAS investigation of ammonia electrooxidation on Pt-Ir and Pt-Ru nanoparticles prepared on engineered catalyst supports

The catalytic activity and surface reactivity of monometallic Pt and bimetallic Pt-Ir and Pt-Ru nanoparticles, supported on two distinct Engineered Catalyst Supports (ECSs), were investigated for the Ammonia Electrooxidation Reaction (AmER) in alkaline media. XRD measurements confirmed alloy formation between Pt-Ir and Pt-Ru nanoparticles, as indicated by the shift of the (111) reflection to higher 2θ values. Cyclic voltammetry, linear sweep voltammetry, and chronoamperometry experiments were conducted to assess the catalytic activity of the Pt, Pt-Ir, and Pt-Ru electrocatalysts. All bimetallic catalysts exhibited lower onset potentials compared to Pt. The differing Tafel slopes between Pt (74 mV dec⁻¹), Pt-Ir (152 mV dec⁻¹), and Pt-Ru (118–197 mV dec⁻¹) suggest that alloying Pt with Ir or Ru alters the reaction mechanisms. Furthermore, the bimetallic Pt-Ir and Pt-Ru catalysts demonstrated greater tolerance for concentrated ammonia solutions relative to Pt. In-situ Polarization Modulation Infrared Reflection Absorption Spectroscopy (PM-IRRAS) provided insights into the formation of N-H species, azide anions (N₃⁻), and N-O compounds. For the Pt-Ru catalyst, an additional peak around ~3600 cm⁻¹ was observed, corresponding to OH⁻ species. The PM-IRRAS results align with the Gerischer–Mauerer mechanism, indicating that partially dehydrogenated ammonia adsorbates act as active intermediates in the oxidation of ammonia over Pt-Ir and Pt-Ru catalysts.

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

First-principles investigation of lanthanides diffusion in HCP zirconium via vacancy-mediated transport

The diffusion of lanthanide fission products plays an important role in the growth of the fuel-cladding chemical interaction (FCCI) region in metallic fuels. The use of a Zr interdiffusion barrier may mitigate the transport of lanthanides from the fuel to the cladding, but the efficacy of such a liner is not yet known. In this paper, the stability and vacancy-mediated diffusion of La, Ce, Pr, and Nd in hexagonal close-packed (HCP) Zr is investigated via density functional theory (DFT) calculations and self-consistent mean field (SCMF) analysis. DFT is used to calculate the formation, binding, and migration energies of vacancies and vacancy-solute pairs. The DFT energetics are used in the KineCluE code to calculate the Onsager transport coefficients. La is found to be the fastest diffusing species in HCP Zr and experiences an almost isotropic diffusion behavior. The other three species (Ce, Pr, and Nd) demonstrate anisotropic diffusion where the diffusion in the basal planes is significantly faster than that along the c-axis. Further, the calculated lanthanide diffusivities in HCP Zr are fitted to an Arrhenius relation and the activation energies and prefactors are reported for the first time. Furthermore, the vacancy drag and the segregation tendencies were analyzed using the calculated off-diagonal transport coefficients. According to our vacancy-mediated diffusion model, lanthanides are expected to be enriched at vacancy sinks at low temperatures, while at high temperatures, lanthanides are depleted at sinks and will preferably diffuse into the bulk. The enrichment/depletion transition temperature depends on the diffusion direction (basal or axial) and hence will be controlled by the grain texture and orientation.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Purification of Degraded Aqueous Piperazine by Ion Exchange and Carbon Treating

Amine purification by ion exchange and carbon treating has been performed for decades to maintain low solvent impurity levels in amines used in natural gas sweetening. The use of these purification methods to remove dissolved metal catalysts and other impurities from amines used in post-combustion carbon capture is less well known. This report presents batch experiments testing the equilibrium adsorption of impurities from degraded piperazine (PZ) onto activated carbon and various ion exchange resins. Compared to cationic exchange, anionic exchange was found to better remove catalytic Fe and UV active species, making it a better candidate for large-scale PZ purification. Carbon treating greatly outperforms anion exchange in the removal of UV active species, but the adsorbents’ difference is less pronounced for the removal of other impurities. Equilibrium adsorption was determined by the UV absorption of the PZ solution at 320 nm, and the results will show this may not be an accurate way to determine equilibrium. Further bench-scale experiments comparing ion exchange and carbon treating are necessary before advancing to larger-scale testing.

amine purification↗
NASA NTRS1988

Numerical simulation of hypersonic flow over an aeroassist flight experiment vehicle

Solutions of wind-tunnel and entry-flight flow around the vehicle are obtained from the Navier-Stokes equations coupled with the chemical species continuity equations if needed. The time-iterative method employs several techniques: shock fitting, chemistry-split ADI and an algebraic grid in conformal spherical-polar space. Sensitivities of the results to numerical parameters and to frozen, equilibrium and finite rate reactions are investigated in the forebody computation. Quantitative results are obtained for the shock layer and the near wake for the entire vehicle corresponding to both ground test and flight conditions. Complex flow characteristics are analyzed on the basis of the complete flowfield over the aerobrake and simplified afterbodies. The method is stable and cost effective, and has yielded shock locations and wall pressure distributions which are in good agreement with wind-tunnel data.

Li, C. P.↗
NASA NTRSDate not supplied

Numerical Simulation of Lean Blowout of Alternative Fuels in 7-element Lean Direct Injector

This research presents the result of numerically simulating 7-element swirl-venturi Lean Direct Injector (SV-LDI) lean blowout (LBO) experiments conducted at NASA Glenn Research Center in May of 2019. After simulating a cold flow case to confirm the pressure drop agrees well with the experiment, additional cases with two different fuels (an average jet fuel and a Gevo alcohol-to-jet fuel) from the National Jet Fuels Combustion Program (NJFCP) were computed to numerically determine the LBO condition. The procedure to approach the LBO follows the method used in the experiment where the air mass flow rate is gradually increased while the fuel supply is maintained. Transient history of global heat release rate as a function of air flow rate is presented, as well as temperature contours at different conditions to give a visual representation of the flame state. The Open National Combustion Code (OpenNCC) used in this research adopted reduced HyChem (Hybrid Chemistry) models along with k-LES turbulence model and a Lagrangian spray model that takes into account droplet internal temperature distribution affected by the shear force on the droplet surface. The transport equations and chemical reaction terms are integrated together to enhance conservation of chemical species that are especially important in the near LBO conditions. After showing the computed range of LBO agrees well with the experimental measurements, time averaged solutions of both fuels at their initial condition and at their limiting condition just before LBO are compared in detail to facilitate the understanding of LBO mechanism.

combustion↗
NASA NTRS2023

Capitol Reef Ecological Conservation: Mapping Vegetation Functional Groups to Inform Invasive Vegetation Management, Ecological Conservation and Restoration in Capitol Reef National Park

Invasive exotic plant (IEP) species have been found within the park boundaries of Capitol Reef National Park (CARE) in Utah. Currently, remotely sensed datasets such as the Rangeland Analysis Platform (RAP) from the United States Department of Agriculture (USDA) have been used to investigate IEP species within the park, but validation of the national RAP program is necessary for informing decisions at a local scale. CARE seeks a remote monitoring solution that can precisely target managerial efforts within the park’s challenging terrain and hard-to-reach locations. To fulfill this objective, we harnessed Landsat 8 Operational Land Imager (OLI) imagery and leveraged Random Forest (RF) modeling to generate classification maps characterizing vegetation functional groups for 2013 and 2022 within the park. Subsequently, the Land Change Modeler (LCM) in Idrisi TerrSet facilitated the production of a predicted classification map for 2033. The team also devised an annual grass probability map to accentuate areas impacted by exotic grasses. A comparative assessment between the RF classification map and the RAP map for 2022 revealed an overall agreement of 47.41%, with disparities primarily arising from differences in bare soil and shrub areas. Significantly, the 2022 RF-generated classification map showcased an impressive overall accuracy of 92.17%. In short, the probability map, the land cover change detection spanning 2013 to 2022, and the forecasting of observed trends into the future aids in the evaluation of invasive plant impacts and facilitation of CARE’s preparedness for potential ecological disturbances. Notably, in comparison to the RAP, the RF classification method generates functional group maps that are more representative of the study area.

Vanchy Li↗
DOE OSTI2025

PFSA-Ionomer Dispersions to Thin-Films: Interplay Between Sidechain Chemistry and Dispersion Solvent

In fuel-cell catalyst layers (CLs), ionomers exists as nanoscale thin-films binding catalyst particles, wherein ionic and gaseous species transport to catalyst sites. To improve film function, ionomers have been designed based on the prototypical perfluorosulfonic acid (PFSA). Since CLs are fabricated through solution-based processing of inks, understanding how ink/dispersion solvent impacts ionomer interactions and function is important to guide future ionomer design. Herein, PFSA and its chemical derivatives are characterized from dispersion (in solvent) to thin film (cast on support) to investigate the impact of solvent water content on ionomer dispersion (structure, acidity) and resulting film properties. Dispersion characterization reveals that secondary aggregate structure depends on the sidechain chemistry, but all PFSA-based ionomers evolve similarly as dispersion becomes water-rich. The differences in aggregation alter ionomer adsorption and thin film morphology. Further, hydrophilic domain spacing and orientation are affected primarily by sidechain chemistry, and modulated by solvent composition. Thin-film conductivity correlates with hydration and nanodomain alignment, signifying the role of morphology in properties. Overall, this work provides insights into how chemistry can be used to coarse-tune structure–property relationships, while processing solvent is for fine-tuning of dispersion-film interplay for controlling CL-ionomer function, which can be translated to other chemistries and ink formulations.

25 ENERGY STORAGE↗
DOE OSTI2024

Elucidating Phosphate and Cadmium Cosorption Mechanisms on Mineral Surfaces with Direct Spectroscopic and Modeling Evidence

The simultaneous sorption of cations and anions at the mineral–water interface can substantially alter their individual sorption characteristics; however, this phenomenon lacks a mechanistic understanding. Our study provides direct spectroscopic and modeling evidence of the molecular cosorption mechanisms of the cadmium ion (Cd 2+ ) and phosphate (P) on goethite and layered manganese (Mn) oxide of birnessite, through in situ attenuated total reflection Fourier-transform infrared (ATR-FTIR), P K-edge X-ray absorption near-edge structure (XANES) spectroscopy, and surface complexation modeling. Phosphate synergistically cosorbed with Cd on goethite predominantly through P-bridged ternary complexes (≡Fe–P–Cd) and electrostatic interactions at wide pH conditions. Likewise, P and Cd exhibited synergistic cosorption on birnessite by forming P-bridged ternary complexes (≡Mn–P–Cd) and weak competitive sorption at the layer edge sites. As pH and Cd loading increased, the surface P species transitioned from a binary complex to a ternary complex and/or Cd 3 (PO 4 ) 2 precipitate for both goethite and birnessite. Compared to that in solution at pH 8, the formation of Cd 3 (PO 4 ) 2 was inhibited by the presence of goethite and birnessite, ascribed to the specific adsorption of P and Cd, more pronounced in birnessite due to the stronger sorption of Cd at its vacant sites. Finally, the discovered cosorption mechanisms of P and Cd have important implications for understanding and predicting their mobility and availability in Cd-contaminated settings.

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

Strong Enhancement of Nanoconfined Water Mobility by a Structure Breaking Salt

For the majority of the water present on earth, the two most important factors influencing its behavior are confinement, in either inorganic or organic matrixes, and the presence of solutes. Here, we investigate the effect of confinement in 3 nm pores on water diffusivity in aqueous solutions with archetypical solutes, a structure making (kosmotrope) NaCl and a structure breaking (chaotrope) KCl, up to 1.0 M in concentration. The water diffusivity in bulk aqueous solutions in such a concentration range is known to decrease very slightly in the presence of NaCl and increase very slightly in the presence of KCl. However, here we observe the water diffusivity in confined H 2 O–KCl increases by a factor of 2 compared to the pure water diffusivity in the same confinement. This unusually strong cumulative effect of confinement and a structure breaking additive may have profound implications for the mobility and transport of aqueous species in nature.

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

Approaches for the Simulation of Coupled Processes in Evolving Fractured Porous Media Enabled by Exascale Computing

Models have historically represented fractured porous media with continuum descriptions that characterize the media using bulk parameters. The impact of small-scale features is not captured in these models, although they may be controlling the performance of subsurface applications. Pore-scale models can simulate processes in small-scale features by representing the pore space geometry explicitly but are computationally expensive for large domains. The alternative multiscale approach entails the combination of pore-scale and continuum-scale descriptions in a single framework. We use Chombo-Crunch, a computational capability that discretizes complex geometries with an adaptive, embedded boundary method to contrast these two approaches. Chombo-Crunch takes advantage of recent computational performance and memory bandwidth improvements resulting from the emergence of exascale computing resources. These combined improvements enable the efficient simulation of reactive transport in fractured media with a high degree of fidelity and the ability to capture the control small-scale processes exert on the overall medium evolution.

42 ENGINEERING↗
DOE OSTI2020

Platform for efficient large-scale storage and analysis of multi-omics data in plant and microbial systems (Final Technical Report)

Genomic variation at the sequence level fundamentally affects the phenotypic state of all organisms at all stages of development, while dynamic processes such as changes in the epigenome (e.g. DNA methylation state) and transcriptome regulate the specific phenotype expressed at any given state of development based upon that genomic variation. In plants, DNA methylation is a particularly important mechanism for both regulating transcriptomic expression and for management of genomic variations that could be deleterious to the organism due to the presence of active retrotransposons in plant genomes. While DNA methylation is heritable, it is also dynamic through a given plant’s development and life cycle, particularly during the development from seed to mature specimen suggesting variations in DNA methylation could be critical regulators of biologically and commercially important phenotypes such as time to flowering; in addition, plant DNA methylation is more complex than that of animals, with methylation of CHG and CHH trinucleotides evident in addition to the better-known CG methylation. The complexity of plant DNA methylation and its interplay with genomic sequence variation, transcriptomics and other epigenomic factors demand a storage and analysis framework that can cope with the complexity both within a single specimen and with analyses that span many individuals and even many species, such as attempts to extend models from model organisms to commercially relevant species. In addition to complexity, the rapid development and proliferation of sequencing technology has led to an explosion of data volume that conventional storage and analysis solutions will likely be unable to cope with in the long run. We proposed to study these with suitable distributed storage and computation and therefore for the application of cloud computing to biological analyses; integrate with existing data sources and compatible with virtually any interface use case, from fully automated shell scripts to notebooks and do all these at scale in this STTR grant.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

Acid-Base Equilibrium and Dielectric Environment Regulate Charge in Supramolecular Nanofibers

Peptide amphiphiles are a class of molecules that can self-assemble into a variety of supramolecular structures, including high-aspect-ratio nanofibers. It is challenging to model and predict the charges in these supramolecular nanofibers because the ionization state of the peptides are not fixed but liable to change due to the acid-base equilibrium that is coupled to the structural organization of the peptide amphiphile molecules. Here, we have developed a theoretical model to describe and predict the amount of charge found on self-assembled peptide amphiphiles as a function of pH and ion concentration. In particular, we computed the amount of charge of peptide amphiphiles nanofibers with the sequence C 16 − V 2 A 2 E 2 . In our theoretical formulation, we consider charge regulation of the carboxylic acid groups, which involves the acid-base chemical equilibrium of the glutamic acid residues and the possibility of ion condensation. The charge regulation is coupled with the local dielectric environment by allowing for a varying dielectric constant that also includes a position-dependent electrostatic solvation energy for the charged species. We find that the charges on the glutamic acid residues of the peptide amphiphile nanofiber are much lower than the same functional group in aqueous solution. There is a strong coupling between the charging via the acid-base equilibrium and the local dielectric environment. Our model predicts a much lower degree of deprotonation for a position-dependent relative dielectric constant compared to a constant dielectric background. Furthermore, the shape and size of the electrostatic potential as well as the counterion distribution are quantitatively and qualitatively different. These results indicate that an accurate model of peptide amphiphile self-assembly must take into account charge regulation of acidic groups through acid–base equilibria and ion condensation, as well as coupling to the local dielectric environment.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗