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 217 records · Page 12

Deciphering the Solvation Structure of Aqueous ZnCl 2 Solutions from X-ray Absorption Spectra Using the Interpretable Graph Neural Network

Machine learning (ML) provides powerful pathways for predicting spectroscopic observables from atomic structures, but its broader impact depends on making model predictions interpretable in terms of physical and chemical principles. Here, we introduce a physics-guided graph neural network (GNN) model that predicts Zn K-edge X-ray spectroscopy (XAS) spectra of aqueous ZnCl 2 solutions. Training data are generated from ab initio XAS calculations on molecular dynamics snapshots obtained using a machine learning interatomic potential. The GNN reproduces experimental spectra across concentrations from dilute (<0.1 m) to highly concentrated (30 m, “water-in-salt”) regimes and scales efficiently to large, disordered liquid systems beyond the reach of conventional ab initio approaches. Gradient-based attribution analysis reveals that the model learns physically meaningful structure-spectrum relationships. Ligand-specific attributions reflect orbital hybridization patterns and the origin of the excitations derived from the density functional theory. Bond-length attributions recover spectral shifts consistent with multiple-scattering theory. Finally, this work bridges data-driven prediction with electronic-structure theory, establishing a general paradigm for interpretable ML that links atomic structure, electronic structure, and spectroscopic observables.

25 ENERGY STORAGE↗

A redox model for NO oxidation, NH 3 oxidation and high temperature standard SCR over Cu-SSZ-13

A kinetic model is developed to predict the influence of temperature and hydrothermal aging on the redox of active Cu sites under standard SCR, NO oxidation and NH 3 oxidation conditions over a practically relevant fully-formulated Cu-SSZ-13 catalyst. NO 2 /N 2 production during NO/NH 3 titration of Cu II sites is utilized to identify rate parameters associated with NO-only RHC (reduction half cycle) and NH 3 -only RHC respectively. Integral N 2 formation during subsequent NO + NH 3 titration is consistent with the production of one NO 2 per two CuII sites reduced during NO-only RHC and one N 2 per six Cu II sites reduced during NH 3 -only RHC. Decreased reduction of Cu II sites by NO/NH 3 upon hydrothermal aging, along with the production of one NO 2 per two CuII sites during NO-only RHC, is accordant with the involvement of proximal ZCuOH and oxygen-bridged dimeric Cu II sites. Oxidation of partially solvated and framework coordinated ZCu (Cu I ) sites occurs in presence of O 2 , does not produce N 2 and can lead to the consumption of Brønsted acid sites. A global OHC kinetic model is developed to predict transient and integral N 2 formation during exposure of CuI sites to a mixture of NO and O 2 . The resulting redox kinetic model quantitatively predicts NO and NH 3 consumption during isothermal transient response Cu redox (TRCR) protocols, along with temperature and age dependent steady-state standard SCR and oxidation conditions. The redox model presented in this work synthesizes recent kinetic, spectroscopic and computational findings to provide a foundational description of active site redox during standard SCR, NO oxidation and NH3 oxidation over Cu-SSZ-13.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Microscopic insights into the solvation of polyethylene glycol chains in water: A machine learning potential approach

Polyethylene glycol (PEG) is a structurally simple, nontoxic, and water-soluble polymer widely utilized in medical and pharmaceutical applications. Notably, when a PEG chain is immersed in water, the surrounding water molecules play a key role in driving conformational changes of this macromolecule. In this study, we explore the solvation behavior of PEG under mechanical strain using molecular dynamics simulations, with an interatomic potential obtained from machine learning. Our focus is on the transition from the favored coil-like conformation to an extended one under external force. Through analyses of radial distribution functions, hydrogen bonding, and solvation dynamics, we uncover how mechanical stretching influences the local hydration environment. Furthermore, we disentangle the enthalpic and entropic contributions to the conformational stability of PEG in water. Surprisingly, our neural network potential model identifies dewetting of PEG C-atoms, and not water H-bonding with PEG O-atoms, as the main enthalpic driving force for the coiling of PEG in water.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Nuclear–Electronic Orbital QM/MM Approach: Geometry Optimizations and Molecular Dynamics

Hybrid quantum mechanical/molecular mechanical (QM/MM) methods allow simulations of chemical reactions in atomistic solvent and heterogeneous environments such as proteins. Herein, the nuclear–electronic orbital (NEO) QM/MM approach is introduced to enable the quantization of specified nuclei, typically protons, in the QM region using a method such as NEO-density functional theory (NEO-DFT). This approach includes proton delocalization, polarization, anharmonicity, and zero-point energy in geometry optimizations and dynamics. Expressions for the energies and analytical gradients associated with the NEO-QM/MM method, as well as the previously developed polarizable continuum model (NEO-PCM), are provided. Geometry optimizations of small organic molecules hydrogen bonded to water in either dielectric continuum solvent or explicit atomistic solvent illustrate that aqueous solvation can strengthen hydrogen-bonding interactions for the systems studied, as indicated by shorter intermolecular distances at the hydrogen-bond interface. We then performed a real-time direct dynamics simulation of a phenol molecule in explicit water using the NEO-QM/MM method. Furthermore, these developments and initial examples provide the foundation for future studies of nuclear–electronic quantum dynamics in complex chemical and biological environments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Zn 2+ -mediated catalysis for fast-charging aqueous Zn-ion batteries

Rechargeable aqueous zinc-ion batteries (AZIBs), renowned for their safety, high energy density and rapid charging, are prime choices for grid-scale energy storage. Historically, ion-shuttling models centring on ion-migration behaviour have dominated explanations for charge/discharge processes in aqueous batteries, like classical ion insertion/extraction and pseudocapacitance mechanisms. However, these models struggle to account for the exceptional performance of AZIBs compared to other aqueous metal-ion batteries. Here, in this study, we present a catalysis model elucidating the Zn 2+ anomaly in aqueous batteries, explaining it through the concept of adsorption in catalysis. Such behaviour can serve the charge/discharge role, predominantly dictated by solvated metal cations and cathode materials. First-principles calculations suggest optimal adsorption/desorption behaviour (water dissociation process) with the Zn 2+ -vanadium nitride (VN) combination. Experimentally, AZIBs implementing VN cathodes demonstrate fast-charging kinetics, showing a capacity of 577.1 mAh g -1 at a current density of 300,000 mA g -1 . The grasp of catalysis steps within AZIBs can drive solutions beyond state-of-the-art fast-charging batteries.

25 ENERGY STORAGE↗

Unexpected electrochemical behavior of an anolyte redoxmer in flow battery electrolytes: solvating cations help to fight against the thermodynamic–kinetic dilemma

Redoxmers are redox-active molecules that can store energy in electrolytes for redox flow batteries (RFBs), and their electrochemical properties are significantly affected by the choice of supporting electrolytes. Herein, we use 2,1,3-benzothiadiazole (BzNSN) as a model system to scrutinize the supporting electrolyte impact. By systemically varying the components of supporting salts, BzNSN not only shows substantial redox potential shifts but also exhibits varying electrochemical stabilities. Specifically, changing the size of cations can effectively alter the coordination between the supporting salt and BzNSN species. From Li + , Na + , K + , to NEt 4 + , the redox potential of BzNSN shifts negatively, from -1.63 V to -1.82 V vs. Ag/Ag + . Molecular dynamics and density functional theory simulations revealed that smaller cations, like Li + , are closer to the charged BzNSN when coordinated, implying stronger coordination, while larger cations, like K + and NEt 4 + , are farther away. Interestingly, the large cation electrolytes also lead to much improved electrochemical stability, evidenced by the extraordinarily enhanced kinetic lifetime from electron paramagnetic resonance measurement. This study demonstrates the first example of tuning an anolyte redoxmer toward a concurrent improvement of lowered redox potentials AND enhanced calendar lives via solvation means, which is usually constrained by the thermodynamic–kinetic relation.

2,1,3-benzothiadiazole↗

Effect of water models on structure and dynamics of lignin in solution

Lignin, a major biomass component, can be an excellent source for different monomers in the polymer industry. However, the complex and heterogeneous structure of lignin poses a significant challenge for designing energy-efficient processes for depolymerization. As many proposed depolymerization processes are solvothermal, it is essential to understand the structure and dynamics of lignin in solution, in particular aqueous solution. Here, we utilize molecular dynamics simulations to understand the effect of water models on the structure and dynamics of different model lignin oligomers (softwood and hardwood) as a function of temperature. We have examined three different water models: TIP3P, TIP4P/Ew, and flexible SPC/Fw. We find that the diffusion constant of lignin oligomers in an aqueous solution differs significantly depending on the water model used. The diffusion constant of lignin in the TIP3P water model is almost twice as large as that in SPC/Fw and TIP4P/Ew. The softwood and hardwood oligomers adopt an extended structure in TIP3P water compared to SPC/Fw and TIP4P/Ew. Given the different levels of sensitivity of transport and structural properties of aqueous lignin on water models, it is important to take these into account when discussing results from a specific water model.

09 BIOMASS FUELS↗

Influence of Charge Fraction on the Phase Behavior of Symmetric Single-Ion Conducting Diblock Copolymers

A series of symmetric poly[(oligo(ethylene glycol) methyl ether methacrylate-co-oligo(ethylene glycol) propyl sodium sulfonate methacrylate)]-block-polystyrene (PsOEGMA-PS) diblock copolymers were synthesized as a model system to probe the effect of charge fraction on the phase behavior of charged-neutral single-ion conducting diblock copolymers. Small-angle X-ray scattering (SAXS) experiments showed that increasing the charge fraction does not alter the ordered phase morphology (lamellar) but increases the order–disorder transition temperature (T ODT ) significantly. Additionally, the effective Flory–Huggins interaction parameter (χ eff ) was found to increase linearly with the charge fraction, similar to the case of conventional salt-doped diblock copolymers. This indicates that the effect of counterion solvation, attributed to the significant mismatch between the dielectric constant of each block, provides the dominant effect in tuning the phase behavior of this charged diblock copolymer. Lastly, we therefore infer that electrostatic cohesion (local charge ordering induced by Coulombic interactions), which is predicted to suppress microphase separation and lead to asymmetric phase diagrams, only plays a minor role in this model system.

36 MATERIALS SCIENCE↗

Machine Learning-Driven Solvent Screening for Biobased 2,3-Butanediol Extraction

Biobased 2,3-butanediol (2,3-BDO) is a valuable biomass-derived chemical due to its versatility in being transformed into a wide variety of products. However, the separation and purification of 2,3-BDO from fermentation broth remain a significant challenge owing to its high boiling point and hydrophilic nature. Herein, we developed a machine learning (ML)-based screening workflow that uses molecular calculations as training data and requires only a small number of experimental measurements for validation to identify alternative solvent candidates for the liquid–liquid extraction (LLE) of 2,3-BDO from aqueous solution. In particular, 130 density functional theory (DFT) calculations with the implicit solvation method not only built a correlation between the computational partition coefficient and the experimental distribution coefficient of 2,3-BDO but also parameterized an Extra-Trees ML model to screen the distribution coefficient for a wider range of 6717 organic solvents. The experimental measurements of only 24 solvents were needed to validate the computational results. A list of 50 prioritized solvents was proposed for 2,3-BDO LLE, and seven additional experimental measurements were conducted to further verify our selected solvents. The impact of the extraction temperature and solvent-to-feed ratio was also investigated for selected solvents in experiments. Furthermore, this work suggested alternative solvents for 2,3-BDO LLE and proposed a versatile workflow that requires fewer experiments and can be applied to a broader range of LLE studies.

Extraction↗

Minimal Experimental Bias on the Hydrogen Bond Greatly Improves Ab Initio Molecular Dynamics Simulations of Water

Experiment Directed Simulations (EDS) is a method within a class of techniques seeking to improve molecular simulations by minimally biasing the system Hamiltonian to reproduce certain experimental observables. In a previous application of EDS to ab initio molecular dynamics (AIMD) simulation based on electronic density functional theory (DFT), the AIMD simulations of water were biased to reproduce its experimentally derived solvation structure. In particular, by solely biasing the O-O pair correlation functions, other structural and dynamical properties that were not biased were improved. In this work, the hypothesis is tested that directly biasing the OH pair correlation (and hence the H-O∙∙∙H hydrogen bonding), will provide an even better improvement of DFT-based water properties in AIMD simulations. The logic behind this hypothesis is that for most electronic DFT descriptions of water the hydrogen bonding is known to be deficient due to anomalous charge transfer and over polarization in the DFT. Using recent advances to the EDS learning algorithm, we thus train a minimal bias on AIMD water that reproduces the O-H radial distribution function derived from the highly accurate MB-pol model of water. Finally, it is then confirmed that biasing the O-H pair correlation alone can lead to improved AIMD water properties, with structural and dynamical properties in even closer to experiment than the previous EDS-AIMD model.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Semiconductor Nanostructures by Scientific Design (Final Technical Report)

We developed and used quantum simulation techniques to predict structural, electronic, optical and thermal properties affecting the scientific design of semiconductor nanostructures. We focused on the study of surfaces and interfaces at the nanoscale, in particular embedding and solvation of semiconductor nanoparticles, and provided a microscopic description of nanostructured materials in realistic environments, directly comparable with experimental conditions. We addressed three major problems: (i) the design of realistic structural models of integrated nanostructures, e.g. quantum dots embedded in solid matrices and nanostructures immersed in fluids; (ii) the development of efficient and accurate calculations to describe the excited state properties of nanostructures, beyond standard Density Functional Theory, so as to understand and predict electronic properties involved in renewable energy applications, for example solar conversion devices. In particular we developed and used methods based on a many body perturbation theory. (iii) the development of calculations of heat transport properties of nanostructured materials, for thermoelectric applications, including predictive calculations of thermal and electrical conductivity and of Seebeck coefficients. Specifically, we developed and used methods based on the Boltzmann transport equation and molecular dynamics.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Probing the self-ionization of liquid water with ab initio deep potential molecular dynamics

The chemical equilibrium between self-ionized and molecular water dictates the acid–base chemistry in aqueous solutions, yet understanding the microscopic mechanisms of water self-ionization remains experimentally and computationally challenging. Herein, Density Functional Theory (DFT)–based deep neural network (DNN) potentials are combined with enhanced sampling techniques and a global acid–base collective variable to perform extensive atomistic simulations of water self-ionization for model systems of increasing size. The explicit inclusion of long-range electrostatic interactions in the DNN potential is found to be crucial to accurately reproduce the DFT free energy profile of solvated water ion pairs in small (64 and 128 H 2 O) cells. The reversible work to separate the hydroxide and hydronium to a distance S is found to converge for simulation cells containing more than 500 H 2 O, and a distance of ∼ 8 Å is the threshold beyond which the work to further separate the two ions becomes approximately zero. The slow convergence of the potential of mean force with system size is related to a restructuring of water and an increase of the local order around the water ions. Calculation of the dissociation equilibrium constant illustrates the key role of long-range electrostatics and entropic effects in the water autoionization process.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Acid-Base Equilibrium and Dielectric Environment Regulate Charge in Supramolecular Nanofibers

Peptide amphiphiles are a class of molecules that can self-assemble into a variety of supramolecular structures, including high-aspect-ratio nanofibers. It is challenging to model and predict the charges in these supramolecular nanofibers because the ionization state of the peptides are not fixed but liable to change due to the acid-base equilibrium that is coupled to the structural organization of the peptide amphiphile molecules. Here, we have developed a theoretical model to describe and predict the amount of charge found on self-assembled peptide amphiphiles as a function of pH and ion concentration. In particular, we computed the amount of charge of peptide amphiphiles nanofibers with the sequence C 16 − V 2 A 2 E 2 . In our theoretical formulation, we consider charge regulation of the carboxylic acid groups, which involves the acid-base chemical equilibrium of the glutamic acid residues and the possibility of ion condensation. The charge regulation is coupled with the local dielectric environment by allowing for a varying dielectric constant that also includes a position-dependent electrostatic solvation energy for the charged species. We find that the charges on the glutamic acid residues of the peptide amphiphile nanofiber are much lower than the same functional group in aqueous solution. There is a strong coupling between the charging via the acid-base equilibrium and the local dielectric environment. Our model predicts a much lower degree of deprotonation for a position-dependent relative dielectric constant compared to a constant dielectric background. Furthermore, the shape and size of the electrostatic potential as well as the counterion distribution are quantitatively and qualitatively different. These results indicate that an accurate model of peptide amphiphile self-assembly must take into account charge regulation of acidic groups through acid–base equilibria and ion condensation, as well as coupling to the local dielectric environment.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Computational studies of Ras and PI3K

Until recently, experimental techniques in molecular cell biology have been the primary means to investigate biological risk upon space radiation. However, computational modeling provides an alternative theoretical approach, which utilizes various computational tools to simulate proteins, nucleotides, and their interactions. In this study, we are focused on using molecular mechanics (MM) and molecular dynamics (MD) to study the mechanism of protein-protein binding and to estimate the binding free energy between proteins. Ras is a key element in a variety of cell processes, and its activation of phosphoinositide 3-kinase (PI3K) is important for survival of transformed cells. Different computational approaches for this particular study are presented to calculate the solvation energies and binding free energies of H-Ras and PI3K. The goal of this study is to establish computational methods to investigate the roles of different proteins played in the cellular responses to space radiation, including modification of protein function through gene mutation, and to support the studies in molecular cell biology and theoretical kinetics models for our risk assessment project.

Ren, Lei↗

Exploring the impact of ions on oxygen K-edge X-ray absorption spectroscopy in NaCl solution using the GW-Bethe-Salpeter-equation approach

X-ray absorption spectroscopy (XAS) is a powerful experimental tool to probe the local structure in materials with the core hole excitations. Here, the oxygen K-edge XAS spectra of the NaCl solution and pure water are computed by using a recently developed GW-Bethe-Salpeter equation approach, based on configurations modeled by path-integral molecular dynamics with the deep-learning technique. The neural network is trained on ab initio data obtained with strongly constrained and appropriately normed density functional theory. The observed changes in the XAS features of the NaCl solution, compared to those of pure water, are in good agreement between experimental and theoretical results. We provided detailed explanations for these spectral changes that occur when NaCl is solvated in pure water. Specifically, the presence of solvating ion pairs leads to localization of electron-hole excitons. As a result, our theoretical XAS results support the theory that the effects of the solvating ions on the H-bond network are mainly confined within the first hydration shell of ions, however beyond the shell the arrangement of water molecules remains to be comparable to that observed in pure water.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Unlocking enhanced gas capture via core scrambling of porous-organic cages

The demand for low-cost, low-energy, and highly selective gas capture and separations is an ongoing driver of porous material development. Porous liquids have been identified as a promising gas separation material by creating permanent porosity in inorganic solvents through inclusion of nanoporous materials that sterically exclude solvent from their internal porosity. Among the nanoporous materials that can be used to form porous liquids, porous-organic cages (POCs) have been one of the most popular due to the inherent tunability of POCs. “Scrambled” POCs with varying functionalities on the POC vertices have been developed and incorporated into porous liquid compositions, increasing their gas adsorption capacity. An unexplored avenue to tailor the properties of porous liquids is through scrambling the functionality of the core of the POC. Here, therefore, we have synthesized a new POC, a CC3-OH derivative with scrambled hydroxides on the core and evaluated the impact on the CO 2 uptake capacity in silicon oil-based porous liquids. Core scrambling of the POC resulted in a twofold increase CO 2 adsorption capacity in the porous liquid, an emergent property that is a dramatic increase beyond a linear combination of the gas adsorption capacity of the neat solvent and the POC. Density functional theory modeling of the CC3 POC and its hydroxide-based derivatives identified that free rotation of the linker hydroxide allowed for forced interaction between the CO 2 molecule and the hydroxide in the pore window. Solvation of the POC may release scrambled core hydroxides from intramolecular bonding with a neighboring imine, allowing for increased gas uptake in the porous liquid over the neat POC. These results identify a key structural relationship of POCs that enables emergent properties in porous liquids and can guide future development of liquid phase gas capture and separation materials for environmental and industrial applications.

Gas capture↗

Binding of Phosphate Species to Ca 2+ and Mg 2+ in Aqueous Solution

Phosphate derivatives and their interaction with metal cations are involved in many important biological phenomena, so an accurate characterization of the phosphate-metal interaction is necessary to properly understand the role of phosphate-metal contacts in mediating biological function. Herein, we improved the standard 12-6 Lennard-Jones (LJ) potential via the usage of the 12-6-4 LJ model, which incorporates ion-induced dipole interactions. Via parameter scanning, we fine-tuned the 12-6-4 LJ polarizability values to obtain accurate absolute binding free energies for the phosphate anions H 2 PO 4 – , HPO 4 2– , PO 4 3– coordinating with Ca 2+ and Mg 2+ . First, we modified the phosphate 12-6-4 LJ parameters to reproduce the solvation free energies of the series of phosphate anions using the thermodynamic integration (TI) method. Then, using the potential mean force (PMF) method, the polarizability of the metal-phosphate interaction was obtained. We show that the free energy profiles of phosphate ions coordinated to Ca 2+ and Mg 2+ generally show similar trends at longer metal–phosphate distances, while the absolute binding energy values increased with deprotonation. The resulting parameters demonstrate the flexibility of the 12-6-4 LJ-type nonbonded model and its usefulness in accurately describing cation–anion interactions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Determining Partial Atomic Charges for Liquid Water: Assessing Electronic Structure and Charge Models

Partial atomic charges provide an intuitive and efficient way to describe the charge distribution and the resulting intermolecular electrostatic interactions in liquid water. Many charge models exist and it is unclear which model provides the best assignment of partial atomic charges in response to the local molecular environment. In this work, we systematically scrutinize various electronic structure methods and charge models (Mulliken, natural population analysis, CHelpG, RESP, Hirshfeld, Iterative Hirshfeld, and Bader) by evaluating their performance in predicting the dipole moments of isolated water, water clusters, and liquid water as well as charge transfer in the water dimer and liquid water. Although none of the seven charge models is capable of fully capturing the dipole moment increase from isolated water (1.85 D) to liquid water (about 2.9 D), the Iterative Hirshfeld method performs best for liquid water, reproducing its experimental average molecular dipole moment, yielding a reasonable amount of intermolecular charge transfer, and showing modest sensitivity to the local water environment. The performance of the charge model is dependent on the choice of the density functional and the quantum treatment of the environment. The computed molecular dipole moment of water generally increases with the percentage of the exact Hartree–Fock exchange in the functional, whereas the amount of charge transfer between molecules decreases. For liquid water, including two full solvation shells of surrounding water molecules (within about 5.5 Å of the central water) in the quantum chemical calculation converges the charges of the central water molecule. Furthermore, our final pragmatic quantum chemical charge-assigning protocol for liquid water is the Iterative Hirshfeld method with M06-HF/aug-cc-pVDZ and a quantum region cutoff radius of 5.5 Å.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗