Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “anion diffusion”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Field-Driven Simulations to Probe the Impact of Ionic Correlations on Solution Transport Coefficients in Binary, Ternary, and Reciprocal Quaternary Aqueous Electrolytes

Ionic correlations play a critical role in governing transport properties of mixed-salt aqueous electrolyte solutions, yet their quantitative characterization remains challenging, particularly for multicomponent aqueous systems. Here, we develop a general nonequilibrium molecular dynamics framework to efficiently compute Onsager transport coefficients and ionic correlations in mixed-salt aqueous solutions. Using field-driven simulations, we obtain accurate Onsager matrices for LiCl/ KCl, KCl/KBr, and LiBr/KCl electrolyte solutions with significantly reduced computational cost relative to equilibrium Green−Kubo methods. The framework enables direct assessment of how attractive cation−anion and repulsive like-ion interactions contribute to conductivity and salt diffusivity across compositions. While static ion pairing exhibits strong composition dependence, dynamic ion correlations remain nearly invariant, leading to constant deviations from Nernst−Einstein predictions. These results highlight the disconnect between static ion association and dynamic transport correlations, and they establish a transferable approach for analyzing ion transport in complex electrolyte environments relevant to separation processes and electrochemical systems.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Dislocation loop evolution in Kr-irradiated ThO 2

Here, the early stage of microstructural evolution of ThO 2 , under krypton irradiation at 600, 800, and 1000°C, was investigated using in situ transmission electron microscopy (TEM). Dislocation loops grew faster, whereas their number density decreased with increasing irradiation temperature. Loop density was found to decrease with ion dose. Interstitial dislocation loops, including Frank loops with Burgers vector of $\textit{a}/3\langle111\rangle$ and perfect loops with Burgers vector of $\textit{a}/2\langle110\rangle$, were determined by traditional TEM and atomic resolution–scanning TEM techniques. Atomistic and mesoscale level modeling are performed to interpret experimental observations. The migration energy barriers of defects in ThO 2 were calculated using density-functional theory. The energetics of different dislocation loop types were studied using molecular dynamics simulations. Loop density and diameter were analyzed using a kinetic rate theory model that considers stoichiometric loop evolution. This analysis reveals that loop growth is governed by the mobility of cation interstitials, whereas loop nucleation is determined by the mobility of anion defects. Lastly, a rate theory model was used to extract the diffusion coefficients of thorium interstitials, oxygen interstitials, and vacancies.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

A simulation-based model studying monoethanolamine and aprotic heterocyclic anion ionic liquid (AHA-IL) mixtures for carbon capture

Capturing of waste CO 2 , particularly from large industrial point sources, is necessary for either permanent storage or CO 2 utilization. Amine solvents used for carbon capture have certain drawbacks, particularly the high energy required by the reboiler for solvent regeneration and degradation of the amine with absorption/regeneration cycling. These challenges may be overcome using ionic liquids (ILs) that can be tailored for operating under typical flue gas stream conditions (8–10% CO 2 , 18–20% H 2 O, 2–3% O 2 , and 67–72% N 2 ). ILs, however, oftentimes suffer from high viscosities, leading to diffusion limitations when compared to amine solvents. This work models the combinations of aprotic heterocyclic anion ILs (AHA-IL) with monoethanolamine (MEA) and studies the carbon capture effectiveness of these mixtures. Furthermore, this study found that all AHA-IL/MEA mixtures exhibited 32 % lower enthalpy of reaction (compared to MEA) and 64 % lower viscosity (compared to AHA-IL), which are desirable characteristics for CO 2 capture solvents.

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

Cations impact radical reaction dynamics in concentrated multicomponent aqueous solutions

Ultraviolet (UV) photolysis of nitrite ions (NO 2 – ) in aqueous solutions produces a suite of radicals, viz., NO · , O – , · OH, and · NO 2 . The O – and NO · radicals are initially formed from the dissociation of photoexcited NO 2 – . The O – radical undergoes reversible proton transfer with water to generate · OH. Both · OH and O – oxidize the NO 2 – to · NO 2 radicals. The reactions of · OH occur at solution diffusion limits, which are influenced by the nature of the dissolved cations and anions. Here, we systematically varied the alkali metal cation, spanning the range from strongly to weakly hydrating ions, and measured the production of NO · , · OH, and · NO 2 radicals during UV photolysis of alkaline nitrite solutions using electron paramagnetic resonance spectroscopy with nitromethane spin trapping. Comparing the data for the different alkali cations revealed that the nature of the cation had a significant effect on production of all three radical species. Radical production was inhibited in solutions with high charge density cations, e.g., lithium, and promoted in solutions containing low charge density cations, e.g., cesium. Through complementary investigations with multinuclear single pulse direct excitation nuclear magnetic resonance (NMR) spectroscopy and pulsed field gradient NMR diffusometry, cation-controlled solution structures and extent of NO 2 – solvation were determined to alter the initial yields of · NO and · OH radicals as well as alter the reactivity of NO 2 – toward · OH, impacting the production of · NO 2 . The implications of these results for the retrieval and processing of low-water, highly alkaline solutions that comprise legacy radioactive waste are discussed.

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

Enhancing CO 2 Transport Across a PEEK‐Ionene Membrane and Water‐Lean Solvent Interface

Abstract Efficient direct air capture (DAC) of CO 2 will require strategies to deal with the relatively low concentration in the atmosphere. One such strategy is to employ the combination of a CO 2 ‐selective membrane coupled with a CO 2 capture solvent acting as a draw solution. Here, the interactions between a leading water‐lean carbon‐capture solvent, a polyether ether ketone (PEEK)‐ionene membrane, CO 2 , and combinations were probed using advanced NMR techniques coupled with advanced simulations. We identify the speciation and dynamics of the solvent, membrane, and CO 2 , presenting spectroscopic evidence of CO 2 diffusion through benzylic regions within the PEEK‐ionene membrane, not spaces in the ionic lattice as expected. Our results demonstrate that water‐lean capture solvents provide a thermodynamic and kinetic funnel to draw CO 2 from the air through the membrane and into the bulk solvent, thus enhancing the performance of the membrane. The reaction between the carbon‐capture solvent and CO 2 produces carbamic acid, disrupting interactions between the imidazolium (Im + ) cations and the bistriflimide anions within the PEEK‐ionene membrane, thereby creating structural changes through which CO 2 can diffuse more readily. Consequently, this restructuring results in CO 2 diffusion at the interface that is faster than CO 2 diffusion in the bulk carbon‐capture solvent.

Walter, Eric D.↗
DOE OSTI2023

Design Strategies of Spinel Oxide Frameworks Enabling Reversible Mg-Ion Intercalation

In this study, reversible Mg 2+ intercalation in metal oxides frameworks is a key enabler for an operational Mg ion battery with high energy density needed for the next generation of energy storage technologies. While functional Mg-ion batteries have been achieved in structures with soft anions (e.g., S 2- and Se 2- ), they do not meet energy density requirements to compete with the current rechargeable lithium-ion batteries due to their low insertion potentials, emphasizing the necessity of finding an oxide-based cathode that operates at high potentials. A leading hypothesis to explain the limited availability of oxide Mg-ion cathodes is the belief that Mg 2+ has sluggish diffusion kinetics in oxides due to strong electrostatic interactions between the Mg 2+ ions and oxide anions in the lattice. From this assessment, it can be hypothesized that such rate limiting kinetic shortcomings can be mitigated by tailoring an oxide framework through creating less stable Mg 2+ - O 2- coordination.

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

Carbon chain abundance in the diffuse interstellar medium

Thanks to the mid-IR sensitivities of the ISO and IRTS orbiting spectrometers it is now possible to search the diffuse interstellar medium for heretofore inaccessible molecular emission. In view of the recent strong case for the presence of C(7-) (Kirkwood et al. 1998, Tulej et al. 1998),and the fact that carbon chains possess prominent infrared active modes in a very clean portion of the interstellar spectrum, we have analyzed the IRTS spectrum of the diffuse interstellar medium for the infrared signatures of these species. Theoretical and experimental infrared band frequencies and absolute intensities of many different carbon chain species are presented. These include cyanopolyynes, neutral and anionic linear carbon molecules, and neutral and ionized, even-numbered, hydrogenated carbon chains. We show that--as a family--these species have abundances in the diffuse ISM on the order of 10(-10) with respect to hydrogen, values consistent with their abundances in dense molecular clouds. Assuming an average length of 10 C atoms per C-chain implies that roughly a millionth of the cosmically available carbon is in the form of carbon chains and that carbon chains can account for a few percent of the visible to near-IR diffuse interstellar band (DIB) total equivalent width (not DIB number).

Carbon/analysis/chemistry↗
NASA NTRS2014

Computational and Experimental Investigation of Li-doped Ionic Liquid Electrolytes: [pyr14][tfsi], [pyr13][fsi], and [EMIM][BF4]

We employ molecular dynamics (MD) simulation and experiment to investigate the structure, thermodynamics, and transport of N-methyl-N-butylpyrrolidinium bis(trifluoromethylsufonyl)imide ([pyr14][TFSI]), N -methyl-N-propylpyrrolidinium bis(fluorosufonyl)imide ([pyr13][FSI]), and 1-ethyl-3-methylimidazolium boron tetrafluoride ([EMIM][BF4]), as a function of Li-salt mole fraction (0.05 xLi+ 0.33) and temperature (298 K T 393 K). Structurally, Li+ is shown to be solvated by three anion neigh- bors in [pyr14][TFSI] and four anion neighbors in both [pyr13][FSI] and [EMIM][BF4], and at all levels of xLi+ we find the presence of lithium aggregates. Pulsed field gradient spin-echo NMR measurements of diffusion and electrochemical impedance spectroscopy measurements of ionic conductivity are made for the neat ionic liquids as well as 0.5 molal solutions of Li-salt in the ionic liquids. Bulk ionic liquid properties (density, diffusion, viscosity, and ionic conductivity) are obtained with MD and show excellent agreement with experiment. While the diffusion exhibits a systematic decrease with increasing xLi+, the contribution of Li+ to ionic conductivity increases until reach- ing a saturation doping level of xLi+ 0.10. Comparatively, the Li+ conductivity of [pyr14][TFSI] is an order of magnitude lower than that of the other liquids, which range between 0.1-0.3 mScm. Our transport results also demonstrate the necessity of long MD simulation runs ( 200 ns) required to converge transport properties at room T. The differences in Li+ transport are reflected in the residence times of Li+ with the anions (Li), which are revealed to be much larger for [pyr14][TFSI] (up to 100 ns at the highest doping levels) than in either [EMIM][BF4] or [pyr13][FSI]. Finally, to comment on the relative kinetics of Li+ transport in each liquid, we find that while the net motion of Li+ with its solvation shell (vehicular) significantly contributes to net diffusion in all liquids, the importance of transport through anion exchange (hopping) increases at high xLi+ and in liquids with large anions.

electrolytes↗
DOE OSTI2022

A reaction mechanism for plasma electrolysis of AgNO 3 forming silver nanoclusters and nanoparticles

In plasma-driven solution electrolysis (PDSE), gas-phase plasma-produced species interact with an electrolytic solution to produce, for example, nanoparticles. An atmospheric pressure plasma jet (APPJ) directed onto a liquid solution containing a metallic salt will promote reduction of metallic ions in solution, generating metallic clusters that nucleate to form nanoparticles. In this article, results from a computational investigation are discussed of a PDSE process in which a radio-frequency APPJ sustained in helium impinges on a silver nitrate solution, resulting in growth of silver nanoparticles. A reaction mechanism was developed and implemented in a global plasma chemistry model to predict nanoparticle growth. To develop the reaction mechanism, density functional theory was used to generate probable silver growth pathways up to Ag 9 . Neutral clusters larger than Ag 9 were classified as nanoparticles. Kinetic reaction rate coefficients for thermodynamically favorable growth pathways were estimated based on an existing, empirically determined base reaction mechanism for smaller Ag particle interactions. These rates were used in conjunction with diffusion-controlled reaction rate coefficients that were calculated for other Ag species. The role of anions in reduction of Agn ions in forming nanoparticles is also discussed. Oxygen containing impurities or admixtures to the helium, air entrainment into the APPJ, and dissociation of saturated water vapor above the solution can produce additional reactive oxygen species in solution, resulting in the production of anions and O 2 – in particular. For a given molarity, delivering a sufficient fluence of reducing species will produce similar nanoparticle densities and sizes for all applied power levels. In conclusion, comparisons are made to alternate models for nanoparticle formation, including charged nanoparticles and use of direct current plasmas.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Cation selectivities in zwitterion grafted nanopores: effect of zwitterion architecture

Selective separation of monovalent cations is a critical challenge in applications such as water purification and lithium recovery from salt brines. Cross-linked zwitterionic amphiphilic copolymer (ZAC-X) membranes have gained attention for their exceptional anion permselectivity, attributed to self-assembled zwitterion-lined nanodomains that interact preferentially with anions according to their hydrated radii r hyd . However, these membranes show minimal selectivity among monovalent cations, despite significant differences in their hydration structures, motivating studies on the underlying mechanisms of cation transport and selectivity in this family of materials. In this study, we conducted molecular dynamics simulations of aqueous salt solutions within zwitterion-functionalized nanopores to elucidate the influence of dipole orientation of the zwitterionic (ZI) ligands on cation diffusivities, partitioning, and permeabilities. To this end, we examined two contrasting ZI ligand architectures: Motif A (surface–cation–anion, S–ZI + –ZI − ) and Motif B (surface–anion–cation, S–ZI − –ZI + ). Our results show that in Motif A, the sulfonate groups of the ZI ligands are localized near the pore center radially, leading to strong electrostatic interactions with small bare cations (Mg 2+ and Li + ). This configuration results in high cation partitioning but low cation diffusion, maintaining solution-diffusion tradeoff typical of functionalized membranes. In contrast, Motif B show that sulfonate groups shift radially toward the mid-region of the pore. This shift, especially for small bare cations, introduces steric constraints that weaken their interactions with the sulfonate groups, thereby enhancing hydration and lowering partitioning, while still maintaining their low self-diffusivity. These findings establish zwitterion dipole orientation as a powerful design lever for tuning cation selectivity in membrane systems and offer molecular-level insights for engineering next-generation ion separation materials.

Morishita, Kazuya [Univ. of Texas, Austin, TX (Uni↗
DOE OSTI2023

Two‐Membrane Dual Non‐Aqueous/Aqueous Electrolyte Flow Cell Operation for Electrochemical Conversion of CO 2 to Methyl Formate

Recently, tandem cathodic reactions have been demonstrated in non-aqueous solvents to couple CO 2 reduction to a secondary reaction to create novel species that are not produced in aqueous CO 2 electrolysis. One reaction that can be performed with high selectivity and durability is the electrochemical conversion of CO 2 to formic acid and in-situ esterification with methanol to produce methyl formate. However, the translation to a high-performance flow electrolyzer is far from trivial, as the non-aqueous catholyte leads to reactor challenges including flooding the gas diffusion electrode. Here, in this study, a two-membrane flow electrolyzer with both anion and cation exchange membranes was used with flowing methanol catholyte and aqueous anolyte. This design prevented methanol from flooding the cathode, which was a pervasive limiting issue for electrolyzers with a single membrane. Methyl formate production at 42.9 % faradaic efficiency was achieved with pure methanol in a flow electrolyzer with stable performance beyond 80 min. However, low-water-content catholyte compositions also led to increased cell resistance and lower operating current densities. Thus, with the present ionomer materials there is a tradeoff between methyl formate selectivity and current density depending on water concentration, highlighting a need for new ionomers tailored for desirable non-aqueous solvents such as methanol.

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

Paradigms of frustration in superionic solid electrolytes

Superionic solid electrolytes have widespread use in energy devices, but the fundamental motivations for fast ion conduction are often elusive. In this Perspective, we draw upon atomistic simulations of a wide range of superionic conductors to illustrate some ways frustration can lower diffusion cation barriers in solids. Based on our studies of halides, oxides, sulfides and hydroborates and a survey of published reports, we classify three types of frustration that create competition between different local atomic preferences, thereby flattening the diffusive energy landscape. These include chemical frustration, which derives from competing factors in the anion–cation interaction; structural frustration, which arises from lattice arrangements that induce site distortion or prevent cation ordering; and dynamical frustration, which is associated with temporary fluctuations in the energy landscape due to anion reorientation or cation reconfiguration. For each class of frustration, we provide detailed simulation analyses of various materials to show how ion mobility is facilitated, resulting in stabilizing factors that are both entropic and enthalpic in origin. We propose the use of these categories as a general construct for classifying frustration in superionic conductors and discuss implications for future development of suitable descriptors and improvement strategies.

25 ENERGY STORAGE↗
DOE OSTI2022

Li 8 MnO 6 : A Novel Cathode Material with Only Anionic Redox

In Li-excess transition-metal-oxide cathode materials, anionic oxygen redox can offer high capacity and high voltages, although peroxo and superoxo species may cause oxygen loss, poor cycling performance, and capacity fading. Previous work showed that undesirable formation of peroxide and superoxide bonds is controlled to some extent by Mn substitution, and the present work uses density functional calculations to examine the reasons for this by studying the anionic redox mechanism Li 8 MnO 6 . This material is obtained by substituting Mn for Sn in Li 8 SnO 6 or for Zr in Li 8 ZrO 6 , and we also compare to previous work on those materials. The calculations predict that Li 8 MnO 6 is stable at room temperature (with a band gap of 3.19 eV as calculated HSE06 and 1.82 eV as calculated with the less reliable with PBE+U), and they elucidate the chemical and structural effects involved in the inhibition of oxygen release in this cathode. Throughout the whole delithiation process, only O 2- ions are oxidized. The directional Mn-O bonds formed from unfilled 3d orbitals effectively inhibit the formation of O-O bonds, and the layered structure is maintained even after removing 3 Li per Li 8 MnO 6 formula unit. The calculated average voltage for removal of 3 Li is 3.69 V by HSE06, and the corresponding capacity is 389 mAh/g. The high voltage of oxygen anionic redox and the high capacity result in a high energy density of 1436 Wh/kg. The Li-ion diffusion barrier for the dominant interlayer diffusion path along the c-axis is 0.57 eV by PBE+U. Finally, these results help us to understand the oxygen redox mechanism in a new lithium-rich Li 8 MnO 6 cathode material and contribute to the design of high-energy-density lithium-ion-battery cathode materials with favorable electrochemical properties based on anionic oxygen redox.

Li8MnO6↗
DOE OSTI2026

Ionic Liquid-Enhanced Interfaces to Boost Reactive C O2 Capture

The addition of ionic liquids (ILs) to a mixture containing a molecular solvent and other ionic species can induce the heterogeneous redistribution of cations and anions at the gas–liquid interface. This nonuniform redistribution of cations and anions driven by the differences in the solvophilicity of ions can improve the thermophysical and interfacial properties of such mixtures, creating a local chemical environment that is conducive to some reactions. In this work, ILs are added to a mixture of potassium hydroxide (KOH) and ethylene glycol (EG), used as a reactive absorbent and electrolyte in the migration-assisted moisture-gradient (MAMG) process for CO 2 capture. Molecular dynamics (MD) simulations are employed to probe into the effects of complex ion–ion and ion–solvent interactions and to examine the chemical composition at the gas–liquid interface. A total of 12 systems are investigated using molecular simulations to identify trends in the performance of IL additives based on the choice of cation, anion, and IL concentration. The cation effects are studied using IL additives based on 1-ethyl-3-methylimidazolium ([EMIM] + ) and 1-butyl-3-methylimidazolium ([BMIM] + ), while the impact of anions is examined using additives based on dicyanamide [DCA] − , triflate [TfO] − , bistriflimide [NTf 2 ] − , and hexafluorophosphate [PF 6 ] − anions, respectively. The influence of the IL concentration is also evaluated at molar concentrations between 1% and 4%. The simulation results indicate that the use of IL additives can affect the physical CO 2 solubility, surface tension, and the localization of CO 2 around the [OH] − ions at the gas–liquid interface. It is also evident that the choice of cations, anions, and IL concentration determines the extent to which the IL additives impact the local physicochemical properties. Physical dissolution, diffusive transport, and interaction with [OH] − are critical intermediate steps toward reactive CO 2 capture using a liquid absorbent. Hence, the improvement in one or more of these properties, aided by IL additives, is expected to improve the overall CO 2 capture performance. Experiments reaffirmed the impact of IL additives on CO 2 capture performance and the sensitivity to the choice of the cation, anion, and concentration of the IL additive.

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

Impact of Nickel Ions on the Oxygen Reduction Reaction Kinetics of Pt and on Oxygen Diffusion through Ionomer Thin Films

The effects of dissolved nickel on the oxygen reduction reaction (ORR) kinetics and oxygen transport properties of perfluorosulfonic acid (PFSA) thin films were investigated using rotating disk electrode (RDE) measurements of ORR on a PFSA-coated platinum electrode. The electrochemical characterization in 0.1 M perchloric acid (HClO 4 ) with and without added Ni 2+ quantitatively measured the impact of ionic interactions between the Ni 2+ cations and sulfonate (SO 3 - ) anions on oxygen transport through the PFSA thin film. Cyclic voltammetry (CV) curves in deaerated electrolyte showed that Ni 2+ cations diffusing through the PFSA thin film interact with the Pt surface altering the hydrogen underpotential deposition and stripping processes and decreasing ORR kinetics. The RDE limiting current results point to reduced permeability of oxygen through PFSA-Ni 2+ compared to PFSA-H + . Furthermore, the results indicate that transition metals leached from Pt alloy catalyst may be detrimental not only to the intrinsic ORR kinetics of the PEFC cathode catalyst through loss of the ORR-enhancing transition metal, but may also inhibit the diffusion of oxygen to the catalytic sites and poison the ORR.

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

Lithium superoxide encapsulated in a benzoquinone anion matrix

Lithium peroxide is the crucial storage material in lithium–air batteries. Understanding the redox properties of this salt is paramount toward improving the performance of this class of batteries. Lithium peroxide, upon exposure to p –benzoquinone (p–C 6 H 4 O 2 ) vapor, develops a deep blue color. This blue powder can be formally described as [Li 2 O 2 ] 0.3 • [LiO 2 ] 0.7 • {Li[p–C 6 H 4 O 2 ]} 0.7 , though spectroscopic characterization indicates a more nuanced structural speciation. Infrared, Raman, electron paramagnetic resonance, diffuse-reflectance ultraviolet-visible and X-ray absorption spectroscopy reveal that the lithium salt of the benzoquinone radical anion forms on the surface of the lithium peroxide, indicating the occurrence of electron and lithium ion transfer in the solid state. As a result, obligate lithium superoxide is formed and encapsulated in a shell of Li[p–C 6 H 4 O 2 ] with a core of Li 2 O 2 . Lithium superoxide has been proposed as a critical intermediate in the charge/discharge cycle of Li–air batteries, but has yet to be isolated, owing to instability. The results reported herein provide a snapshot of lithium peroxide/superoxide chemistry in the solid state with redox mediation.

25 ENERGY STORAGE↗
DOE OSTI2022

Extended kinetic lattice grand canonical Monte Carlo simulation method for transport of multicomponent ion mixtures through a model nanopore system

Abstract An extended version of the original kinetic lattice grand canonical Monte Carlo simulation method combined with mean field theory (KLGCMC/MF) ( J. Chem. Phys. 2007 , 127, 024706) is presented for the study of transport of multicomponent ion mixtures through a model nanopore. Comparison of the extended KLGCMC/MF (eKLGCMC/MF) simulation results with Poisson–Nernst–Planck (PNP) calculations is also made to confirm the validity of the extended simulation approach. Unlike the original version of KLGCMC/MF simulation method that treats only a binary ionic solution with one cation and one anion species, this extended version can deal with a system that includes ternary ion mixtures. A diffusion probability algorithm is also added to the extended version of the simulation method to describe the inhomogeneous diffusivity of ions that is often observed in the ion permeation through nanopores. Both Legendre and Chebyshev polynomials of the second kind were tested as a basis set for the basis set expansion (BSE) method with which to calculate the reaction field energy in the eKLGCMC/MF simulation. It turned out that the Legendre polynomials perform better than the Chebyshev polynomials, and as a result, the Legendre polynomials were implemented in the current version of eKLGCMC/MF simulation algorithm. The presented eKLGCMC/MF simulation method with new features finds its potential applications in nanopore systems where the correlation between ion species with the same sign of charges plays a key role such as oscillating ion currents or anomalous mole fraction effects.

Choi, Inhyeok↗