Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Enthalpy”

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

Determining the Reaction Kinetics and Thermodynamics of a Diels–Alder Network Using Dynamic Gel Criteria

We undertook a detailed rheological investigation to evaluate the kinetic parameters of the forward and reverse Diels–Alder (DA) reactions of a model network cross-linked using a furan prepolymer and a common aromatic bismaleimide. At high temperature where the Winter–Chambon’s criterion of frequency-independence was more applicable, a multiwave technique permitted van’t Hoff analysis and calculation of the reaction thermodynamic parameters, specifically the enthalpy and entropy of the reaction: ΔH° = –38.3 ± 5.2 kJ mol –1 and ΔS° = –94.3 ± 13.4 J mol –1 . At mild temperatures where the G'–G" crossover point is experimentally convenient to measure gelation, isothermal tests were used to obtain reasonable fDA kinetic parameters from Eyring analysis such as the apparent activation enthalpy and entropy of ΔH$^{‡}_{fDA}$ = 76.8 ± 6.9 kJ mol –1 and ΔS$^{‡}_{fDA}$= –82.8 ± 22.2 J mol –1 K –1 . Comparable rheokinetic methods include cross-linking density measurements and stress relaxation tests to calculate effective kinetics, whereas the critical gel conversion was consistently applied here. As a result, rate data are fitted with the Arrhenius equation for comparison purposes and the Eyring equation to demonstrate its broader utility.

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

Bottlebrush Block Copolymer Thin Films

The self-assembly of block copolymers (BCPs) in thin films is governed by interfacial interactions (enthalpy) between the blocks and interfaces (at both the substrate and the surface). Advantage can also be taken of entropy to control the polymer orientation. Here, we synthesized a series of bottlebrush block copolymers (BBCPs), with deuterated polystyrene (DPS) as the core block and poly­(solketal acrylate) (PSA) as the corona block, where the backbone length (N BB) and grafting densities (GDs) were varied. The hydrophobic PSA block was converted to a hydrophilic poly­(glyceryl acrylate) (PGA) block by solid-state hydrolysis, bringing the BBCP from a disordered state into an ordered state with a lamellar microdomain morphology. The orientation of the morphology as a function of distance from the interfaces was systematically investigated by atomic force microscopy (AFM), interference microscopy, grazing-incidence small-angle X-ray scattering (GISAXS), and grazing-incidence small angle neutron scattering (GISANS). For N BB = 1 (diblock BCP), a mixed lamellar orientation was found. For N BB = 5 (star-like BBCP), a vertical orientation of the lamellar microdomains was found dominantly at the interfaces that propagated into the bulk of the film due to the entropic penalty associated with constraints on the junction points. Increasing the N BB to 50 (rod-like BBCP) and 100 (worm-like BBCP) yielded a random orientation, as the side chains became overcrowded when aligned in the same direction. By reducing the grafting density (GD) of a N BB = 50 BBCP from 100% to 75% and 50% by the incorporation of spacers into the backbone, the steric hindrance of side chains is reduced markedly, facilitating a vertical orientation of the lamellar microdomains. However, the vertical orientation becomes more random within the bulk of the thin film for different BBCP architectures. By modifying the substrate with a layer of PS, the preferential interactions between the core block of the BBCP and the interfaces induce a reconfiguration of the BBCP, drawing the core block out to contact substrate directly driven by enthalpy. Under these conditions, N BB = 5 showed a predominantly parallel orientation of the lamellar microdomains with the corona block minimizing contact with the substrate. Higher N BB values of 50 and 100 showed a more random orientation of the lamellar microdomains at different depths on the PS-modified substrate. Overall, thin-film morphology was regulated by tuning the configuration of the BBCP, varying the entropic contribution to the microdomain orientation.

Chen, Zhan↗
DOE OSTI2022

Variable-Temperature Scattering and Spectroscopy Characterizations for Temperature-Dependent Solution Assembly of PffBT4T-Based Conjugated Polymers

The solution structure of conjugated polymers (CPs) from which the films are cast is critical for tailoring the thin-film morphology thus device performance. Here, we used multimodal variable-temperature scattering and spectroscopy tools to fully quantify the solution assembly of poly[(5,6-difluoro-2,1,3-benzothiadiazol-4,7-diyl)-alt-(3,3'''-dialkyl-2,2';5',2'';5'',2'''-quaterthiophen-5,5'''-diyl)] (PffBT4T) polymers with varying side-chain lengths at different assembly temperatures. The conformational and aggregation behaviors for PffBT4T-based CPs were found to be very sensitive to both temperature and side chain length using ultraviolet–visible (UV–vis) spectroscopy, nuclear magnetic resonance (NMR) spectroscopy, dynamic light scattering (DLS), and small-angle neutron scattering (SANS). We found that with slightly increasing side chain length from 2-octyldodecyl (C8C12) to 2-nonyltridecyl (C9C13), PffBT4T-based CPs show a significant decrease in aggregation-to-dissolved chain transition temperature (10 °C), degree of aggregation, enthalpy change of aggregation, and size of the aggregates in solution. At room temperature, PffBT4T polymer strongly aggregated to form fabric structure with the film thickness of a few nanometers in thickness and hundreds of nanometers in length, as probed by atomic force microscopy (AFM), transmission electronic microscopy (TEM), and dynamic light scattering (DLS). At the elevated temperature above the aggregation-to-dissolved chain transition temperature, PffBT4T is fully dissolved and adopts a semiflexible coil conformation with the persistence length of 3.1 nm for PffBT4T-C8C12 and a slightly increased persistence length of 3.4 nm for PffBT4T-C9C13, according to temperature-dependent SANS measurements. Longer side chains of PffBT4T-C9C13 also lead to less aggregation enthalpy gain compared with PffBT4T-C8C12. This work provides a solution structure manipulating strategy of CPs and thus will inspire the molecular design and processing protocols of CPs toward higher performance electronic devices.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Reactive Species and Reaction Pathways for the Oxidative Cleavage of 4-Octene and Oleic Acid with H 2 O 2 over Tungsten Oxide Catalysts

The oxidative cleavage of carbon-carbon double bonds (C=C) in alkenes and fatty acids produces aldehydes and acids valued as chemical intermediates. Solid tungsten oxide catalysts are low cost, non-toxic, and selective for the oxidative cleavage of C=C bonds with hydrogen peroxide (H 2 O 2 ) and are, therefore, a promising option for continuous processes. Despite the relevance of these materials, the elementary steps involved and their sensitivity to the form of W sites present on surfaces have not been described. In this work, we combine in situ spectroscopy and rate measurements to identify significant steps in the reaction, the reactive species present on catalysts, and examine differences between the kinetics of this reaction on isolated W atoms grafted to alumina and on those exposed on crystalline WO 3 nanoparticles. Raman spectroscopy shows that W-peroxo complexes (W-(η 2 -O 2 ) formed from H 2 O 2 ) react in a kinetically relevant step to produce epoxides, which undergo hydrolysis at protic surface sites. Subsequently, the CH 3 CN solvent deprotonates diols to form alpha hydroxy ketones that react to form aldehydes and water following nucleophilic attack of H 2 O 2 . Turnover rates for oxidative cleavage, determined by in situ site titrations, on WO x -Al 2 O 3 are 75% greater than those on WO 3 at standard conditions. These differences reflect activation enthalpies (ΔH ‡ ) for oxidative cleavage of 4-octene that are much lower for the isolated WO x sites (36±3 and 60±6 kJ∙mol -1 for WO x -Al 2 O 3 and WO 3 , respectively) and correlate strongly with the difference between enthalpies of adsorption for epoxyoctane (ΔH ads,epox ), which resembles the transition state for epoxidation. The WO x -Al 2 O 3 catalysts mediate oxidative cleavage of oleic acid with H 2 O 2 following a mechanism comparable to that for 4-octene. The WO 3 materials, however, form only the epoxide and do not cleave the C-C bond or produce aldehydes and acids. These differences reflect the distinct site requirements for these reaction pathways and indicate that acid sites required for diol formation are strongly inhibited by oleic acids and epoxides on WO 3 whereas the Al 2 O 3 support provides sites competent for this reaction and increase the yield of the oxidative cleavage products.

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

Demonstrating the Critical Role of Solvation in Supported Ti and Nb Epoxidation Catalysts via Vapor-Phase Kinetics

Catalytic oxidation of hydrocarbons with hydrogen peroxide (H 2 O 2 ) has been of the utmost importance for several decades. The vast majority of studies have been performed in the condensed phase, even though condensed phases introduce complex solvent effects and can promote the leaching of active sites. In response, we have built a custom reactor system to understand H 2 O 2 activation and selective oxidation in the vapor-phase. In this report, we study the epoxidation of cyclohexene with H 2 O 2 over four Lewis-acidic metal oxide catalysts: Ti and Nb grafted on SiO 2 and on the Zr based metal–organic framework, NU-1000. The M–SiO2 materials are highly selective to the formation of epoxides and diols, as they can be in the condensed phase, while the NU-1000 based materials are far more prone to overoxidation to CO 2 , which appears to be connected to their strong reactant adsorption. Apparent activation energies are calculated for all materials when operating in the same kinetic regime, and the heats of cyclohexene adsorption into their pores are then used to directly compare intrinsic enthalpies of activation in the vapor vs condensed phase for the M–SiO 2 catalysts. Nb–SiO 2 catalysts exhibit similar intrinsic enthalpies of activation in the vapor and condensed phases, whereas the condensed phase transition state in Ti–SiO 2 is 24 kJ/mol lower in energy than that of the same material in the vapor phase. These experiments establish another methodology for understanding the various roles of solvent in selective oxidation reactions and studying these reactions under conditions that differ significantly from the thousands of prior studies in the condensed phase.

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

Entropy Stabilization Effects and Ion Migration in 3D “Hollow” Halide Perovskites

A recently discovered new family of 3D halide perovskites with the general formula (A) 1–x (en) x (Pb) 1–0.7x (X) 3–0.4x (A = MA, FA; X = Br, I; MA = methylammonium, FA = formamidinium, en = ethylenediammonium) is referred to as “hollow” perovskites owing to extensive Pb and X vacancies created on incorporation of en cations in the 3D network. The “hollow” motif allows fine tuning of optical, electronic, and transport properties and bestowing good environmental stability proportional to en loading. To shed light on the origin of the apparent stability of these materials, we performed detailed thermochemical studies, using room temperature solution calorimetry combined with density functional theory simulations on three different families of “hollow” perovskites namely en/FAPbI 3 , en/MAPbI 3 , and en/FAPbBr 3 . We found that the bromide perovskites are more energetically stable compared to iodide perovskites in the FA-based hollow compounds, as shown by the measured enthalpies of formation and the calculated formation energies. The least stable FAPbI 3 gains stability on incorporation of the en cation, whereas FAPbBr 3 becomes less stable with en loading. This behavior is attributed to the difference in the 3D cage size in the bromide and iodide perovskites. Configurational entropy, which arises from randomly distributed cation and anion vacancies, plays a significant role in stabilizing these “hollow” perovskite structures despite small differences in their formation enthalpies. With the increased vacancy defect population, we have also examined halide ion migration in the FA-based “hollow” perovskites and found that the migration energy barriers become smaller with the increasing en content.

36 MATERIALS SCIENCE↗
DOE OSTI2026

A High-Fidelity Molecular Model of the Cu(111) Repeating Unit

Dynamic processes at surfaces are central to heterogeneous catalysis, but their atomistic mechanism(s) can prove difficult to elucidate due to variations in material structure and the corresponding impact on reactivity. Moreover, disparities between reaction conditions and those employed for spectroscopic characterization at surfaces can inhibit detailed understanding of catalysis-relevant chemistries. Herein, we substantiate the so-called “cluster-surface” analogy by leveraging a low-valent tricopper architecture ( 1 ) as a model system for small molecule activation at Cu(111). Two reaction classes are explored: the adsorption of carbon monoxide (CO) and the dissociative adsorption of dihydrogen (H 2 ). These processes serve as an ideal testbed to compare the reactivity of a molecular cluster ( 1 ) to that of a heterogeneous surface, as both reactions have empirical data from measurements performed on crystalline Cu(111). Cluster 1 reversibly binds CO. Variable temperature NMR analysis with 13 CO reveals a favorable enthalpy but large negative entropy (−5.1 kcal × mol –1 and −22.9 cal × mol –1 × K –1 , respectively) for CO binding, affording a process that is marginally endergonic at room temperature (ΔG ads (298.15 K) = 1.7 ± 0.5 kcal × mol –1 ). Similarly, analogous to a Cu(111) surface, 1 is shown to oxidatively add (chemisorb) H 2 . Kinetic parameters were determined for this process and the activation enthalpy (8.4 ± 0.5 kcal × mol –1 ) closely mirrors that established for H 2 binding at the Cu(111) facet (6.0 to 12.4 kcal × mol –1 ). Together, these results showcase that a trinuclear cluster can reproduce the small molecule binding and activation energetics of a bulk crystalline surface, setting the stage for studying less-defined surface processes in an atomically precise molecular setting.

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

Computational Discovery of Ultralow Thermal Conductivity in the Energy-Degenerate Polymorphic Crystal Family A 2 M 2 M’Q 4

Crystalline materials, characterized by their well-defined lattices, typically exhibit a unique global thermodynamic minimum for a specific composition. However, in this study, we discover a quaternary chalcogenide family, A 2 M 2 M’Q 4 (A: alkali metals; M: coinage metal; M’: transition or group-IVA metals; Q: chalcogens), that exhibits pervasive energy (near-)degeneracy. For a given composition, multiple structurally distinct polymorphs exist within a formation enthalpy window of only a few milli-electron volts per atom. We quantify this inherent structural flexibility using a dedicated descriptor, σ f : the standard deviation of formation enthalpies among degenerate (meta)stable polymorphs. The consistently low σf observed across the A 2 M 2 M’Q 4 family signifies a characteristically shallow and frustrated potential energy landscape, which drives pronounced lattice anharmonicity, marking these materials as prime candidates for ultralow lattice thermal conductivity (κ L ). Employing an advanced high-throughput computational framework that integrates thermodynamics, lattice dynamics, and thermal conductivity calculations, we screen 1215 A 2 M 2 M’Q 4 compounds, identifying 30 stable candidates with κ L < 0.5 W m –1 K –1 at 300 K. Among them, Rb 2 Ag 2 SnTe 4 and Rb 2 Au 2 HfTe 4 , two representatives from the IVA and TM subgroups, are predicted to show ultralow room-temperature κ L of 0.174 W m –1 K –1 and 0.295 W m –1 K –1 , respectively. A systematic analysis suggests that the nonbonding and antibonding states induced by “dual rattlers” are the origin of low thermal conductivity in these compounds. Our results position the A 2 M 2 M’Q 4 family as a rich source of intrinsic thermal insulators and suggest that polymorphic energy degeneracy may serve as a valuable signpost for identifying crystalline families with potential anharmonicity.

cations↗
DOE OSTI2025

Gamma Radiation Effects in Amorphous and Crystalline Polyaryletherketones

This work investigates the effect of gamma radiation on chain scission and crosslinking in polyaryletherketone (PAEK) thermoplastic filaments, namely, polyetheretherketone (PEEK) and polyetherketoneketone (PEKK) that is either semicrystalline (PEKK-c) or highly amorphous (PEKK-a). The PAEK thermoplastic filaments were exposed to total ionizing doses (TID) of gamma radiation from 1 to 8 MGy in the Gamma Irradiation Facility at Sandia National Laboratories. A combination of thermal and FT-IR spectral analysis of the irradiated specimens was used to study the chain scission and crosslinking behavior between PEEK, PEKK-c, and PEKK-a. Crystallinity decreased with increased TID and melt enthalpy decreased. Despite clear indications of changes to melt and crystallization behavior, any change to the glass-transition temperature was minor until a post-irradiation thermal annealing occurred, which caused an increase in crystallinity. Gamma radiation impacted PEKK-a more due to the lower degree of initial crystallinity, where a reduction in melt enthalpy, crystallinity, and crystallization kinetics was observed. Despite these effects, the ability to achieve a higher degree of crystallization (X c < 20%) by annealing was not affected. The results presented here provide important insight into the performance of a class of commercialized PEKK thermoplastics exposed to a harsh radiation environment.

Crystallinity↗
DOE OSTI2024

The role of charge in microdroplet redox chemistry

In charged water microdroplets, which occur in nature or in the lab upon ultrasonication or in electrospray processes, the thermodynamics for reactive chemistry can be dramatically altered relative to the bulk phase. Here, we provide a theoretical basis for the observation of accelerated chemistry by simulating water droplets of increasing charge imbalance to create redox agents such as hydroxyl and hydrogen radicals and solvated electrons. We compute the hydration enthalpy of OH - and H + that controls the electron transfer process, and the corresponding changes in vertical ionization energy and vertical electron affinity of the ions, to create OH • and H • reactive species. We find that at ~ 20 - 50% of the Rayleigh limit of droplet charge the hydration enthalpy of both OH - and H + have decreased by >50 kcal/mol such that electron transfer becomes thermodynamically favorable, in correspondence with the more favorable vertical electron affinity of H + and the lowered vertical ionization energy of OH - . We provide scaling arguments that show that the nanoscale calculations and conclusions extend to the experimental microdroplet length scale. The relevance of the droplet charge for chemical reactivity is illustrated for the formation of H 2 O 2 , and has clear implications for other redox reactions observed to occur with enhanced rates in microdroplets.

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

Engineering metal oxidation using epitaxial strain

The oxides of platinum group metals are promising for future electronics and spintronics due to the delicate interplay of spin-orbit coupling and electron correlation energies. However, their synthesis as thin films remains challenging due to their low vapour pressures and low oxidation potentials. Here we show how epitaxial strain can be used as a control knob to enhance metal oxidation. Using Ir as an example, we demonstrate the use of epitaxial strain in engineering its oxidation chemistry, enabling phase-pure Ir or IrO 2 films despite using identical growth conditions. The observations are explained using a density-functional-theory-based modified formation enthalpy framework, which highlights the important role of metal-substrate epitaxial strain in governing the oxide formation enthalpy. We also validate the generality of this principle by demonstrating epitaxial strain effect on Ru oxidation. The IrO 2 films studied in our work further revealed quantum oscillations, attesting to the excellent film quality. The epitaxial strain approach we present could enable growth of oxide films of hard-to-oxidize elements using strain engineering.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A theoretical investigation of the effect of Ga alloying on thermodynamic stability, electronic-structure, and oxidation resistance of Ti 2 AlC MAX phase

We present a systematic investigation of thermodynamic stability, phase-reaction, and chemical activity of Al containing disordered Ti 2 (Al-Ga)C MAX phases using machine-learning driven high-throughput framework to understand the oxidation resistance behavior with increasing temperature and exposure to static oxygen. The A-site (at Al) disordering of Ti 2 AlC with Ga shows significant change in the chemical activity of Al with increasing temperature and exposure to static oxygen, which is expected to enable surface segregation of Al, thereby, the formation of Al 2 O 3 and improved oxidation resistance. We performed in-depth convex hull analysis of ternary Ti-Al-C, Ti-Ga-C, and Ti-Al-Ga-C based MAX phase, and provide detailed contribution arising from electronic, chemical and vibrational entropies. The thermodynamic analysis shows change in the Gibbs formation enthalpy (ΔG form ) at higher temperatures, which implies an interplay of temperature-dependent enthalpy and entropic contributions in oxidation resistance Ga doped Ti 2 AlC MAX phases. A detailed electronic structure and chemical bonding analysis using crystal orbital Hamilton population method reveal the origin of change in phases stability and in oxidation resistance in disorder Ti2(Al 1-x Ga x )C MAX phases. Our electronic structure analysis correlate well with the change in oxidation resistance of Ga doped MAX phases. We believe our study provides a useful guideline to understand to role of alloying on electronic, thermodynamic, and oxidation related mechanisms of bulk MAX phases, which can work as a precursor to understand oxidation behavior of twodimensional MAX phases, i.e., MXenes (transition metal carbides, carbonitrides and nitrides).

36 MATERIALS SCIENCE↗
DOE OSTI2024

Semi-empirical model for Henry’s law constant of noble gases in molten salts

Henry’s law constant, which describes the proportionality of dissolved gas to partial pressure of free gas in liquid–gas equilibrium systems, can also be applied to mass transport applications. In this work, we investigated an approach for determining the solubility of noble gases in a molten salt liquid utilizing the equilibrium concept of Henry’s gas constant. Henry’s gas constant is described as a mathematical function dependent on the van der Waals radius of the noble gas and the temperature of the molten salt. The alteration in Gibbs free energy encompasses contributions from both surface and volume energies. Enthalpy and entropy are deduced from these surface and volume energies in the Gibbs free energy formulation. A comparative analysis was conducted between the conventional method and our proposed model. Moreover, useful chemical properties can be determined from examination of surface and volume energies. Our findings provide an accurate and general theory of Gibbs free energy that can be validated experimentally based on the model proposed herein. This work unifies the prediction of Henry gas constant and subsequently the entropy and enthalpy calculation for noble gases in a molten salt solution to a single functional form using van der Waals radius of the gas and temperature of the system. This functional form is then used to perform a multiple regression method to find two parameters corresponding to the surface energy and volume energy. These two parameters are consistent between all combinations of noble gas and molten salt.

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

Defect graph neural networks for materials discovery in high-temperature clean-energy applications

We present a graph neural network approach that fully automates the prediction of defect formation enthalpies for any crystallographic site from the ideal crystal structure, without the need to create defected atomic structure models as input. Here we used density functional theory reference data for vacancy defects in oxides, to train a defect graph neural network (dGNN) model that replaces the density functional theory supercell relaxations otherwise required for each symmetrically unique crystal site. Interfaced with thermodynamic calculations of reduction entropies and associated free energies, the dGNN model is applied to the screening of oxides in the Materials Project database, connecting the zero-kelvin defect enthalpies to high-temperature process conditions relevant for solar thermochemical hydrogen production and other energy applications. The dGNN approach is applicable to arbitrary structures with an accuracy limited principally by the amount and diversity of the training data, and it is generalizable to other defect types and advanced graph convolution architectures. In conclusion, it will help to tackle future materials discovery problems in clean energy and beyond.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Conductivity, structure, and thermodynamics of Y 2 Ti 2 O 7 –Y 3 NbO 7 solid solutions

The defect fluorite yttrium niobate Y 3 NbO 7 and pyrochlore yttrium titanate Y 2 Ti 2 O 7 solid solutions have been synthesized via a solid state synthesis route. The resulting stoichiometry of the oxides is Y 2+x Ti 2-2x Nb x O 7 , where x = 0 to x = 1. All of the samples were single-phase; however, for those with a predominant fluorite phase, a small amount of additional pyrochlore phase was detected. The volume of the solid solution unit cells linearly increases with increase in yttrium niobate content. The water uptake increases with (x) and the protonic defect concentration reaches almost 4.5 × 10 -3 mol mol -1 at 300 °C. The calculated enthalpy of formation from oxides suggests strong stability for all of the compositions, with the values of enthalpy ranging from -84.6 to -114.3 kJ mol -1 . The total conductivity does not have a visible dependence on Y 3 NbO 7 content. For each compound, the total conductivity is higher in wet air. Interestingly, for samples where x < 0.5, the ratio of conductivity in hydrogen to air increases with increasing temperature, while for x > 0.5, the trend is the opposite.

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

Boundary layer description of directional polymer crystallisation

Nearly fifty years ago Lovinger and Gryte suggested that the directional crystallization of a polymer was analogous to the quiescent isothermal crystallization experiment but at a supercooling where the crystal growth velocity was equal to the velocity of the moving front. Our experiments showed that this equivalence holds in a detailed manner at low directional velocities. To understand the underlying physics of these situations, we modeled the motion of a crystallization front in a liquid where the left side boundary is suddenly lowered below the melting point (Stefan's problem) but with the modification that the crystallization kinetics follow a version of the Avrami model. Our numerical results surprisingly showed that the results of the polymer analog track with the Stefan results which were derived for a simple liquid that crystallizes completely at its melting point; in particular, the position of the crystal growth-front evolved with time exactly as in the Stefan problem. The numerical solution also showed that the temperature in the immediate vicinity of the growth-front decreased with increasing front velocity, which is in line with Lovinger and Gryte's ansatz. To provide a clear theoretical understanding of these numerical results we derive a boundary layer solution to the governing coupled differential equations of the polymer problem. The analytical results are in agreement with our observations from experiments and numerical computations but show that this equivalence between the small molecule and polymer analog only holds in the limit where the crystallization enthalpy is much larger than the rate at which heat is conducted away in the polymer. In particular, in the context of the temperature profile, the enthalpy generated by the crystallisation process which is spread out over a narrow spatial region can be approximated as a point source whose location and temperature correspond to the Lovinger–Gryte ansatz.

36 MATERIALS SCIENCE↗
DOE OSTI2023

An extensive assessment of the performance of pairwise and many-body interaction potentials in reproducing ab initio benchmark binding energies for water clusters n = 2–25

We assess the performance of 7 pairwise additive (TIP3P, TIP4P, TIP4P-ice, TIP5P, OPC, SPC, SPC/E) and 8 families of many-body potentials (q-AQUA, HIPPO, AMOEBA, EFP, TTM, WHBB, MB-pol, MB-UCB) in reproducing high-level ab initio benchmark values, CCSD(T) or MP2 at the complete basis set (CBS) limit for the binding energy and the many-body expansion (MBE) of water clusters n = 2–11, 16–17, 20, 25. By including a large range of cluster sizes having dissimilar hydrogen bonding networks, we obtain an understanding of how these potentials perform for different hydrogen bonding arrangements that are mostly outside of their parameterization range. While it is appropriate to compare the results of ab initio based many-body potentials directly to the electronic binding energies (D e 's), the pairwise additive ones are compared to the enthalpies at T = 298 K, ΔH(298 K), as the latter class of force fields are parametrized to reproduce enthalpies (implicitly accounting for zero-point energy corrections) rather than binding energies. We find that all pairwise additive potentials considered overestimate the reference ΔH values for the n = 2–25 clusters by >13%. For the water dimer (n = 2) in particular, the errors are in the range 83–119% for the pairwise additive potentials studied since these are based on an effective rather than the true 2-body interaction specifically designed as a means of partially accounting for the missing many-body terms. This stronger 2-body interaction is achieved by an enhanced monomer dipole moment that mimics its increase from the gas phase monomer to the condensed phase value. Indeed, for cluster sizes n ≥ 4 the percent deviations become slightly smaller (albeit all exceeding 13%). In contrast, we find that the many-body potentials perform more accurately in reproducing the electronic binding energies (D e 's) throughout the entire cluster range (n = 2–25), all reproducing the ab initio benchmark binding energies within ±7% of the respective CBS values. We further assess the ability of a subset of the many-body potentials (MB-UCB, q-AQUA, MB-pol, and TTM2.1-F) to also reproduce the magnitude of the ab initio many-body energy terms for water cluster sizes n = 7, 10, 16 and 17. The potentials show an overall good agreement with the available benchmark values. However, we identify characteristic differences upon comparing the many-body terms at both the ab initio-optimized geometries and the respective potential-optimized geometries to the reference ab initio values. Additionally, by applying this analysis to a wide range of cluster sizes, trends in the MBE of the potentials with increasing cluster size can be identified. Finally, in an attempt to draw a parallel between the pairwise additive and many-body potentials, we report the analysis of the individual molecular dipole moments for water clusters with 1 to ~4 solvation shells with the TTM2.1-F potential. We find that the internally solvated water molecules have in general a larger molecular dipole moment ranging from 2.6–3.0 D. This justifies the use of an enhanced, with respect to the gas-phase value, molecular dipole moment for the pairwise additive potentials, which is intended to fold in the many body terms into an effective (enhanced) pairwise interaction through the choice of the charges. These results have important implications for the development of future generations of efficient, transferable, and highly accurate classical interaction potentials in both the pairwise additive and many-body categories.

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

Energetics of formation and stability in high pressure steam of barium lanthanide cobaltite double perovskites

Here, this study concerns energetics of formation and the stability in high water partial pressure of BaLnCo 2 O 6-δ , (Ln = La, Pr, Nd, and Gd) (BLnC) and BaGd 1-x La x Co 2 O 6-δ , where x = 0.2, 0.5, and 0.7 (BGLC) double perovskite cobaltites. Those materials are extensively studied due to their potential applications as a positrode in electrochemical devices. Therefore, their stability under such conditions is a key issue. All investigated materials are thermodynamically stable relative to binary oxides and exhibit strongly exothermic enthalpies of formation. Moreover, BaGd 0.3 La 0.7 Co 2 O 6-δ and BaGd 0.8 La 0.2 Co 2 O 6-δ remain the main perovskite structure up to 3 bars of water vapor at 400 °C. At higher steam pressure, reaching 10 bar at 300 °C, the partial decomposition to constituent oxides and hydroxides was observed. The BGLC compounds exhibit higher negative formation enthalpies in comparison to single-Ln compositions, which does not translate into higher chemical stability under high steam pressures since the BLnC series retained the main perovskite structure at higher temperatures as well as in higher water vapor pressures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗