SEARCH · Engineering Papers
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.
Stimuli-responsive “solvate-sponge”-(DMF) 3 NaClO 4 exhibits linear chains of DMF–Na + ions with ClO 4 − anions in the interstitial space. At increased pressure or temperature, DMF is expelled (reversibly), resulting in a new stoichiometry-(DMF) 2 NaClO 4 .
Evolution of a lithium cation solvation complex in low concentration electrolytes (LCEs), high concentration electrolytes (HCEs), and localized high concentration electrolytes (LHCEs) reveals competition of transport, desolvation, and deposition events.
Correction for ‘The solvation structure, transport properties and reduction behavior of carbonate-based electrolytes of lithium-ion batteries' by Tingzheng Hou et al. , Chem. Sci. , 2021, 12 , 14740–14751, https://doi.org/10.1039/D1SC04265C.
Empirical measurements of solution vapor pressure of ternary acetonitrile (MeCN) H 2 O–NaCl–MeCN mixtures were recorded, with NaCl concentrations ranging from zero to the saturation limit, and MeCN concentrations ranging from zero to an absolute mole fraction of 0.64. After accounting for speciation, the variability of the Henry's law coefficient at vapor–liquid equilibrium (VLE) of MeCN ternary mixtures decreased from 107% to 5.1%. Solute speciation was modeled using a mass action solution model that incorporates solute solvation and ion-pairing phenomena. Two empirically determined equilibrium constants corresponding to solute dissociation and ion pairing were utilized for each solute. When speciation effects were considered, the solid–liquid equilibrium of H 2 O–NaCl–MeCN mixtures appear to be governed by a simple saturation equilibrium constant that is consistent with the binary H 2 O–NaCl saturation coefficient. Further, our results indicate that the precipitation of NaCl in the MeCN ternary mixtures was not governed by changes in the dielectric constant. Our model indicates that the compositions of the salt-induced liquid–liquid equilibrium (LLE) boundary of the H 2 O–NaCl–MeCN mixture correspond to the binary plateau activity of MeCN, a range of concentrations over which the activity remains largely invariant in the binary water–MeCN system. Broader comparisons with other ternary miscible organic solvent (MOS) mixtures suggest that salt-induced liquid–liquid equilibrium exists if: (1) the solution displays a positive deviation from the ideal limits governed by Raoult's law; and (2) the minimum of the mixing free energy profile for the binary water-MOS system is organic-rich. Finally, this work is one of the first applications of speciation-based solution models to a ternary system, and the first that includes an organic solute.
To scrutinise the impact of electric fields on the structure and vibrations of biomolecules in the presence of water, we study the sequential solvation of lithium diglycine up to three water molecules with cryogenic infrared action spectroscopy.
Based on THz-spectra we quantify hydrophilic solvation and correlate with mixing enthalpy.
Tuning the solvation structure of lithium ions via electrolyte engineering has proven effective for lithium metal (Li) anodes. Further advancement that bypasses the trial-and-error practice relies on the establishment of molecular design principles. Expanding the scope of our previous work on solvent fluorination, we report here an alternative design principle for non-fluorinated solvents, which potentially have reduced cost, environmental impact, and toxicity. By studying non-fluorinated ethers systematically, we found that the short-chain acetals favor the [gauche, gauche] molecular conformation due to hyperconjugation, which leads to weakened monodentate coordination with Li + . The dimethoxymethane electrolyte showed fast activation to >99% coulombic efficiency (CE) and high ionic conductivity of 8.03 mS cm -1 . The electrolyte performance was demonstrated in anode-free Cu$∥$LFP pouch cells at current densities up to 4 mA cm -2 (70 to 100 cycles) and thin-Li$∥$high-loading-LFP coin cells (200–300 cycles). Overall, we demonstrated and rationalized the improvement in Li metal cyclability by the acetal structure compared to ethylene glycol ethers. We expect further improvement in performance by tuning the acetal structure.
Virtual agents and cloud computing have enabled chemists to easily access automated simulations of explicitly solvated molecules.
This study demonstrates comprehensive approaches for predicting solvation structures in Ca-ion battery electrolytes by integrating ab initio calculations and machine learning force fields.
The combination of machine learning (ML) models with chemistry-related tasks requires the description of molecular structures in a machine-readable way. The nature of these so-called molecular descriptors has a direct and major impact on the performance of ML models and remains an open problem in the field. Structural descriptors like SMILES strings or molecular graphs lack size-independence and can be memory intensive. Machine-learned descriptors can be of low dimensionality and constant size but lack physical significance and human interpretability. Sigma profiles, which are unnormalized histograms of the surface charge distributions of solvated molecules, combine physical significance with low dimensionality and size-independence, making them a suitable candidate for a universal molecular descriptor. However, their widespread adoption in ML applications requires open access to sigma profile generation, which is currently not available. This work details the development of OpenSPGen – an open-source tool for generating sigma profiles. Also presented are studies on the effect of different settings on the efficacy of the generated sigma profiles at predicting thermophysical material properties when used as inputs to a Gaussian process as a simple surrogate ML model. We find that a higher level of theory does not translate to more accurate results. We also provide further recommendations for sigma profile calculation and use in ML models.
Microscopic linear liquid free-streams offer a means of gently delivering biological samples into a probe beam in vacuum while maintaining the sample species in a fully solvated state. By employing gas dynamic forces to form the microscopic liquid stream (as opposed to a conventional solid-walled convergent nozzle), liquid free-streams down to 300 nm diameter have been generated. Such 'Gas Dynamic Virtual Nozzles' (GDVN) are ideally suited to injecting complex biological species into an X-ray Free Electron Laser (XFEL) to determine the structure of the biological species via Serial Femtosecond Crystallography (SFX). GDVN injector technology developed for this purpose is described.
Infrared spectra in the C–H stretch region are reported for the allyl (CH 2 CHCH 2 ) and allyl peroxy (CH 2 =CH–CH 2 OO·) radicals solvated in superfluid helium nanodroplets. Nine bands in the spectrum of the allyl radical have resolved rotational substructure. We have assigned three of these to the ν 1 (a 1 ), ν 3 (a 1 ), and ν 13 (b 2 ) C–H stretch bands and four others to the ν 14 /(ν 15 +2ν 11 ) (b 2 ) and ν 2 /(ν 4 +2ν 11 ) (a 1 ) Fermi dyads, and an unassigned resonant polyad is observed in the vicinity of the ν 1 band. Experimental coupling constants associated with Fermi dyads are consistent with quartic force constants obtained from density functional theory computations. The peroxy radical was formed within the He droplet via the reaction between allyl and O 2 following the sequential pick-up of the reactants. Five stable conformers are predicted for the allyl peroxy radical, and a computed two-dimensional potential surface for rotation about the CC–OO and CC–CO bonds reveals multiple isomerization barriers greater than ≈300 cm –1 . Furthermore, the C–H stretch infrared spectrum is consistent with the presence of a single conformer following the allyl + O 2 reaction within helium droplets.
The aqueous hydrogen molecule is studied with molecular dynamics simulations at ambient temperature and pressure conditions, using a newly developed flexible and polarizable H{sub 2} molecule model. The design and implementation of this model, compatible with an existing flexible and polarizable force field for water, is presented in detail. The structure of the hydration layer suggests that first-shell water molecules accommodate the H{sub 2} molecule without major structural distortions and two-dimensional, radial-angular distribution functions indicate that as opposed to strictly tangential, the orientation of these water molecules is such that the solute is solvated with one of the free electron pairs of H{sub 2}O. The calculated self-diffusion coefficient of H{sub 2}(aq) agrees very well with experimental results and the time dependence of mean square displacement suggests the presence of caging on a time scale corresponding to hydrogen bond network vibrations in liquid water. Orientational correlation function of H{sub 2} experiences an extremely short-scale decay, making the H{sub 2}–H{sub 2}O interaction potential essentially isotropic by virtue of rotational averaging. The inclusion of explicit polarizability in the model allows for the calculation of Raman spectra that agree very well with available experimental data on H{sub 2}(aq) under differing pressure conditions, including accurate reproduction of the experimentally noted trends with solute pressure or concentration.
In this work we achieve three milestones: (1) we present a subsystem DFT method capable of running ab-initio molecular dynamics simulations accurately and efficiently. (2) In order to rid the simulations of inter-molecular self-interaction error, we exploit the ability of semilocal frozen density embedding formulation of subsystem DFT to represent the total electron density as a sum of localized subsystem electron densities that are constrained to integrate to a preset, constant number of electrons; the success of the method relies on the fact that employed semilocal nonadditive kinetic energy functionals effectively cancel out errors in semilocal exchange–correlation potentials that are linked to static correlation effects and self-interaction. (3) We demonstrate this concept by simulating liquid water and solvated OH{sup •} radical. While the bulk of our simulations have been performed on a periodic box containing 64 independent water molecules for 52 ps, we also simulated a box containing 256 water molecules for 22 ps. The results show that, provided one employs an accurate nonadditive kinetic energy functional, the dynamics of liquid water and OH{sup •} radical are in semiquantitative agreement with experimental results or higher-level electronic structure calculations. Our assessments are based upon comparisons of radial and angular distribution functions as well as the diffusion coefficient of the liquid.
A dielectric response theory of solvation beyond the conventional Born model for polar fluids is presented. The dielectric response of a polar fluid is described by a Born response mode and a linear combination of Debye-Hückel-like response modes that capture the nonlocal response of polar fluids. The Born mode is characterized by a bulk dielectric constant, while a Debye-Hückel mode is characterized by its corresponding Debye screening length. Both the bulk dielectric constant and the Debye screening lengths are determined from the bulk dielectric function of the polar fluid. The linear combination coefficients of the response modes are evaluated in a self-consistent way and can be used to evaluate the electrostatic contribution to the thermodynamic properties of a polar fluid. In conclusion, our theory is applied to a dipolar hard sphere fluid as well as interaction site models of polar fluids such as water, where the electrostatic contribution to their thermodynamic properties can be obtained accurately.
Lignin as a potential renewable source of biofuels, chemicals, and other value-added products has gained much attention. However, the complexity of lignin structure poses a significant challenge for developing efficient valorization techniques. As most processes involve solvothermal conditions to minimize energy cost, lignin depolymerization is governed by reaction conditions (temperature and pressure), and solvents. In this work, binding of β -O-4 linkage consisting lignin dimers on MWW 2-dimensional (2D) zeolite is investigated using periodic density functional theory (DFT). Furthermore, the effect of different terminated surfaces (H:OH % = 100:0; 50:50; 0:100 %), different temperatures (323, 353, 373 K), and different solvents (water and methanol) on the binding modes is quantified. Here, our work shows that in the gas phase the binding strength increases 10 to 15 kcal/mol upon increasing the number of hydroxyl groups on the surface. Also, the phenolic dimer binds more strongly than the non-phenolic dimer, and the binding strength of model compounds increases in the presence of solvent. Analysis of structural changes in the presence of the solvent reveals that aromatic rings are parallel to the zeolite surface and primary interaction with zeolite is through the hydroxyl groups near the β -O-4 linkage. Furthermore, while the solvation energy decreases with increasing temperature the opposite trend is observed for the binding energy with the surface.
Proper statistical mechanics understanding of nanoparticle solvation processes requires an accurate description of the molecular structure of the solvent. Achieving this goal with standard molecular dynamics (MD) simulation methods is challenging due to large length scales. An alternative approach to this problem can be formulated using classical density functional theory (cDFT), where a full configurational description of the positions of all the atoms is replaced by collective atomic site densities in the molecule. Using an example of the negatively charged silica-like system in an aqueous polar environment represented by a two-site water model, we demonstrate here that cDFT can reproduce MD data at a fraction of the computational cost. An important implication of this result is the ability to understand how the solvent molecular features may affect the system’s properties at the macroscopic scale. A concrete example highlighted in this work is the analysis of nanoparticle interactions with sizes of up to 100 nm in diameter.
Heparan sulfates are found in all animal tissues and have essential roles in living systems. This family of biomacromolecules modulates binding to calcium ions (Ca 2+ ) in low free energy reactions that influence biochemical processes from cell signaling and anticoagulant efficacy to biomineralization. Despite their ubiquity, the thermodynamic basis for how heparans and similarly functionalized biomolecules regulate Ca 2+ interactions is not yet established. Using heparosan (Control) and heparins with different positions of sulfate groups, we quantify how SO 3 − and COO − content and SO 3 − position modulate Ca 2+ binding by isothermal titration calorimetry. The free energy of all heparin-Ca2+ interactions (ΔG rxn ) is dominated by entropic contributions due to favorable water release from polar, hydrophilic groups. Heparin with both sulfate esters (O-SO 3 − ) and sulfamides (N-SO 3 − ) has the strongest binding to Ca 2+ compared to heparosan and to heparin with only O-SO 3 − groups (~3X). By linking Ca 2+ binding thermodynamics to measurements of the interfacial energy for calcite (CaCO 3 ) crystallization onto polysaccharides, we show molecule-specific differences in nucleation rate can be explained by differences in water structuring during Ca 2+ interactions. A large entropic term (-TΔS rxn ) upon Ca 2+ –polysaccharide binding correlates with high interfacial energy to CaCO 3 nucleation. Combining our measurements with literature values indicates many Ca 2+ –polysaccharide interactions have a shared thermodynamic signature. The resulting enthalpy–entropy compensation relationship suggests these interactions are generally dominated by water restructuring involving few conformational changes, distinct from Ca 2+ –protein binding. Our findings quantify the thermodynamic origins of heparin-specific interactions with Ca 2+ and demonstrate the contributions of solvation and functional group position during biomacromolecule-mediated ion regulation.