Search NASA⌕ Search

SEARCH · Search NASA

Results for “Solvation models”

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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Effect of Solvent Composition on Non-DLVO Forces and Oriented Attachment of Zinc Oxide Nanoparticles

Oriented attachment (OA) occurs when nanoparticles in solution align their crystallographic axes prior to colliding and subsequently fuse into single crystals. Traditional colloidal theories such as DLVO provide a framework for evaluating OA but fail to capture key particle interactions due to the atomistic details of both the crystal structure and the interfacial solution structure. Using zinc oxide as a model system, we investigated the effect of the solvent on short-ranged and long-ranged particle interactions and the resulting OA mechanism. In situ TEM imaging showed that ZnO nanocrystals in toluene undergo long-range attraction comparable to 1kT at separations of 10 nm and 3kT near particle contact. These observations were rationalized by considering non-DLVO interactions, namely dipole-dipole forces and torques between the polar ZnO nanocrystals. Langevin dynamics simulations showed stronger interactions in toluene compared to methanol solvents, consistent with the experimental results. Concurrently, we performed atomic force microscopy measurements using ZnO-coated probes for the short-ranged interaction. Our data provided valuable insights into another type of non-DLVO interaction, namely the repulsive solvation force. Specifically, the solvation force was stronger in water compared to ethanol and methanol, due to the stronger hydrogen bonding and denser packing of water molecules at the interface. In conclusion, our results highlight the importance of non-DLVO forces in a general framework for understanding and predicting particle aggregation and attachment.

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

Adsorption and solvation of ethanol at the water liquid-vapor interface: a molecular dynamics study

The free energy profiles of methanol and ethanol at the water liquid-vapor interface at 310K were calculated using molecular dynamics computer simulations. Both alcohols exhibit a pronounced free energy minimum at the interface and, therefore, have positive adsorption at this interface. The surface excess was computed from the Gibbs adsorption isotherm and was found to be in good agreement with experimental results. Neither compound exhibits a free energy barrier between the bulk and the surface adsorbed state. Scattering calculations of ethanol molecules from a gas phase thermal distribution indicate that the mass accommodation coefficient is 0.98, and the molecules become thermalized within 10 ps of striking the interface. It was determined that the formation of the solvation structure around the ethanol molecule at the interface is not the rate-determining step in its uptake into water droplets. The motion of an ethanol molecule in a water lamella was followed for 30 ns. The time evolution of the probability distribution of finding an ethanol molecule that was initially located at the interface is very well described by the diffusion equation on the free energy surface.

NASA Discipline Exobiology↗
NASA NTRS1996

Mechanism of unassisted ion transport across membrane bilayers

To establish how charged species move from water to the nonpolar membrane interior and to determine the energetic and structural effects accompanying this process, we performed molecular dynamics simulations of the transport of Na+ and Cl- across a lipid bilayer located between two water lamellae. The total length of molecular dynamics trajectories generated for each ion was 10 ns. Our simulations demonstrate that permeation of ions into the membrane is accompanied by the formation of deep, asymmetric thinning defects in the bilayer, whereby polar lipid head groups and water penetrate the nonpolar membrane interior. Once the ion crosses the midplane of the bilayer the deformation "switches sides"; the initial defect slowly relaxes, and a defect forms in the outgoing side of the bilayer. As a result, the ion remains well solvated during the process; the total number of oxygen atoms from water and lipid head groups in the first solvation shell remains constant. A similar membrane deformation is formed when the ion is instantaneously inserted into the interior of the bilayer. The formation of defects considerably lowers the free energy barrier to transfer of the ion across the bilayer and, consequently, increases the permeabilities of the membrane to ions, compared to the rigid, planar structure, by approximately 14 orders of magnitude. Our results have implications for drug delivery using liposomes and peptide insertion into membranes.

NASA Discipline Exobiology↗
DOE OSTI2020

Integrating solvation shell structure in experimentally driven molecular dynamics using x-ray solution scattering data

In the past few decades, prediction of macromolecular structures beyond the native conformation has been aided by the development of molecular dynamics (MD) protocols aimed at exploration of the energetic landscape of proteins. Yet, the computed structures do not always agree with experimental observables, calling for further development of the MD strategies to bring the computations and experiments closer together. Here, we report a scalable, efficient MD simulation approach that incorporates an x-ray solution scattering signal as a driving force for the conformational search of stable structural configurations outside of the native basin. We further demonstrate the importance of inclusion of the hydration layer effect for a precise description of the processes involving large changes in the solvent exposed area, such as unfolding. Utilization of the graphics processing unit allows for an efficient all-atom calculation of scattering patterns on-the-fly, even for large biomolecules, resulting in a speed-up of the calculation of the associated driving force. The utility of the methodology is demonstrated on two model protein systems, the structural transition of lysine-, arginine-, ornithine-binding protein and the folding of deca-alanine. We discuss how the present approach will aid in the interpretation of dynamical scattering experiments on protein folding and association.

Hsu, Darren J.↗
NASA NTRS1990

Cavities in molecular liquids and the theory of hydrophobic solubilities

Thermal configurational data on neat liquids are used to obtain the work of formation of hard spherical cavities of atomic size in six molecular solvents: n-hexane, n-dodecane, n-undecyl alcohol, chloroform, carbon tetrachloride, and water. These results are used to test a recent suggestion that the differences between nonaqueous solvents and liquid water in solvation of inert gases are not principally due to the hydrogen-bonded structure of liquid water but rather to the comparatively small size of the water molecule. The frequencies of occurrence of cavities in liquid water can be meaningfully distinguished from those in the organic solvents. Liquid water has a larger fractional free volume, but that free volume is distributed in smaller packets. With respect to cavity work, water is compared to a solvent of the same molecular density and composed of hard spheres of the same size as the water molecule. That comparison indicates that the hard-sphere liquid finds more ways to configure its free volume in order to accommodate an atomic solute of substantial size and thus, would be more favorable solvent for inert gases. The scaled particle model of inert gas solubility in liquid water predicts cavity works 20% below the numerical data for TIP4P water at 300 K and 1.0 g/cm3 for cavity radii near 2.0 angstroms. It is argued that the sign of this difference is just the sign that ought to be expected and that the magnitude of this difference measures structural differences between water and the directly comparable hard-sphere liquid. In conjunction with previous data, these results indicate that atomic sized cavities should be considered submacroscopic.

NASA Discipline Exobiology↗
DOE OSTI2024

Computing equilibrium free energies through a nonequilibrium quench

Many methods to accelerate sampling of molecular configurations are based on the idea that temperature can be used to accelerate rare transitions. These methods typically compute equilibrium properties at a target temperature using reweighting or through Monte Carlo exchanges between replicas at higher temperatures. A recent paper [G. M. Rotskoff and E. Vanden-Eijnden, Phys. Rev. Lett. 122, 150602 (2019)] demonstrated that accurate equilibrium densities of states can also be computed through a nonequilibrium “quench” process, where sampling is performed at a higher temperature to encourage rapid mixing and then quenched to lower energy states with dissipative dynamics. Here, we provide an implementation of the quench dynamics in LAMMPS and evaluate a new formulation of nonequilibrium estimators for the computation of partition functions or free energy surfaces (FESs) of molecular systems. We show that the method is exact for a minimal model of N-independent harmonic springs and use these analytical results to develop heuristics for the amount of quenching required to obtain accurate sampling. We then test the quench approach on alanine dipeptide, where we show that it gives an FES that is accurate near the most stable configurations using the quench approach but disagrees with a reference umbrella sampling calculation in high FE regions. We then show that combining quenching with umbrella sampling allows the efficient calculation of the free energy in all regions. Moreover, by using this combined scheme, we obtain the FES across a range of temperatures at no additional cost, making it much more efficient than standard umbrella sampling if this information is required. Finally, we discuss how this approach can be extended to solute tempering and demonstrate that it is highly accurate for the case of solvated alanine dipeptide without any additional modifications.

Chemistry↗
DOE OSTI2022

Small angle neutron scattering and lipidomic analysis of a native, trimeric PSI-SMALP from a thermophilic cyanobacteria

The use of styrene-maleic acid copolymers (SMAs) to produce membrane protein-containing nanodiscs without the initial detergent isolation has gained significant interest over the last decade. We have previously shown that a Photosystem I SMALP from the thermophilic cyanobacterium, Thermosynechococcus elongatus (PSI-SMALP), has much more rapid energy transfer and charge separation in vitro than detergent isolated PSI complexes. In this study, we have utilized small-angle neutron scattering (SANS) to better understand the geometry of these SMALPs. These techniques allow us to investigate the size and shape of these particles in their fully solvated state. Further, the particle's proteolipid core and detergent shell or copolymer belt can be interrogated separately using contrast variation, a capability unique to SANS. Here we report the dimensions of the Thermosynechococcus elongatus PSI-SMALP containing a PSI trimer. At ~1.5 MDa, PSI-SMALP is the largest SMALP to be isolated; our lipidomic analysis indicates it contains ~1300 lipids/per trimeric particle, >40-fold more than the PSI-DDM particle and > 100 fold more than identified in the 1JB0 crystal structure. Interestingly, the lipid composition to the PSI trimer in the PSI-SMALP differs significantly from bulk thylakoid composition, being enriched ~50 % in the anionic sulfolipid, SQDG. Finally, utilizing the contrast match point for the SMA 1440 copolymer, we also can observe the ~1 nm SMA copolymer belt surrounding this SMALP for the first time, consistent with most models of SMA organization.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

3-D Geological Modeling for Numerical Flow Simulation Studies of Gas Hydrate Reservoirs at the Kuparuk State 7-11-12 Pad in the Prudhoe Bay Unit on the Alaska North Slope

Accurate reservoir evaluation requires reliable three-dimensional (3-D) geological models. Here, this study conducted 3-D geological modeling for numerical flow simulation of the B1 sand gas hydrate reservoir at the Kuparuk State 7-11-12 pad, Prudhoe Bay Unit, Alaska North Slope. The model integrates well logs, core, and seismic data to address spatial heterogeneity in geological structures and reservoir properties. Two modeling types were performed: structural framework modeling and petrophysical property modeling. For structural framework modeling, seismic data and well log markers were used to reproduce subsurface structures characterized by a normal fault system. A volume-based modeling algorithm and stair-stepping grid were applied. The resulting 3-D model comprised 2,640,000 grid cells across 264 layers, including seven fault grids. For petrophysical property modeling, total porosity was initially modeled using sequential Gaussian simulation with collocated cokriging. To reproduce the upward coarsening of the B1 sand, upscaled log-derived total porosity and a three-dimensional (3-D) trend depicting total porosity variation were used as primary and secondary data, respectively. Gas hydrate saturation distribution was modeled similarly, with secondary data from estimated porosity distribution and seismic-derived acoustic impedance map enhancing accuracy. Results indicate higher gas hydrate saturation in the upper part of the B1 sand and areas with higher acoustic impedance. Intrinsic permeability was modeled from the total porosity and clay-bound water volume, and effective permeability was derived from the gas hydrate saturation and intrinsic permeability distributions based on the “Tokyo model”. Effective permeability distributions were influenced by the total porosity, gas hydrate saturation, and intrinsic permeability. Within the same layer, higher gas hydrate saturation leads to decreased effective permeability. In total, 100 sets of multiple scenarios were prepared, providing input data for dynamic flow simulations to evaluate the effects of lateral heterogeneity in reservoir properties and the hydraulic characteristics of faults on production behavior for preassessment before the long-term production test.

58 GEOSCIENCES↗
DOE OSTI2022

Machine learning enabled quantification of the hydrogen bonds inside the polyelectrolyte brush layer probed using all-atom molecular dynamics simulations

The configuration of densely grafted charged polyelectrolyte (PE) brushes is strongly dictated by the properties and behavior of the counterions that screen the PE brush charges and the solvent molecules (typically water) that solvate the brush molecules and these screening counterions. Only recently, efforts have been made to study the PE brushes atomistically, thereby shedding light on the properties of brush-supported ions and water molecules. However, even for such efforts, there are limitations associated with using a generic definition to estimate certain properties of water and ions inside the brush layer. For example, water–water hydrogen bonds (HBs) will behave differently for locations outside and inside the brush layer, given the fact that the densely closely grafted PE brush molecules create a soft nanoconfinement where the water connectivity becomes highly disrupted: therefore, using the same definition to quantify the HBs inside and outside the brush layer will be unwise. In this paper, we address this limitation by employing an unsupervised machine learning (ML) approach to predict the water–water hydrogen bonding inside a cationic PE brush layer modeled using all-atom molecular dynamics (MD) simulations. Here, the ML method, which relies on a clustering approach and uses the equilibrium coordinates of the water molecules (obtained from the all-atom MD simulations) as the input, is capable of identifying the structural modification of water–water HBs (revealed through appropriate clustering of the data) inside the PE brush layer induced soft nanoconfinement. Such capabilities would not have been possible by using a generic definition of the HBs. Our calculations lead to four key findings: (1) the clusters formed inside and outside the brush layer are structurally similar; (2) the margin of the cluster is shorter inside the PE brush layer confirming the possible disruption of the HBs inside the PE brush layer; (3) the average “hydrogen–acceptor-oxygen–donor-oxygen” angle that defines the HB is reduced for the HBs formed inside the brush layer; (4) the use of the generic definition (definition usable for characterizing the HBs in brush-free bulk) leads to an overprediction of the number of HBs formed inside the PE brush layer.

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

Selecting the Optimal Fluorinated Ether Co-Solvent for Lithium Metal Batteries

To guide the selection of a suitable fluorinated ether (FE) co -solvent for lithium metal batteries, it is crucial to understand the relationship between the organic structures of the FEs and the electrochemical performance of an FE-containing electrolyte. In this work, 1,1,2,2-tetrafluoro-1-(2,2,2-trifluoroethoxy)ethane (FEE), 1,1,2,2-tetrafluoro-3-(1,1,2,2-tetrafluoroethoxy)-propane (TTE), and 1,2-bis(1,1,2,2-tetrafluoroethoxy)ethane (OFDEE) were chosen as representative FE co-solvents because of their distinct structural properties. The structure-activity relationship between the FEs and the electrochemical performance of Li||LiNi 0.6 Mn 0.2 Co 0.2 O 2 (Li||NMC622) cells was correlated and quantified by Fourier-transform infrared and multi-dimensional nuclear magnetic resonance techniques. Sand's model was also employed to assess the extent of lithium dendrite formation in the cells using various FE electrolytes. The cycling performance of Li||NMC622 cells using different FE co-solvents follows the order FEE > TTE > OFDEE. Further, since the direct measurement of Sand's time is difficult, we introduced relative Sand's time to probe the diffusion behavior of each electrolyte, and the results showed that the best performance was obtained in the electrolyte with the longest relative Sand's time. Moreover, the lithium metal cell using the electrolyte with FEE co-solvent showed similar capacity retention compared with the baseline electrolyte at room temperature, but it demonstrated significantly improved low-temperature performance. The results indicate that FEE is a promising co-solvent candidate for improving the low-temperature performance of lithium metal batteries because it possesses not only non-solvating behavior but also very low viscosity and non-flammability. The advanced electrolyte LiPF 6 -FEC-DMC-FEE enables very stable cycling of lithium metal batteries at various temperatures.

25 ENERGY STORAGE↗
DOE OSTI2019

Influence of Electronic Polarization on the Spectral Density

Accurate spectral densities are necessary for computing realistic exciton dynamics and nonlinear optical spectra of chromophores in condensed-phase environments, including multichromophore pigment–protein systems. However, due to the significant computational cost of computing spectral densities from first principles, requiring many thousands of excited-state calculations, most simulations of realistic systems rely on treating the environment as fixed-point charges. In this work, using a number of representative systems ranging from solvated chromophores to the photoactive yellow protein (PYP), we demonstrate that the quantum mechanical (QM) electronic polarization of the environment is key to obtaining accurate spectral densities and line shapes within the cumulant framework. We show that the QM environment can enhance or depress the coupling of fast chromophore degrees of freedom to the energy gap, altering the electronic–vibrational coupling and the resulting vibronic progressions in the absorption spectrum. In analyzing the physical origin of peaks in the spectral density, we identify vibrational modes that couple the electron and the hole as being particularly sensitive to the QM screening of the environment. For PYP, we reveal the need for careful determination of the appropriate QM region to obtain reliable spectral densities. Our findings suggest that the QM polarization of the environment can be crucial not just for excitation energies but also for electronic–vibrational coupling in complex systems with implications for the correct modeling of linear and nonlinear optical spectroscopy in the condensed phase as well as energy transfer in pigment–protein complexes.

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

Revealing the roles of the solid–electrolyte interphase in designing stable, fast-charging, low-temperature Li-ion batteries

Designing the solid–electrolyte interphase (SEI) is critical for stable, fast-charging, low-temperature Li-ion batteries. Fostering a “fluorinated interphase,” SEI enriched with LiF, has become a popular design strategy. Although LiF possesses low Li-ion conductivity, many studies have reported favorable battery performance with fluorinated SEIs. Such a contradiction suggests that optimizing SEI must extend beyond chemical composition design to consider spatial distributions of different chemical species. In this work, we demonstrate that the impact of a fluorinated SEI on battery performance should be evaluated on a case-by-case basis. Sufficiently passivating the anode surface without impeding Li-ion transport is key. We reveal that a fluorinated SEI containing excessive and dense LiF severely impedes Li-ion transport. In contrast, a fluorinated SEI with well-dispersed LiF (i.e., small LiF aggregates well mixed with other SEI components) is advantageous, presumably due to the enhanced Li-ion transport across heterointerfaces between LiF and other SEI components. An electrolyte, 1 M LiPF 6 in 2-methyl tetrahydrofuran (2MeTHF), yields a fluorinated SEI with dispersed LiF. This electrolyte allows anodes of graphite, μSi/graphite composite, and pure Si to all deliver a stable Coulombic efficiency of 99.9% and excellent rate capability at low temperatures. Pouch cells containing layered cathodes also demonstrate impressive cycling stability over 1,000 cycles and exceptional rate capability down to −20 °C. Through experiments and theoretical modeling, we have identified a balanced SEI-based approach that achieves stable, fast-charging, low-temperature Li-ion batteries.

25 ENERGY STORAGE↗
DOE OSTI2020

Electron Paramagnetic Resonance Characteristics and Crystal Structure of a Tutton Salt Analogue: Copper-Doped Cadmium Creatininium Sulfate

Electron paramagnetic resonance and crystallographic studies on copper-doped cadmium creatininium sulfate (CdCrnS) were undertaken to study the characteristics of a copper-hexahydrate complex in an organic analogue of Tutton's salt. X-ray diffraction experiments determined the crystal structure of CdCrnS at both 100 and 298 K. CdCrnS, like Tutton salt, crystallizes in the monoclinic space group P2 1 /n. The unit cell contains two cadmium hexahydrate complexes, four creatininium ions, four sulfates, and four additional solvation waters. Both crystallography and EPR find that the doped copper replaces the cadmium in the structure. Single-crystal EPR measurements at room temperature determined the g and copper hyperfine (A Cu ) tensors (principal values: g = 2.437, 2.134, and 2.080 and A Cu = -327, -84.8, and 7.33 MHz). EPR spectra of the powder at room temperature gave g = 2.448, 2.125, and 2.085 and A Cu = -315, -75.0, and 35.0 MHz and at 110 K gave g = 2.462, 2.116, and 2.077 and A Cu = -340, -30.0, and 35.0 MHz. The room-temperature tensors are close to the "rigid lattice limit" values found in copper-doped Tutton salts but with a higher gmin and weaker A Cu x coupling than average. A small but measurable temperature dependency of the tensors indicated the presence of a dynamic Jahn-Teller (JT) effect. In addition, the EPR line width changed dramatically with temperature, which is like that found in all copper-doped Tutton crystals. Utilizing the model of Silver-Getz for the g-value variation gave an estimate for the energy difference (δ 12 = 640 cm -1 ) between the ground and next highest JT configurations. An empirical correlation appears to exist between δ 12 and g min and A Cu x for the copper hexahydrates studied in similar crystals. Finally, this suggests a relationship between the amount of unpaired spins in the copper d-orbital x lobe and the gap between wells of the adiabatic potential surface.

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

Ferroelectric Large Polarons and Defect Tolerance in Multi-Component Lead Halide Perovskites

Solvation plays a pivotal role in chemistry and biology. A solid-state analogy of solvation is polaron formation, but the magnitude of Coulomb screening is typically an order-of-magnitude weaker than that of solvation in aqueous solutions. Here the PI aims to explore a new class of polarons, the ferroelectric large polaron, which allows efficient Coulomb screening of an electron or hole by extended ordering of dipoles from symmetry-broken unit cells. This kind of local ordering is reflected in the ferroelectric-like THz dielectric responses of lead halide perovskites (LHPs) and may be partially responsible for their exceptional optoelectronic performances. A charge carrier may be localized to and/or induce the formation of nanoscale domain boundaries of locally ordered dipoles in a ferroelectric or paraelectric material, i.e., crystal structure with polar unit cells or polar fluctuations. The ability to form ferroelectric large polarons can result in the efficient screening of charge carriers from scattering with other charge carriers, with charged defects, and with longitudinal optical phonons, thus contributing to enhanced optoelectronic properties. During the past funding period, the PI has explored efficient charge carrier screening in three-dimensional (3D) LHP crystals and has developed the ferroelectric large polaron model to explain a range of carrier properties in these materials. The PI has also developed a new experimental tool, two-dimensional optical Kerr effect (2D-OKE), which is particularly powerful in probing photophysical properties with exquisite energy resolution at or near the bandgap. During the next funding periods, the PI aims to establish the applicability of the ferroelectric large polaron proposal as a general principle for the design/search of defect tolerant semiconductors in optoelectronics. Preliminary experiments on model semiconductor systems beyond LHPs have demonstrated the feasibility of the proposed research.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Using Machine Learning to Greatly Accelerate Path Integral Ab Initio Molecular Dynamics

Ab initio molecular dynamics (AIMD) has become one of the most popular and robust approaches for modeling complicated chemical, liquid, and material systems. However, the formidable computational cost often limits its widespread application in simulations of the largest-scale systems. Furthermore, the situation becomes even more severe in cases where the hydrogen nuclei may be better described as quantized particles using a path integral representation. Here, we present a computational approach that combines machine learning with recent advances in path integral contraction schemes, and we achieve a 2 orders of magnitude acceleration over direct path integral AIMD simulation while at the same time maintaining its accuracy.

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

How the choice of exchange–correlation functional affects DFT-based simulations of the hydrated electron

Hydrated electrons are anionic species that are formed when an excess electron is introduced into liquid water. Building an understanding of how hydrated electrons behave in solution has been a long-standing effort of simulation methods, of which density functional theory (DFT) has come to the fore in recent years. The ability of DFT to model the reactive chemistry of hydrated electrons is an attractive advantage over semi-classical methodologies; however, relatively few density functional approximations (DFAs) have been used for the hydrated electron simulations presented in the literature. Here, we simulate hydrated electron systems using a series of exchange–correlation (XC) functionals spanning Jacob’s ladder. We calculate a variety of experimental and other observables of the hydrated electron and compare the XC functional dependence for each quantity. We find that the formation of a stable localized hydrated electron is not necessarily limited to hybrid XC functionals and that some hybrid functionals produce delocalized hydrated electrons or electrons that react with the surrounding water at an unphysically fast rate. Here, we further characterize how different DFAs impact the solvent structure and predicted spectroscopy of the hydrated electron, considering several methods for calculating the hydrated electron’s absorption spectrum for the best comparison between structures generated using different density functionals. None of the dozen or so DFAs that we investigated are able to correctly predict the hydrated electron’s spectroscopy, vertical detachment energy, or molar solvation volume.

Ab-initio molecular dynamics↗
DOE OSTI2023

Using Metadynamics to Reveal Extractant Conformational Free Energy Landscapes

Understanding the impact of extractant functionalization on metal-binding energetics in liquid-liquid extraction is essential to guide the development of better separation processes. Traditionally, computational extractant design uses electronic structure calculations on metal-ligand clusters to determine the metal-binding energy of the lowest energy state. Although highly accurate, this approach does not account for all of the relevant physics encountered under experimental conditions. Such methodologies often neglect entropic contributions such as temperature effects and ligand flexibility, in addition to approximating solvent-extractant interactions with implicit solvent models. In this study, we use classical molecular dynamics simulations with an advanced sampling method, metadynamics, to map out extractant molecule conformational free energies in the condensed phase. Here we generate the complete conformational landscape in solution for a family of bidentate malonamide-based extractants with different functionalizations of the headgroup and the side chains. In particular, we show how such alkyl functionalization reshapes the free energy landscape, affecting the free energy penalty of organizing the extractant into the cis-like metal-binding conformation from the trans-like conformation of the free extractant in solution. Specifically, functionalizing alkyl tails to the center of the headgroup has a greater influence on increasing molecular rigidity and disfavoring the binding conformation than functionalizing side chains. These findings are consistent with trends in metal-binding energetics based on experimentally reported distribution ratios. We also consider a different bidentate extractant molecule, carbamoylmethylphosphine oxide, and show how the choice of solvent can further reshape the conformational energetic landscape. This study demonstrates the feasibility of using molecular dynamics simulations with advanced sampling techniques to investigate extractant conformational energetics in solution, which, more broadly, will enable extractant design that accounts for entropic effects and explicit solvation.

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

Diversity-oriented synthesis of polymer membranes with ion solvation cages

Microporous polymers feature shape-persistent free volume elements (FVEs), which are permeated by small molecules and ions when used as membranes for chemical separations, water purification, fuel cells, and batteries. It remains a significant challenge to identify FVEs with analyte specificity, due to difficulties in generating microporous polymer libraries with sufficient diversity for screening their properties. Here, we describe a diversity-oriented synthetic (DOS) strategy for microporous polymer membranes from which we identified those whose FVEs serve as solid solvation cages for lithium ions (Li + ). Furthermore, key elements of our strategy included diversification of bis(catechol)-type monomers via multi-component Mannich reactions to introduce Li+-coordinating functionality within individual FVEs, topology-enforcing polymerizations for generating macromolecular skeletal diversity for networking FVEs into different pore architectures, and several classes of on-polymer reactions for diversifying pore geometries and dielectric properties. Lead candidate polymer membranes featuring explicit ion solvation cages exhibited both higher ionic conductivity and higher cation transference number than control membranes where FVEs were aspecific, which indicates conventional bounds for membrane permeability and selectivity for ion transport can be overcome.4 These advantages are tied to enhanced Li + partitioning from the electrolyte when the cages are present, higher diffusion barriers for anions within the pores, and network-enforced restrictions on the number of solvent molecules bound to Li+ by comparison to the bulk electrolyte, which reduces the effective mass of the working ion. Such membranes show promise as anode-stabilizing interlayers in high-voltage lithium-metal batteries for electric mobility.

42 ENGINEERING↗