Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “extrapolation”

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 613 records · Page 34

Exploring toroidal anvil profiles for larger sample volumes above 4 Mbar

Abstract With the advent of toroidal and double-stage diamond anvil cells (DACs), pressures between 4 and 10 Mbar can be achieved under static compression, however, the ability to explore diverse sample assemblies is limited on these micron-scale anvils. Adapting the toroidal DAC to support larger sample volumes offers expanded capabilities in physics, chemistry, and planetary science: including, characterizing materials in soft pressure media to multi-megabar pressures, synthesizing novel phases, and probing planetary assemblages at the interior pressures and temperatures of super-Earths and sub-Neptunes. Here we have continued the exploration of larger toroidal DAC profiles by iteratively testing various torus and shoulder depths with central culet diameters in the 30–50 µm range. We present a 30 µm culet profile that reached a maximum pressure of 414(1) GPa based on a Pt scale. The 300 K equations of state fit to our P–V data collected on gold and rhenium are compatible with extrapolated hydrostatic equations of state within 1% up to 4 Mbar. This work validates the performance of these large-culet toroidal anvils to > 4 Mbar and provides a promising foundation to develop toroidal DACs for diverse sample loading and laser heating.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Data-driven emulation of modal aerosol microphysics via neural operator-based modeling

The complexity and the small characteristic scales of aerosol microphysical processes pose a big challenge for accurate and efficient Earth system simulations at regional and global scales. In this work, we construct and evaluate a surrogate model: the aerosol deep operator network (ADON), a physics-inspired dual-net architecture for emulating the aerosol microphysics parameterization suite in the version 2 of the Energy Earth System Model (E3SMv2). The current version of the surrogate model is trained on a dataset comprising 9.8 million samples obtained from a global E3SMv2 simulation with the horizontal resolution of about one degree under cloud-free conditions. Incorporating domain spatial and temporal coordinates, as well as principle components extracted from training data, the dual-net surrogate model effectively captures the intricate representations of aerosol and the relationship with atmospheric state variables, achieving an R-squared score over $$95.7\%$$ for all the lognormal aerosol modes in the extrapolated regime. The validated model provides feature importance of input variables and their impact on the predictive capacity of the surrogate model in relation to the E3SM. The computational cost of online inference time deployed on CPUs and GPUs with lower precisions highlights ADON’s efficiency and potential in robust predictive modeling for large-scale Earth system computations.

Bai, Zhe↗

Investigating property-porosity relationships for micro-architected lattice structures

Micro-architected structures are increasingly valued for their light weight and tunable mechanical properties; this class of material includes sheet structures like the Schöen gyroid and beam lattices like the octet. For design purposes, it is critical to understand how to tune the relative density (RD) to obtain desired mechanical properties. This study investigates the mechanical response of Ti-6Al-4V gyroid structures, spanning a broad range of RDs (0.03–0.90), unit cell sizes (1–4 cm), and sheet thicknesses (0.2–7.8 mm). Results demonstrate that the classical Gibson-Ashby power law scaling between RD and modulus and yield stress does not adequately capture the response over a wide range of RDs, nor does it extrapolate correctly to the fully dense solid. Analytical models, in concert with experimental results and high-fidelity finite element calculations, show that the deviation from Gibson-Ashby reflects a transition from structure-dominated to material-dominated behavior. Distinctions are drawn between scaling relationships and property-porosity models, and an analytical model is proposed to better capture the evolution of mechanical properties with RD. These observations are mirrored in other architected topologies, like the octet, and highlight the importance of understanding the relationship between mechanical properties and geometry for design purposes.

Materials science↗

A hybrid calorimetry-simulation model of mixing enthalpy for molten salt

Calorimetric determination of enthalpies of mixing (ΔH mix ) in multicomponent molten salts is often interpreted using empirical models that lack physically meaningful parameters. However, for improving pyrochemical separation of spent nuclear fuel, where lanthanides are major fission products and critical elements, a deeper thermodynamic understanding of the link between excess thermodynamic properties and solvation structure is critically needed. In this work, we implement a hybrid and physics-informed framework, MIVM+Calorimetry+AIMD, which integrates experimentally measured ΔH mix (via high temperature drop calorimetry) with solvation structures from ab initio molecular dynamics (AIMD). This approach is demonstrated using LaCl 3 mixed with eutectic LiCl-KCl (58 mol% – 42 mol%) at 873 K and 1133 K. MIVM-derived parameters enable extrapolation of excess Gibbs energy and La 3+ activity across compositions. In contrast, direct ΔH mix predictions from AIMD and polarizable ion model simulations deviate significantly. By incorporating experimentally benchmarked solvation structures into an interpretable thermodynamic model, the MIVM+Calorimetry+AIMD formalism achieves higher accuracy and generalizable method for studying molten salts, offering a robust path for understanding and optimizing molten salt chemistry relevant to nuclear fuel cycles and separation science.

Goncharov, Vitaliy G. [Washington State Univ., Pul↗

Learning continuous models for continuous physics

Abstract Dynamical systems that evolve continuously over time are ubiquitous throughout science and engineering. Machine learning (ML) provides data-driven approaches to model and predict the dynamics of such systems. A core issue with this approach is that ML models are typically trained on discrete data, using ML methodologies that are not aware of underlying continuity properties. This results in models that often do not capture any underlying continuous dynamics—either of the system of interest, or indeed of any related system. To address this challenge, we develop a convergence test based on numerical analysis theory. Our test verifies whether a model has learned a function that accurately approximates an underlying continuous dynamics. Models that fail this test fail to capture relevant dynamics, rendering them of limited utility for many scientific prediction tasks; while models that pass this test enable both better interpolation and better extrapolation in multiple ways. Our results illustrate how principled numerical analysis methods can be coupled with existing ML training/testing methodologies to validate models for science and engineering applications.

97 MATHEMATICS AND COMPUTING↗

The quest for superheavy elements and the limit of the periodic table

The borders of the periodic table of the elements and of the chart of nuclides are not set in stone. The desire to explore the properties of atoms and their nuclei in a regime of very large numbers of electrons, protons and neutrons has motivated new experimental facilities to create new elements and nuclides at the limits of atomic number and mass. But the small production rates and short lifetimes of superheavy nuclei and their atoms mean that ‘atom-at-a-time’ studies are the only experimental way to probe them. The physical and chemical data obtained so far, augmented by theoretical calculations, indicate significant deviations from extrapolations from lighter elements and isotopes. This situation raises the following question: how much further can one push the limits of the periodic table? Here, in this Review, we describe the major challenges in the field of the superheavy elements and speculate about future directions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The design space of E(3)-equivariant atom-centred interatomic potentials

Abstract Molecular dynamics simulation is an important tool in computational materials science and chemistry, and in the past decade it has been revolutionized by machine learning. This rapid progress in machine learning interatomic potentials has produced a number of new architectures in just the past few years. Particularly notable among these are the atomic cluster expansion, which unified many of the earlier ideas around atom-density-based descriptors, and Neural Equivariant Interatomic Potentials (NequIP), a message-passing neural network with equivariant features that exhibited state-of-the-art accuracy at the time. Here we construct a mathematical framework that unifies these models: atomic cluster expansion is extended and recast as one layer of a multi-layer architecture, while the linearized version of NequIP is understood as a particular sparsification of a much larger polynomial model. Our framework also provides a practical tool for systematically probing different choices in this unified design space. An ablation study of NequIP, via a set of experiments looking at in- and out-of-domain accuracy and smooth extrapolation very far from the training data, sheds some light on which design choices are critical to achieving high accuracy. A much-simplified version of NequIP, which we call BOTnet (for body-ordered tensor network), has an interpretable architecture and maintains its accuracy on benchmark datasets.

Computer Science↗

Advancing molecular machine learning representations with stereoelectronics-infused molecular graphs

Molecular representation is a critical element in our understanding of the physical world and the foundation for modern molecular machine learning. Previous molecular machine learning models have used strings, fingerprints, global features and simple molecular graphs that are inherently information-sparse representations. However, as the complexity of prediction tasks increases, the molecular representation needs to encode higher fidelity information. This work introduces a new approach to infusing quantum-chemical-rich information into molecular graphs via stereoelectronic effects, enhancing expressivity and interpretability. Learning to predict the stereoelectronics-infused representation with a tailored double graph neural network workflow enables its application to any downstream molecular machine learning task without expensive quantum-chemical calculations. We show that the explicit addition of stereoelectronic information substantially improves the performance of message-passing two-dimensional machine learning models for molecular property prediction. We show that the learned representations trained on small molecules can accurately extrapolate to much larger molecular structures, yielding chemical insight into orbital interactions for previously intractable systems, such as entire proteins, opening new avenues of molecular design. Finally, we have developed a web application (simg.cheme.cmu.edu) where users can rapidly explore stereoelectronic information for their own molecular systems.

Boiko, Daniil A↗

Differentiable modelling to unify machine learning and physical models for geosciences

Process-based modelling offers interpretability and physical consistency in many domains of geosciences but struggles to leverage large datasets efficiently. Machine-learning methods, especially deep networks, have strong predictive skills yet are unable to answer specific scientific questions. Here, in this Perspective, we explore differentiable modelling as a pathway to dissolve the perceived barrier between process-based modelling and machine learning in the geosciences and demonstrate its potential with examples from hydrological modelling. ‘Differentiable’ refers to accurately and efficiently calculating gradients with respect to model variables or parameters, enabling the discovery of high-dimensional unknown relationships. Differentiable modelling involves connecting (flexible amounts of) prior physical knowledge to neural networks, pushing the boundary of physics-informed machine learning. It offers better interpretability, generalizability, and extrapolation capabilities than purely data-driven machine learning, achieving a similar level of accuracy while requiring less training data. Additionally, the performance and efficiency of differentiable models scale well with increasing data volumes. Under data-scarce scenarios, differentiable models have outperformed machine-learning models in producing short-term dynamics and decadal-scale trends owing to the imposed physical constraints. Differentiable modelling approaches are primed to enable geoscientists to ask questions, test hypotheses, and discover unrecognized physical relationships. Future work should address computational challenges, reduce uncertainty, and verify the physical significance of outputs.

58 GEOSCIENCES↗

A Permanently Porous Chalcogen-Bonded Organic Framework

The nature of connectivity between constituent atomic or molecular building blocks is fundamental in shaping the properties and functionality of materials. The extrapolation of emergent interatomic interactions to enable functional materials has driven transformative technological advancements. However, the bonding interactions used in material design have been largely static since the emergence of dynamic covalent chemistry ~30 years ago. Here we demonstrate that non-covalent chalcogen bonding (Ch-bonding) is a distinct mode of interatomic connectivity for constructing functional materials by design. This is established by leveraging self-complementary assembly of 1,2,5-telluradiazole moieties to construct a honeycomb-type permanently porous Ch-bonded organic framework, assembled and stabilized solely through non-covalent Te...N contacts. Empirical and computational studies of electronic structure, structural healing and lattice dynamics highlight the pi-type electronic communication, controlled assembly and modulated lattice dynamics in Trip3Tez-I arising directly from the unique nature of the Te...N Ch-bonding that holds substantial implications for next generation crystalline semiconductors. In addition to introducing a distinct class of permanently porous frameworks, this work establishes Ch-bonding as a programmable molecular tool for constructing functional materials with distinct properties.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Destructive reverse bias pinning in perovskite/silicon tandem solar modules caused by perovskite hysteresis under dynamic shading

In this work, we demonstrate how perovskite hysteresis can result in permanent reductions in power output in perovskite/silicon tandem modules—including irreversible hotspot-induced damage—from only brief periods of shading. We show that reverse bias events in which a perovskite cell is biased above a threshold voltage—which in this work we find to be as low as –1.1 V—produces a temporary reduction in power output that is of sufficient magnitude to keep the cell pinned in reverse bias after the shading event ends. As a hysteretic phenomena, this crucial failure mode may be overlooked by static models of perovskite-based solar cells. Higher reverse bias voltages exacerbate the temporary reduction in short-circuit photocurrent, which is also sensitive to the level of illumination under reverse bias. Numerical device modelling demonstrates that this effect is consistent with our understanding of perovskite hysteresis as a consequence of mobile ion-electron coupling controlling rates of non-radiative recombination over time. Measurements of the dynamic response of single-junction perovskite cells are extrapolated to two-terminal and four-terminal perovskite/silicon tandem module modelling. We validate these models with measurements from an equivalent electronic circuit that represents a two-terminal perovskite–silicon tandem mini module. Two module-level solutions are discussed that address this issue, which includes increasing the number of bypass diodes and choosing better suited silicon bottom cells with higher shunter resistance in two-terminal tandem modules.

14 SOLAR ENERGY↗

Electron leakage through heterogeneous LiF on lithium–metal battery anodes

The solid–electrolyte interphase (SEI) that forms on lithium ion battery (LIB) anodes prevents degradation-causing transfer of electrons to the electrolyte. Grain boundaries (GBs) between different SEI components, like LiF, have been suggested to accelerate Li+ transport. However, using the non-equilibrium Green's function technique with density functional theory (NEGF-DFT), we find that GBs enhance electron tunneling in thin LiF films by 1–2 orders of magnitude, depending on the bias. Extrapolating to thicker films using the Wentzel–Kramers–Brillouin (WKB) method emphasizes that safer batteries require passivation of GBs in the SEI.

25 ENERGY STORAGE↗

Thermodynamic and kinetic properties of layered-CaCo2O4 for the Ca-ion batteries: a systematic first-principles study

One of the more promising directions in multivalent energy storage is systems based on Ca ion intercalation due to the potential for high voltage and capacity. A major challenge for enabling such a battery is to find cathode materials capable of fast ionic diffusion and reversible insertion of Ca ions. Here, on the basis of first-principles calculations, we have demonstrated that layered CaCo 2 O 4 exhibits favorable thermodynamic and kinetic properties that should enable topotactic Ca ion intercalation reactions. The P3-type layered Ca x Co 2 O 4 (0 < x < 1) with either of space groups of P 1 or P 2 1 / m are stable at multiple Ca concentrations and show a smooth voltage plateau higher than 3 V up to X = 0.5. The energy barriers of the single Ca ion migration are as low as 0.36 eV and 0.27 eV at the dilute and high vacancy concentration limits, respectively. Therefore, although varying the vacancy environments of the diffusing atom influences the migration barriers, they do not exceed 0.6 eV. Stochastic analysis of Ca hopping events performed by ab initio molecular dynamics (AIMD) simulation has shown that the migration barriers are lower than 0.32 eV. Therefore, the Ca diffusivity at room temperature extrapolated from the AIMD results is comparable to Li diffusivity (>10 -10 cm 2 s -1 ) in conventional Li cathode materials, suggesting the feasibility of layered Ca x Co 2 O 4 as multivalent cathode materials. Finally, the structural factors that enable fast diffusion are discussed.

25 ENERGY STORAGE↗

Experimental and computational investigation of the bond energy of thorium dicarbonyl cation and theoretical elucidation of its isomerization mechanism to the thermodynamically most stable isomer, thorium oxide ketenylidene cation, OTh + CCO

Collision-induced dissociation (CID) of [Th,2C,2O] + with Xe is performed using a guided ion beam tandem mass spectrometer (GIBMS). The only products observed are ThCO + and Th + by sequential loss of CO ligands. The experimental findings and theoretical calculations support that the structure of [Th,2C,2O] + is the bent homoleptic thorium dicarbonyl cation, Th + (CO) 2 , having quartet spin, which is both thermodynamically and kinetically stable enough in the gas phase to be observed in our GIBMS instrument. Analysis of the kinetic energy-dependent cross sections for this CID reaction yields the first experimental determination of the bond dissociation energy (BDE) of (CO)Th + –CO at 0 K as 1.05 ± 0.09 eV. A theoretical BDE calculated at the CCSD(T) level with cc-pVXZ (X = T and Q) basis sets and a complete basis set (CBS) extrapolation is in very good agreement with the experimental result. Although the doublet spin bent thorium oxide ketenylidene cation, OTh + CCO, is calculated to be the most thermodynamically stable structure, it is not observed in our experiment where [Th,2C,2O] + is formed by association of Th + and CO in a direct current discharge flow tube (DC/FT) ion source. Potential energy profiles of both quartet and doublet spin are constructed to elucidate the isomerization mechanism of Th + (CO) 2 to OTh + CCO. The failure to observe OTh + CCO is attributed to a barrier associated with C–C bond formation, which makes OTh + CCO kinetically inaccessible under our experimental conditions. Furthermore, chemical bonding patterns in low-lying states of linear and bent Th + (CO) 2 and OTh + CCO isomers are also investigated.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Transfer learning of memory kernels for transferable coarse-graining of polymer dynamics

The present work concerns the transferability of coarse-grained (CG) modeling in reproducing the dynamic properties of the reference atomistic systems across a range of parameters. In particular, we focus on implicit-solvent CG modeling of polymer solutions. The CG model is based on the generalized Langevin equation, where the memory kernel plays the critical role in determining the dynamics in all time scales. Thus, we propose methods for transfer learning of memory kernels. The key ingredient of our methods is Gaussian process regression. By integration with the model order reduction via proper orthogonal decomposition and the active learning technique, the transfer learning can be practically efficient and requires minimum training data. Through two example polymer solution systems, we demonstrate the accuracy and efficiency of the proposed transfer learning methods in the construction of transferable memory kernels. The transferability allows for out-of-sample predictions, even in the extrapolated domain of parameters. Built on the transferable memory kernels, the CG models can reproduce the dynamic properties of polymers in all time scales at different thermodynamic conditions (such as temperature and solvent viscosity) and for different systems with varying concentrations and lengths of polymers.

Ma, Zhan↗

Understanding multi-scale battery degradation with a macro-to-nano zoom through its hierarchy

Lithium-ion batteries (LIBs) feature structural and chemical complexities across a broad range of length scales. It is the hierarchy of the battery structure that determines its functionality. An in-depth understanding of the battery function, degradation, and failure mechanisms requires a thorough and systematic investigation from structural, chemical, mechanical, and dynamic perspectives. Here we present a macro-to-nano zoom through the hierarchy of a commercial 18650 type LIB using a suite of state-of-the-art X-ray microscopy techniques. Damage, deformation, and heterogeneity at different length scales are visualized and are associated with different degradation phenomena and mechanisms. Our results highlight the importance of the mechanical properties of the cathode material, which could impact both the immediate and the long-term cell behaviors significantly. While this study focuses on a commercial lithium cell with a standard configuration, our findings can be extrapolated and are applicable to the development of next-generation energy storage technology, e.g. solid-state batteries.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Sub 20 cm -1 computational prediction of the CH bond energy – a case of systematic error in computational thermochemistry

Here, the bond dissociation energy of methylidyne, D 0 (CH), is studied using an improved version of the High-Accuracy Extrapolated ab initio Thermochemistry (HEAT) approach as well as the Feller–Peterson–Dixon (FPD) model chemistry. These calculations, which include basis sets up to nonuple (aug-cc-pCV9Z) quality, are expected to be capable of providing results substantially more accurate than the ca. 1 kJ mol -1 level that is characteristic of standard high-accuracy protocols for computational thermochemistry. The calculated 0 K CH bond energy (27 954 ± 15 cm -1 for HEAT and 27 956 ± 15 cm -1 for FPD), along with equivalent treatments of the CH ionization energy and the CH + dissociation energy (85 829 ± 15 cm -1 and 32 946 ± 15 cm -1 , respectively), were compared to the existing benchmarks from Active Thermochemical Tables (ATcT), uncovering an unexpected difference for D 0 (CH). This has prompted a detailed reexamination of the provenance of the corresponding ATcT benchmark, allowing the discovery and subsequent correction of a systematic error present in several published high-level calculations, ultimately yielding an amended ATcT benchmark for D 0 (CH). Finally, the current theoretical results were added to the ATcT Thermochemical Network, producing refined ATcT estimates of 27 957.3 ± 6.0 cm -1 for D0(CH), 32 946.7 ± 0.6 cm -1 for D 0 (CH + ), and 85 831.0 ± 6.0 cm -1 for IE(CH).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Theoretical study of the CO 2 –O 2 van der Waals complex: potential energy surface and applications

A four-dimensional-potential energy surface (4D-PES) of the atmospherically relevant carbon dioxide–oxygen molecule (CO 2 –O 2 ) van der Waals complex is mapped using the ab initio explicitly correlated coupled cluster method with single, double, and perturbative triple excitations (UCCSD(T)-F12b), and extrapolation to the complete basis set (CBS) limit using the cc-pVTZ-F12/cc-pVQZ-F12 bases and the l –3 formula. An analytic representation of the 4D-PES was fitted using the method of interpolating moving least squares (IMLS). These calculations predict that the most stable configuration of CO 2 –O 2 complex corresponds to a planar slipped-parallel structure with a binding energy of V ~ –243 cm –1 . Another isomer is found on the PES, corresponding to a non-planar cross-shaped structure, with V ~ –218 cm –1 . The transition structure connecting the two minima is found at V ~ –211 cm –1 . We also performed comparisons with some CO 2 –X van der Waals complexes. Moreover, we provide a SAPT analysis of this molecular system. Then, we discuss the complexation induced shifts of CO 2 and O 2 . Afterwards, this new 4D-PES is employed to compute the second virial coefficient including temperature dependence. A comparison between quantities obtained in our calculations and those from experiments found close agreement attesting to the high quality of the PES and to the importance of considering a full description of the anisotropic potential for the derivation of thermophysical properties of CO 2 –O 2 mixtures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗