Search NASA⌕ Search

SEARCH · Search NASA

Results for “Binary mixtures”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Statistics of Wide Pre-main-sequence Binaries in the Orion OB1 Association

Statistics of low-mass pre-main-sequence binaries in the Orion OB1 association with separations ranging from 0."6 to 20" (220 to 7400 au at 370 pc) are studied using images from the VISTA Orion mini survey and astrometry from Gaia. The input sample based on the CVSO catalog contains 1137 stars of K and M spectral types (masses between 0.3 and 0.9 M ⊙ ), 1021 of which are considered to be association members. There are 135 physical binary companions to these stars with mass ratios above ~0.13. The average companion fraction is 0.09 ± 0.01 over 1.2 decades in separation, slightly less than, but still consistent with, the field. We found a difference between the Ori OB1a and OB1b groups, the latter being richer in binaries by a factor of 1.6 ± 0.3. No overall dependence of the wide-binary frequency on the observed underlying stellar density is found, although in the Ori OB1a off-cloud population, these binaries seem to avoid dense clusters. The multiplicity rates in Ori OB1 and in sparse regions like Taurus differ significantly, hinting that binaries in the field may originate from a mixture of diverse populations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Mixed Tin-Titanium Oxides by Atomic Layer Deposition on Planar Substrates: Physical and Electronic Structure

Ternary atomic layer deposition of Sn-doped TiO x (Sn:TiO x ) and Ti-doped SnO x (Ti:SnO x ) was performed using tetrakisdimethylamidotitanium (IV) and tin(IV) with water at 150 °C on (1 0 0) Si and quartz substrates. The physical structure of the films was investigated using scanning electron, atomic force, and Kelvin probe force microscopies as well as spectroscopic ellipsometry, Raman spectroscopy, and X-ray diffraction. All as-deposited films are amorphous and possessed flat and conformal topographies with low surface roughness (<1 nm). Annealed Sn:TiO x and Ti:SnO x crystallize into the tetragonal anatase TiO 2 and cassiterite SnO 2 structures, respectively, without large variations to the film topography. The chemical composition of the as-deposited and annealed films was investigated through X-ray photoelectron spectroscopy. While dopants are homogenously mixed, the ternary ALD to form Sn:TiO x and Ti:SnO x surprisingly does not follow the rule of mixtures, despite similarities in both the growth rates of binary SnO x and TiO x and the precursors used. This finding is also supported by spectroscopic ellipsometry. Here, we postulate the deviation from the rule of mixtures is due to non-ideal surface chemistry of the binary ALD of SnO x . The optical bandgap of all films was quantified by the Tauc analysis and the valence band position was investigated by valence band X-ray photoelectron spectroscopy. Intra-bandgap states above the valence band maximum were observed for as-deposited Ti:SnO x and Sn:TiO x and ascribed to 5 s states from Sn 2+ ions. The presence of these intra-bandgap states induces enhanced visible absorption and bathochromic shifts to the apparent bandgap for all as-deposited films. Annealing the films removes the intra-bandgap states and increases the optical bandgap. Simplified band diagrams created using the data available demonstrated that the conduction and valence band positions did not vary between the positions of pure TiO 2 or SnO 2 with monotonic increases to the dopant incorporation.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Inferring topological transitions in pattern-forming processes with self-supervised learning

Abstract The identification of transitions in pattern-forming processes are critical to understand and fabricate microstructurally precise materials in many application domains. While supervised methods can be useful to identify transition regimes, they need labels, which require prior knowledge of order parameters or relevant microstructures describing these transitions. Instead, we develop a self-supervised, neural-network-based approach that does not require predefined labels about microstructure classes to predict process parameters from observed microstructures. We show that assessing the difficulty of solving this inverse problem can be used to uncover microstructural transitions. We demonstrate our approach by automatically discovering microstructural transitions in two distinct pattern-forming processes: the spinodal decomposition of a two-phase mixture and the formation of binary-alloy microstructures during physical vapor deposition of thin films. This approach opens a path forward for discovering unseen or hard-to-discern transitions and ultimately controlling complex pattern-forming processes.

Abram, Marcin↗
DOE OSTI2021

Kinetically Controlled Self-Assembly of Binary Polymer-Grafted Nanocrystals into Ordered Superstructures via Solvent Vapor Annealing

Polymer-inorganic nanocomposites based on polymer-grafted nanocrystals (PGNCs) are enabling technologically relevant applications owing to their unique physical, chemical, and mechanical properties. While diverse PGNC superstructures have been realized through evaporation-driven self-assembly, this approach presents multifaceted challenges in experimentally probing and controlling assembly kinetics. Here, we report a kinetically controlled assembly of binary superstructures from a homogeneous disordered PGNC mixture utilizing solvent vapor annealing (SVA). Using a NaZn 13 -type superstructure as a model system, we demonstrate that varying the solvent vapor pressure during SVA allows for exquisite control of the rate and extent of PGNC assembly, providing access to nearly complete kinetic pathways of binary PGNC crystallization. Characterization of kinetically arrested intermediates reveals that assembly follows a multistep crystallization pathway involving spinodal-like preordering of PGNCs prior to NaZn 13 nucleation. Our work opens up new avenues for the synthesis of multicomponent PGNC superstructures exhibiting multifunctionalities and emergent properties through a thorough understanding of kinetic pathways.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2021

Direct numerical simulation of single-species and binary-species boundary layers at high pressure

Direct numerical simulations of single-species and binary-species temporal boundary layers at high pressure are performed. The main objective is to investigate flow physics in a binary-species boundary layer at high pressure, with special attention to mass diffusion. The working fluids are nitrogen in the single-species cases, and a mixture of nitrogen and methane in the binary-species cases. An investigation of mean profiles shows that velocity and temperature profiles have steep gradients near the wall, whereas the mass fraction profiles do not have a steep gradient. This result indicates that the similarity law between velocity, temperature and mass fraction does not hold true in the binary-species boundary layer. The comparison of turbulent fluctuations shows that the qualitative characteristics of velocity and temperature are similar to each other. In contrast, profiles of mass fraction fluctuations are largely different from those of velocity and temperature. This result indicates that mass diffusion is not similar to momentum and thermal diffusions in binary-species boundary layers.

Toki, Takahiko↗
DOE OSTI2022

Factorization of Binary Matrices: Rank Relations, Uniqueness and Model Selection of Boolean Decomposition

The application of binary matrices are numerous. Representing a matrix as a mixture of a small collection of latent vectors via low-rank decomposition is often seen as an advantageous method to interpret and analyze data. In this work, we examine the factorizations of binary matrices using standard arithmetic (real and nonnegative) and logical operations (Boolean and $\mathbb{Z}$ 2 ). We examine the relationships between the different ranks, and discuss when factorization is unique. In particular, we characterize when a Boolean factorization X = W$\land$H has a unique W, a unique H (for a fixed W), and when both W and H are unique, given a rank constraint. We introduce a method for robust Boolean model selection, called BMFk, and show on numerical examples that BMFk not only accurately determines the correct number of Boolean latent features but reconstruct the pre-determined factors accurately.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Simplified Universal Equations for Ionic Conductivity and Transference Number

Nernst-Einstein equation can provide a reasonable estimate of the ionic conductivity of dilute solutions. For concentrated solutions, alternate methods such as Green–Kubo relations and Einstein relations are more suitable to account for ion-ion interactions. Such computations can be expensive for multicomponent systems. Simplified mathematical expressions like the Nernst-Einstein equation do not exist for concentrated multicomponent mixtures. Newman’s treatment of multicomponent concentrated solutions yields a conductivity relation in terms of species concentration and Onsager phenomenological coefficients. However, the estimation of these phenomenological coefficients is not straightforward. Here, mathematical formulations that relate the phenomenological coefficients with the friction coefficients are developed, leading to simplified, ready-to-use expressions of conductivity and transference numbers that can be used for a wide range of ionic mixtures. This approach involves spectral decomposition of the matrix of Onsager phenomenological coefficients. The general analytical expressions for conductivity and transference number are simplified for binary electrolytes, and numerical solutions are provided for ternary and quaternary mixtures with ion dissociation.

Electrochemistry↗
DOE OSTI2025

SpecDis: Value Added Distance Catalog for 4 Million Stars from DESI Year-1 Data

We present the SpecDis value-added stellar distance catalog accompanying DESI Data Release 1. SpecDis trains a feed-forward neural network (NN) with Gaia parallaxes and gets the distance estimates. To build up an unbiased training sample, we do not apply selections on parallax error or signal-to-noise (S/N) of the stellar spectra, and instead, we incorporate parallax error into the loss function. Moreover, we employ principal component analysis to reduce the noise and dimensionality of stellar spectra. Validated by independent external samples of member stars with precise distances from globular clusters, dwarf galaxies, stellar streams, combined with blue horizontal branch stars, we demonstrate that our distance measurements show no significant bias up to 100 kpc, and are much more precise than Gaia parallax beyond 7 kpc. The median distance uncertainties are 23%, 19%, 11%, and 7% for S/N < 20, 20 ≤ S/N < 60, 60 ≤ S/N < 100, and S/N ≥ 100. Selecting stars with ${\mathrm{log}}\,g\lt 3.8$ and distance uncertainties smaller than 25%, we have more than 74,000 giant candidates within 50 kpc of the Galactic center and 1500 candidates beyond this distance. Additionally, we develop a Gaussian mixture model to identify unresolvable equal-mass binaries by modeling the discrepancy between the NN-predicted and the geometric absolute magnitudes from Gaia parallaxes and identify 120,000 equal-mass binary candidates. Our final catalog provides distances and distance uncertainties for >4 million stars, offering a valuable resource for Galactic astronomy.

astronomy data analysis↗
DOE OSTI2024

Generic Behavior of Ultrastability and Anisotropic Molecular Packing in Codeposited Organic Semiconductor Glass Mixtures

Vapor-deposited glass mixtures of organic semiconductors commonly serve as active layers in organic electronic devices, whose lifetime and performance are strongly influenced by the stability and structure of these mixed glasses. Here, we study the stability and anisotropic molecular packing of six co-deposited organic semiconductor glass mixtures with 50:50 weight ratio, by differential scanning calorimetry and spectroscopic ellipsometry. We also find that all six binary systems exhibit high kinetic stability and significantly reduced enthalpy relative to the corresponding liquid-cooled glassy mixtures (ultrastable behavior), even for systems where the glass transition temperatures of the components differ by more than 90 K. Furthermore, we demonstrate that the birefringence of a co-deposited glass mixture, a measure of its anisotropic packing, can be predicted from the birefringence of glasses of the two pure components. These results for stability and structure are expected to be applicable to other co-deposited organic semiconductor glass mixtures, so long as the two components mix well in the glass and individually can form ultrastable glasses. Therefore, our findings are significant for designing novel electronic devices with enhanced device lifetime and increased operational efficiency.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Undercooling minimization in ultrasound coupled DTA measurements of molten salts

Molten salts are ionic liquids that are used for the electrolytic pyroprocessing of metals and as heat transfer fluids in very high temperature processes. Recently, halide-based molten salt reactors (MSR) have gained momentum for high density and environmentally responsible electricity generation. The function of these reactors and their fuel cycle depend on a knowledge of the halide salt’s thermodynamic properties [1]. Therefore, high accuracy phase equilibria of MSR relevant base halide and actinide containing salts are needed. Modern thermal analysis of the phase transitions of halide salts is usually done during heating at relatively high scan rates with commercial devices. Normally such measurements are adequate for pure or pseudo-binary salts. However, as the number of components in the mixture increases, accurately resolving the liquidus becomes increasingly difficult. Reversing the scanning mode greatly increases the sensitivity of phase transition measurements but can decrease accuracy due to undercooling [2] — a common occurrence in molten halide systems [3]. This work presents the development of a differential thermal analysis (DTA) cell intended for use in radiological gloveboxes. Non-contact ultrasonic agitation of the halide salt is implemented to limit kinetic limitations on crystallization during cooling to minimize undercooling [4]. Measurements on halide salts are also presented to elucidate the effect of non-contact mixing on undercooling. Reference [1] S. Boyd and C. Taylor, “3 - Chemical fundamentals and applications of molten salts,” in Molten Salt Reactors and Thorium Energy, T. J. Dolan, Ed., Woodhead Publishing, 2017, pp. 29–91. doi: 10.1016/B978-0-08-101126-3.00003-8. [2] K. Nitsch, A. Cihlár, and M. Rodová, “Molten state and supercooling of lead halides,” J. Cryst. Growth, vol. 264, no. 1, pp. 492–498, Mar. 2004, doi: 10.1016/j.jcrysgro.2004.01.011. [3] L. Rycerz, “Practical remarks concerning phase diagrams determination on the basis of differential scanning calorimetry measurements,” J. Therm. Anal. Calorim., vol. 113, no. 1, pp. 231–238, Jul. 2013, doi: 10.1007/s10973-013-3097-0. [4] Md. H. Zahir, S. A. Mohamed, R. Saidur, and F. A. Al-Sulaiman, “Supercooling of phase-change materials and the techniques used to mitigate the phenomenon,” Appl. Energy, vol. 240, pp. 793–817, Apr. 2019, doi: 10.1016/j.apenergy.2019.02.045.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Bridging confined phase behavior of CH 4 -CO 2 binary systems across scales

Phase behavior of confined fluids may deviate significantly from that of the bulk fluid due to the fluid-wall interactions being a significant portion of all intermolecular interactions under confinement. Despite recent advancements in understanding confined phase behavior of pure fluids, confined phase behavior of mixtures remains an understudied topic. In this work, we examine the confined phase behavior of a CH 4 -CO 2 binary system by combining Monte Carlo (MC) simulations, a cubic equation of state (EoS), and the lattice Boltzmann method (LBM). First, the effects of confinement on density and phase distribution in nano-size pores are established using Gibbs Ensemble MC calculations, which produce precise results of liquid and vapor confined pressures and account for the modification of the phase change location. By comparing the phase envelopes of bulk and confined mixtures at a fixed temperature, here it is observed that the phase envelopes shrink with reductions in pore size. Based on this observation, we extend a modified Peng-Robinson EoS, which was originally developed for pure fluids under confinement, to mixtures via van-der-Waals-type mixing rules and by accounting for shifts in the critical properties of confined CH 4 -CO 2 . The resulting phase envelopes are in good agreement with the MC data. In addition, a local density model is used in combination with the confined EoS to calculate adsorption isotherms of CH 4 -CO 2 mixtures and to characterize the behavior of confined matter in nanopores. Finally, we incorporate this EoS in a multicomponent multiphase LBM that uses a pseudopotential model to represent intermolecular forces. This workflow utilizes multiscale simulation techniques to bridge the behavior of multicomponent systems across scales and to shed light on the confined phase behavior of CH 4 -CO 2 binary systems.

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

Nonane and Hexanol Adsorption in the Lamellar Phase of a Nonionic Surfactant: Molecular Simulations and Comparison to Ideal Adsorbed Solution Theory

Adsorption of n-nonane/1-hexanol (C9/C6OH) mixtures into the lamellar phase formed by a 50/50 w/w triethylene glycol mono-n-decyl ether (C10E3)/water system was studied using configurational-bias Monte Carlo simulations in the osmotic Gibbs ensemble. The interactions were described by the Shinoda–Devane–Klein coarse-grained force field. Prior simulations probing single-component adsorption indicated that C9 molecules preferentially load near the center of the bilayer increasing the bilayer thickness, whereas C6OH molecules are more likely to be found near the interface of the polar and non-polar moieties swelling the bilayer in the lateral dimension. Here, we extend this work to binary C9/C6OH adsorption to probe whether the difference in the spatial preferences may lead to a synergistic effect and enhanced loadings for the mixture. Comparing loading trends and the thermodynamics of binary adsorption to unary adsorption reveals that C9–C9 interactions lead to the largest enhancement, whereas C9–C6OH and C6OH–C6OH interactions are less favorable for this bilayer system. As a result, ideal adsorbed solution theory yields satisfactory predictions of the binary loading.

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

Sorption and oxidation of Co(II) at the surface of birnessite: Impacts of aqueous Mn(II)

The sorption of aqueous Mn(II) (1 mM) and Co(II) (50 and 200 μM) onto hexagonal birnessite (0.1 g L –1 ) was studied under anoxic conditions at pH 6.5 and 7.5 in binary and ternary experiments using batch kinetic experiments and XRD, ATR-FTIR, and Co K-edge EXAFS analyses. In the binary systems, sorption of Co(II) was accompanied by partial oxidation to Co(III) yielding a mixture of corner-sharing Co(II) and edge-sharing Co(III) complexes at both pH values, while Mn(II)-birnessite interaction resulted in coordination of Mn(II/III) at layer vacancy sites at pH 6.5, and led to reductive transformation of birnessite into feitknechtite at pH 7.5. In the ternary systems, strong mutual interferences between Co(II) and Mn(II) co-sorbates reduced the rate and extent of sorption relative to the binary experiments. The introduction of Mn(II) aq into the high-Co systems (200 μM) halted the slow sorption of Co(II) that was attributed to the incorporation of Co(III) into layer vacancies, while Co(II) aq prevented Mn(II)-driven transformation of birnessite into feitknechtite at pH 7.5. In the low-Co system (50 μM), reductive transformation of birnessite by Mn(II) aq at pH 7.5 produced Co(II)-substituted feitknechtite, a conversion that was accompanied by reduction of sorbed Co(III) to Co(II) which was partially released to solution. The strong effects of co-sorption are attributed to the similarity in sorption mechanisms of Co(II) and Mn(II), which both sorb as a mixture of di- and trivalent species. Finally, the results of this work demonstrate that aqueous Mn(II) may significantly affect the reactivity of phyllomanganate sorbents towards dissolved Co(II), and therefore impact the speciation and solubility of this trace metal in anoxic and suboxic geochemical environments.

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

Effect of MgSO 4 addition on alkali sulphates induced hot corrosion of a β-NiAl coating

The effect of MgSO 4 addition on alkali sulphates induced hot corrosion of a β-NiAl coating was studied by performing a series of tests with deposits of Na 2 SO 4 , Na 2 SO 4 -20 %K 2 SO 4 and Na 2 SO 4 -12 %K 2 SO 4 -35 %MgSO 4 (mol.% for all) at 700 °C in air with addition of SO 2 . For elucidating the observed corrosion phenomena, various salt-oxide/NiAl powder mixtures were exposed to the same hot corrosion test conditions. Two pseudo-binary phase diagrams of the Na, K, Mg-sulphate systems were calculated using an in-house developed thermodynamic database. Addition of MgSO 4 to Na 2 SO 4 -20 %K 2 SO 4 deposit salt inhibited the hot corrosion attack during 24 h exposure, which is related to consumption of aggressive K 2 SO 4 by the formation of solid K 2 Mg 2 (SO 4 ) 3 . However, extending the exposure with Na 2 SO 4 -12 %K 2 SO 4 -35 %MgSO 4 deposit to 100 h resulted in severe corrosion attack. The latter observation is explained by a reaction between K 2 Mg 2 (SO 4 ) 3 and NiAl causing formation of a Mg, Ni and Al containing spinel accompanied by release of the aggressive K 2 SO 4 into the liquid phase.

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

Intergranular oxidation and oxygen transport in Ni-20Cr during exposure to model dry and wet oxidation environments

Selective oxidation is a critical component of the proposed internal oxidation mechanism of stress corrosion cracking (SCC) for Ni-based alloys in primary water reactor (PWR) operation. To isolate and study this process, the binary Ni-20Cr was exposed to pressurized water and a CO/CO2 gas mixture. In both environments, the oxygen potential was kept under the Ni/NiO line to stimulate penetrative intergranular attack. A combination of electron microscopy and atom probe tomography was used to study the nanometer-scale details of the passivation and penetrative intergranular oxidation processes at high-energy grain boundaries. Oxygen transport towards the terminating oxidation front is elucidated with dedicated usage of oxygen isotopic tracer exchange experiments to support classical theories of internal oxidation, revealing preferred oxygen transport paths were visualized with sub-nanometer resolution.

Weiser, Martin↗
DOE OSTI2023

Property Measurements of NaCl-UCl 3 and NaCl-KCl-UCl 3 Molten Salts (Rev.1)

Thermochemical and thermophysical property values of several salt compositions of interest are needed by molten salt reactor (MSR) developers to design, license, and operate their reactors. Thermochemical and thermophysical properties being measured at Argonne include thermal transitions, phase behavior, heat capacity, density, volumetric thermal expansion of the liquid phase, thermal diffusivity, thermal conductivity, and viscosity. Several properties of eutectic compositions in the ternary NaCl-KCl-UCl 3 and binary NaCl-UCl 3 systems that may be used by MSR developers as fuel bearing salts are being measured. A 65.8 mol % NaCl–34.2 mol % UCl 3 mixture and a near-eutectic mixture of 50.9 mol % NaCl–24.4 mol % KCl–24.7 mol % UCl 3 were synthesized and the thermochemical properties of the mixtures were measured by using differential scanning calorimetry (DSC). The measured transition temperatures were compared to transition temperatures predicted using two models. A thermodynamic model of the binary NaCl-UCl 3 system was constructed using data in the Molten Salt Thermal Properties Database–Thermochemical Version 2.0 (MSTDB-TC V2.0). A ternary NaClKCl-UCl 3 system model constructed at Argonne and was described in a previous report. These comparisons can be used to validate the models. Thermophysical property values of molten salts are needed to model how salt retains and transfers heat in an MSR system. These property values are essential to the entire MSR design because molten salt is used as both the fuel and the coolant material in a salt fueled reactor. Heat capacity of the synthesized NaCl-UCl 3 and NaCl-KCl-UCl 3 salt mixtures was measured by using DSC and thermal diffusivity was measured by using laser flash analysis (LFA) at temperatures spanning the typical operating range of an MSR.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Local Structural Disorder in Metavanadates MV 2 O 6 (M = Zn and Cu) Synthesized by the Deep Eutectic Solvent Route: Photoactive Oxides with Oxygen Vacancies

Metavanadates MV 2 O 6–δ (M = Zn and Cu) are synthesized by using a deep eutectic solvent (DES), a mixture of hydrogen bond donor and acceptor, as a reaction medium. Dissolution of stable binary metal oxides in a DES followed by a heat treatment yields phase-pure vanadates. According to in situ powder X-ray diffraction, ternary phases (α-Zn 2 V 2 O 7 and β-Cu 2–x V 2 O 7 , x ~ 0.27) are intermediates in the reaction pathway taken. Identifying a polymorphic phase transformation temperature for CuV 2 O 6 as well as the narrow temperature range between formation and decomposition for both metavanadates allows for tackling the challenge of the synthesis of these materials. The oxygen vacancy introduced by the DES route is accompanied by the formation of reduced V 4+ and Cu + in the oxide matrix, based on X-ray photoelectron spectroscopy. These oxygen vacancies modify the vibrational modes in the corresponding Raman spectra and are also responsible for broad optical absorptions in the 1.8–1.1 eV range. The optical band gaps of the materials are found at 1.8 eV (CuV 2 O 6 ) and 2.2 eV (ZnV 2 O 6 ), approximately 0.1–0.3 eV below the values reported in the literature. The reduced band gaps and sub-band gap photon absorption are key features of the oxygen-deficient metavanadates. Surface photovoltage spectroscopy reveals that all the synthesized vanadates are n-type materials with electrons as the majority charge carriers, and photoelectrochemical measurements confirm photoanodic currents for methanol oxidation. Furthermore, these results provide insight into the synthesis and structure–property relationship of the metavanadates for their potential use as photoanodes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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↗