Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Dilute solution”

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

In Situ Imaging of Faujasite Surface Growth Reveals Unique Pathways of Zeolite Crystallization

Zeolite crystallization occurs by complex processes involving a variety of possible mechanisms. The sol gel media used to prepare zeolites leads to heterogeneous mixtures of solution and solid states with diverse solute species. At later stages of zeolite synthesis when growth occurs predominantly from solution, classical two-dimensional nucleation and spreading of layers on crystal surfaces via the addition of soluble species is the dominant pathway. At earlier stages, these processes occur in parallel with nonclassical pathways involving crystallization by particle attachment (CPA). The relative roles of solution- and solid-state species in zeolite crystallization have been a subject of debate. Here, in this work, we investigate the growth mechanism of a commercially relevant zeolite, faujasite (FAU). In situ atomic force microscopy (AFM) measurements reveal that supernatant solutions extracted from a conventional FAU synthesis at various times do not result in growth, indicating that FAU growth predominantly occurs from the solid state through a disorder-to-order transition of amorphous precursors. Elemental analysis shows that supernatant solutions are significantly more siliceous than both the original growth mixture and the FAU zeolite product; however, in situ AFM studies using a dilute clear solution with a lower Si/Al ratio revealed three-dimensional growth of surfaces that is distinct from layer-by-layer and CPA pathways. This unique mechanism of growth differs from those observed in studies of other zeolites. Given that relatively few zeolite frameworks have been the subject of mechanistic investigation by in situ techniques, these observations of FAU crystallization raise the question whether its growth pathway is characteristic of other zeolite structures.

36 MATERIALS SCIENCE↗

Chemostatic concentration–discharge behaviour observed in a headwater catchment underlain with discontinuous permafrost

Concentration–discharge dynamics were evaluated in a small (~ 2.25 km2) headwater catchment underlain with discontinuous permafrost on the Seward Peninsula of western Alaska. A large storm, during which 48 mm of rain fell over a 24-h period, enabled the evaluation of solute concentration–discharge response to a sizeable hydrological event, while water stable isotopes enabled an appraisal of the contributions of event water. Under normal catchment conditions, chemostatic behaviour was observed for solutes typically derived from mineral weathering (e.g. calcium, magnesium, sodium and silica). The chemostatic behaviour observed for most solutes under normal catchment conditions indicated that catchment storage and residence times are sufficiently long for many solute generating reactions to approach equilibrium. Following the storm however, most solutes exhibited dilutive and highly variable behaviour. This likely indicated the exceedance of a discharge threshold where chemostatic behaviour could no longer be maintained for most solutes. Dissolved organic carbon and silica were the only solutes monitored to exhibit chemostatic behaviour during all time periods.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Membrane-based solvent extraction for the recovery of rare earths from phosphate mining process streams

This study reports on the capture of rare earth elements (REEs) from phosphate industry process streams, including phosphoric acid (PA) sludge and phosphogypsum (PG), using a membrane solvent extraction (MSX) process. While MSX has been proven effective for a relatively concentrated feed, its effectiveness for dilute REEs solutions remains unexplored. Investigated PA-sludge and PG particles contain total REEs concentrations of ∼1100 and ∼320 ppm, respectively. Acid leaching, implemented to dissolve the REEs, significantly dilutes the REEs concentration to ∼210 ppm for PA-sludge leachate and ∼60 ppm for PG leachate. These low concentrations, compounded by the higher levels of non-REE ions and radioactive species, uranium (U) and thorium (Th), poses challenges to the MSX process. Here, we demonstrated that N,N,N′,N′-tetraoctyl-diglycolamide (TODGA) selectively binds REEs from a >3 M nitric-acid leachate while effectively rejecting U and Th. Concentrations of light REEs in strip solution were doubled compared to the feed, while heavy REEs were preferentially extracted. Furthermore, >99% purity gypsum, free of U and Th, was precipitated during the acid leaching process, aiding separation by removing significant amounts of non-REEs species (e.g., calcium) prior to the MSX process. Molecular simulations support the experimental data, suggesting preferential separation of heavy over light REEs. Based on these results, a cost-effective integrated process including pretreatment, acid leaching, MSX, and wastewater treatment is proposed for the co-recovery of REEs, phosphoric acid, gypsum, and U. This study shows MSX as a technically and economically feasible process for the recovery of REEs from low-concentration process streams, offering advantages over conventional solvent extraction.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The effect of lignin degradation products on the generation of pseudo-lignin during dilute acid pretreatment

Pseudo-lignin is an insoluble material that tends to deposit on acidic pretreated fiber surface and negatively impacts the biological conversion of biomass to valued added products. Employing a laboratory generated acidic hydrolysate solution from dilute acid pretreatment of bamboo residues, we sought to seek an improved understanding of pseudo-lignin formation from fragmented carbohydrates. In addition, we also performed experiments to investigate how lignin model compounds interact with carbohydrates en route to pseudo-lignin’s formation. The yields of pseudo-lignin generated from xylose (CXPL) were higher than those generated from glucose (CGPL). Gas chromatography-Mass spectrometry and Foloin-Ciocalteau analysis suggested that there were more aromatic compounds in hydrolyzate during CGPL’s formation. The presence of lignin phenolic model compounds had both a positive and negative effect upon quantitative yields of pseudo-lignin. Gel permeation chromatography and NMR analysis ( 31 P and 2D-HSQC NMR) of the pseudo-lignin samples generated in the presence and absence of 4-hydroxybenzoic acid (4-HA), employed as a lignin model compound, revealed drastic differences in molecular weight, changes in hydroxyl group content, and hydrocarbon bonds, suggesting significant differences in the pseudo-lignin reaction pathways when certain lignin model compounds are present. Finally, our findings demonstrate the effects of both hemicellulose carbohydrate makeup and lignin model compound impact on the formation of pseudo-lignin and its chemical properties, suggesting a potential way to minimize pseudo-lignin formation.

59 BASIC BIOLOGICAL SCIENCES↗

SANS quantification of bound water in water-soluble polymers across multiple concentration regimes

Contrast-variation small-angle neutron scattering (CV-SANS) is a widely used technique for quantifying hydration water in soft matter systems, but it is predominantly applied in the dilute regime or for systems with a well-defined structure factor. Here, CV-SANS was used to quantify the number of hydration water molecules associating with three water-soluble polymers with different critical solution temperatures and types of water–solute interactions in dilute, semidilute, and concentrated solution through the exploration of novel methods of data fitting and analysis. Multiple SANS fitting workflows with varying levels of model assumptions were evaluated and compared to give insight into SANS model selection. These fitting pathways ranged from general, model-free algorithms to more standard form and structure factor fitting. In addition, Monte Carlo bootstrapping was evaluated as a method to estimate parameter uncertainty through simulation of technical replicates. The most robust fitting workflow for dilute solutions was found to be form factor fitting without CV-SANS ( i.e. polymer in 100% D 2 O). For semidilute and concentrated solutions, while the model-free approach can be mathematically defined for CV-SANS data, the addition of a structure factor imposes physical constraints on the optimization problem, suggesting that the optimal fitting pathway should include appropriate form and structure factor models. The measured hydration numbers were consistent with the number of tightly bound water molecules associated with each monomer unit, and the concentration dependence of the hydration number was largely governed by the chemistry-specific interactions between water and polymer. Polymers with weaker water–polymer interactions ( i.e. those with fewer hydration water molecules) were found to have more bound water at higher concentrations than those with stronger water–polymer interactions due to the increase in the number of forced water–polymer contacts in the concentrated system. This SANS-based method to count hydration water molecules can be applied to polymers in any concentration regime, which will lead to improved understanding of water–polymer interactions and their impact on materials design.

36 MATERIALS SCIENCE↗

Linear Viscoelastic Properties of Putative Cyclic Polymers Synthesized by Reversible Radical Recombination Polymerization (R3P)

Linear viscoelastic properties in both melt and solution states are reported for a series of poly(3,6-dioxa-1,8-octanedithiol) (polyDODT) made by reversible radical recombination polymerization (R 3 P) under conditions designed to produce linear (LDODT), cyclic (RDODT) and linear-cyclic mixtures (LRDODT). PolyDODT is amorphous (T g < -50 °C) and highly flexible (entanglement molecular weight M e,lin ≈ 1850 g/mol for LDODT). PolyDODT’s low T g and low M e,lin enable characterization over a wide dynamic range and a wide range of dimensionless weight-average molecular weight Z w = M w /M e,lin . Measurements at temperatures from -57 °C to 100 °C provide up to 18 deades of reduced frequency, which is necessary to characterize RDODT melts with Z w from 23 to 300. The two highest molecular weight polymers in the present RDODT series have such high M w (406k and 556k g/mol) that mass spectrometry, NMR spectroscopy, and even chemical assays for chain ends are unable to rule out up to 2 mol% linear contaminant. By studying the samples in solution (using dilution to reduce Z w ) we could compare their dynamics with those of previously established high-purity polystyrene (PS) rings (limited to Z w ≤ 13.6). RDODT solutions with Z w < 15 (concentrations <5wt% for RDODT-406k and 556k) have dynamic moduli G^* that accord with LCCC-purified PS rings in terms of the frequency dependence (including the absence of a plateau), the progression of shapes of G* as a function of Z w , and the linear scaling of their zero-shear viscosity η0 with M w . The shape of G* as a function of Z w for solutions of RDODT-406k and -556k also accord with lower M w RDODT melts (which have ≤ 1.3mol% linear contaminant). Thus, measurement of the linear viscoelastic properties of appropriate concentrations of high M w (>200k g/mol) putative cyclic polymers, in which linear chains evade spectroscopic detection, may provide an alternative means (though not fully proven) of validation of sample purity. When Z w > 15 (including all seven RDODT melts and eight of their solutions), G* has a rubbery plateau. This suggests that the onset of entanglement-like behavior in rings requires 4-5-fold greater Z w than is required for linear chains. Further, the plateau moduli of RDODT samples are indistinguishable from G N o of the corresponding LDODT (melt or matched-concentration solutions). In entangled linear polymers, the observation that G N o is independent of Z w follows from limitations on lateral fluctuations due to neighboring chains becoming independent of position along a given chain. The present results for RDODT suggest that this holds for sufficiently long endless chains, too. While the RDODT have the same G N o as entangled LDODT, when Z w > 60 the terminal relaxation, if reached at all, of RDODT extends to orders of magnitude lower frequency than an entangled linear polymer of the same Z w . Consequently, the viscosity of RDODT with Z w > 60 increases with Z w much more strongly than the 3.4-power observed for entangled linear polymers. Lastly, these novel polymers, with disulfide-linked backbone and broad relaxation time distribution may prove important in relation to biodegradable elastomers and materials with exceptional low-frequency dissipation, extending at least 12 decades below the onset of the rubbery plateau.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Unexpected Conformational Behavior of Poly(poly(ethylene glycol) methacrylate)-Poly(propylene carbonate)-Poly(poly(ethylene glycol) methacrylate) (PPEGMA-PPC-PPEGMA) Amphiphilic Block Copolymers in Micellar Solution and at the Air-Water Interface

Self-assemblies (bulk micelles and Langmuir monolayers) of a new amphiphilic block copolymer, poly(poly(ethylene glycol) methacrylate)-poly(propylene carbonate)-poly(poly(ethylene glycol) methacrylate) (PPEGMA-PPC-PPEGMA), were investigated in dilute aqueous solution and at the air-water interface. An important feature that distinguishes PPEGMA-PPC-PPEGMA from conventional linear PEG-based polymer surfactants is that in the hydrophilic PPEGMA block the PEG moieties exist as side chains attached to the poly(methacrylate) (PMA) backbone. This study was performed using three PPEGMA-PPC-PPEGMA samples having an identical PPC molecular weight (5.6 kDa) and different PPEGMA molecular weights (7.3, 2.8 and 2.1 kDa on either side) (named as “G7C6G7”, “G3C5G3”, and “G2C6G2”, respectively). Cryo-TEM images revealed that G2C6G2 micelles form large clusters, whereas G7C6G7 micelles exist as unimers; the behavior of G3C6G3 is intermediate. The average core surface distance between two adjacent G2C6G2 micelles in a cluster was estimated to be only about 3.7 nm, which suggests that in the micelle corona layer only the PEG side chains (but not the PMA backbone) exist as hydrated brush chains. A calculation showed that at the small micelle separation distance observed, the van der Waals attractive forces become significant, and are thus likely responsible for the cluster formation. Steady shear rheometry measurements confirmed that G2C6G2 micelles are indeed only weakly bound in clusters, and these micelle clusters are thus prone to breakup even under mild steady shear. The micelle breakup force estimated from the shear viscosity data is in reasonable agreement with an estimate obtained from the theoretical van der Waals interaction potential. Langmuir monolayers formed at the air-water interface by PPEGMA-PPC-PPEGMA were also characterized; surface pressure-area measurements were performed. Detailed features of the surface pressure-area isotherms further supported that only the PEG side chains are hydrated in the subphase, while the PMA backbone remains unhydrated and situated on the water surface.

Lee, Jaewon↗

FY24 Task 4: Studies of Phosphate and Fluoride Solubility for Dissolution of High Phosphate Tank Waste

Knowledge gaps have been identified in phosphate solubility in almost every single- and multi-component system, and particularly for aluminate in phosphate/hydroxide, where uncertainties in predicted values of aluminum solubility are greater than 50%. Solubility data of multicomponent, aqueous electrolytes containing sodium hydroxide, sodium fluoride, sodium phosphate, sodium nitrite, sodium nitrate, and dissolved gibbsite are also sparse. For example, for solutions of (i) sodium nitrate and sodium phosphate, data are only available at 30 and 50 °C, and for (ii) dissolved gibbsite in sodium hydroxide and sodium phosphate, only two data points are available at 20 and 40 °C. There is no data available on mixtures of sodium hydroxide, sodium phosphate, sodium fluoride, and dissolved gibbsite. These knowledge gaps were identified in a technical review of waste solubility data and the impact of dilution on solution stabilities and will be addressed in this work to predict aluminum hydroxide and sodium phosphate solubility in multicomponent electrolytes upon dilution, and upon variation of temperature. Results will provide the technical basis to develop accurate models for (i) gibbsite solubility and mass transfer of aluminum between solid and liquid forms following the sluicing of sludge and saltcake with water; and (ii) further blending of these suspensions with bismuth from bismuth phosphate waste, and zirconium and uranium left from the fuel decladding. This work will be essential to developing a disposition path for retrieval solutions from bismuth phosphate wastes. This effort would also support sludge washing to further reduce phosphate concentration if that process were to be added back to the flowsheet in the future.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Adsorption Equilibrium, Kinetics, and Column Breakthrough Data for Aqueous Solutions of Binary-Acid and Ternary-Acid Mixtures of Acetic Acid, Butyric Acid, and Lactic Acid on IRN-78 Ion-Exchange Resin at Initial pH Levels of ∼3–7 and at 25–55 °C

This work is part of an effort to develop thermophysical property data and models supporting adsorptive process development for organic acid separation from a dilute aqueous solution of fermentation broth. It presents systematic experimental measurements for aqueous-phase adsorption equilibrium, kinetics, and column breakthrough for three binary-acid aqueous mixtures (acetic acid + lactic acid, butyric acid + lactic acid, and acetic acid + butyric acid) and one ternary-acid aqueous mixture (acetic acid + butyric acid + lactic acid) on Amberlite IRN-78 ion-exchange resin. The equilibrium measurements covered broad ranges of the initial acid concentration (100–400 mmol/L), initial pH (∼3–7), and temperature (25–55 °C). The equilibrium data for the binary-acid and ternary-acid aqueous mixtures indicate selective adsorption of lactic acid at initial pH levels of ∼3 and 4, while acetic acid and butyric acid are selectively adsorbed at initial pH levels of 5, 6, and 7. The subsequent kinetics and column breakthrough experiments were performed at 200 mmol/L with equimolar ratios, an initial pH of 6, and 25 °C. In conclusion, the measurements provide essential data sets for rigorous thermodynamic modeling and process simulation of adsorptive separation processes for organic acid separation from fermentation broth.

Adsorption↗

Molar Absorptivities of U(VI), U(IV), and Pu(III) in Nitric Acid Solutions of Various Concentrations Relevant to Developing Nuclear Fuel Recycling Flowsheets

Testing of a co-decontamination (CoDCon) tributyl phosphate (TBP) based solvent extraction flowsheet is being performed to optimize conditions for achieving a target U/Pu ratio in the Pu-containing product. The flowsheet is designed to directly produce a mixture of U and Pu with a U/Pu mass ratio of 7/3. The system is monitored in real time using optical spectroscopy, which allows for control of the U/Pu ratio in the aqueous U/Pu nitrate product. On-line spectrophotometric monitoring of the aqueous stream at the reductive stripping step involves determination of U(VI), U(IV) and Pu(III) concentrations and relies on molar absorptivities of these species in nitric acid at a number of wavelengths corresponding to their peaks’ maxima. The magnitudes of these molar absorptivities are a sensitive function of nitric acid concentration and have to be determined as precisely as possible using an off-line spectrophotometry under well controlled conditions. This step (called training set acquisition) is typically performed not in a high contamination area of radiological glovebox where flow-through cells and other on-line equipment are installed but in a less aggressive environment of a radiological fume hood with much better control of stock solutions quality, dilution factors, optical absorbance scale calibration, etc. This is followed by the calibration transfer between the off-line instrument and the instrument deployed during extraction and stripping runs. This paper reports values of molar absorptivities of U(VI), U(IV), and Pu(III) in nitric acid solution of 0.5 M to 4 M concentration and compares them with available technical literature data. Results of off-line spectrophotometric analysis of selected aqueous grab samples from one of CoDCon runs show satisfactory agreement with total uranium and plutonium concentrations determined in the same samples by ICP-MS.

Sinkov, Sergey I.↗

Utilizing integrated neutron diffraction and elastoplastic self-consistent crystal plasticity model to quantitatively assess the strengthening mechanism in Al–12.5Ce and Al–12.5Ce–0.4Mg alloys

An integrated in-situ neutron diffraction and elastic plastic self-consistent crystal plasticity (EPSC-CP) modeling scheme is performed on a binary Al–12Ce alloy and a ternary Al–12Ce–0.4Mg alloys. Using this scheme, the constitutive parameters, i.e. elastic constants and slip system parameters of individual phases can be calibrated which can be used in microstructure-based CP models to predict materials performance. From this study, it is shown that the elastic constants of Al 11 Ce 3 intermetallics calculated from density function theory calculation in the literature are rather accurate. When applied to the EPSC-CP model, the lattice strains of both the binary and ternary alloys are correctly predicted as compared with experiments, and large lattice strain differences between Al (100) plane and Al 11 Ce 3 (010) plane are demonstrated. The slip system parameters calibrated by the scheme shows that the addition of 0.4 wt% Mg in the alloy has little influence on the critical resolved shear stress of initial dislocation glide in the Al matrix which caused plastic yield in the material. This can be explained by the very dilute Mg solute content in the Al solid solution, causing large spacing of Al–Mg lattice misfit sites and little impact on resistance of dislocation glide at initial yield. The 0.4 wt% Mg addition, on the other hand, has a large influence on the hardening term in the slip system parameters, indicating those Al–Mg misfit sites do help dislocation accumulation during the deformation. The impact of dilute Mg addition on the Al slip system parameters is also reflected in the flow behavior of the ternary alloy: little impact on the yield stress, but a large impact on working hardening and tensile strength of the materials which is consistent with the literature.

36 MATERIALS SCIENCE↗

Selective recovery of rare earth elements with ligand-functionalized polymers in fixed-bed adsorption columns

Rare earth elements (REE) are a group of valuable metals with growing demand and broad applications. Mineral ores, the traditional sources of REE, require significant capital investment and their refinement has been a source of environmental contamination. Industrial fluids and natural REE-bearing liquids are potential alternative sources for these metals. This work investigated the performance and REE selectivity of polymer resin beads functionalized with N,N-bis(phosponomethyl)glycine (BPG) for extraction of REE from saline solutions in fixed-bed adsorption columns. Competitive batch adsorption experiments were conducted with various metals (Nd, Gd, Ho, Al, Fe, Co, Ni, Ba, Pb, Th, and U) and the BPG-functionalized resins were up to 137 times more selective for REE than aminated resins. In column experiments, the BPG-functionalized resins preferentially adsorbed heavier metals and REE were strongly retained in the functionalized column, taking 270 times longer than the amine column to reach 10% breakthrough and 128 times longer to reach 50% breakthrough. REE bound to the BPG-functionalized resins were recovered with a dilute HNO3 solution, yielding REE concentrations up to 236 times higher than the influent feedstock. This work provides new insight into the operational performance of novel functionalized adsorbents for recovery of REE from saline fluids.

42 ENGINEERING↗

Structural Characterization of Disulfide-Linked p53-Derived Peptide Dimers

Disulfide bonds provide a convenient method for chemoselective alteration of peptide and protein structure and function. We previously reported that mild oxidation of a p53-derived bisthiol peptide (CTFANLWRLLAQNC) under dilute non-denaturing conditions led to unexpected disulfide-linked dimers as the exclusive product. The dimers were antiparallel, significantly α-helical, resistant to protease degradation, and easily reduced back to the original bisthiol peptide. Here, in this work, we examine the intrinsic factors influencing peptide dimerization using a combination of amino acid substitution, circular dichroism (CD) spectroscopy, and X-ray crystallography. CD analysis of peptide variants suggests critical roles for Leu6 and Leu10 in the formation of stable disulfide-linked dimers. The 1.0 Å resolution crystal structure of the peptide dimer supports these data, revealing a leucine-rich LxxLL dimer interface with canonical knobs-into-holes packing. Two levels of higher-order oligomerization are also observed in the crystal: an antiparallel “dimer of dimers” mediated by Phe3 and Trp7 residues in the asymmetric unit and a tetramer of dimers mediated by Trp7 and Leu10. In CD spectra of Trp-containing peptide variants, minima at 227 nm provide evidence for the dimer of dimers in dilute aqueous solution. Importantly, and in contrast to the original dimer model, the canonical leucine-rich core and robust dimerization of most peptide variants suggests a tunable molecular architecture to target various proteins and evaluate how folding and oligomerization impact various properties, such as cell permeability.

36 MATERIALS SCIENCE↗

Variation of the Passive Film on Compositionally Concentrated Dual-Phase Al 0.3 Cr 0.5 Fe 2 Mn 0.25 Mo 0.15 Ni 1.5 Ti 0.3 and Implications for Corrosion

The passive film on a dual-phase Al 0.3 Cr 0.5 Fe 2 Mn 0.25 Mo 0.15 Ni 1.5 Ti 0.3 FCC + Heusler (L2 1 ) compositionally concentrated alloy formed during extended exposure to an applied potential in the passive range in dilute chloride solution was characterized. Each phase, with its own distinct composition of passivating elements, formed unique passive films separated by a heterophase interface. High-resolution, surface sensitive characterization enabled chemical analysis of the passive film formed over individual phases. The film formed over the L2 1 phase had a higher concentration of Al, Ni, and Ti, while the film formed over FCC phase was of similar thickness but contained comparatively higher Cr, Fe, and Mo concentrations, consistent with the differences in bulk microstructure composition. The passive film was continuous across phase boundaries and the distribution of passivating elements (Al, Cr, and Ti) indicated both phases were independently passivated. Spatially resolved analysis of the surface chemistry of the dual-phase CCA revealed that the cation with the highest composition in passive film formed on the FCC phase was Cr (52.4 at. pct) and for the L2 1 phase was Ti (53.1 at. pct) despite the bulk concentration of each element being below 20 at. pct in their respective phases. Al, Cr, and Ti were enriched in both phases within the passive film relative to their respective bulk compositions. In parallel studies, single-phase alloys with compositions representative of the FCC and L2 1 phases were synthesized to evaluate the corrosion behavior of each phase in isolation. The corrosion behavior of the dual-phase alloy showed passivity evidenced by a pitting potential of 0.615 V SCE in 0.01 M NaCl. The pitting potential and other electrochemical parameters suggested a combination of behaviors of both single-phase samples, suggesting that the global corrosion behavior may be represented by a composite theory applied to phases, their area fractions, and interphase length. However, the interphase in the dual-phase CCA was a local corrosion initiation site and may limit localized corrosion protectiveness. The alloy design implications for optimization of second phase structure and morphology are discussed.

36 MATERIALS SCIENCE↗

Effects of cluster expansion on the locations of phase transition boundary as a first step to quantify uncertainty in first principles statistical mechanics framework

Predicting phase diagrams from first principle calculations eliminates the need of tedious experimental trials and errors. Fully automating first principle phase diagram calculations without any sort of human intervention has been a long daunting task and troubling scientific communities for decades. This grand problem remains not fully resolved, largely due to the vastly high-dimensional parameter space associated with density functional theory, cluster expansion, lattice Monte Carlo, and the substantial uncertainty propagating through a set of complex simulations. As a first step to tackle this grand problem, we reported a first demonstration of how sensitive phase boundary locations can be to various cluster expansion fittings and input DFT training data. To the best knowledge of the authors, this study reported the first ever attempt to quantify uncertainty in first principles statistical mechanics framework. In addition, a semi-automated phase transition detection algorithm has been devised in this paper to deal with the associated statistical errors and uncertainties from Monte Carlo method and its predecessor DFT calculations and cluster expansions. This algorithm has been applied in a classical cluster expansion Mg-Cd binary alloy system to detect phase transitions at various locations of phase diagram using different cluster expansions to demonstrate its predictive power and quantify uncertainties. The results suggested that using Chebyshev basis function shifted the transition locations toward the dilute solid solution phase and expanded the phase stability of concentrated ordered phases to wider composition range. The addition of perturbed defect configurations into training data set lowered the order-disorder transition temperature to be away from true transition temperature, suggesting the transition nature is indeed configurational disorder dominated rather than defect assisted. Finally, we found that the weighting has negligible effect on the transition locations, except for the case of using Chebyshev basis function to fit all non-weighted configurations that can be susceptible to Monte Carlo sampling hysteresis.

36 MATERIALS SCIENCE↗

Modeling dissociation of ionic liquids with electrolyte NRTL model

Dissociation of ionic liquids (ILs) is critical to electrochemical process performance by affecting ionic conductivity. Prior studies of thermodynamic modeling either assumed no dissociation or complete dissociation that were not consistent with experimental observations. This work presented an eNRTL model for three aqueous imidazolium-based IL systems: [Emim][EtSO 4 ], [Emim][TFA], and [Emim][TfO]. The model accurately correlates the IL dissociation and VLE data in the entire concentration range from pure IL to infinite dilution aqueous solution. The system non-ideality is dominated by the dissociation chemistry. Dissociation behavior is driven by strong short-range interaction of the molecular IL-water pair relative to the dissociated IL-water pair, consistent with prior studies indicating that molecular IL and water can form an energetically favorable complex. The interaction energy parameters calculated in this work qualitatively agrees with the hydrogen bonding ability of the anions previously reported. This research not only provides a model that well describes the IL dissociation but also demonstrates the physical significance of eNRTL model.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Influence of Long-Range Forces on the Transition States and Dynamics of NaCl Ion-Pair Dissociation in Water

We study NaCl ion-pair dissociation in a dilute aqueous solution using computer simulations both for the full system with long-range Coulomb interactions and for a well-chosen reference system with short-range intermolecular interactions. Analyzing results using concepts from Local Molecular Field (LMF) theory and the recently proposed AI-based analysis tool “State predictive information bottleneck” (SPIB), we show that the system with short-range interactions can accurately reproduce the transition rate for the dissociation process, the dynamics for moving between the underlying metastable states, and the transition state ensemble. Contributions from long-range interactions can be largely neglected for these processes because long-range forces from the direct interionic Coulomb interactions are almost completely canceled (>90%) by those from solvent interactions over the length scale where the transition takes place. Thus, for this important monovalent ion-pair system, short-range forces alone are able to capture detailed consequences of the collective solvent motion, allowing the use of physically suggestive and computationally efficient short-range models for the dissociation event. Here, we believe that the framework here should be applicable to disentangling mechanisms for more complex processes such as multivalent ion disassociation, where previous work has suggested that long-range contributions may be more important.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Lineshape Distortions in Internal Reflection Two-Dimensional Infrared Spectroscopy: Tuning across the Critical Angle

Reflection mode two-dimensional infrared spectroscopy (R-2DIR) has recently emerged as a tool that expands the utility of ultrafast IR spectroscopy toward a broader class of materials. The impact of experimental configurations on the potential distortions of the transient reflectance (TR) spectra has not been fully explored, particularly in the vicinity of the critical angle θ c and through the crossover from total internal reflection to partial reflection. Here we study the impact on the spectral lineshape of a dilute bulk solution as θ c is varied across the incident angle by tuning the refractive index of the solvent. We demonstrate the significance of several distortions, including the appearance of phase twisted lineshapes and apparent changes in the spectral inhomogeneity, and show how these distortions impact the interpretation of the TR and R-2DIR spectroscopies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗