Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Solvation”

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 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Localized High-Concentration Electrolytes Get More Localized

Liquid electrolytes in batteries are typically treated as macroscopically homogeneous ionic transport media, despite having complex chemical composition and atomistic solvation structures, thus leaving a knowledge gap of microstructural characteristics. A promising electrolyte system design has recently been developed with the title of a localized high-concentration electrolyte (LHCE). These electrolytes are a summation of three components – a charge-carrying salt, a solvent that solvates the salt, and a diluent that is sparingly interactive with the salt in terms of ion solvation. This creates the perception of a high-concentration electrolyte (HCE) with high anion-cation interactions that supports a salt-derived, stable solid-electrolyte interphase, while an LHCE improves bulk electrolyte properties, such as ionic conductivity and viscosity, when compared to an HCE. Much trial by error design of this LHCE concept has been supported by simulations to confirm ionic interactions, cycling data to confirm cell stability, and post-mortem analysis to provide evidence of mechanisms. However, a systematic design to extend the capabilities of LHCEs is still to be desired. Here, we analyzed a specific LHCE mixture to dictate the best viable mixture that extends the localization of solvation structures while maintaining viable cell-level capabilities. We seek to observe and understand solvent-to-salt interactions that impact electrolyte performance. These interactions range from single-ion solvation, multidentate interactions (i.e., solvation-driven clustering), to more widespread network-level interactions. With the support of both experimental and computational analysis of different electrolyte features, we have begun a methodology of formulating the most ideal LHCE mixtures for improving the cyclability of high-capacity battery chemistries.

25 ENERGY STORAGE↗
DOE OSTI2016

Rate theory of solvent exchange and kinetics of Li{sup +} − BF{sub 4}{sup −}/PF{sub 6}{sup −} ion pairs in acetonitrile

In this paper, we describe our efforts to apply rate theories in studies of solvent exchange around Li{sup +} and the kinetics of ion pairings in lithium-ion batteries (LIBs). We report one of the first computer simulations of the exchange dynamics around solvated Li{sup +} in acetonitrile (ACN), which is a common solvent used in LIBs. We also provide details of the ion-pairing kinetics of Li{sup +}-[BF{sub 4}] and Li{sup +}-[PF{sub 6}] in ACN. Using our polarizable force-field models and employing classical rate theories of chemical reactions, we examine the ACN exchange process between the first and second solvation shells around Li{sup +}. We calculate exchange rates using transition state theory and weighted them with the transmission coefficients determined by the reactive flux, Impey, Madden, and McDonald approaches, and Grote-Hynes theory. We found the relaxation times changed from 180 ps to 4600 ps and from 30 ps to 280 ps for Li{sup +}-[BF{sub 4}] and Li{sup +}-[PF{sub 6}] ion pairs, respectively. These results confirm that the solvent response to the kinetics of ion pairing is significant. Our results also show that, in addition to affecting the free energy of solvation into ACN, the anion type also should significantly influence the kinetics of ion pairing. These results will increase our understanding of the thermodynamic and kinetic properties of LIB systems.

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

Structure and Dynamics of NaCl Ion Pairing in Solutions of Water and Methanol

Ion pairing can have profound effects upon the ionic strength of electrolyte solutions but is poorly understood in solutions containing more than one solvent. In this paper, a combined density functional theory and molecular dynamics approach is used to examine the effect of both methanol concentration and interionic distance upon the structure and dynamics within successive solvation shells of Na + and Cl – in water/methanol binary solutions. The structure and dynamics of the first and second solvation shells were studied along a reaction coordinate associated with ion pair formation using potential of mean force simulations. The lifetimes of the solvent–solvent hydrogen bonds become perturbed when the second solvation shells of the ions begin to interact. In contrast, the structural properties within the first and second solvation shells of the ions were found to be largely independent of both methanol concentration and interionic distance until a contact ion pair is formed. Thus, as the ions are brought together, the effect of the opposing ion manifests itself in the solvation dynamics before any structural changes are observed. Lastly, as anticipated based upon the decreased dielectric constant of the binary solution, ion pair formation becomes energetically more favorable as the concentration of methanol increases.

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

The hydrophobic effect in a simple isotropic water-like model: Monte Carlo study

Using Monte Carlo computer simulations, we show that a simple isotropic water-like model with two characteristic lengths can reproduce the hydrophobic effect and the solvation properties of small and large non-polar solutes. Influence of temperature, pressure, and solute size on the thermodynamic properties of apolar solute solvation in a water model was systematically studied, showing two different solvation regimes. Small particles can fit into the cavities around the solvent particles, inducing additional order in the system and lowering the overall entropy. Large particles force the solvent to disrupt their network, increasing the entropy of the system. At low temperatures, the ordering effect of small solutes is very pronounced. Above the cross-over temperature, which strongly depends on the solute size, the entropy change becomes strictly positive. Pressure dependence was also investigated, showing a “cross-over pressure” where the entropy and enthalpy of solvation are the lowest. These results suggest two fundamentally different solvation mechanisms, as observed experimentally in water and computationally in various water-like models.

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

Consequences of Intrapore Liquids on Reactivity, Selectivity, and Stability for Aldol Condensation Reactions on Anatase TiO 2 Catalysts

Abstract This study provides evidence and mechanistic interpretations for the significant consequences of intrapore non‐polar liquids on acetone aldol condensation turnover rates, selectivity to primary dimer products, and catalyst stability for reactions at Lewis acid‐base site pairs on TiO 2 surfaces. These non‐polar liquids confer such benefits through the preferential stabilization of transition states (TS) for adsorption (“entry”) and desorption (“exit”) steps, which place their respective reactants or products within a solvating outer sphere environment. The extent to which non‐polar fluids ( n ‐heptane) form an intrapore liquid phase within TiO 2 voids was obtained from N 2 uptakes using established formalisms that consider the different molal volume, surface tension, and volatility between N 2 and n ‐heptane. Acetone condensation rates are limited by C−H activation, an “entry” step that forms bound prop‐1‐en‐2‐olates via a TS stabilized by intrapore liquids, leading to higher aldol condensation turnover rates as n ‐heptane pressure increases and active TiO 2 surfaces become increasingly immersed within a non‐polar liquid phase. These liquids solvate the late TS structures that mediate the desorption of primary C 6 condensation products even more effectively than those involved in prop‐1‐en‐2‐olate formation or in nucleophilic attack events that later form C−C bonds. Such preferential solvation favors desorption over C−C coupling events, thus disfavoring the formation of larger oligomers that become stranded at active sites, thus leading to much slower deactivation. Moreover, solvation by non‐polar liquids also leads to C 6 alkanones as the sole products formed in a single surface sojourn. These effects of a non‐polar dense phase circumvent the inherent stability, reactivity, and selectivity hurdles that have precluded practical aldol condensation catalysis on Lewis acid‐base pairs at oxide surfaces; these consequences are demonstrated here for TiO 2 catalysts, acetone aldol condensation reactions, and n ‐heptane as the non‐polar liquid but through strategies, concepts, and mechanistic features that extend to other systems. More generally, these observations and their mechanistic origins demonstrate how a contacting liquid preferentially solvates TS structures for elementary steps that involve either reactants arriving from or products entering into an outer sphere environment that contains a dense non‐polar phase.

Kadam, Shashikant A.↗
DOE OSTI2022

Film formation mechanisms in mixed-dimensional 2D/3D halide perovskite films revealed by in situ grazing-incidence wide-angle X-ray scattering

We report mixed dimensionality 2D/3D hybrid halide perovskites retain the stability of 2D perovskites (formula (A')2(A)n-1PbnI3n+1) and long diffusion lengths of the 3D materials (formula AMX3), thereby affording solar cells with extended stability as well as state of the art efficiencies approaching that of the 3D materials. These films are made by spin coating precursor solutions with arbitrarily large average layer thickness n ($\langle$n$\rangle$ > 7) to give films with both 2D and 3D phases present. While the 2D and 3D perovskite film formation mechanisms have been studied previously, little is understood about composite 2D/3D film formation. We used in situ grazing-incidence wide-angle scattering (GIWAXS) with synchrotron radiation to fully characterize films fabricated from precursor solutions with stoichiometries of (BA) 2 (MA) n-1 Pb n I 3n+1 ($\langle$n$\rangle$ = 3, 4, 5, 7, “12,” “50,” and ∞ (MAPbI 3 )). Four different mechanisms are seen depending on the stoichiometry in the precursor solution. When $\langle$n$\rangle$ = 3, oriented perovskite grows directly. For $\langle$n$\rangle$ = 4-5, the needle-shaped solvate (MA) 2 Pb 3 I 8 ·2DMF forms concomitantly with the oriented perovskite which leads to 2D/3D composite films. For $\langle$n$\rangle$ ≥ 7, this solvate dominates and unoriented perovskite grows in secondarily. For 3D perovskite, the solvate is the only phase formed. Kinetic analysis shows faster and earlier growth of the solvate with increasing $\langle$n$\rangle$. Using this knowledge, we show that the undesirable solvate mechanism can be reduced by addition of excess ammonium iodide in the precursor solution to give highly oriented films.

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

Continuous recovery of high-purity Kraft lignin from black liquor via simultaneous, liquid-phase acidification and purification

Here, a key advantage of the Sequential Liquid-Lignin Recovery and Purification (SLRP) process for the recovery of lignin from Kraft black liquor is the precipitation of a solvated, “liquid-lignin” phase (vs. the typical solid) during the CO 2 acidification step. Here, we have discovered that this lignin can also be maintained in a solvated state during the H 2 SO 4 acidification step if aqueous acetic acid from the Aqueous Lignin Purification with Hot Agents (ALPHA) process is incorporated into the H 2 SO 4 solution. As the liquid-lignin is converted from its salt form at pH = 9.5 to its acid form at pH = 2.5, the lignin initially solvated with water becomes solvated and thus liquefied with the organic–water solution. Maintaining the lignin in a solvated liquid state has two advantages: (1) the entire SLRP + ALPHA process, from starting black liquor to final lignin product, can be operated continuously; and (2) mass transfer of metals impurities out of the liquid-lignin phase is significantly enhanced, resulting in both a cleaner Kraft lignin product (~2000 ppm Na and 0.5-1.0 % ash) and reduced water washing.

09 BIOMASS FUELS↗
DOE OSTI2019

In Situ Liquid Secondary Ion Mass Spectrometry: A Surprisingly Soft Ionization Process for Investigation of Halide Ion Hydration

The understanding of ion solvation phenomenon is of significance due to their influences on many important chemical, biological and environmental processes. Mass spectrometry (MS) based methods have been used to investigate this topic with molecular insights. As ion-solvent interactions are weak, ionization processes should be as soft as possible in order to retain solvation structures. An in situ liquid secondary ion MS (SIMS) approach developed in our group has been recently utilized in investigations of Li ion solvation in non-aqueous solution and detected a series of solvated Li ions. As traditionally SIMS has long been recognized as a hard ionization process with strong damage occurring at the sputtering interface, it is very interesting to study further how soft in situ liquid SIMS can be. In this work, we used halide ion hydration as an example to compare the ionization performance of in situ liquid SIMS approach with regular electrospray ionization MS (ESI-MS). Results show that although ESI has been recognized as a soft ionization method, nearly no solvated halide ions were detected by ESI-MS analysis, which acquired only strong signals of salt ion clusters. In contrast, in liquid SIMS spectra, a series of obvious hydrated halide ion compositions could be observed. We further evaluated the hydration numbers of halide ions and revealed the effects of the ion size, charge density and polarizability on the hydration phenomenon. Our findings demonstrated that the in situ liquid SIMS approach is surprisingly soft, and expected to have very broad applications on investigation of various ion-solvent interactions and many other interesting chemical processes (e.g., the initial nucleation of nanoparticle formation) in liquid environment.

Zhang, Yanyan↗
DOE OSTI2021

Photoelectron Spectroscopy and Theoretical Study on Monosolvated Cyanate Analogue Clusters ECX-· Sol (ECX- = NCSe-, AsCSe- and AsCS-; Sol = H2O, CH3CN)

Six monosolvated cyanate analogue clusters ECX-?Sol (ECX– = NCSe-, AsCSe-, and AsCS-; Sol = H2O and CH3CN) were investigated using negative ion photoelectron spectroscopy (NIPES). NIPES experiments show that these clusters possess overall similar spectra compared to their respective isolated ECX- anions but shift to higher electron binding energy with solvent CH3CN stabilizing the excess electrons slightly more than H2O. For the ECX-·H2O series, vertical detachment energies and their increments relative to the bare species are measured to be 3.700 / 0.370, 3.085 / 0.425, 3.085 / 0.430 eV for NCSe-, AsCSe- and AsCS-, respectively, while the corresponding values in ECX-·CH3CN series are 3.835 / 0.505, 3.145 / 0.475, and 3.135 / 0.480 eV. Ab initio electronic structure calculations indicate the excess charges located at the terminal N and Se atoms in NCSe- and migrated to the central C atom in AsCSe- and AsCS-. For NCSe-, the solvation is driven by the interactions with the two negatively charged terminal ends, while for AsCSe- and AsCS-, the solvation revolves around the interactions with the central C atom, where all the excess negative charge is concentrated. Two nearly degenerate isomers for NCSe-·H2O are identified, one forming a single strong N…H-O hydrogen bond (HB) and the other featuring a bidentate HB with two hydroxyl H atoms pointing to N and Se ends. In contrast, the negative central C atom in AsCSe-/AsCS- allows the formation of a bifurcated HB with H2O. Similar effects are observed for the acetonitrile case, in which the three H atoms of methyl group interact with the two negatively charged terminal ends in NCSe-, while prefer binding to the central negative carbon atom in AsCSe-/AsCS-. The different binding motifs derived in this work may suggest distinctly different solvation properties in NCSe- versus AsCSe-/AsCS-, with the former anion leading to asymmetric solvation at the N end of the solute, while the latter species creating more ‘isotropic’ solvation around the central C equatorial plane. This work was supported by U.S. Department of Energy (DOE), Office of Science, Office of Basic Energy Sciences, Division of Chemical Sciences, Geosciences, and Biosciences, and performed using EMSL, a national scientific user facility sponsored by DOE’s Office of Biological and Environmental Research and located at Pacific Northwest National Laboratory, which is operated by Battelle Memorial Institute for the DOE. The theoretical calculations were conducted on EMSL’s “Cascade” Supercomputer.

Yuan, Qinqin↗
DOE OSTI2024

CH 3 Radical Generation in Microplasmas for Up-Conversion of Methane

The conversion of methane, CH 4 , into higher value chemicals using low temperature plasmas is challenged by both improving efficiency and selectivity. One path towards selectivity is capturing plasma produced methyl radicals, CH 3 , in a solvent for aqueous processing. Due to the rapid reactions of methyl radicals in the gas phase, the transport distance from production of the CH 3 to its solvation should be short, which then motivates the use of microplasmas. The generation of CH 3 in Ar/CH 4 /H 2 O plasmas produced in nanosecond pulsed dielectric barrier discharge microplasmas is discussed using results from a computational investigation. The microplasma is sustained in the channel of a microfluidic chip in which the solvent flows along one wall or in droplets. CH 3 is primarily produced by electron-impact of and dissociative excitation transfer to CH 4 , as well as CH 2 reacting with CH 4 . CH 3 is rapidly consumed to form C 2 H 6 which, in spite of being subject to these same dissociative processes, accumulates over time, as do other stable products including C 3 H 8 and CH 3 OH. The gas mixture and electrical properties were varied to assess their effects on CH 3 production. CH 3 production is largest with 5% CH 4 in the Ar/CH 4 /H 2 O mixture due to an optimal balance of electron-impact dissociation, which increases with CH 4 percentage, and dissociative excitation transfer and CH 2 reacting with CH 4 , which decrease with CH 4 percentage. Design parameters of the microchannels were also investigated. Increasing the permittivity of the dielectrics in contact with the plasma increased the ionization wave intensity which increased CH 3 production. Increased energy deposition per pulse generally increased CH 3 production as does lengthening pulse length up to a certain point. The arrangement of the solvent flow in the microchannel can also affect the CH 3 density and fluence to the solvent. The fluence of CH 3 to the liquid solvent is increased if the liquid is immersed in the plasma as a droplet or is a layer on the wall where the ionization wave terminates. The solvation dynamics of CH 3 with varying numbers of droplets was also examined. Here, the maximum density of solvated methyl radicals CH 3aq occurs with a large number of droplets in the plasma. However, the solvated CH 3aq density can rapidly decrease due to desolvation, emphasizing the need to quickly react the solvated species in the solvent.

03 NATURAL GAS↗
DOE OSTI2019

Effect of Mixed-Solvent Environments on the Selectivity of Acid-Catalyzed Dehydration Reactions

The composition of the liquid phase can alter the rates of individual reaction steps and thus alter the selectivity of acid-catalyzed reactions, but these solvent effects are difficult to anticipate for design purposes. Herein, we report the kinetics and selectivity of Brønsted acid-catalyzed 1,2-propanediol dehydration in pure water and in aqueous mixtures of the polar aprotic cosolvents γ-valerolactone, 1,4-dioxane, tetrahydrofuran, N-methyl-2-pyrrolidone, tetramethylene sulfoxide, and dimethyl sulfoxide at 433 K. We find that the major product of 1,2-propanediol dehydration is propanal in most mixed-solvent environments with selectivities between 1 and 68 mol %. In contrast, 1,2-propanediol dehydration in aqueous mixtures of dimethyl sulfoxide affords acetone as the major product with up to 48% selectivity with minimal propanal formation. We use classical molecular dynamics simulations to probe these solvent effects by computing the difference between the solvation free energies of 1,2-propanediol and propanal in aqueous mixtures of polar aprotic cosolvents and in pure water. We find that the difference in the solvation free energies is correlated with the rates of propanal formation in all mixed-solvent environments, indicating that the solvent-mediated stabilization of the product state relative to the reactant state translates to increased selectivity toward the same product. Similar agreement between simulated solvation free energies and experimental reaction rates/selectivities is established for the acid-catalyzed dehydration of cis- and trans-1,2-cyclohexanediol and 1,3-cyclohexanediol. Finally, analysis of the solvation environment around 1,2-propanediol shows that dimethyl sulfoxide uniquely competes against water to solvate reactive hydroxyl groups, which causes a change in reaction mechanism in this solvent system that leads to the formation of acetone rather than propanal. Here, these results represent a step toward the computationally efficient screening of solvent systems for acid-catalyzed, liquid-phase processes.

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

Applications of Polarizable Continuum Models To Determine Accurate Solution-Phase Thermochemical Values Across a Broad Range of Cation Charge – The Case of U(III–VI)

Contributing factors to the solution-phase correction to the free energy of the molecular clusters U(H 2 O) n 3+/4+ and UO 2 (H 2 O) m 1+/2+ ( n = 8, 9, 30, 41, 77; m = 4, 5, 30, 41, 77) have been examined as a function of cavity type in the integrated-equation-formalism-protocol (IEF) and SMD polarizable continuum models (PCMs). It is observed that the free energy correction, G corr , does not smoothly converge to zero as the number of explicitly solvating water molecules approaches the bulk limit, and the convergence behavior varies significantly with cavity and model. The rates of convergence of the gas-phase hydration energy, Δ G hyd , wherein the bare metal ion is inserted into a molecular water cluster and Δ G corr for the reaction exhibit wide variations as a function of ion charge, cavity, and model. This is the likely source of previously reported discrepancies in predicted free energies of solvation for metal ions when using different PCM cavities and/or models. The cancellation of errors in Δ G hyd and Δ G corr is optimal for clusters consisting of only a second solvation shell of explicit water molecules ( n = m = 30). Finally, the UFF cavity within IEF, in particular, exhibits the most consistent cancellation of errors when using a molecular cluster consisting of a second shell of solvating water for all oxidation states of uranium, leading to accurate free energies of solvation Δ G solv for these species.

chemistry↗
DOE OSTI2022

A Sobering Examination of the Feasibility of Aqueous Aluminum Batteries

Recent reports on aqueous aluminum (Al) batteries boast encouraging metrics including high energy density, nonflammability, and low raw material cost. However, the poorly understood solvation environment around the trivalent aluminum ions (Al3+) raises uncertainty regarding their feasibility. Furthermore, the corrosiveness of the acidic solutions and degree of hydration of aluminum trifluoromethanesulfonate (Al(OTf)3) salt is routinely mischaracterized. This study provides a rigorous examination of aqueous Al electrolytes, with the first compelling evidence for a dynamic octahedral solvation structure around Al3+ dominated by labile water and OH-, without Al-OTf contact ion pairs, at high concentrations. This solvation behavior and the concomitant electrostatic hydrolysis differs greatly from previously reported water-in-salt electrolytes, and occurs due to the high charge density of the Lewis acidic Al3+. Nuclear magnetic resonance spectroscopy and other physicochemical measurements quantitatively reveal how species activity evolves with concentration. The high proton activity observed in transport and electrochemical measurements constitutes a prohibitive barrier to the feasibility of rechargeable aluminum batteries with aqueous electrolyte.

Pastel, Glenn R.↗
DOE OSTI2022

Molecular engineering of fluoroether electrolytes for lithium metal batteries

Fluoroether solvents are promising electrolyte candidates for high-energy-density lithium metal batteries, where high ionic conductivity and oxidative stability are important metrics for design of new systems. Recent experiments have shown that these performance metrics, particularly stability, can be tuned by changing the fraction of ether and fluorine content. However, little is known about how different molecular architectures influence the underlying ion transport mechanisms and conductivity. Here, we use all-atom molecular dynamics simulations to elucidate the ion transport and solvation characteristics of fluoroether chains of varying length, and having different ether segment and fluorine terminal group contents. The design rules that emerge from this effort are that solvent size determines lithium-ion transport kinetics, solvation structure, and solvation energy. In particular, the mechanism for lithium-ion transport is found to shift from ion hopping between solvation sites located in different fluoroether chains in short-chain solvents, to ion–solvent co-diffusion in long-chain solvents, indicating that an optimum exists for molecules of intermediate length, where hopping is possible but solvent diffusion is fast. Consistent with these findings, our experimental measurements reveal a non-monotonic behavior of the effects of solvent size on lithium-ion conductivity, with a maximum occurring for medium-length solvent chains. Finally, a key design principle for achieving high ionic conductivity is that a trade-off is required between relying on shorter fluoroether chains having high self-diffusivity, and relying on longer chains that increase the stability of local solvation shells.

25 ENERGY STORAGE↗
DOE OSTI2026

Integrated machine learning-molecular dynamics framework for electrolyte property prediction

Electrochemical stability windows determine the operating range of battery electrolytes, yet accurate prediction remains challenging because stability emerges from statistical ensembles of local solvation environments rather than single ground-state molecular structures. Traditional density functional theory calculations on energy-minimized clusters cannot capture the thermal variations in local coordination environments and geometries that govern decomposition, while SMILES-based machine learning methods lack explicit representation of three-dimensional solvation structure and ion pairing. Here, we introduce a structure-aware machine learning framework that predicts frontier orbital energies (HOMO and LUMO) directly from molecular dynamics-sampled solvation configurations, achieving sub-0.6 eV accuracy at computational costs 3–4 orders of magnitude lower than first-principles methods. Across twelve representative battery electrolytes, we demonstrate that solvent-separated and contact ion pairs exhibit strong size- and local chemistry dependent electronic stability, with variations in coordination shifts of HOMO or LUMO level by 2–3 eV, and that extended solvation structure and partially desolvated environment further modulate stability by up to 3 eV. By encoding the statistical nature of electrochemical failure through ensemble sampling of explicit solvation geometries, our approach enables high-throughput screening and rational design of next-generation battery electrolytes with mechanistic understanding of structure–property relationships.

Energy - Storage↗
DOE OSTI2018

Solvophobic and solvophilic contributions in the water-to-aqueous guanidinium chloride transfer free energy of model peptides

We report on the solvation free energy of two different conformations (helix and extended) of two different peptides (deca-alanine and deca-glycine) in two different solvents (water and aqueous guanidinium chloride, GdmCl). The free energies are gained using the quasichemical organization of the potential distribution theorem, an approach that naturally provides the repulsive (solvophobic or cavity) and attractive (solvophilic) contributions to solvation. The solvophilic contribution is further parsed into a chemistry contribution arising from solute interaction with the solvent in the first solvation shell and a long-range contribution arising from non-specific interactions between the solute and the solvent beyond the first solvation shell. The cavity contribution is obtained for two different envelopes, Σ SE , which theory helps identify as the solvent excluded volume, and Σ G , a larger envelope beyond which solute-solvent interactions are Gaussian. The Σ SE envelope is independent of the solvent, as expected on the basis of the insensitivity to the solvent type of the distance of closest approach between protein heavy atoms and solvent heavy atoms, but contrary to the intuition based on treating solvent constituents as spheres of some effective radii. For both envelopes, the cavity contribution in water is proportional to the surface area of the envelope. The same does not hold for GdmCl(aq), revealing the limitation of using molecular area to assess solvation energetics. The Σ G -cavity contribution predicts that GdmCl(aq) should favor the more compact state, contrary to the role of GdmCl in unfolding proteins. The chemistry contribution attenuates this effect, but still the net local (chemistry plus Σ G -packing) contribution is inadequate in capturing the role of GdmCl. With the inclusion of the long-range contribution, which is dominated by van der Waals interaction, aqueous GdmCl favors the extended conformation over the compact conformation. Our result emphasizes the importance of weak, but attractive, long-range dispersion interactions in protein solution thermodynamics.

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

Probing Chromophore Energetics and Couplings for Singlet Fission in Solar Cell Applications: Final technical report

The overall goal of this project is to provide a molecular-level understanding of the mechanisms and requirements for efficient exciton multiplication through singlet fission processes. Such understanding would open the paths for improving alternative sources of energy including an improved efficiency of dye-sensitized solar cells. Using high-resolution anion photoelectron spectroscopy, electronic structures of chromophores are explored by probing directly the singlet and triplet states involved in singlet fission, determining precisely their energetics, vibrational structures and couplings. To carry out these experiments, we built a flexible-source instrument that allows for cryogenic cooling of precursor anions, providing improved spectral resolutions. Moreover, our velocity map imaging spectrometer includes novel design considerations that provide better focus with less stringent parameters. Using this instrument, we studied the three-ring polycyclic aromatic hydrocarbon molecules, anthracene and fluoranthene, the latter can be considered as a covalently-bonded benzyl and naphthyl system. Starting from the radical anion ground state, we accessed the singlet and triplet states at the same time, and the results revealed intriguing vibronic coupling signatures in both electronic states, which may possibly be related to singlet fission mechanisms. We further studied aromatic oxide molecules, phenoxy, 1-naphthoxy, and 2-naphthoxy, whose large molecular dipole moments supported highly delocalized dipole bound states. An intrinsic property of singlet fission is the delocalization of excitation to nearby chromophore molecules. In an effort to extend our studies to solvated chromophore clusters, we development new instrumental capabilities that allowed for formation of weakly bound clusters around any ion by utilizing liquid-nitrogen cooled octupole ion traps. In characterizing this new approach for gas phase cluster formation, we studied the structures of a variety of solvated clusters, from solvated ionic liquid cation, where intermolecular interactions are dominated by electrostatic interactions, to solvated peptides, where hydrogen-bonding with solvent molecules can directly alter the peptide structure. Overall, we found that the formation of weakly-bound clusters in a cryogenic ion trap can be a viable approach for accessing structures that better represent those found in condensed phase. In the process, we also implemented a straightforward double-resonance setup and established a protocol for isotopic substitution using cryogenic ion traps, both of which aided in decongesting crowded experimental spectra that often accompany these larger solvated systems.

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

Quantitative Prediction of Uranium Speciation and Amidoxime Binding in Seawater from Advanced Simulation Techniques

The goal of this project was to quantitatively predict stability of uranium species and thermochemistry of uranium binding with amidoxime ligands in realistic seawater conditions. The specific aims of the project were: (i) determine the stable species of uranium and its solvation structure in realistic seawater conditions; also speciation and solvation of uranium's major competing species in seawater, vanadium; (ii) accurately predict the free energies (and thereby stability constants) of binding between uranium and amidoxime; (iii) accurately predict the reaction enthalpies and entropies of binding between uranium and amidoxime as well as between vanadium and amidoxime. Main results and conclusions are as follows. First-principles molecular dynamic simulations found that the structure of the Ca 2 UO 2 (CO 3 ) 3 complex is very stable and that one Ca ion binds to the center UO 2 (CO 3 ) 3 4- anion stronger than the other Ca ion. This finding suggests that using time-resolved EXAFS spectra may confirm the asymmetry in binding of the two Ca ions in the aqueous Ca 2 UO 2 (CO 3 ) 3 complex. To consider the common ions, solvation of the Ca 2 UO 2 (CO 3 ) 3 complex in seawater was simulated by classical molecular dynamics simulation. It was found that the structure of the Ca 2 UO 2 (CO 3 ) 3 complex is very stable in the model seawater. A Na+ ion was found to be closely associated with the Ca 2 UO 2 (CO 3 ) 3 complex by indirectly interacting with one axial oxygen atom of the UO 2 group bridged by a water molecule. In addition, the Na + ion interacts closely with one Ca 2+ ion than the other. The present simulations revealed the key role of common ions such as Na + in impacting the solvation, structure, and apparent charge of the Ca 2 UO 2 (CO 3 ) 3 complex in seawater. The sequential processes of UO 2 2+ binding with the three CO 3 2- groups and the two Ca 2+ ions in pure water and in 0.1 M NaCl were simulated by classical molecular dynamics with both the non-constrained brute-force approach and umbrella sampling. The simulated free energies show excellent agreement with the experiment. The coordination structures of UO 2 (CO 3 ) 3 4- and CaUO 2 (CO 3 ) 3 2- were found to be significantly affected by the presence of Na + ions, leading to a monodentate binding of a carbonate group to U and a water molecule entering the first coordination shell of U. To shed light on the binding between uranyl and glutardiamidoxime (H 2 B), an important model and state-of-the-art ligand for seawater uranium extraction, a suite of computational methods was employed. From molecular dynamics (MD) simulations, it was found that the binding configurations of B 2- with uranyl favor the twofold distorted η 2 binding between the oximate ends (C=N-O-) and U, while HB - prefers the chelating mode for the oximate end with the neutral end being solvated by water. The free energies of sequential ligand binding to form UO 2 B, [UO 2 B 2 ] 2- , and [UO 2 (HB)B] - were simulated with umbrella sampling and very good agreement with the experimental values was achieved, which corroborates the structural insights into the binding mode. State-of-the-art polymeric sorbents employ both amidoximate and carboxylate groups on the side chains to achieve optimal U uptake and selectivity, so we simulated the binding of a model amidoximate–carboxylate bifunctional ligand with uranyl. Classical MD and free-energy simulations in 0.5 M NaCl showed that the carboxylate group binds first to uranyl, leading to a loose intermediate state, and then, the amidoximate group binds, resulting in a more stable and tight chelate state. Binding of the second bifunctional ligand follows a similar process, and the two ligands prefer a trans configuration around the uranyl group. The simulated free energies indicate that the two bifunctional ligands bind with uranyl 14 kcal/mol stronger than the two ligands with only amidoximate groups, confirming an important synergy between amidoximate and carboxylate groups in binding uranyl.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗