Multiple linear regression and thermodynamic fluctuations are equivalent for computing thermodynamic derivatives from molecular simulation
Not Available
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.
Not Available
Dendrite formation during electrodeposition while charging lithium metal batteries compromises their safety. Although high-shear-modulus ( G s ) solid-ion conductors (SICs) have been prioritized to resolve the pressure-driven instabilities that lead to dendrite propagation and cell shorting, it is unclear whether these or alternatives are needed to guide uniform lithium electrodeposition, which is intrinsically density-driven. Here, we show that SICs can be designed within a universal chemomechanical paradigm to access either pressure-driven dendrite-blocking or density-driven dendrite-suppressing properties, but not both. This dichotomy reflects the competing influence of the SIC’s mechanical properties and the partial molar volume of Li + $(V_{\mathrm{Li}^+})$ relative to those of the lithium anode ( G Li and V Li ) on plating outcomes. Within this paradigm, we explore SICs in a previously unrecognized dendrite-suppressing regime that are concomitantly ‘soft’, as is typical of polymer electrolytes, but feature an atypically low $(V_{\mathrm{Li}^+})$ that is more reminiscent of ‘hard’ ceramics. Li plating (1 mA cm -2 ; T = 20 °C) mediated by these SICs is uniform, as revealed using synchrotron hard X-ray microtomography. As a result, cell cycle life is extended, even when assembled with thin Li anodes (~30 µm) and either high-voltage NMC-622 cathodes (1.44 mAh cm -2 ) or high-capacity sulfur cathodes (3.02 mAh cm -2 ).
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.
Not Available
When exposed to fire, the integrity of cement-based materials is governed by thermally-induced changes in the mechanical properties of their binding phase, i.e., the calcium–silicate–hydrate (C–S–H) gel. However, the effect of temperature on the structure, density, and mechanical properties of C–S–H remains only partially known. Here, based on reactive molecular dynamics simulations, we reveal the nature of thermally-induced damage in C–S–H gels. In general, we show that, at the atomic scale, exposure to high temperature results in partial dehydration, volumetric shrinkage, disordering, and stiffening in the C–S–H grains. However, we show that the thermal response of C–S–H strongly depends on its chemical composition, wherein C–S–H systems associated with lower Ca/Si molar ratios are able to undergo higher temperatures before amorphization. Based on these results, we demonstrate that the stiffness of C–S–H gels (i.e., including porosity—as probed by nanoindentation) is governed by a competition between the stiffening of the grains and the decrease in packing density—wherein the latter eventually become predominant.
Thermodynamic properties of Ca–Pb alloys are investigated by electromotive force (emf) measurements to determine equilibrium cell potentials and phase properties for their application in energy storage systems such as liquid metal batteries. Using the electrochemical cell Ca(s) | CaF 2 (s) | Ca(in Pb) at 700–1060 K, cell emf is measured for thirteen Ca–Pb alloys at mole fractions, x Ca = 0.06–0.80. At 873 K, the equilibrium potentials of liquid Ca–Pb alloys are 0.57–0.62 V versus Ca and the activity values are as low as a Ca = 6.2 × 10 -8 at x Ca = 0.06. In addition, the emf values as a function of temperature provide partial molar quantities (entropy and enthalpy) as well as phase transitions which are corroborated by determining transition temperatures and phase constituents using differential scanning calorimetry (DSC) and powder X-ray diffraction (XRD). This study establishes the fundamental data necessary for the design of Pb-containing liquid metal electrodes through the integration of electrochemical, thermal, and structural properties of Ca–Pb electrodes.
The coupled transport of charge and heat provide fundamental insights into the microscopic thermodynamics and kinetics of materials. We describe a sensitive ac differential resistance bridge that enables measurements of the temperature difference on two sides of a coin cell with a resolution of better than 10 μK. We use this temperature difference metrology to determine the ionic Peltier coefficients of symmetric Li-ion electrochemical cells as a function of Li salt concentration, solvent composition, electrode material, and temperature. The Peltier coefficients Π are negative, i.e., heat flows in the direction opposite to the drift of Li ions in the applied electric field, large, –Π > 30 kJ mol –1 , and increase with increasing temperature at T > 300 K. The Peltier coefficient is approximately constant on time scales that span the characteristic time for mass diffusion across the thickness of the electrolyte, suggesting that heat of transport plays a minor role in comparison to the changes in partial molar entropy of Li at the interface between the electrode and electrolyte. Here, our work demonstrates a new platform for studying the non-equilibrium thermodynamics of electrochemical cells and provides a window into the transport properties of electrochemical materials through measurements of temperature differences and heat currents that complement traditional measurements of voltages and charge currents.
In this work, a classical force field based on the General Amber Force Field (GAFF) was refined for the simulation of choline chloride (ChCl) and ethylene glycol (EG) mixtures over a wide composition range by scaling the partial charges and van der Waals parameters. The scaling factors were derived by fitting the simulation results to only eight experimental density and viscosity data points of pure EG, and ChCl/EG mixtures at 1:2, 1:4, and 1:6 molar ratios. Using the refined force field, properties essential for electrochemical applications such as density, viscosity, self-diffusion coefficient, and ionic conductivity were calculated, and excellent agreement to experimental results was found even for compositions and temperatures not used in the fitting procedure. In addition, new experimental data for density, viscosity, and ionic conductivity are reported as a function of temperature and composition for this mixture. To the best of our knowledge, this is the first classical force field developed for the study of ChCl/EG mixtures over a composition range that includes the eutectic point. Using the new model, the liquid dynamics was studied in terms of ionic conductivity. It was found that the dynamics in ChCl/EG mixtures with ChCl mole fraction higher than 20% is similar to that of ionic liquids, high temperature molten salts, and highly concentrated water-in-salt electrolytes.
Advancing circularity through effective chemical recycling strategies is essential for enabling the sustainable use of emerging biodegradable polyesters in film applications. In this context, poly(pentylene adipate-co-terephthalate) with a 40/60 adipate-to-terephthalate ratio (PPeAT60), a significantly stiffer and partially biobased copolyester recently reported by our team, offers improved mechanical performance and processability compared to commercial poly(butylene adipate-co-terephthalate), while providing a platform to evaluate advanced recycling approaches for film applications. Here, this study describes depolymerization of PPeAT60 through amine-catalyzed methanolysis, followed by the repolymerization of the recovered monomers. The resulting repolymerized material exhibited a molar mass distribution, chemical structure, and thermal stability comparable to the virgin polymer, confirming the technical feasibility of the methanolysis-based recycling route. Both polymers displayed nearly identical melting, glass transition, and decomposition temperatures, crystallinity, and modulus, demonstrating full structural and mechanical reproducibility after recycling. Variations in crystallization kinetics, elongation at break, and barrier properties likely arise from trace impurities (<1%) introduced during processing acting as nucleating and stress-inducing agents. Overall, these results demonstrate that PPeAT60 retains its key properties after chemical recycling, establishing methanolysis as a promising platform for recycling of this partially biobased polymer that may be able to replace both conventional polyolefins and PBAT in flexible packaging applications.
A fundamental analysis of multicomponent gas transport models was performed in application to the oxygen electrodes of solid oxide cells. It is common practice to neglect the effect of pressure gradients within oxygen electrodes, even though a net molar flux at the electrolyte surface implies that a pressure gradient must exist. The influence of both Darcy velocity and Knudsen flux are considered in the context of ordinary (Fickian) diffusion, the dusty gas model, and the binary friction model. Comparisons between the models and different sets of assumptions are made via parametric studies on operating load, oxygen partial pressure, microstructural properties, and electrode thickness. Results show that the pressure gradient will have a significant impact on the oxygen concentration distribution and therefore the concentration overpotential. In electrolysis mode, pressure increases up to 1 atm are predicted, indicating that pressure at the electrode/electrolyte interface could contribute to electrode delamination. Additionally, it is found that Darcy's law is insufficient for calculating the pressure distribution without accounting for the flux due to Knudsen diffusion. Additionally, it is found that for the range of properties typical of oxygen electrodes, there is a negligibly small difference between the dusty gas model and binary friction model from a practical standpoint.
Thermodynamic properties of binary Gd-Bi alloys (mole fraction, x Gd = 0.02–0.35) were determined using electromotive force (emf) measurements in molten LiCl-KCl-GdCl 3 electrolyte in complement with structural, microstructural, and thermal characterization. X-ray diffractometry (XRD) and scanning electron microscopy (SEM) identified the phase constituents of Gd-Bi alloys (x Gd = 0.02–0.35) as Bi matrix and GdBi compound. According to differential scanning calorimetry (DSC) measurements up to 1250 K, a eutectic transition [L = Bi + GdBi] was observed at 543 K but the reported peritectic transition [L + GdBi = GdBi 2 ] was not detected, suggesting the absence of the GdBi 2 compound. Based on structural, microstructural, and thermal characterization of Gd-Bi alloys (x Gd = 0.02–0.35), the GdBi 2 compound was not detected, and its thermodynamic stability is questionable. The emf of Gd-Bi alloy (x Gd = 0.16) at 700–1048 K relative to Gd(s) was measured by electrodepositing pure Gd metal at 25 K increments. Further, the emf values of Gd-Bi alloys were determined via coulometric titration of Gd into liquid Bi using a two-phase [L + GdBi] alloy as a reference electrode at 773–973 K. The solubility of Gd in liquid Bi was also estimated to be 0.50 mol% (773 K), 0.99 mol% (873 K), and 2.21 mol% (973 K) and the excess partial molar Gibbs energy (Δ$\mathrm{\bar{G}^{ex}_{Gd}}$) of liquid Gd-Bi alloys was as low as 65 kJ mol –1 , indicating strong chemical interactions between Gd and Bi.
One hundred years ago Hermann Staudinger was strongly criticized by his scientific peers for his macromolecular hypothesis, but today it is hard to imagine a world without polymers. His hypothesis described polymers as macromolecules composed of large numbers of structural units connected by covalent bonds. In the 1990s the concept of supramolecular polymers emerged in the scientific literature as discrete entities of large molar mass comparable to that of classical polymers but built through non-covalent bonds among monomers. Supramolecular polymers exist in biological systems, and potentially blend the physical properties of covalent polymers with unique features such as high degrees of internal order within the polymeric structure, defined shapes, and novel dynamics. This trend article provides a summary of seminal contributions in supramolecular polymerization and provides recent examples from the Stupp laboratory to demonstrate the potential applications of an exciting class of materials composed fully or partially of supramolecular polymers. Finally, in closing, we provide our perspective on future opportunities provided by this field at the onset of a second century of polymers. It is our objective here to demonstrate that this second century could be as prosperous, if not more so, than the preceding one.
Glymes have been extensively studied as solvents for Li-battery electrolytes, most recently in equimolar mixtures with lithium bis(trifluoromethanesulfonyl)imide (LiTFSI), due to their ability to form stable solvates. However, directly quantifying free and coordinated glyme molecules in the liquid state has been challenging due to several experimental limitations. Here, in this work, new vibrational probes are demonstrated for studying the solvation structures of diglyme and triglyme in LiTFSI electrolytes. These IR probes make use of an amine group to report the solvation state of glymes at salt-to-solvent molar ratios ranging from 1:5 to 1:10. Characterization of the thermodynamic properties of the solvent exchange occurring in the first solvation shell of lithium ions (Li + ) showed an equilibrium constant for these probes close to unity at room temperature. This result demonstrates that the probes exhibit a similar solvation behavior to their glyme analogue. Concentration dependence studies also revealed a lack of significant amounts of contact ion pairs at the studied concentrations. Moreover, the first solvation shell of Li + appears to be formed by two partially chelating glyme molecules, establishing that even triglyme with multiple chelation sites does not fully coordinate the cation. Complementary molecular dynamics (MD) simulations agree with the experimental results and suggest that at these concentrations, TFSI – predominantly forms solvent-separated ion pairs. However, the simulations do not properly capture the partial solvation structure of the glyme molecules in the solvation shell of Li + as derived from the experiments.
The quest for achieving high coercivity in Sm(Fe,Co,Ti) 12 alloys, despite their inherent strong magnetocrystalline anisotropy, has posed significant challenges. Recently, (Sm,Zr)(Fe,Co,Ti) 12 monocrystalline particles have exhibited coercivity μ 0 H c > 1.2 T, showcasing promising prospects and significant potential for both manufacturing and research endeavors. This study delves into the structural evolution of (Sm,Zr)(Fe,Co,Ti) 12 (1:12) alloy particles made via the calciothermic reduction-diffusion synthesis process as influenced by the molar ratios of Ca atoms to O 2- ions (Ca/O), annealing time and annealing temperature. Critical insight that informs conditions to optimize the magnetic response is gained via systematic experimentation and advanced electron microscopy. Complex structural features, including core-shell morphologies and intricate multiphase compositions within individual particles, are unveiled. An optimal Ca/O ratio of 1.30 produces particles with a coercivity up to μ 0 H c = 1.63 T, while higher Ca/O ratios induce the formation of a Sm-rich TbCu 7 -type (1: 7 ) phase, which only partially transforms into the desired 1:12 phase during annealing. Persistent remnants of the 1:7 phase locally impact atomic structure, particle morphology, and coercivity. Furthermore, these findings underscore the complex interplay between synthesis parameters, resulting structures, and magnetic properties, informing the design and optimization of high-performance permanent magnets comprised of the (1:12) compound.
The Martian mantle is considered to have a higher Fe/Mg ratio than the Earth's mantle. Ringwoodite, γ-(Mg,Fe) 2 SiO 4 , is likely the dominant polymorph of olivine in the core-mantle boundary (CMB) region of Mars. We synthesized anhydrous iron-rich ringwoodite with molar Mg/(Mg + Fe) = 0.44 and determined its thermal equation of state up to 35 GPa and 750 K by synchrotron X-ray diffraction. Using a third order Birch-Murnaghan equation of state, we obtain K T0 = 182 (3) GPa, K' = 4.6 (2), and α 0 = 3.18 (6) × 10 -5 K -1 . Using these results and an updated mineralogical model with an iron-rich composition of Mg/(Mg + Fe) = 0.75 for the Martian mantle, we estimate ~1900 K for the temperature of the D1000 seismic discontinuity inside Mars. The resulting adiabat predicts a warm aerotherm, which could explain the presence of partial melt at the CMB of Mars recently detected with seismic data from the 2019 InSight mission.
Highlights: • Effect of the CaO/FeO molar ratio in alkali-activated Fe-rich slags is unraveled. • Low CaO/FeO molar ratios enhance the early kinetics and compressive strength. • Reaction extent and strength further increase with higher CaO/FeO ratio after 7 days. • Fe{sup 2+}-trioctahedral layers of the binder are hydroxylated and in a silicate phase. • The final strength is determined by the influence of Ca on the reaction products. In the search for CO{sub 2} mitigating cement technologies, alkali-activated Fe-rich slags present a possible alternative. The influence of the chemical composition of the synthetic Fe-silicate slag on the reaction mechanisms is assessed by varying the CaO/FeO molar ratio. The alkali-activated Fe-rich slags consist of trioctahedral layers that are shown to be hydroxylated and in contact with silicates, in a similar way to phyllosilicate minerals. The formation of this phyllosilicate-like structure is hindered during the first days for samples with high CaO/FeO ratios, because of the partial incorporation of Ca in the trioctahedral layers. At later ages, the samples with higher CaO/FeO ratios gain in reaction extent, to exceed the reaction extent of low CaO/FeO ratios at 28 days and beyond. The increase in compressive strength is even more pronounced than the effect on the reaction extent, which underlines the importance of Ca on the inherent strength of the binding phases.