Determining Ion Activity Coefficients in Ion-Exchange Membranes with Machine Learning and Molecular Dynamics Simulations
Not Available
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
SEARCH · Search NASA
Results for “Activity coefficient”
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.
Nonrandom two‐liquid activity coefficient model with association theory
Abstract Considered one of the most versatile and widely used classical thermodynamic models to correlate phase equilibria behavior of nonideal systems, the nonrandom two‐liquid (NRTL) theory does not explicitly account for specific chemical associations such as hydrogen bonding. This deficiency has been recognized as the cause for unsatisfactory representation of association systems like methanol–alkanes binaries. This work presents a practical approach to integrate Wertheim's perturbation theory for association contribution with the classical NRTL model. Specifically, the association contribution is calculated from pre‐determined molecule‐specific association strengths while the physical interaction contribution is captured with NRTL binary interaction parameters. The resulting association NRTL model correlates fluid phase equilibria of association systems with few adjustable parameters and offers improved predictive capability for higher order systems.
Hao, Yifan↗
Salt Activity Coefficient and Chain Statistics in Poly(ethylene oxide)Based Electrolytes
Explore the source record for details and available documents.
Fang, Chao↗
Effect of Dielectric Saturation on Ion Activity Coefficients in Ion Exchange Membranes
Not Available
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Effect of Silicon on Activity Coefficients of P, Bl, CD, SN, and AG in Liquid Fe-Si, and Implications for Differentiation and Core Formation
Cores of differentiated bodies (Earth, Mars, Mercury, Moon, Vesta) contain light elements such as S, C, Si, and O. We have previously measured small effects of Si on metal-silicate partitioning of Ni and Co [1,2], and larger effects for Mo, Ge, Sb, As [2]. The effect of Si on many siderophile elements could be an important, and as yet unquantified, influence on the core-mantle partitioning of SE. Here we report new experiments designed to quantify the effect of Si on the partitioning of Bi, Cd, Sn, Ag, and P between metal and silicate melt. The results will be applied to Earth, Mars, Mercury, Moon, and Vesta, for which we have excellent constraints on the mantle Bi, Cd, Sn, Ag, and P concentrations from mantle and/or basalt samples.
Righter, K.↗
Effect of Silicon on Activity Coefficients of P, Bi, Cd, Sn, and Ag in Liquid Fe-Si, and Implications for Core Formation
Cores of differentiated bodies (Earth, Mars, Mercury, Moon, Vesta) contain light elements such as S, C, Si, and O. We have previously measured small effects of Si on Ni and Co, and larger effects on Mo, Ge, Sb, As metal/silicate partitioning. The effect of Si on metal-silicate partitioning has been quantified for many siderophile elements, but there are a few key elements for which the effects are not yet quantified. Here we report new experiments designed to quantify the effect of Si on the partitioning of Bi, Cd, Sn, Ag, and P between metal and silicate melt. The results will be applied to Earth, Mars, Moon, and Vesta, for which we have good constraints on the mantle Bi, Cd, Sn, Ag, and P concentrations from mantle and/or basalt samples.
Righter, K.↗
Effect of Pressure and Oxidation State on Pnictogen Metal-Silicate Partitioning
The partitioning of trace elements between metal and silicate melts serves as the foundation for understanding the differentiation of a planetary body into a metallic core and silicate mantle. Element activity influences metal-silicate partitioning behavior. Activity coefficients are directly dependent on composition and temperature, and can be indirectly dependent on oxygen fugacity and pressure. Distinguishing the effect of pressure from other variables on the activity coefficients and partitioning is important for understanding the chemical evolution of different planetary bodies during differentiation. In this study, we investigated the influence of pressure on the activity coefficients of Cu, Mo, Pd, Pt, As, Sb, and Bi in Fe-Si metallic liquids. All of these elements exhibit moderate to high activity coefficients in Fe-Si liquids at low pressure, which significantly controls their metal-silicate partitioning behavior. Identifying whether this strong dependence persists at higher pressures is critical to modeling and understanding the chemical consequences of core formation. New experiments at 10 GPa were used to derive activity coefficients for these metals which can be compared to activity coefficients determined at 1 GPa. Experiments were conducted at 10 GPa and 2373K using a 10/5 assembly in the 880-ton multi anvil press at NASA Johnson Space Center. The standard 10/5 COMPRES assembly was slightly modified to accommodate a sample capsule machined from single-crystal MgO that minimizes melt percolation out of the sample volume during the experiment. Experiment starting materials were comprised of 70 wt.% Knippa basalt and 30 wt.% metal. The metal mixture (~85 wt.% Fe) was created by adding the elements of interest (Cu, Mo, Pd, Pt, As, Sb, Bi) to Fe metal powder. Varied amounts of Si metal (0-10 wt.% Si) were added to the metal-silicate mixtures to generate a systematic series of starting materials. For each element, an epsilon interaction parameter in Fe-Si liquid was derived from the results of our 0-10 wt.% Si metal series. To investigate whether pressure influences the trace element activity coefficients in the 1-10 GPa pressure range, we compared our results at 10 GPa to those at 1-4 GPa [1-3]. Our results can also be directly compared to interaction parameters for Au, P, V, Mn, Ga, Zn, Cd, Sn, W, Pb, and Nb previously determined at 10 GPa and 2373K following the same methods [4]. Combined, this suite of interaction parameters will directly inform metal- silicate partitioning between 1 and 10 GPa, and assess whether these values can be extrapolated to modeling differentiation processes at pressures >10 GPa. [1] Righter et al. (2018) GCA 232, 101-123 [2] Righter et al. (2019) MaPS 54, 1379-1394 [3] Steenstra et al. (2020) Icarus, 335, 113391 [4] Righter et al. (2020) Geochem. Persp. Let. 15, 44-49
Kelsey Prissel↗
Modeling ethanol/water adsorption in all-silica zeolites using the real adsorbed solution theory
A comprehensive set of single-component and binary isotherms were collected for ethanol/water adsorption into the siliceous forms of 185 known zeolites using grand-canonical Monte Carlo simulations. Using these data, a systematic analysis of ideal/real adsorbed-solution theory (IAST/RAST) was conducted and activity coefficients were derived for ethanol/water mixtures adsorbed in different zeolites based on RAST. It was found that activity coefficients of ethanol are close to unity while activity coefficients of water are larger in most zeolites, indicating a positive excess free energy of the mixture. This observation can be attributed to water/ethanol interactions being less favorable than water/water interactions in the single-component adsorption of water at comparable loadings. The deviation from ideal behavior can be highly structure-dependent but no clear correlation with pore diameters was identified. Furthermore, our analysis also demonstrates the following: (1) accurate unary isotherms in the low-loading regime are critical for obtaining physically sensible activity coefficients; (2) the global regression scheme to solve for activity model parameters performs better than fitting activity models to activity coefficients calculated locally at each binary state point; and (3) including the dependence on adsorption potential offers only a minor benefit for describing binary adsorption at the lowest fugacities. Finally, the Margules activity model was found incapable of capturing the non-ideal adsorption behavior over the entire range of fugacities and compositions in all zeolites, but for conditions typical of solution-phase adsorption, RAST predictions using zeolite-specific or even bulk Margules parameters provide an improved description compared to IAST.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Screening of metal complexes and organic solvents using the COSMOSAC-LANL model to enhance the energy density in a non-aqueous redox flow cell: an insight into the solubility
In this paper, we have proposed a first-principles methodology to screen transition metal complexes against a particular organic solvent and organic solvents against a particular transition metal complex based on their solubility information without the knowledge of heat of fusion and melting temperature. The energy density of a non-aqueous redox flow cell directly depends on the solubility of the redox active species in the non-aqueous medium. Here, we have used the “COSMOSAC-LANL” activity coefficient model which is based on first-principles COSMO calculations where the microscopic information is passed to the macroscopic world via a dielectric continuum solvation model, followed by a post-statistical thermodynamic treatment of the self-consistent properties of the solute particle to calculate the solubility. To model the activity coefficient at infinite dilution for the binary mixtures, a 3-suffix Margules (3sM) function is introduced for the quantitative estimation of the asymmetric interactions and, for the combinatorial term, the Staverman–Guggenheim (SG) form is used. The new activity coefficient model is separately called the “LANL” activity coefficient model. The metal complex and the organic solvent have been treated as a simple binary mixture. The present model has been applied to a set of 14 different organic solvents and 16 different transition metal complexes. Using the new LANL activity coefficient model in combination with the ADF-COSMOSAC-2013 model, we have shown how one can improve the solubility of a transition metal complex in an organic solvent. We applied our model to screen 84 binary mixtures to predict the compatible pair of redox active species and organic solvent to increase the energy density. The solvation mechanism of the transition metal complexes in the organic solvents was obtained using the new model. The results have been compared with the experimental and theoretical results where they are available.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Generalization of thermodynamic Langmuir isotherm for mixed‐gas adsorption equilibria
Abstract This work presents a comprehensive thermodynamic model for both pure component isotherms and mixed‐gas adsorption equilibria. A generalization of thermodynamic Langmuir isotherm, the proposed model assumes competitive adsorption of multiple adsorbates on adsorbent surface for mixed‐gas adsorption equilibria, and it applies an area‐based adsorption nonrandom two‐liquid activity coefficient model in the activity coefficient calculations for the adsorbate phase. The resulting generalized Langmuir isotherm properly captures both surface loading dependence and adsorbate phase composition dependence for mixed‐gas adsorption equilibria. The model is validated with accurate representations of gas adsorption equilibrium data for varieties of unary, binary, and ternary gas systems. The model results are further compared with those calculated from extended Langmuir isotherm and Ideal Adsorbed Solution Theory.
Hamid, Usman↗
Thermodynamic modeling of adsorption at the liquid-solid interface
Adsorptive separation techniques are significantly energy efficient in comparison to conventional thermal separation techniques such as distillation. Despite extensive research and development activities undertaken for mixed-gas adsorption, the use of adsorption techniques for the separation of multicomponent liquid mixtures is still limited. A major barrier is the lack of accurate adsorption thermodynamic models, which form the scientific foundation of process simulation of such systems, making the translation to industrial scale challenging. In this work, we have rigorously computed the surface excess of adsorption for six binary liquid mixtures on silica gel at 303 K using the frameworks of the generalized Langmuir isotherm (gL) and the adsorbed solution theory (AST). The six binary liquid mixtures were formed by the pair-wise combinations of four components: benzene, 1,2-dichloroethane, cyclohexane, and n-heptane. We have based our calculations by considering simultaneous equilibria of three phases: saturated vapor phase, bulk liquid phase, and adsorbed phase. Further, the composition of the corresponding saturated vapor phase was estimated by the Nonrandom Two-Liquid activity coefficient model and experimental vapor-liquid equilibria data. The activity coefficients of the adsorbed phase, the central issue of multicomponent adsorption thermodynamics, were calculated using the adsorption Nonrandom Two-Liquid activity coefficient model. Devoid of simplifying assumptions, gL and AST provide rigorous thermodynamic frameworks for adsorption equilibria of multicomponent liquid mixtures.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Effect of Metal Chloride Impurities on Equilibrium Potential of Fe/FeCl 2 in Eutectic LiCl-KCl
In this study, interaction of dilute metal chloride solutes in molten eutectic LiCl-KCl was studied via equilibrium potential measurements. Open circuit potential (OCP) was measured between a pure iron rod and a Ag/AgCl reference electrode with concentration of FeCl 2 varying from 0.2–1 mol% at 530 °C. The Nernst equation was used to calculate activity coefficient of FeCl 2 based on OCP. Measurements were made in LiCl-KCl-FeCl 2 with and without additions of 0.5 mol% CsCl, NaCl, and LaCl 3 . The activity coefficient of FeCl 2 was found to be unaffected by these small concentrations of CsCl, LaCl 3 , and NaCl. However, there was observed evidence that the activity coefficient of FeCl 2 , above 0.8 mol% Fe, may be affected by the presence of LaCl 3 at higher concentrations beginning around 0.5 mol% La. Over the full range of conditions tested, the activity coefficient of FeCl 2 in molten LiCl-KCl ranged from 2.95 × 10 −5 to 1.01 × 10 −4 .
Electrochemistry↗
Machine Learning the COSMO Model for Predicting Thermodynamics of Electrolyte Mixtures
Bottom-up design of electrolyte mixtures for battery systems requires predicting macro thermodynamic properties from molecular constituents. For instance, molten salt electrolyte batteries require conditions far above room temperature to operate. Therefore, discovering mixtures with increasingly lower eutectic melting points is desirable. A model that can approximate chemical activity is a valuable tool to search through the vast compositional design space. Machine learning can predict properties of materials such as vibrational free energies, electronic energy gaps, and thermal conductivities. Moreover, they can learn physical models such as interatomic potentials. The COSMO-SAC model uses theory and empirical parameterization to predict liquid-vapor and liquid-solid properties using first-principles calculations. However, obtaining activity coefficients required for parameterizing the COSMO-SAC model is costly and limited to a select chemical space. In this work, we explored if machine learning methods could improve the COSMO-SAC model and bridge density functional theory calculations to liquid phase thermodynamic properties. Our data-driven approach uses existing databases for sigma-profiles of organic solvents and reconciles their methodological differences via ensemble averaging. First, an optimal machine learning model is constructed for each dataset. Our machine learning algorithms use the sigma-profile as an input feature to predict binary mixtures' activity coefficients using multi-output regression. Each dataset uses different choices of functionals, methods, and basis sets. Therefore, our ensemble model attempts to predict corrected activity coefficients given the combination of all the model outputs. The activity coefficients used for training are generated using the COSMO-SAC model. This approach enables the extraction of meaningful information from the existing datasets to improve the COSMO-SAC model for obtaining thermodynamic properties of electrolyte mixtures. With the liquid phase activities, we can identify electrolyte mixtures that meet desired phase equilibria conditions.
Thermodynamics↗
Impact of phase state and non-ideal mixing on equilibration timescales of secondary organic aerosol partitioning
Abstract. Evidence has accumulated that secondary organic aerosols (SOAs) exhibit complex morphologies with multiple phases that can adopt amorphous semisolid or glassy phase states. However, experimental analysis and numerical modeling on the formation and evolution of SOA still often employ equilibrium partitioning with an ideal mixing assumption in the particle phase. Here we apply the kinetic multilayer model of gas–particle partitioning (KM-GAP) to simulate condensation of semi-volatile species into a core–shell phase-separated particle to evaluate equilibration timescales of SOA partitioning. By varying bulk diffusivity and the activity coefficient of the condensing species in the shell, we probe the complex interplay of mass transfer kinetics and the thermodynamics of partitioning. We found that the interplay of non-ideality and phase state can impact SOA partitioning kinetics significantly. The effect of non-ideality on SOA partitioning is slight for liquid particles but becomes prominent in semisolid or solid particles. If the condensing species is miscible with a low activity coefficient in the viscous shell phase, the particle can reach equilibrium with the gas phase long before the dissolution of concentration gradients in the particle bulk. For the condensation of immiscible species with a high activity coefficient in the semisolid shell, the mass concentration in the shell may become higher or overshoot its equilibrium concentration due to slow bulk diffusion through the viscous shell for excess mass to be transferred to the core phase. Equilibration timescales are shorter for the condensation of lower-volatility species into semisolid shell; as the volatility increases, re-evaporation becomes significant as desorption is faster for volatile species than bulk diffusion in a semisolid matrix, leading to an increase in equilibration timescale. We also show that the equilibration timescale is longer in an open system relative to a closed system especially for partitioning of miscible species; hence, caution should be exercised when interpreting and extrapolating closed-system chamber experimental results to atmosphere conditions. Our results provide a possible explanation for discrepancies between experimental observations of fast particle–particle mixing and predictions of long mixing timescales in viscous particles and provide useful insights into description and treatment of SOA in aerosol models.
54 ENVIRONMENTAL SCIENCES↗
Guiding Principles for Geochemical/Thermodynamic Model Development and Validation in Nuclear Waste Disposal: A Close Examination of Recent Thermodynamic Models for H + —Nd 3+ —NO 3 - (—Oxalate) Systems
Development of a defensible source-term model (STM), usually a thermodynamical model for radionuclide solubility calculations, is critical to a performance assessment (PA) of a geologic repository for nuclear waste disposal. Such a model is generally subjected to rigorous regulatory scrutiny. In this article, we highlight key guiding principles for STM model development and validation in nuclear waste management. We illustrate these principles by closely examining three recently developed thermodynamic models with the Pitzer formulism for aqueous H + —Nd 3+ —NO 3 - (—oxalate) systems in a reverse alphabetical order of the authors: the XW model developed by Xiong and Wang, the OWC model developed by Oakes et al., and the GLC model developed by Guignot et al., among which the XW model deals with trace activity coefficients for Nd(III), while the OWC and GLC models are for concentrated Nd(NO 3 ) 3 electrolyte solutions. The principles highlighted include the following: (1) Principle 1. Validation against independent experimental data: A model should be validated against experimental data or field observations that have not been used in the original model parameterization. We tested the XW model against multiple independent experimental data sets including electromotive force (EMF), solubility, water vapor, and water activity measurements. The results show that the XW model is accurate and valid for its intended use for predicting trace activity coefficients and therefore Nd solubility in repository environments. (2) Principle 2. Testing for relevant and sensitive variables: Solution pH is such a variable for an STM and easily acquirable. All three models are checked for their ability to predict pH conditions in Nd(NO 3 ) 3 electrolyte solutions. The OWC model fails to provide a reasonable estimate for solution pH conditions, thus casting serious doubt on its validity for a source-term calculation. In contrast, both the XW and GLC models predict close-to-neutral pH values, in agreement with experimental measurements. (3) Principle 3. Honoring physical constraints: Upon close examination, it is found that the Nd(III)-NO 3 association schema in the OWC model suffers from two shortcomings. Firstly, its second stepwise stability constant for Nd(NO 3 ) 2+ (log K 2 ) is much higher than the first stepwise stability constant for NdNO 3 2+ (log K 1 ), thus violating the general rule of (log K 2 –log K 1 ) < 0, or $\frac{K1}{K2}$>1. Secondly, the OWC model predicts abnormally high activity coefficients for Nd(NO 3 ) 2 + (up to ~900) as the concentration increases. (4) Principle 4. Minimizing degrees of freedom for model fitting: The OWC model with nine fitted parameters is compared with the GLC model with five fitted parameters, as both models apply to the concentrated region for Nd(NO 3 ) 3 electrolyte solutions. The latter appears superior to the former because the latter can fit osmotic coefficient data equally well with fewer model parameters. The work presented here thus illustrates the salient points of geochemical model development, selection, and validation in nuclear waste management.
12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
Studies of Protein Solution Properties Using Osmotic Pressure Measurements
Examination of the protein crystallization process involves investigation of the liquid and solid state and a protein's properties in these states. Liquid state studies such as protein self association in solution by light scattering methods or other methods have been used to examine a protein Is properties and therefore its crystallization process and conditions. Likewise can osmotic pressure data be used to examine protein properties and various published osmotic pressure studies were examined by us to correlate osmotic pressure to protein solution properties. The solution behavior of serum albumin, alpha - chymotrypsin, beta - lactoglobulin and ovalbumin was examined over a range of temperatures, pH values and different salt types and concentrations. Using virial expansion and a local composition model the non ideal solution behavior in form of the activity coefficients (thermodynamic) was described for the systems. This protein activity coefficient data was related to a protein's solubility behavior and this process and the results will be presented.
Agena, S.↗
Thermodynamic modeling of aqueous acetic acid, butyric acid, and lactic acid solutions
Based on the activity coefficient – fugacity coefficient approach, a rigorous thermodynamic modeling study is presented for accurate correlation of vapor-liquid equilibrium data of aqueous solutions of acetic acid (293 to 391 K), butyric acid (325 to 436 K), lactic acid (378 to 409 K), and acetic acid + butyric acid binary mixture (358 to 421 K). In addition, the pH data of the three aqueous, single carboxylic acid solutions were measured at 298 to 328 K and successfully correlated. Given that these aqueous carboxylic acid solutions exhibit various degrees of association behavior in both vapor and liquid phases, the thermodynamic models considered for this study include the Redlich-Kwong equation of state (RK-EoS) and the Hayden-O’Connell equation of state (HOC-EoS) for the vapor phase fugacity coefficients and the electrolyte non-random two-liquid model (eNRTL) and the association electrolyte non-random two-liquid model (AeNRTL) for the liquid phase activity coefficients. The combination of the HOC-EoS for the vapor phase and the AeNRTL model for the liquid phase is found to provide the best correlation results, consistent with the fact that the HOC-EoS and the AeNRTL model explicitly account for association behaviors in the vapor phase and the liquid phase, respectively.
09 BIOMASS FUELS↗
Partitioning of Ni/2+/ between basaltic and synthetic melts and olivines An experimental study
Olivine/liquid distribution coefficients for Ni in basaltic and simple synthetic melts are found in the 1400-1070 C temperature range. The results suggest that: (1) the Ni distribution coefficients and NiO activity coefficients are strongly temperature dependent, (2) the Ni distribution coefficients are functions of melt composition and structure at a constant temperature although the dependence is small for most basaltic magmas, (3) differences in melt composition do not significantly influence generally 'basaltic' melts, (4) the Ni distribution coefficients may be used as geothermometers to predict the olivine crystallization temperatures of basaltic magmas to within plus or minus 50 C, (5) the Ni distribution coefficients are greater than one at all natural magmatic temperatures, and (6) the Ni distribution coefficients do not significantly vary as functions of the Ni content of the melt to at least one mole percent.
Leeman, W. P.↗