Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “entropy functional”

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 631 records · Page 35
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Enhancement and maximum in the isobaric specific-heat capacity measurements of deeply supercooled water using ultrafast calorimetry

Knowledge of the temperature dependence of the isobaric specific heat (C p ) upon deep supercooling can give insights regarding the anomalous properties of water. If a maximum in C p exists at a specific temperature, as in the isothermal compressibility, it would further validate the liquid–liquid critical point model that can explain the anomalous increase in thermodynamic response functions. The challenge is that the relevant temperature range falls in the region where ice crystallization becomes rapid, which has previously excluded experiments. Here, we have utilized a methodology of ultrafast calorimetry by determining the temperature jump from femtosecond X-ray pulses after heating with an infrared laser pulse and with a sufficiently long time delay between the pulses to allow measurements at constant pressure. Evaporative cooling of ~15-µm diameter droplets in vacuum enabled us to reach a temperature down to ~228 K with a small fraction of the droplets remaining unfrozen. We observed a sharp increase in C p , from 88 J/mol/K at 244 K to about 218 J/mol/K at 229 K where a maximum is seen. The C p maximum is at a similar temperature as the maxima of the isothermal compressibility and correlation length. From the C p measurement, we estimated the excess entropy and self-diffusion coefficient of water and these properties decrease rapidly below 235 K.

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

Towards accurate prediction of configurational disorder properties in materials using graph neural networks

Abstract The prediction of configurational disorder properties, such as configurational entropy and order-disorder phase transition temperature, of compound materials relies on efficient and accurate evaluations of configurational energies. Previous cluster expansion methods are not applicable to configurationally-complex material systems, including those with atomic distortions and long-range orders. In this work, we propose to leverage the versatile expressive capabilities of graph neural networks (GNNs) for efficient evaluations of configurational energies and present a workflow combining attention-based GNNs and Monte Carlo simulations to calculate the disorder properties. Using the dataset of face-centered tetragonal gold copper without and with local atomic distortions as an example, we demonstrate that the proposed data-driven framework enables the prediction of phase transition temperatures close to experimental values. We also elucidate that the variance of the energy deviations among configurations controls the prediction accuracy of disorder properties and can be used as the target loss function when training and selecting the GNN models. The work serves as a fundamental step toward a data-driven paradigm for the accelerated design of configurationally-complex functional material systems.

Chemistry↗
DOE OSTI2023

Symmetry-resolved entanglement entropy, spectra & boundary conformal field theory

We perform a comprehensive analysis of the symmetry-resolved (SR) entanglement entropy (EE) for one single interval in the ground state of a 1 + 1D conformal field theory (CFT), that is invariant under an arbitrary finite or compact Lie group, G. We utilize the boundary CFT approach to study the total EE, which enables us to find the universal leading order behavior of the SREE and its first correction, which explicitly depends on the irreducible representation under consideration and breaks the equipartition of entanglement. We present two distinct schemes to carry out these computations. The first relies on the evaluation of the charged moments of the reduced density matrix. This involves studying the action of the defect-line, that generates the symmetry, on the boundary states of the theory. This perspective also paves the way for discussing the infeasibility of studying symmetry resolution when an anomalous symmetry is present. The second scheme draws a parallel between the SREE and the partition function of an orbifold CFT. This approach allows for the direct computation of the SREE without the need to use charged moments. From this standpoint, the infeasibility of defining the symmetry-resolved EE for an anomalous symmetry arises from the obstruction to gauging. Finally, we derive the symmetry-resolved entanglement spectra for a CFT invariant under a finite symmetry group. We revisit a similar problem for CFT with compact Lie group, explicitly deriving an improved formula for U(1) resolved entanglement spectra. Using the Tauberian formalism, we can estimate the aforementioned EE spectra rigorously by proving an optimal lower and upper bound on the same. In the abelian case, we perform numerical checks on the bound and find perfect agreement.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Active learning of ternary alloy structures and energies

Abstract Machine learning models with uncertainty quantification have recently emerged as attractive tools to accelerate the navigation of catalyst design spaces in a data-efficient manner. Here, we combine active learning with a dropout graph convolutional network (dGCN) as a surrogate model to explore the complex materials space of high-entropy alloys (HEAs). We train the dGCN on the formation energies of disordered binary alloy structures in the Pd-Pt-Sn ternary alloy system and improve predictions on ternary structures by performing reduced optimization of the formation free energy, the target property that determines HEA stability, over ensembles of ternary structures constructed based on two coordinate systems: (a) a physics-informed ternary composition space, and (b) data-driven coordinates discovered by the Diffusion Maps manifold learning scheme. Both reduced optimization techniques improve predictions of the formation free energy in the ternary alloy space with a significantly reduced number of DFT calculations compared to a high-fidelity model. The physics-based scheme converges to the target property in a manner akin to a depth-first strategy, whereas the data-driven scheme appears more akin to a breadth-first approach. Both sampling schemes, coupled with our acquisition function, successfully exploit a database of DFT-calculated binary alloy structures and energies, augmented with a relatively small number of ternary alloy calculations, to identify stable ternary HEA compositions and structures. This generalized framework can be extended to incorporate more complex bulk and surface structural motifs, and the results demonstrate that significant dimensionality reduction is possible in thermodynamic sampling problems when suitable active learning schemes are employed.

Chemistry↗
DOE OSTI2024

Cosmic neutrino decoupling and its observable imprints: insights from entropic-dual transport

Abstract Very different processes characterize the decoupling of neutrinos to form the cosmic neutrino background (CνB) and the much later decoupling of photons from thermal equilibrium to form the cosmic microwave background (CMB). The CνB emerges from the fuzzy, energy-dependent neutrinosphere and encodes the physics operating in the early universe in the temperature rangeT∼ 10 MeV toT∼ 10 keV. This is the epoch where beyond Standard Model (BSM) physics, especially in the neutrino sector, may be influential in setting the light element abundances, the necessarily distorted fossil neutrino energy spectra, and other light particle energy density contributions. Here we use techniques honed in extensive CMB studies to analyze the CνB as calculated in detailed neutrino energy transport and nuclear reaction simulations of the protracted weak decoupling and primordial nucleosynthesis epochs. Our moment method, relative entropy, and differential visibility approach can leverage future high precision CMB and light element primordial abundance measurements to provide new insights into the CνB and any BSM physics it encodes. We demonstrate that the evolution of the energy spectrum of the CνB throughout the weak decoupling epoch is accurately captured in the Standard Model by only three parameters per species, a non-trivial conclusion given the deviation from thermal equilibrium and the impact of the decrease of electron-positron pairs. Furthermore, we can interpret each of the three parameters as physical characteristics of a non-equilibrium system. Though the treatment presented here makes some simplifying assumptions including ignoring neutrino flavor oscillations, the success of our compact description within the Standard Model motivates its use also in BSM scenarios. We further demonstrate how observations of primordial light element abundances can be used to place constraints on the CνB energy spectrum, deriving response functions that can be applied for general deviations from a thermal spectrum. Combined with the description of those deviations that we develop here, our methods provide a convenient and powerful framework to constrain the impact of BSM physics on the CνB.

Astronomy & Astrophysics↗
DOE OSTI2022

Multi-task graph neural networks for simultaneous prediction of global and atomic properties in ferromagnetic systems *

Abstract We introduce a multi-tasking graph convolutional neural network, HydraGNN, to simultaneously predict both global and atomic physical properties and demonstrate with ferromagnetic materials. We train HydraGNN on an open-source ab initio density functional theory (DFT) dataset for iron-platinum with a fixed body centered tetragonal lattice structure and fixed volume to simultaneously predict the mixing enthalpy (a global feature of the system), the atomic charge transfer, and the atomic magnetic moment across configurations that span the entire compositional range. By taking advantage of underlying physical correlations between material properties, multi-task learning (MTL) with HydraGNN provides effective training even with modest amounts of data. Moreover, this is achieved with just one architecture instead of three, as required by single-task learning (STL). The first convolutional layers of the HydraGNN architecture are shared by all learning tasks and extract features common to all material properties. The following layers discriminate the features of the different properties, the results of which are fed to the separate heads of the final layer to produce predictions. Numerical results show that HydraGNN effectively captures the relation between the configurational entropy and the material properties over the entire compositional range. Overall, the accuracy of simultaneous MTL predictions is comparable to the accuracy of the STL predictions. In addition, the computational cost of training HydraGNN for MTL is much lower than the original DFT calculations and also lower than training separate STL models for each property.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

First-principles calculation of lattice distortion, electronic structure, and bonding properties of GeTe-based and PbSe-based high-entropy chalcogenides

The massive amount of wasted heat energy from industry has pushed the development of thermoelectric (TE) materials that directly convert heat into electricity to a new level of concern. Recently, multicomponent alloys such as GeTe-based and PbSe-based high-entropy (HE) chalcogenides have attracted a great deal of attention due to their potential application as TE materials. The nature of the interatomic bonding, lattice distortion (LD), and the electronic structure in this class of materials is not fully understood. Herein, we report a comprehensive computational investigation of nine GeTe-based HE alloys with eight metallic elements (Ag, Pb, Sb, Bi, Cu, Cd, Mn, and Sn) with large supercells of 1080 atoms each; seven PbSe-based HE solid solutions: Pb0.99−ySb0.012SnySe1−2xTexSx (x = 0.1, 0.2, 0.25, 0.3, 0.35, 0.4, 0.45, with y = 0) with supercells of 1000 atoms each; and five Pb0.99−ySb0.012SnySe1−2xTexSx (y = 0.05, 0.1, 0.15, 0.2, 0.25 with x = 0.25) solid solutions with supercells of 1000 atoms each. All these HE models are theoretically investigated for the first time. The electronic structure, interatomic bonding, charge transfer, and lattice distortion (LD) are investigated by first-principles calculations based on density functional theory. Multicomponent HE alloys can cause a significant LD, which affects their mechanical, thermal, and TE properties. The calculations for the GeTe-based HE chalcogenides showed that they are semiconductors with a narrow bandgap, except for m8, which has a semi-metallic characteristic, and this makes them good candidates for TE applications. For most of these models, the Fermi level shifts upward and locates deeply in the conduction bands, resulting in the enhancement of the electrical conductivity (σ). The bonding properties showed that most bonds in m5 are more dispersed, indicating highest LD and lower lattice thermal conductivity. For PbSe-based HE solid solutions, the LD calculations showed that the models Pb0.99Sb0.012Se0.5Te0.25S0.25 and Pb0.89Sb0.012Sn0.1Se0.5Te0.25S0.25 have the higher LD, and thus a lower lattice thermal conductivity. Such investigations are in high demand since it enables us to design new HE chalcogenides for TE applications. We use the novel concept of total bond order density as a single quantum mechanical metric to characterize the internal cohesion of these HE alloys and correlate with calculated properties, especially the mechanical properties. This work provides a solid database for HE chalcogenides and a road map for many potential applications. Moreover, the computational procedure we developed can be used to design new HE chalcogenides for specific TE applications.

Physics↗
DOE OSTI2022

Characterizing the effect of incrementally increasing dry bulb temperature on linear and nonlinear measures of heart rate variability in nonpregnant, mid-gestation, and late-gestation sows

Abstract Characterizing the sow physiological response to an increased heat load is essential for effective heat stress mitigation. The study objective was to characterize the effects of a 400-min heating episode on sow heart rate variability (HRV) at different reproductive stages. HRV is a commonly used noninvasive proxy measure of autonomic function. Twenty-seven sows were enrolled in the study according to their gestation stage at time of selection: 1) nonpregnant (NP; n = 7), 2) mid-gestation (MID; 57.3 ± 11.8 d gestation; n = 11), and 3) late-gestation (LATE; 98.8 ± 4.9 d gestation; n = 8). The HRV data utilized in the study were collected from each pig as the dry bulb temperature in the room increased incrementally from 19.84 ± 2.15 °C to 35.54 ± 0.43 °C (range: 17.1–37.5 °C) over a 400-min period. After data collection, one 5-min set of continuous heart rate data were identified per pig for each of nine temperature intervals (19–20.99, 21–22.99, 23–24.99, 25–26.99, 27–28.99, 29–30.99, 31–32.99, 33–34.99, and 35–36.99 °C). Mean inter-beat interval length (RR), standard deviation of r-r intervals (SDNN), root mean square of successive differences (RMSSD), high frequency spectral power (HF), sample entropy (SampEn), short-term detrended fluctuation analysis (DFAα1), and three measures (%REC, DET, LMEAN) derived from recurrence quantification analysis were calculated for each data set. All data were analyzed using the PROC GLIMMIX procedure in SAS 9.4. Overall, LATE sows exhibited lower RR than NP sows (P < 0.01). The standard deviation of r–r intervals and RMSSD differed between each group (P < 0.01), with LATE sows exhibiting the lowest SDNN and RMSSD and NP sows exhibiting the greatest SDNN and RMSSD. Late-gestation sows exhibited lower HF than both MID and NP sows (P < 0.0001), greater DFA values than NP sows (P = 0.05), and greater DET compared to MID sows (P = 0.001). Late-gestation also sows exhibited greater %REC and LMEAN compared to MID (P < 0.01) and NP sows (all P < 0.01). In conclusion, LATE sows exhibited indicators of greater autonomic stress throughout the heating period compared to MID and NP sows. However, temperature by treatment interactions were not detected as dry bulb increased. Future studies are needed to fully elucidate the effect of gestational stage and increasing dry bulb temperature on sow HRV.

Agriculture↗
DOE OSTI2020

Crystal structure and physical properties of Yb 2 In and Eu 2–x Yb x In alloys

While binary R E 2 In , where R E = rare earth , have been reported a few decades ago, recent investigations revealed intriguing new physical insights. For instance, the discovery of a nearly ideal first-order ferromagnetic transition in Eu 2 In calls for further exploration of structures and properties of R E 2 In , in particular for the least-documented R E = Eu and Yb cases. In this work, we investigate Eu 2 – x Yb x In pseudobinaries with nominal values of x = 0.25 , 0.5, 0.75, 1, 1.5, 2 by powder x-ray diffraction (including as function of temperature from 100 to 375 K for Yb 2 In ), magnetization (5–300 K), as well as electrical resistivity (5–300 K) and calorimetric (2–150 K) measurements for Yb 2 In . Compared to other RE , Yb or Eu always raise challenging questions linked to their valence states. From average atomic volume, Yb is anticipated to be divalent in Yb 2 In , at least between 100 and 375 K, which is in line with the absence of 4 f magnetism. In agreement with x-ray diffraction and magnetization data, the resistivity of Yb 2 In is rather featureless and typical of a metal. Establishing Yb 2 In as a nonmagnetic isostructural reference for Eu 2 In allows one to use its heat capacity to revisit that of the latter, and get experimental insights into the exceptional magnetocaloric effect of the compound with Eu. In particular, we show that a third of the total magnetic entropy ( S m ≈ 35.6 J mo l – 1 K – 1 at T = 100 K ) is concentrated in a 3 K temperature window around the T C of Eu 2 In . Starting from the ferromagnetic compound Eu 2 In [ T C = 55.2 ( 5 ) K ] , we show that Yb substitutions in Eu 2 – x Yb x In lead to a decrease in both the Curie temperature [ T C = 41 ( 2 ) and 32(2) K for x = 0.25 and 0.5] and magnetic saturation, while weakening the first-order character of the transition as x increases. A significant isothermal entropy change of 5.1 ( 4 ) J mo l – 1 K – 1 for Δ B = 2 T is found at 44 K in Eu 1.75 Yb 0.25 In , demonstrating that the giant magnetocaloric effect of Eu 2 In can be tuned to lower temperatures by Yb substitutions.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Disorder-induced magnetoelastic behaviors of MnTexSbyBi1-x-y alloys

This dataset contains input and output files from density functional theory (DFT) simulations used to study the disorder-induced magnetoelastic behaviors of MnTexSbyBi1-x-y (0 ≤ x + y ≤ 1) alloys and their binary end members MnTe, MnSb, and MnBi. The alloys adopt the hexagonal NiAs-type (nickeline) structure and span ternary (MnTexSb1-x, MnTexBi1-x, MnBixSb1-x), and quaternary compositions across the full MnTe–MnSb–MnBi composition triangle. For each alloy composition, the dataset provides DFT calculations in three magnetic configurations: A-type antiferromagnetic (AFM), C-type AFM, and ferromagnetic (FM). Every magnetic configuration folder contains the fully relaxed crystal structure (CONTCAR), VASP input parameters (INCAR), and the main VASP output file (OUTCAR), from which total electronic energies, Mn magnetic moments, lattice parameters, and percent volume changes between magnetic states are extracted. These data are used to construct compositional phase diagrams, evaluate thermodynamic stability (formability), and map magnetoelastic responses across the alloy space. For A-type AFM and FM configurations, additional data are provided as follows: (i) FORCE_CONSTANTS and thermal_properties.yaml files at the top level of A-type_AFM/ and FM/ folders — present only for compositions marked with an asterisk (*) in Table I of the main text. These are derived from Phonopy finite-displacement calculations on full disordered 128-atom supercells and provide vibrational free energy, entropy (Svib)contribution from explicit disorder calculations. (Table I of the associated main manuscript) (ii) A VCA/ subfolder within A-type_AFM/ and FM/, containing FORCE_CONSTANTS and thermal_properties.yaml from Virtual Crystal Approximation phonon calculations (without spin-orbit coupling). VCA data are available for all compositions and are used to estimate vibrational contributions to the Gibbs free energy across the full composition space. (iii) A SOC/ subfolder containing CONTCAR, INCAR, and OUTCAR from spin-orbit coupling calculations, providing relativistic corrections to electronic energies and lattice parameters (Tables S2–S3 of the SM, and Table I of the main manuscript). (iv) A SOC/VCA/ subfolder containing FORCE_CONSTANTS and thermal_properties.yaml from VCA phonon calculations performed within the SOC framework, combining relativistic and vibrational thermodynamic corrections. The computed properties are used to map the AFM–FM magnetic crossover near MnTe0.75Sb0.25, demonstrate disorder- and spin-induced phonon broadening, identify a semiconductor-to-metal crossover, and quantify the pronounced magnetoelastic volume response near the magnetic phase boundary.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Tensile deformation behavior of twist grain boundaries in CoCrFeMnNi high entropy alloy bicrystals

High entropy alloys (HEA) are a class of materials that consist of multiple elemental species in similar concentrations. The use of elements in far from dilute concentrations introduces a multi-dimensional composition design space by which the properties of metallic systems can be tailored. While the mechanical behavior of HEAs has been the subject of active research recently, the role of grain boundaries (GBs) in their deformation behavior remains poorly understood. Motivated by recent experiments on HEAs demonstrating that GBs act as nucleation sites for deformation twins, herein, we leverage atomistic simulations to construct a series of equiatomic CoCrFeMnNi HEA bicrystals with $\langle 110 \rangle$ and $\langle 111 \rangle$ symmetric twist GBs and examine their tensile behavior and underlying deformation mechanisms at 77 K. Simulation results reveal that plastic deformation proceeds by the nucleation of partial dislocations from GBs, which then grow with further loading by bowing into the bulk crystals leaving behind stacking faults. Variations in the nucleation stress exist as function of GB character, defined in this work by the twist angle. Our results provide future avenues to explore GBs as a microstructure design tool to develop HEAs with tailored properties.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Accurate prediction of short-range order and its effect on thermodynamic, structural, and electronic properties of disordered alloys: exemplified in archetypical Cu 3 Au

Electronic-structure methods based on density-functional theory (DFT) were used to quantify the effect of chemical short-range order (SRO) on thermodynamic, structural, and electronic properties of archetypal face-centered-cubic (fcc) Cu3Au alloy. We showed that SRO can be tuned to alter bonding and lattice dynamics (i.e., phonons) and detail how these properties are changed with SRO. Thermodynamically favorable SRO significantly improved the phase stability of fcc Cu3Au from -0.0343 eV-atom -1 to –0.0682 eV-atom -1 . We used our DFT-based linear-response theory to predict SRO and its electronic origin, and accurately estimate the observed transition temperature, ordering instability (L1 2 ), and Warren-Cowley SRO parameters, in agreement with experiments. The accurate prediction of real-space SRO gives an edge over computationally and resource intensive approaches such as monte-carlo methods or experiments, which will enable large scale molecular dynamic simulations by providing supercells with optimized SRO. Here we also analyzed phonon dispersion and estimated the vibrational entropy change (from 9kB at 300 K to 6kB at 100 K) in fcc Cu3Au. We established from SRO analysis that exclusion of chemical interactions may lead to a skewed view of true properties in chemically complex alloys. The first-principles methods described in this work are generally applicable to any arbitrary solid-solution alloys, including multi-principal-element alloys, therefore, holds promise for designing technologically useful materials.

36 MATERIALS SCIENCE↗
DOE OSTI2023

A deep neural network regressor for phase constitution estimation in the high entropy alloy system Al-Co-Cr-Fe-Mn-Nb-Ni

High Entropy Alloys (HEAs) are composed of more than one principal element and constitute a major paradigm in metals research. The HEA space is vast and an exhaustive exploration is improbable. Therefore, a thorough estimation of the phases present in the HEA is of paramount importance for alloy design. Machine Learning presents a feasible and non-expensive method for predicting possible new HEAs on-the-fly. A deep neural network (DNN) model for the elemental system of: Mn, Ni, Fe, Al, Cr, Nb, and Co is developed using a dataset generated by high-throughput computational thermodynamic calculations using Thermo-Calc. The features list used for the neural network is developed based on literature and freely available databases. A feature significance analysis matches the reported HEAs phase constitution trends on elemental properties and further expands it by providing so far-overlooked features. The final regressor has a coefficient of determination ( r 2 ) greater than 0.96 for identifying the most recurrent phases and the functionality is tested by running optimization tasks that simulate those required in alloy design. The DNN developed constitutes an example of an emulator that can be used in fast, real-time materials discovery/design tasks.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Phase equilibrium of liquid water and hexagonal ice from enhanced sampling molecular dynamics simulations

We study the phase equilibrium between liquid water and ice Ih modeled by the TIP4P/Ice interatomic potential using enhanced sampling molecular dynamics simulations. Our approach is based on the calculation of ice Ih-liquid free energy differences from simulations that visit reversibly both phases. The reversible interconversion is achieved by introducing a static bias potential as a function of an order parameter. The order parameter was tailored to crystallize the hexagonal diamond structure of oxygen in ice Ih. Furthermore, we analyze the effect of the system size on the ice Ih-liquid free energy differences, and we obtain a melting temperature of 270 K in the thermodynamic limit. This result is in agreement with estimates from thermodynamic integration (272 K) and coexistence simulations (270 K). Since the order parameter does not include information about the coordinates of the protons, the spontaneously formed solid configurations contain proton disorder as expected for ice Ih.

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

Ab Initio Molecular Dynamics Study of Pd Nucleation on a γ-Al 2 O 3 Surface

We perform ab initio molecular dynamics (AIMD) free energy simulations to investigate the stability and dynamics of dispersed Pd atoms and subnanometer clusters on γ-Al 2 O 3 . We consider nucleation for varying Pd loadings on dry γ-Al 2 O 3 (100), dry γ-Al 2 O 3 (110) and wet γ-Al 2 O 3 (110) and examine temperature, entropy and water coverage effects on the stability and dynamics of Pd clusters. At low Pd loading (1.07 Pd atoms/nm 2 on the (100) and 0.75 Pd/nm 2 on the (110)), association of the Pd atoms is thermodynamically unfavorable on the dry surfaces and almost thermoneutral on the wet (110) owing to water blocking of the AlIII defect sites. In agreement with past studies, the Pd atoms are quite mobile but the diffusion free energy barriers do not increase with the water coverage, as previously suggested, but rather remain similar in magnitude owing to translational entropy at finite temperatures. At these Pd loadings, oxidized Pd(II) atoms on the wet (110) surface are dispersed and immobile. At higher Pd loading (1.12 and 1.5 Pd atoms/nm 2 ) on the (110) surface, Pd 3 and Pd 4 clusters are thermodynamically very stable and form fast. At the 1.12 Pd atoms/nm 2 loading, the formation of the Pd 3 cluster is two orders of magnitude slower when the surface is hydroxylated as the system becomes somewhat trapped before the dinuclear cluster Pd 2 collides with a Pd atom.

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

Unravelling the mechanism of pressure induced polyamorphic transition in an inorganic molecular glass

The atomic structure of a germanium doped phosphorous selenide glass of composition Ge 2.8 P 57.7 Se 39.5 is determined as a function of pressure from ambient to 24 GPa using Monte-Carlo simulations constrained by high energy x-ray scattering data. The ambient pressure structure consists primarily of P 4 Se 3 molecules and planar edge shared phosphorus rings, reminiscent of those found in red phosphorous as well as a small fraction of locally clustered corner-sharing GeSe 4 tetrahedra. This low-density amorphous phase transforms into a high-density amorphous phase at ~6.3 GPa. The high-pressure phase is characterized by an extended network structure. The polyamorphic transformation between these two phases involves opening of the P 3 ring at the base of the P 4 Se 3 molecules and subsequent reaction with red phosphorus type moieties to produce a cross linked structure. The compression mechanism of the low-density phase involves increased molecular packing, whereas that of the high pressure phase involves an increase in the nearest-neighbor coordination number while the bond angle distributions broaden and shift to smaller angles. The entropy and volume changes associated with this polyamorphic transformation are positive and negative, respectively, and consequently the corresponding Clapeyron slope for this transition would be negative. This result has far reaching implications in our current understanding of the thermodynamics of polyamorphic transitions in glasses and glass-forming liquids.

42 ENGINEERING↗
DOE OSTI2020

Temperature and Solvent Effects on H 2 Splitting and Hydricity: Ramifications on CO 2 Hydrogenation by a Rhenium Pincer Catalyst

The catalytic hydrogenation of carbon dioxide holds immense promise for applications in sustainable fuel synthesis and hydrogen storage. Mechanistic studies that connect thermodynamic parameters with the kinetics of catalysis can provide new understanding and guide predictive design of improved catalysts. Reported here are thermochemical and kinetic analyses of a new pincer-ligated rhenium complex ( tBu POCOP)Re(CO) 2 ( tBu POCOP = 2,6-bis(di-tert-butylphosphinito)phenyl) that catalyzes CO 2 hydrogenation to formate with faster rates at lower temperature. Because the catalyst follows the prototypical “outer sphere” hydrogenation mechanism, comprehensive studies of temperature and solvent effects on the H 2 splitting and hydride transfer steps are expected to be relevant to many other catalysts. Strikingly large entropy associated with cleavage of H2 results in a strong temperature dependence on the concentration of [( tBu POCOP)Re(CO) 2 H] – present during catalysis, which is further impacted by changing the solvent from toluene to tetrahydrofuran to acetonitrile. New methods for determining the hydricity of metal hydrides and formate at temperatures other than 298 K were developed, providing insight into how temperature can influence the favorability of hydride transfer during catalysis. These thermochemical insights guided the selection of conditions for CO 2 hydrogenation to formate with high activity (up to 364 h –1 at 1 atm or 3330 h –1 at 20 atm of 1:1 H 2 CO 2 ). In cases where hydride transfer is the highest individual kinetic barrier, entropic contributions to outer sphere H 2 splitting lead to a unique temperature dependence: catalytic activity increases as temperature decreases in tetrahydrofuran (200-fold increase upon cooling from 50 to 0 °C) and toluene (4-fold increase upon cooling from 100 to 50 °C). Furthermore, the ramifications on catalyst structure-function relationships are discussed, including comparisons between “outer sphere” mechanisms and metal–ligand cooperation mechanisms.

14 SOLAR ENERGY↗