The role of surfactant force field on the properties of liquid/liquid interfaces
Not Available
Servis, Michael J.↗
SEARCH · Search NASA
Results for “classical molecular dynamics”
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.
Diffusion in undoped and Cr-doped amorphous UO 2
UO 2 fuel pellets are often doped with chromium oxide to obtain favourable properties such as higher density, improved thermal stability, large grain sizes, improved pellet-clad interaction margins, and increased fission gas retention during transients. Chromium has a low solubility limit in UO 2 , with past experimental work reporting solubility limits ranging between 0.004 to 0.06 wt.% Cr. Due to its low solubility, segregation of Cr ions to the grain boundary may occur. Further, the complexity of these boundaries may be high as observed in other ceramics resulting in disordered or amorphous regions along the boundary, affecting a range of material and operational properties of the fuel pellet. To assess these disordered regions, in this work we study amorphous undoped and Cr doped UO 2 systems (containing 10–50 at.% Cr 3+ ) that have been modelled using classical molecular dynamics methods incorporating Cr 3+ into the well-used CRG potential library. Diffusion coefficients, pre-exponential factors, and activation energies for diffusion were computed for oxygen ions, assessing the impact of structure and extrinsic species on migration. Oxygen diffusion was observed to be much faster in the undoped amorphous system compared to its crystalline counterpart. Oxygen diffusion in doped systems decreased with increasing Cr concentration, highlighting the importance of additives to retain fission products and other migratory species.
11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
Sulfonated polybenzimidazole membrane with graphene oxide additive for 2,3-butanediol/water separation: A molecular simulation
Membrane separation for 2,3-butanediol (2,3-BDO) recovery from fermentation broth is highly valued for sustainable and renewable processes, but it requires efficient membrane materials. Here, this work evaluates the sulfonated polybenzimidazole (sPBI) and its graphene oxide (GO) doped composite membrane for separating 2,3-BDO and water via atomistic simulations. Density functional theory calculations are applied to identify various forms of sPBI structures and quantify their binding interactions with 2,3-BDO and water. Classical molecular dynamic simulations are used to evaluate the structural changes, diffusivity, and selectivity of 2,3-BDO and water in different sPBI models, GO surfaces, and GO-doped sPBI composite models. Our results suggest that sPBI slightly increases the crystallinity of the membrane structures, enhances the adsorption strength for both 2,3-BDO and water, and improves the water/2,3-BDO selectivity by 2–3 times. The GO surfaces display a maximum selectivity at a surface coverage of 0.1–0.15 for both hydroxyl and epoxy surface groups. The addition of GO flakes to sPBI creates new interaction sites for 2,3-BDO and water at the interface of sPBI and GO, and the water/2,3-BDO selectivity of GO-doped sPBI models is further increased up to 3 times. This work illustrates how the integrated addition of sPBI and GO flakes offers a promising approach to selective separation of 2,3-BDO and water, providing theoretical guidance for polybenzimidazole-based membranes in the potential application of 2,3-BDO recovery.
2,3-butanediol↗
Crowded electrolytes containing redoxmers in different states of charge: Solution structure, properties, and fundamental limits on energy density
Nonaqueous redox flow batteries use liquid electrolytes containing redox-active organic molecules (redoxmers) as their energy storage medium. To maximize energy density, the redoxmer concentration needs to be maximized while maintaining low viscosity and high ionic conductivity. During charge, a redoxmer molecule pairs with an ion in the electrolyte while another ion migrates across the membrane to maintain electric neutrality. In a crowded electrolyte, this reconstitution changes physical and chemical properties of the solution. To explore these behaviors, a phenothiazine redoxmer fully miscible with acetonitrile was used, and electrochemical charge was mimicked by chemical oxidation. The solutions were examined using small-angle X-ray scattering, nuclear magnetic resonance, and conductometry and modeled using classical molecular dynamics. Overall, our study indicates that physical and structural properties of redoxmer solutions in both states of charge make it exceedingly difficult to increase the redoxmer concentrations over 2 M at any temperature without compromising dynamic properties of such solutions. The cause for this limitation is proximity to a gel-like regime in which fluidity, diffusivity, and ionic conductivity exponentially decrease with increasing concentration. This tendency is compounded by non-Arrhenius behavior of the electrolyte: a small increase in the concentration outruns gains in fluidity and conductivity at a higher temperature. Thus the properties of crowded electrolytes generally make it impossible to operate when gel-like behavior sets in. Pushing the redoxmer concentration to 2.5-3 M might be possible for small redoxmer molecules, but it would require the use of ionic liquid electrolytes at 340-360 K. (C) 2021 Elsevier B.V. All rights reserved.
25 ENERGY STORAGE↗
Predicting thermophysical properties of molten salts in the MgCl 2 -NaCl-KCl-LiCl system with a shell-model potential
Ternary eutectic salts composed of MgCl 2 , NaCl, and KCl, referred to as MNK salts, have recently emerged as promising candidates as high-temperature heat transfer fluids and thermal energy storage media. Here we performed classical molecular dynamics (MD) simulations to predict the densities, specific heat capacities, viscosities, and ionic self-diffusivities for MNK salts over a wide temperature range. The impact of LiCl additive on their thermophysical properties was also investigated. To capture the electronic polarization of Cl anions by neighboring cations, we developed a novel shell-model potential using the force-matching method and a dataset of ab initio calculated interatomic forces. Our extensive MD simulations predict structure and properties for pure salts and binary/ternary salt mixtures in the MgCl 2 -NaCl-KCl-LiCl system in overall good agreement with available experimental and theoretical data, which corroborates the accuracy and reliability of our developed potential.
36 MATERIALS SCIENCE↗
Unraveling the solvation phenomena of Ln 3+ cations in room-temperature ionic liquids: A computational study
Rare-earth elements (REEs), classified as critical materials, are difficult to separate due to their similar chemical properties and slight differences in ionic radius. Room-temperature ionic liquids (RTILs) have gained significant attention for REE separation because of their unique physicochemical properties and environmental advantages. In this study, we have investigated the solvation mechanisms of two RE cations (Ln 3+ ), Nd 3+ (light REE) and Yb 3+ (heavy REE), in two 1-butyl-3-methylimidazolium ([BMIM] + )-based RTILs using bis(trifluoromethylsulfonyl)imide ([NTf 2 ] − ) and acetate ([OAc] − ) as the respective anions, employing classical molecular dynamics (MD) simulations and density functional theory (DFT) techniques. From MD simulations, it is evident that Ln 3+ cations are primarily solvated by RTIL anions in the first solvation shell, while [BMIM] + cations form the second solvation shell with their numbers depending on the size and composition of the first solvation shell. The neutralizing [NO 3 − ] counterions are mainly solvated by [BMIM] + cations in their first solvation shell. Relative free energy change of solvation calculations using thermodynamic integration (TI) method indicate that Yb 3+ is more strongly solvated than Nd 3+ in both RTILs, a trend further supported by DFT calculations. While both methods predict consistent qualitative behavior, differences in models and energy evaluations lead to variations in absolute values. Overall, this study provides a comprehensive understanding of the solvation behavior of Ln 3+ cations in RTILs, demonstrating a stronger solvation preference for the heavy Ln 3+ cation (Yb 3+ ). In conclusion, these findings have implications for the design of RTIL-based separation processes for REEs.
Ash, Tamalika [Ames Lab., and Iowa State Univ., Am↗
PLAT domain protein 1 (PLAT1/PLAFP) binds to the Arabidopsis thaliana plasma membrane and inserts a lipid
Robust agricultural yields depend on the plant's ability to fix carbon amid variable environmental conditions. Over seasonal and diurnal cycles, the plant must constantly adjust its metabolism according to available resources or external stressors. The metabolic changes that a plant undergoes in response to stress are well understood, but the long-distance signaling mechanisms that facilitate communication throughout the plant are less studied. The phloem is considered the predominant conduit for the bidirectional transport of these signals in the form of metabolites, nucleic acids, proteins, and lipids. Lipid trafficking through the phloem in particular attracted our attention due to its reliance on soluble lipid-binding proteins (LBP) that generate and solubilize otherwise membrane-associated lipids. The Phloem Lipid-Associated Family Protein (PLAFP) from Arabidopsis thaliana is generated in response to abiotic stress as is its lipid-ligand phosphatidic acid (PA). PLAFP is proposed to transport PA through the phloem in response to drought stress. To understand the interactions between PLAFP and PA, nearly 100 independent systems comprised of the protein and one PA, or a plasma membrane containing varying amounts of PA, were simulated using atomistic classical molecular dynamics methods. In these simulations, PLAFP is found to bind to plant plasma membrane models independent of the PA concentration. When bound to the membrane, PLAFP adopts a binding pose where W41 and R82 penetrate the membrane surface and anchor PLAFP. This triggers a separation of the two loop regions containing W41 and R82. Subsequent simulations indicate that PA insert into the β -sandwich of PLAFP, driven by interactions with multiple amino acids besides the W41 and R82 identified during the insertion process. Finally, fine-tuning the protein-membrane and protein-PA interface by mutating a selection of these amino acids may facilitate engineering plant signaling processes by modulating the binding response.
59 BASIC BIOLOGICAL SCIENCES↗
Diffusion in intact secondary cell wall models of plants at different equilibrium moisture content
Secondary plant cell walls are composed of carbohydrate and lignin polymers, and collectively represent a significant renewable resource. Leveraging these resources depends in part on a mechanistic understanding for diffusive processes within plant cell walls. Common wood protection treatments and biomass conversion processes to create biorefinery feedstocks feature ion or solvent diffusion within the cell wall. X-ray fluorescence microscopy experiments have determined that ionic diffusion rates are dependent on cell wall hydration as well as the ionic species through non-linear relationships. In this work, we use classical molecular dynamics simulations to map the diffusion behavior of different plant cell wall components (cellulose, hemicellulose, lignin), ions (Na + , K + , Cu2 + , Cl - ) and water within a model for an intact plant cell wall at various hydration states (3-30 wt% water). From these simulations, we analyze the contacts between different plant cell wall components with each other and their interaction with the ions. Generally, diffusion increases with increasing hydration, with lignin and hemicellulose components increasing diffusion by an order of magnitude over the tested hydration range. Ion diffusion depends on charge. Positively charged cations preferentially interact with hemicellulose components, which include negatively charged carboxylates. As a result, positive ions diffuse more slowly than negatively charged ions. Measured diffusion coefficients are largely observed to best fit piecewise linear trends, with an inflection point between 10 and 15% hydration. These observations shed light onto the molecular mechanisms for diffusive processes within secondary plant cell walls at atomic resolution.
09 BIOMASS FUELS↗
Building Durable Multimetallic Electrocatalysts from Intermetallic Seeds
When combined with earth-abundant metals, Pt-based alloy nanoparticles (NPs) can be cost-effective electrocatalysts. However, these NPs can experience leaching of non-noble-metal components under harsh electrocatalytic conditions. In this work, the Skrabalak group has demonstrated a novel NP construct in which Pt-based random alloy surfaces are stabilized against non-noble-metal leaching by their deposition onto intermetallic seeds. These core@shell NPs are highly durable electrocatalysts, with the ability to tune catalytic performance by the core@shell architecture, surface alloy composition, and NP shape. This versatility was demonstrated in a model system in which random alloy (ra-) PtM surfaces were deposited onto ordered intermetallic (i-) PdCu seeds using seed-mediated co-reduction (SMCR). In the initial demonstration, ra-PtCu shells were deposited on i-PdCu seeds, with these core@shell NPs exhibiting higher specific and mass activities for the oxygen reduction reaction (ORR) when compared to similarly sized ra-PtCu NPs. These NPs also showed outstanding durability, maintaining ~85% in specific activity after 5000 cycles. Characterization of the NPs after use revealed minimal loss of Cu. The activity enhancement was attributed to the strained surface that arises from the lattice mismatch between the intermetallic core and random alloy surface. The outstanding durability was attributed to the ordered structure of the intermetallic core. The origin of this durability enhancement was investigated by classical molecular dynamics simulations, where Pt atoms were found to have a lower potential energy when deposited on an intermetallic core than when deposited on a random alloy core. Also, ordering of Cu atoms at the core@shell interface appears to enhance the overall binding between the core and the shell materials. Inspired by this initial demonstration, SMCR has been used to achieve shells of different random alloy compositions, PtM (M = Ni, Co, Cu, or Fe). This advance is significant because ligand effects vary as a function of PtM identity and Pt/M ratio. These features also influence the degree of surface strain imparted from the lattice mismatch between the core and shell materials. Like the initial demonstration, standout features of these core@shell NPs were high durability and resistance to non-noble metal leaching. Moving forward, efforts have been directed toward integrating shape-control to this core@shell NP construct. This integration is motivated by the shape-dependent catalytic performance of NPs derived from the selective expression of specific facets. Considering the initial i-PdCu@ra-PtCu system, NPs with a cubic shape have been achieved by judicious selection of capping ligands during SMCR. Evaluation of these NPs as catalysts for the electrooxidation of formic acid found that the nanocubic shape enhances catalytic performance compared to similar core@shell NPs with a spherical morphology. We envision that SMCR can be applied to other NP systems to achieve highly durable catalysts as the syntheses of monodisperse and shape-controlled intermetallic seeds are advanced. This Account highlights the role of intermetallic cores in providing more durable electrocatalysts. More broadly, the versatility of SMCR is highlighted as a route to integrate architecture, alloy surfaces, and shape within one NP system, and how this achievement is inspiring new high-performance and robust catalysts is discussed.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Improving the Accuracy of Atomistic Simulations of the Electrochemical Interface
Atomistic simulation of the electrochemical double layer is an ambitious undertaking, requiring quantum mechanical description of electrons, phase space sampling of liquid electrolytes, and equilibration of electrolytes over nanosecond time scales. All models of electrochemistry make different trade-offs in the approximation of electrons and atomic configurations, from the extremes of classical molecular dynamics of a complete interface with point-charge atoms to correlated electronic structure methods of a single electrode configuration with no dynamics or electrolyte. Here, we review the spectrum of simulation techniques suitable for electrochemistry, focusing on the key approximations and accuracy considerations for each technique. Additionally, we discuss promising approaches, such as enhanced sampling techniques for atomic configurations and computationally efficient beyond density functional theory (DFT) electronic methods, that will push electrochemical simulations beyond the present frontier.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Exploring the Role of Inorganic and Organic Interfaces on CO 2 and CH 4 Partitioning: Case Study of Silica, Illite, Calcite, and Kerogen Nanopores on Gas Adsorption and Nanoscale Transport Behaviors
The adsorption, partitioning, and diffusion of CO 2 and CH 4 in organic–inorganic pores composed of kerogen with silica, illite, or calcite are studied using density functional theory and classical molecular dynamics simulations. The adsorption and partitioning behavior of CO 2 and CH 4 molecules is found to be a function of chemistry of the solid interface, pore size, and surface area. CO 2 molecules are preferentially adsorbed compared to CH 4 molecules on silica, illite, calcite, and kerogen surfaces due to significant contributions of electrostatic interactions with the atoms of the associated surfaces. CO 2 and CH 4 molecules display higher adsorption energy on calcite and silica compared to those on illite and kerogen because of enhanced interactions with the positively charged calcium ions on the calcite surface and the hydroxyl group (-OH) on the silica surface. Enhanced internal nanopores and the availability of adsorption active sites in illite and kerogen matrices aid in higher partitioning of the gas molecules into these pores. CO 2 has a significant influence on the clustering and swelling of kerogen fragments compared to CH 4 . Higher CO 2 adsorption on inorganic and kerogen surfaces results in lower overall self-diffusivities compared to CH 4 . This study illustrates surface and pore size controls on the adsorption, partitioning, and self-diffusivities of CO 2 and CH 4 at conditions relevant for subsurface energy applications. The chemistry of the surfaces along with the availability of internal pores and active adsorption sites influences the adsorption behavior of gases. Our results suggest that CO 2 stored in depleted gas reservoirs may preferentially inhabit nanoconfined pores.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Structure and Shape of Surface-Mediated Assembly of Surfactants
Achieving controls on the self-assembly of cationic and block copolymer micelles has significant implications for advancing novel material systems for energy, environmental, and biological applications and for tuning fluid flow behavior in subsurface geologic environments. While it is known that the aggregation of micelles is influenced by the composition of the surrounding fluid, the influence of a solid interface has been less explored. In this study, we probe the organization of cetyltrimethylammonium bromide (CTAB) micelles in the absence and presence of the Pluronic P123 block copolymer and quartz substrate using transmission and grazing-incidence small-angle X-ray scattering measurements and classical molecular dynamics (MD) simulations. In the absence of the quartz interface, CTAB with and without P123 molecules assemble as ellipsoid core–shell micelles. Densification of the core on adding P123 is noted from the decrease in the core radius from 15.9 to 14.6 Å. Here, the presence of a quartz substrate causes the micelles to elongate, which is noted by the emergence of a power-law slope in the low q region (<0.02 Å –1 ). Moreover, in the presence of both P123 and the quartz substrate, the micelle shape changes from ellipsoid core–shell to cylindrical core–shell, and a significantly higher number of aggregates are formed. The higher number of aggregates, and faster aggregation kinetics are linked to the organization of the solvent structures as noted from the MD simulations. These findings demonstrate the effects of adding a block copolymer and including a quartz substrate on the self-assembly of cationic (CTAB) micelles, which can be used to inform the design of functional material systems that harness surfactant chemistries.
42 ENGINEERING↗
Investigating the Accuracy of Water Models through the Van Hove Correlation Function
In this study, we present molecular-simulation-based calculations of the Van Hove correlation function (VHF) of water using multiple modeling approaches: classical molecular dynamics with simple three-site nonpolarizable models, with a polarizable model, and with a reactive force field; density functional tight-binding molecular dynamics; and ab initio molecular dynamics. Due to the many orders of magnitude difference in the computational cost of these approaches, we investigate how small and short the simulations can be while still yielding sufficiently accurate and interpretable results for the VHF. We investigate the accuracy of the different models by comparing them to recently published inelastic X-ray scattering measurements of the VHF. We find that all of the models exhibit qualitative agreement with the experiments, and in some models and for some properties, the agreement is quantitative. This work lays the foundation for future simulation approaches to calculating the VHF for aqueous solutions in bulk and under nanoconfinement.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Capturing Deuteration Effects in a Molecular Mechanics Force Field: Deuterated THF and the THF–Water Miscibility Gap
Deuteration is a common chemical modification used in conjunction with experiments such as neutron scattering, NMR, and Fourier-transform infrared for the study of molecular systems. Under the Born–Oppenheimer (BO) approximation, while the underlying potential energy surface remains unchanged by isotopic substitutions, isotopic substitution still alters intramolecular vibrations, which in turn may alter intermolecular interactions. Molecular mechanics (MM) force fields used in classical molecular dynamics simulations are assumed to represent local approximations of the BO potential energy surfaces, and hence, MD simulations using simple isotopic mass substitutions should capture BO-compatible isotope effects. However, standard MM force-field parameterizations do not directly fit to the local harmonic quantum mechanical (QM) Hessian that describes the BO surface, but rather to QM normal-modes and/or mass-dependent internal-coordinate derived distortion energies. In this study, using tetrahydrofuran (THF)–water mixtures as our model system, we show that not only does a simple mass-substitution approach fail to capture an experimentally characterized deuteration effect (the loss of the closed-loop miscibility gap associated with the complete deuteration of THF) but also it is necessary to generate new MM force-field parameters that correctly describe isotopic dependent vibrations to capture the experimental deuteration effect. We show that the origin of this failure is a result of using mass-dependent features to fit the THF MM force field, which unintentionally biases the bonded terms of the force field to represent only the isotopologue used during the original force-field parameterization. In addition, we make use of our isotopologue-corrected force field for D 8 THF to examine the molecular origins of the isotope-dependent loss of the THF–water miscibility gap.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Computational Predictions of the Hydrolysis of 2,4,6-Trinitrotoluene (TNT) and 2,4-Dinitroanisole (DNAN)
Hydrolysis is a common transformation reaction that can affect the environmental fate of many organic compounds. In this study, three proposed mechanisms of alkaline hydrolysis of 2,4,6-trinitrotoluene (TNT) and 2,4-dinitroaniline (DNAN) were investigated with plane-wave density functional theory (DFT) combined with ab initio and classical molecular dynamics (AIMD/MM) free energy simulations, Gaussian basis set DFT calculations, and correlated molecular orbital theory calculations. Most of the computations in this study were carried out using the Arrows web based tools. For each mechanism: Meisenheimer complex formation, nucleophilic aromatic substitution, and proton abstraction reaction energies and activation barriers were calculated for the reaction at each relevant site. For TNT, it was found that the most kinetically favorable first hydrolysis steps involve Meisenheimer complex formation by attachment of OH$^{-}$ at the C1 and C3 arene carbons and proton abstraction from the methyl group. The nucleophilic aromatic substitution reactions at the C2 and C4 arene carbons were found to be thermodynamically favorable. However, the calculated activation barriers were slightly lower than previous studies, but still found to be $\Delta$G$^{\ddagger} \approx 18$ kcal/mol using PBEo AIMD/MM free energy simulations, suggesting the reactions are not kinetically significant. For DNAN, the barriers of nucleophilic aromatic substitution were even greater ($\Delta$G$^{\ddagger} > 29$ kcal/mol PBEo AIMD/MM). Further, the most favorable hydrolysis reaction for DNAN was found to be a two-step process in which the hydroxyl first attacks the C1 carbon to form a Meisenheimer complex at the C1 arene carbon C1-(OCH$_3$)OH$^{-}$, subsequently, the methoxy anion (-OCH$_3$) at the C1 arene carbon dissociates and the proton shuttles from the C1-OH to the dissociated methoxy group, resulting in methanol and an aryloxy anion.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Using Computationally-Determined Properties for Machine Learning Prediction of Self-Diffusion Coefficients in Pure Liquids
The ability to predict transport properties of liquids quickly and accurately will greatly improve our understanding of fluid properties both in bulk and complex mixtures, as well as in confined environments. Such information could then be used in the design of materials and processes for applications ranging from energy production and storage to manufacturing processes. As a first step, we consider the use of machine learning (ML) methods to predict the diffusion properties of pure liquids. Recent results have shown that Artificial Neural Networks (ANNs) can effectively predict the diffusion of pure compounds based on the use of experimental properties as the model inputs. In the current study, a similar ANN approach is applied to modeling diffusion of pure liquids using fluid properties obtained exclusively from molecular simulations. A diverse set of 102 pure liquids is considered, ranging from small polar molecules (e.g., water) to large nonpolar molecules (e.g., octane). Self-diffusion coefficients were obtained from classical molecular dynamics (MD) simulations. Since nearly all the molecules are organic compounds, a general set of force field parameters for organic molecules was used. The MD methods are validated by comparing physical and thermodynamic properties with experiment. Computational input features for the ANN include physical properties obtained from the MD simulations as well as molecular properties from quantum calculations of individual molecules. Furthermore, fluid properties describing the local liquid structure were obtained from center of mass radial distribution functions (COM-RDFs). Feature sensitivity analysis revealed that isothermal compressibility, heat of vaporization, and the thermal expansion coefficient were the most impactful properties used as input for the ANN model to predict the MD simulated self-diffusion coefficients. The MD-based ANN successfully predicts the MD self-diffusion coefficients with only a subset (2 to 3) of the available computationally determined input features required. A separate ANN model was developed using literature experimental self-diffusion coefficients as model targets. Although this second ML model was not as successful due to a limited number of data points, a good correlation is still observed between experimental and ML predicted self-diffusion coefficients.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Empirically Optimized One-Electron Pseudopotential for the Hydrated Electron: A Proof-of-Concept Study
Mixed quantum-classical molecular dynamics simulations have been important tools for studying the hydrated electron. They generally use a one-electron pseudopotential to describe the interactions of an electron with the water molecules. Furthermore, this approximation shows both the strength and weakness of the approach. On the one hand, it enables extensive statistical sampling and large system sizes that are not possible with more accurate ab initio molecular dynamics methods. On the other hand, there has (justifiably) been much debate about the ability of pseudopotentials to accurately and quantitatively describe the hydrated electron properties. These pseudopotentials have largely been derived by fitting them to ab initio calculations of an electron interacting with a single water molecule. In this paper, we present a proof-of-concept demonstration of an alternative approach in which the pseudopotential parameters are determined by optimizing them to reproduce key experimental properties. Specifically, we develop a new pseudopotential, using the existing TBOpt model as a starting point, which correctly describes the hydrated electron vertical detachment energy and radius of gyration. In addition to these properties, this empirically optimized model displays a significantly modified solvation structure, which improves, for example, the prediction of the partial molar volume.
38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
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↗