Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Predictive Data analytics”

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 1,027 records · Page 57

Data-Driven Kinetic Reaction Networks for Separation Chemistry

Understanding complex, multistep chemical reactions at the molecular level is a major challenge whose solution would greatly benefit the design and optimization of numerous chemical processes. The separation of rare-earth (4f) and actinide (5f) elements is an example where improving our chemical understanding is important for designing and optimizing new chemistries, even with a limited number of observations. Here, in this work, we leverage data-driven artificial intelligence and machine-learning approaches to develop kinetic reaction networks that describe the liquid–liquid extraction mechanism of uranium using N,N-di-2-ethylhexyl-isobutyramide (DEHiBA). Specifically, we compare and contrast the properties of two classes of models: (1) purely data-driven models that are regularized using chemistry-agnostic, L1 regression and (2) chemistry-informed models that are regularized using relative reaction energies provided by quantum mechanical calculations. We observe that purely data-driven models are unbiased, simple, and accurate in their predictions of experimental measurements when provided with sufficient data but are difficult to fully constrain and interpret. In contrast, chemistry-informed models exhibit significantly improved chemical interpretability and consistency, providing a detailed description of the separation process while achieving high accuracy through ensemble averaging. Overall, the dominant species predicted to be extracted into the organic phase is UO 2 (NO 3 ) 2 (DEHiBA) 2 , agreeing with experimental slope analysis, thermodynamic modeling, EXAFS, and crystal structures. This work demonstrates that leveraging the fundamental structure of the problem can lead to efficient learning schemes that provide both accurate predictions and chemical insights at a low computational cost.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Data-driven predictions of complex organic mixture permeation in polymer membranes

Membrane-based organic solvent separations are rapidly emerging as a promising class of technologies for enhancing the energy efficiency of existing separation and purification systems. Polymeric membranes have shown promise in the fractionation or splitting of complex mixtures of organic molecules such as crude oil. Determining the separation performance of a polymer membrane when challenged with a complex mixture has thus far occurred in an ad hoc manner, and methods to predict the performance based on mixture composition and polymer chemistry are unavailable. Here, we combine physics-informed machine learning algorithms (ML) and mass transport simulations to create an integrated predictive model for the separation of complex mixtures containing up to 400 components via any arbitrary linear polymer membrane. We experimentally demonstrate the effectiveness of the model by predicting the separation of two crude oils within 6-7% of the measurements. Integration of ML predictors of diffusion and sorption properties of molecules with transport simulators enables for the rapid screening of polymer membranes prior to physical experimentation for the separation of complex liquid mixtures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Probing Pore Size and Connectivity in Porous Silicas Using 13C MAS NMR Spectroscopy of Supercritical Methane

Measuring the pore size and pore-size distributions and exploring the fluid-exchange dynamics between different types of pores in porous materials remains a significant experimental challenge but is critical to understanding catalysis, chromatography, nutrient cycling, and a whole range of geochemical phenomena, including shale gas and tight gas extraction. Here, we present the results of 1D 13 C NMR and 2D exchange spectroscopy (EXSY) NMR investigations of a porous silica using supercritical methane (scCH 4 ) as a direct probe of pore size and fluid exchange between pore types. The results show that the 13 C chemical shift of scCH 4 adsorbed in nanometer-scale silica pores becomes more negative with increasing pore diameter, in agreement with trends reported for gas hydrates, zeolites, MOFs, and clays and other microporous (<10 nm pores) geochemical materials. These 13 C chemical shifts follow a natural log–linear trend with pore size in the vacuum-dried porous silicas studied here, allowing one to predict pore size based on the 13 C chemical shift in dry silica nanopores. The EXSY and 1D 13 C NMR results both show that CH 4 exchanges between pore and bulk fluid environments over rate scales from 0.01 to 1 kHz in vacuum-dried samples and that exchange dynamics must be considered when interpreting 13 C NMR of pore-adsorbed CH 4 . Introducing H 2 O to the system causes the 13 C chemical shift to become more negative with increasing H 2 O content as a result of H 2 O preferentially filling small pores, forcing CH 4 to occupy larger pores that are better connected to the bulk environment. Likewise, the 1D 13 C data show a decreased distribution of exchange rates between pore-adsorbed and bulk CH 4 in the presence of H 2 O.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Carbon-Carbon Aeroshell Chemical Reaction Model in an Open RPS in Martian Environment

A new RPS is being designed to be open to the vacuum of space to improve its thermal efficiency during deep space missions. This design, however, could lead to some potential complications if it is ever used in the low-pressure gas environment of mission destinations like Mars. This report documents the development of a chemical kinetics reaction model of the FWPF carbon-carbon aeroshell in the Martian gaseous environment. This model uses an “infinite sink” assumption where the RPS is open to the Martian environment with open exchange of flowing gases, so the gas chemistry over the aeroshell does not change with the reaction and represents the worst case scenario of the chemical reaction of the FWPF in the Martian environment. The model combines the kinetics equation of both oxidation reactions that would occur on Mars: the Langmuir-Hinshelwood kinetics of the Boudouard reaction of CO 2 with the chemical kinetics of O 2 oxidation of the carbon-carbon. The FWPF chemical reaction model was benchmarked using the Opila’s experimental chemical reaction test data (1) in a simulated Martian environment to increase the validity of the model’s predicted mass loss rate. The output of the benchmarked chemical reaction model showed a trade-off between the maximum temperature of the aeroshell to have less than 10% reaction at the mission life. For example, an RPS with a five-year intended mission life on Mars has a maximum aeroshell temperature of 835 K / 562 °C, increasing the mission life to ten years reduces the maximum operating temperature to 810 K / 527 °C. If the aeroshell operating temperature in the new RPS on the surface of Mars is 960 Kelvin, the model calculates that the life is limited to 0.22 years / 81 days. Increasing the aeroshell temperature to 1343 Kelvin, further reduces the life to 0.375 days / 9 hrs. This illustrates the thermally activated nature of the chemical kinetics model in limiting the life of the aeroshell.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Carbon-Carbon Aeroshell Chemical Reaction Model in an Open RPS in Martian Environment

A new RPS is being designed to be open to the vacuum of space to improve its thermal efficiency during deep space missions. This design, however, could lead to some potential complications if it is ever used in the low-pressure gas environment of mission destinations like Mars. This report documents the development of a chemical kinetics reaction model of the FWPF carbon-carbon aeroshell in the Martian gaseous environment. This model uses an “infinite sink” assumption where the RPS is open to the Martian environment with open exchange of flowing gases, so the gas chemistry over the aeroshell does not change with the reaction and represents the worst case scenario of the chemical reaction of the FWPF in the Martian environment. The model combines the kinetics equation of both oxidation reactions that would occur on Mars: the Langmuir-Hinshelwood kinetics of the Boudouard reaction of CO2 with the chemical kinetics of O2 oxidation of the carbon-carbon. The FWPF chemical reaction model was benchmarked using the Opila’s experimental chemical reaction test data (1) in a simulated Martian environment to increase the validity of the model’s predicted mass loss rate. The output of the benchmarked chemical reaction model showed a trade-off between the maximum temperature of the aeroshell to have less than 10% reaction at the mission life. For example, an RPS with a five-year intended mission life on Mars has a maximum aeroshell temperature of 835 K / 562 °C, increasing the mission life to ten years reduces the maximum operating temperature to 810 K / 527 °C. If the aeroshell operating temperature in the new RPS on the surface of Mars is 960 Kelvin, the model calculates that the life is limited to 0.22 years / 81 days. Increasing the aeroshell temperature to 1343 Kelvin, further reduces the life to 0.375 days / 9 hrs. This illustrates the thermally activated nature of the chemical kinetics model in limiting the life of the aeroshell.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Multipartite entanglement in the one-dimensional spin- 1 2 Heisenberg antiferromagnet

Multipartite entanglement refers to the simultaneous entanglement between multiple subsystems of a many-body quantum system. While multipartite entanglement can be difficult to quantify analytically, it is known that it can be witnessed through the quantum Fisher information (QFI), a quantity that can also be related to dynamical Kubo response functions. In this work, we first show that the finite temperature QFI can generally be expressed in terms of a static structure factor of the system, plus a correction that vanishes as T→0. This implies that the static structure factor witnesses multipartite entanglement near quantum critical points at temperatures below a characteristic energy scale of the system. Therefore, in systems with a known static structure factor, we can deduce finite temperature scaling of multipartite entanglement and low temperature entanglement depth without knowledge of the full dynamical response function of the system. This is particularly useful to study 1D quantum critical systems in which sub-power-law divergences can dominate entanglement growth, where the conventional scaling theory of the QFI breaks down. Furthermore, the 1D spin- 1 2 antiferromagnetic Heisenberg model is an important example of such a system, and we show that multipartite entanglement in the Heisenberg chain diverges nontrivially as ~ln(1/T) 3/2 . We verify these predictions with calculations of the QFI using conformal field theory and matrix product state simulations. Finally, we discuss the implications of our results for experiments to probe entanglement in quantum materials, comparing to neutron scattering data in KCuF 3 , a material well described by the Heisenberg chain.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Crash Risks Evaluation of Urban Expressways: A Case Study in Shanghai

We report that proactive traffic safety management systems can reduce crashes by identifying crash precursors, evaluating real-time crash risks, and implementing suitable interventions. The basic prerequisite for developing such a system is to propose a reliable crash risk evaluation model that takes real-time traffic flow data as input. Previous studies have primarily focused on real-time crash prediction using some statistical or machine-learning methods. However, further quantitative evaluation and classification of crash risks have been ignored. In this study, we conduct a systematic crash risk evaluation workflow, including crash risk prediction, crash risk quantification, and crash risk classification. Specifically, the crash risk prediction using an extended logit model is proposed, from which CAS, CSD, UAS, DAS, DTV are identified to be contributing factors of crash risks. Then a crash risk quantification model based on the parameter evaluation of the extended logit model is developed. The crash risks of urban expressways and their spatial-temporal evolution trends are quantified. Finally, the crash risks are classified into high crash risk level, moderate crash risk level, and low crash risk level by the k-means cluster algorithm. Then the threshold boundaries of different crash risk levels are determined. The research results provide a proactive guidance for traffic safety management of urban expressways.

33 ADVANCED PROPULSION SYSTEMS↗

Identifying crystal structures and chemical reactions at the interface of stanene on Bi 2 Te 3

Synthesizing monolayers and heterostructures is an enabling approach to extract new physical phenomena from bulk materials. Among the structures amenable to this approach is stanene, which is a monolayer of tin, similar to graphene, and has been predicted to host one-dimensional topological states at its edges. Stanene can be tuned by decorating with different adatoms, which makes it a promising platform on which to engineer topological devices. Here, we deposit Sn on Bi 2 Te 3 and characterize the growth using anomalous synchrotron x-ray scattering and x-ray photoelectron spectroscopy (XPS). X-ray diffraction data reveal the formation of epitaxial Sn-based structures, along with penetration of Sn into the Bi 2 Te 3 , with Sn intercalating between the upper 10 Bi 2 Te 3 quintuple layers. Additionally, XPS data show deposited Sn reacting to form SnTe and Bi at the Bi 2 Te 3 surface. The calculated heat of reaction for Sn and Bi2Te3 is consistent with an exothermic reaction to SnTe and Bi. Using thermodynamic calculations as a guide, we identify several candidate substrates that can stabilize the stanene phase.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Molecular Modes Elucidate the Nuclear Magnetic Resonance Relaxation of Viscous Fluids

The Bloembergen, Purcell, and Pound (BPP) theory of nuclear magnetic resonance (NMR) relaxation in fluids dating back to 1948 continues to be the linchpin in interpreting NMR relaxation data in applications ranging from characterizing fluids in porous media to medical imaging (MRI). The BPP theory is founded on assuming molecules are hard spheres with 1 H– 1 H dipole pairs reorienting randomly; assumptions that are severe in light of modern understanding of liquids. Nevertheless, it is intriguing to this day that the BPP theory was consistent with the original experimental data for glycerol, a hydrogen-bonding molecular fluid for which the hard-sphere-rigid-dipole assumption is inapplicable. To better understand this incongruity, atomistic molecular simulations are used to compute 1 H NMR T 1 relaxation dispersion (i.e., frequency dependence) in two contrasting cases: glycerol, and a (non hydrogen-bonding) viscosity standard. At high viscosities, simulations predict distinct functional forms of T 1 for glycerol compared to the viscosity standard, in agreement with modern measurements, yet both in contrast to BPP theory. The cause of these departures from BPP theory is elucidated, without assuming any relaxation models and without any free parameters, by decomposing the simulated T 1 response into dynamic molecular modes for both intramolecular and intermolecular interactions. The decomposition into dynamic molecular modes provides an alternative framework to understand the physics of NMR relaxation for viscous fluids.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A system identification approach for non-intrusive reduced order modeling of radiation-induced photocurrents

In this study, development of compact photocurrent models is currently dominated by analytical techniques that rely on physical assumptions to render the governing equations solvable in a closed form. Violation of these assumptions can reduce the accuracy of the models and/or limit their scope. In this paper we show that system identification of nonlinear state-space systems can serve as an alternative numerical basis for non-intrusive reduced order modeling of photocurrent effects. To that end we develop a compact gray box photocurrent model (GBPM) by using a state-space representation with a low-dimensional latent state equation that mimics a mathematical model for the response of an idealized class of devices to ionizing radiation. In so doing we obtain a model that learns the dynamics of a quantity of interest directly from its measurements without requiring snapshots of the internal device state or its discretized model, and can be inferred from very small data sets. To demonstrate the approach we train the GBPM using a small experimental data set for a Z5236 Zener diode and a small synthetic data set obtained by simulating a synthetic pn-junction device. We then compare the GBPMs with black box models trained on the same data and show that performance of the latter is limited by the size of the data set, while the former are able to achieve excellent performance in both the reproductive and the predictive regimes.

97 MATHEMATICS AND COMPUTING↗

Discovering equations that govern experimental materials stability under environmental stress using scientific machine learning

Abstract While machine learning (ML) in experimental research has demonstrated impressive predictive capabilities, extracting fungible knowledge representations from experimental data remains an elusive task. In this manuscript, we use ML to infer the underlying differential equation (DE) from experimental data of degrading organic-inorganic methylammonium lead iodide (MAPI) perovskite thin films under environmental stressors (elevated temperature, humidity, and light). Using a sparse regression algorithm, we find that the underlying DE governing MAPI degradation across a broad temperature range of 35 to 85 °C is described minimally by a second-order polynomial. This DE corresponds to the Verhulst logistic function, which describes reaction kinetics analogous to self-propagating reactions. We examine the robustness of our conclusions to experimental variance and Gaussian noise and describe the experimental limits within which this methodology can be applied. Our study highlights the promise and challenges associated with ML-aided scientific discovery by demonstrating its application in experimental chemical and materials systems.

36 MATERIALS SCIENCE↗

Disorderly Conduct of Benzamide IV: Crystallographic and Computational Analysis of High Entropy Polymorphs of Small Molecules

Benzamide, a simple derivative of benzoic acid and a common intermediate of pharmaceutical compounds, was reported to form two polymorphs in 1832, but the single crystal structure of the more stable form was not solved until 1959. Nearly 50 years later, the second form was characterized by powder diffraction, followed shortly thereafter by characterization of a third form, a polytype of the most thermodynamically stable Form I. These two new forms, Forms II and III, are metastable. Herein, we describe a fourth polymorph, Form IV, discovered by melt crystallization concurrently with its crystallization under confinement at small length scales (<10 nm), where it is stable indefinitely. Form III exists under confinement in larger pores, and melting point data for different pore sizes corroborate the existence of Form IV below 10 nm. Form IV is highly disordered, precluding indexing of powder diffraction data other than hk0 reflections. Nonetheless, a combination of powder X-ray diffraction and computational crystal structure prediction reveals that Form IV contains a 2D motif resembling that of Form II, but with longrange order in the third dimension masked by ubiquitous stacking faults. This approach relies on distilling a large number of candidate structures to a few possible disorder models based on benzamide tetrads that organize in 2D parquet-like tiles, with organization along the third dimension, that can be modeled with various stacking fault configurations having distinct intermolecular interactions and translations in the dimension orthogonal to the tiling planes. These observations reveal a bewildering crystallographic complexity for such a simple molecule. Nonetheless, the approach described herein demonstrates that challenging structures that may be abandoned prematurely because of poor crystallinity, twinning, or disorder can be solved.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

LDRD Abbreviated report: High-Order General-Discrete-Ordinates Method Enabling Efficient Deterministic Transport in Hydrodynamic Simulations

Deterministic transport simulations for national-security and energy applications often operate in high-dimensional phase-space, where accuracy and cost both become major challenges. A common numerical artifact in such problems is the “ray-effect,” which appears as unphysical streaks. Beyond misinterpretation, these artifacts can contaminate tightly coupled physics, such as fluid dynamics, radiation-hydrodynamics, and laser-plasma interactions, eroding the predictive capability of entire multiphysics workflows. Our objective was to make high-dimension studies practical on modern hardware while mitigating the ray-effect without relying on prohibitively expensive sampling approaches such as Monte Carlo methods. We developed the Generic Discretization Library (GenDiL), a Graphics Processing Unit (GPU)-first framework that uses high-order Discontinuous Galerkin (DG) methods and matrix-free algorithms to reduce memory usage and improve computational efficiency, critical for phase-space simulations. GenDiL supports phase-space adaptivity in both mesh size and polynomial order (hp-adaptivity) to place resolution only where it is needed. A central capability is Local Dimensional Refinement (LDR), which couples lower-dimension continuum models to higher-dimension kinetic models through stable and conservative interfaces, so that high-fidelity physics is applied only in regions where it is essential. Building on the GenDiL framework, we developed the General SN (GSN) family of algorithms as a true generalization of the polar SN approach (discrete ordinates, often denoted SN). Rather than tying discrete ordinates to a specific polar change of coordinates, GSN formulates transport on an arbitrary change of coordinates chosen to reduce ray-effect. We studied two complementary variants: an analytic variant, where the coordinate map is prescribed in advance by a closed-form function; and a data-driven variant, where a quantity of interest, such as the net flux, guides the coordinate system. GenDiL provides the library infrastructure for efficient GPU execution, but the GSN concept is algorithmic and independent of any one library. Across representative high-dimension tests, including non-symmetric solutions, both variants delivered strong ray-effect mitigation at practical cost, moving four- to six-dimensional analysis toward repeatable, routine studies.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Invariant Molecular Representations for Heterogeneous Catalysis

Catalyst screening is a critical step in the discovery and development of heterogeneous catalysts, which are vital for a wide range of chemical processes. In recent years, computational catalyst screening, primarily through density functional theory (DFT), has gained significant attention as a method for identifying promising catalysts. However, the computation of adsorption energies for all likely chemical intermediates present in complex surface chemistries is computationally intensive and costly due to the expensive nature of these calculations and the intrinsic idiosyncrasies of the methods or data sets used. This study introduces a novel machine learning (ML) method to learn adsorption energies from multiple DFT functionals by using invariant molecular representations (IMRs). To do this, we first extract molecular fingerprints for the reaction intermediates and later use a Siamese-neural-network-based training strategy to learn invariant molecular representations or the IMR across all available functionals. Our Siamese network-based representations demonstrate superior performance in predicting adsorption energies compared with other molecular representations. Notably, when considering mean absolute values of adsorption energies as 0.43 eV (PBE-D3), 0.46 eV (BEEF-vdW), 0.81 eV (RPBE), and 0.37 eV (scan+rVV10), our IMR method has achieved the lowest mean absolute errors (MAEs) of 0.18 0.10, 0.16, and 0.18 eV, respectively. These results emphasize the superior predictive capacity of our Siamese network-based representations. The empirical findings in this study illuminate the efficacy, robustness, and dependability of our proposed ML paradigm in predicting adsorption energies, specifically for propane dehydrogenation on a platinum catalyst surface.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Deep Learning Approaches for Predicting the Surface Tension of Ionic Liquids

Ionic liquids (ILs) are a novel class of solvents that have attracted significant attention due to their unique and tunable properties. Among their physiochemical characteristics, surface tension plays a critical role in various industrial applications including electrolytes, heat transfer fluids, and separation processes. However, because of the exploratory nature of IL design and the vast combinatorial space of possible anion–cation pairs, the experimental determination of these properties is often impractical, being both time-consuming and costly. To overcome these challenges, computational approaches are increasingly employed to develop accurate predictive models that can accelerate IL discovery and design. In this study, we present two deep learning (DL) models for predicting the surface tension of ILs across a broad temperature range at a constant pressure. The models use simplified molecular input line entry system, SMILES, representations of ILs to extract molecular features as inputs. Both DL models demonstrate excellent agreement with experimental data, achieving an R 2 value of 0.990 and a root-mean-square error of 0.792 mN/m. In conclusion, these results offer valuable insights for the rapid screening and rational design of ILs with tailored surface tension values.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Baseflow Identification via Explainable AI With Kolmogorov‐Arnold Networks

Abstract Hydrological models often involve constitutive laws that may not be optimal in every application. We propose to replace such laws with the Kolmogorov‐Arnold networks (KANs), a class of neural networks designed to identify symbolic expressions. We demonstrate KAN's potential on the problem of baseflow identification, a notoriously challenging task plagued by significant uncertainty. KAN‐derived functional dependencies of the baseflow components on the aridity index outperform their original counterparts; they demonstrate that water availability, rather than potential evapotranspiration, drives baseflow by constraining actual evapotranspiration under arid conditions. On a test set, they increase the Nash‐Sutcliffe efficiency (NSE) by 65%, decrease the root mean squared error by 29%, and increase the Kling‐Gupta efficiency by 34%. This superior performance is achieved while reducing the number of fitting parameters from three to two. Next, we use data from 378 catchments across the continental United States to refine the water‐balance equation at the mean‐annual scale. The KAN‐derived equations based on the refined water balance outperform both the current aridity index model, with up to a 105% increase in NSE, and the KAN‐derived equations based on the original water balance. While the performance of our model and tree‐based machine learning methods is similar, KANs offer the advantage of simplicity and transparency and require no specific software or computational tools. This case study focuses on the aridity index formulation, but the approach is flexible and transferable to other hydrological processes. Plain Language Summary Equations used in hydrologic model are often suboptimal, resulting in reduced prediction accuracy and efficiency. We implemented Kolmogorov‐Arnold networks (KAN), a machine learning algorithm for deriving symbolic formulations, to estimate groundwater recharge and showed that it outperforms an existing state‐of‐the‐art semi‐empirical formulation. In hydrology, Nash‐Sutcliffe efficiency (NSE), root mean squared error (RMSE), and Kling‐Gupta efficiency (KGE) are commonly used to evaluate model performance. Higher NSE and KGE values indicate better performance, while lower RMSE values are preferable. Our results show that NSE increased by 71%, RMSE decreased by 32%, and KGE improved by 25%. In addition, KAN identifies an optimal functional form and can be used to derive new analytical formulas using the prior knowledge. The KAN‐inspired equation outperformed the original formulation and reduced the fitting parameters. Furthermore, we refined the water‐balance equation at the mean‐annual scale and showed that, based on the new water‐balance equation, KAN can derive new formulations that are superior to the original aridity index formulations (up to 105% increase in NSE) and KAN‐derived equations based on the original water balance. These findings highlight the significant potential of KAN to advance the scientific understanding of a wide range of hydrologic processes. Key Points Kolmogorov‐Arnold networks (KANs) enhance interpretability of machine‐learned hydrological models KAN‐derived symbolic formulations outperform state‐of‐the‐art semi‐empirical aridity indices KAN‐identified functional form yields an analytical index with fewer fitting parameters and improved performance

baseflow↗

Palladium at high pressure and high temperature: A combined experimental and theoretical study

Palladium is one of the most important technological materials, yet its phase diagram remains poorly understood. At ambient conditions, its solid phase is face-centered cubic (fcc). However, another solid phase of Pd, body-centered cubic (bcc), was very recently predicted in two independent theoretical studies to occur at high pressures and temperatures. In this work, we report an experimental study on the room-temperature equation of state (EOS) of Pd to a pressure of 80 GPa, as well as a theoretical study on the phase diagram of Pd including both fcc-Pd and bcc-Pd. Our theoretical approach consists in ab initio quantum molecular dynamics (QMD) simulations based on the Z methodology which combines both direct Z method for the simulation of melting curves and inverse Z method for the calculation of solid–solid phase transition boundaries. We obtain the melting curves of both fcc-Pd and bcc-Pd and an equation for the fcc–bcc solid–solid phase transition boundary as well as the thermal EOS of Pd which is in agreement with experimental data and QMD simulations. We uncover the presence of another solid phase of Pd on its phase diagram, namely, random hexagonal close-packed (rhcp), and estimate the location of the rhcp-bcc solid–solid phase transition boundary and the rhcp–fcc–bcc triple point. We also discuss the topological similarity of the phase diagrams of palladium and silver, the neighbor of Pd in the periodic table. We argue that Pd is a reliable standard for shock-compression studies and present the analytic model of its principal Hugoniot in a wide pressure range.

36 MATERIALS SCIENCE↗

Palladium at high pressure and high temperature: A combined experimental and theoretical study

Palladium is one of the most important technological materials, yet its phase diagram remains poorly understood. At ambient conditions, its solid phase is face-centered cubic (fcc). However, another solid phase of Pd, body-centered cubic (bcc), was very recently predicted in two independent theoretical studies to occur at high pressures and temperatures. In this work, we report an experimental study on the room-temperature equation of state (EOS) of Pd to a pressure of 80 GPa, as well as a theoretical study on the phase diagram of Pd including both fcc-Pd and bcc-Pd. Our theoretical approach consists in ab initio quantum molecular dynamics (QMD) simulations based on the Z methodology which combines both direct Z method for the simulation of melting curves and inverse Z method for the calculation of solid–solid phase transition boundaries. We obtain the melting curves of both fcc-Pd and bcc-Pd and an equation for the fcc–bcc solid– solid phase transition boundary as well as the thermal EOS of Pd which is in agreement with experimental data and QMD simulations. We uncover the presence of another solid phase of Pd on its phase diagram, namely, random hexagonal close-packed (rhcp), and estimate the location of the rhcp-bcc solid–solid phase transition boundary and the rhcp–fcc–bcc triple point. We also discuss the topological similarity of the phase diagrams of palladium and silver, the neighbor of Pd in the periodic table. We argue that Pd is a reliable standard for shock-compression studies and present the analytic model of its principal Hugoniot in a wide pressure range.

36 MATERIALS SCIENCE↗