Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Physics - Plasma physics , Mathematics and Computing, Physics”

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 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

A theoretical study of the upper bound of surface elevation variance in the Phillips initial stage during wind-wave generation

The resonance mechanism in the initial of wind-wave generation proposed by Phillips is a foundation of wind-wave generation theory, but a precise theoretical quantification of wave energy growth in this initial stage has not been obtained yet after more than six decades of research. In this study, we aim to address this knowledge gap by developing an analytical approach based on a novel complex analysis method to theoretically investigate the temporal evolution of the wave energy in the Phillips initial stage. We quantitatively derive and analyse the growth behaviour of the surface wave energy and obtain an analytical solution for its upper bound. Our result highlights the crucial effects of surface tension. Because the phase velocity of gravity–capillary waves has a minimal value at a critical wavenumber, gravity–capillary waves and gravity waves (which neglect surface tension) exhibit distinct resonance curve properties and wave energy growth behaviours. For gravity waves, the resonance curve extends indefinitely; for gravity–capillary waves, it either forms a finite-length curve or does not exist, depending on the wind speed. The leading-order term of the upper-bound solution of the energy of gravity waves increases linearly over time, while for gravity–capillary waves, the term increases linearly over time under strong wind conditions but remains finite under weak wind conditions. This theoretical study provides an analytical framework for the generation of wind-waves in the Phillips initial stage, which may inspire further theoretical, numerical and experimental research.

58 GEOSCIENCES↗
DOE OSTI2025

Confusion-Driven Machine Learning of Structural Phases of a Flexible, Magnetic Stockmayer Polymer

We use a semisupervised, neural-network-based machine learning technique, the confusion method, to investigate structural transitions in magnetic polymers, which we model as chains of magnetic colloidal nanoparticles characterized by dipole–dipole and Lennard-Jones interactions. As input for the neural network, we use the particle positions and magnetic dipole moments of equilibrium polymer configurations, which we generate via replica-exchange Wang–Landau simulations. We demonstrate that by measuring the classification accuracy of neural networks, we can effectively identify transition points between multiple structural phases without any prior knowledge of their existence or location. We corroborate our findings by investigating relevant conventional order parameters. Our study furthermore examines previously unexplored low-temperature regions of the phase diagram, where we find new structural transitions between highly ordered helicoidal polymer configurations.

36 MATERIALS SCIENCE↗
DOE OSTI2025

First-Principles-Based Study of the Decomposition of Phenol and Hydroquinone on Pt(111) Combined with Quantitative Information from XPS Spectra to Address the Impact of Coverage and Number of Hydroxyl Functional Groups

A combined first-principles-based and experimental X-ray photoelectron spectroscopy approach was used to investigate the thermal decomposition of two model biofuel compounds, phenol and hydroquinone, on Pt(111) at both low and high coverages. The DFT-based approach yields adsorption geometries and energies, activation barriers and core-level binding energy shifts for C 1s and O 1s. Increasing the coverage in the theoretical model leads to slight shifts in core-level binding energies─toward higher values for C 1s and lower values for O 1s. It also alters the energy profiles of the decomposition reaction pathway, resulting in weaker adsorption energies and changes in both reaction and activation barriers. At low temperatures, we observe a multilayer for phenol and hydroquinone upon adsorption, with desorption occurring at 200 and 270 K, respectively. Following desorption of the multilayer, decomposition proceeds via initial O–H bond scission, followed by two parallel pathways involving either C–H or C–C bond scission, whereby in the case of phenol C–H bond scission occurs first. Here, we further provide characteristic core level binding energies by theoretical calculations that are subsequently used in experimental analyses, establishing a reference database for key spectra of phenolic functionalities applicable to a range of catalytic reactions.

09 BIOMASS FUELS↗
DOE OSTI2025

How Solvent Mediates the Catalytic Fate of 2-Methoxyphenol in Hydrodeoxygenation Catalysis?

Kinetic and isotopic assessments combined with theoretical calculations shed light on the mechanistic differences in terms of the number of hydrogen addition steps and the electronic charges of reactive hydrogen in the catalytic sequence leading to the Ar−OCH 3 scission of 2-methoxyphenol occurring at various interfaces of vapor−Ru, cyclohexane−Ru, and water−Ru. At vapor−Ru and cyclohexane−Ru interfaces, the initial hydrogen adatom (H*) additions on the phenyl ring disrupt its aromaticity, as a step required to occur before kinetically relevant Ar−OCH 3 scission on an ensemble of Ru sites. The number of H* addition events prerequisite to the formation of the kinetically relevant Ar−OCH 3 scission transition state, however, is larger at the cyclohexane−Ru interface than those at the vapor−Ru and water−Ru interfaces. The preferential solvation of 2-methoxyphenol derived reaction precursors and H* weakens their adsorption on Ru sites, which leads to a decrease in the activation barrier for the third H* addition step and in turn a larger number of H* additions. At the water−Ru interface, water as a polar protic solvent promotes the Ar−OCH 3 scission by enabling a mechanistically distinct catalytic route an adjacent water molecule from the solvent layer assists with the intramolecular proton shuttling from the hydroxyl group (−OH) to the vicinal methoxy group (−OCH 3 ) of 2-methoxyphenol, evolving a highly charged transition state that is further stabilized by the water solvent layer, thereby leading to Ar−OCH 3 scission rate constants that are much higher than those at the vapor−Ru interface and, when comparing at the same H 2 chemical potential, higher Ar−OCH 3 scission turnovers. Furthermore, these distinct mechanistic details of the Ar−OCH 3 scission at various reaction interfaces, inferred from experiment and theory, illustrate how a solvent either solvates reaction intermediates and transition states or directly participates in the formation of kinetically relevant transition states, leading to changes in the free energy landscape, thus rewriting the catalytic fates in hydrodeoxygenation reactions.

09 BIOMASS FUELS↗
DOE OSTI2025

Tuning Catalytic Reactivity via Wetting Control through Oxygen Vacancies: Ru Clusters on Anatase TiO 2 and CeO 2 Supports

The shape of supported metal particles regulates their catalytic reactivity and is determined by the degree of wetting between the metal particle and the support surface. Flattened particles that wet support surfaces were reported in various catalytic systems, particularly in the subnanometer size regime. Such consequential metal–support wetting phenomena are poorly understood, and methods to study them on powder catalysts under realistic conditions are lacking. Here, we investigate the size-dependent wetting behaviors of Ru particles on two reducible-oxide supports, anatase TiO 2 (TiO 2 -A) and CeO 2 , under reducing catalytic conditions. X-ray absorption spectroscopy (XAS), low-energy ion scattering (LEIS), and density functional theory (DFT) are combined to determine the shape of Ru particles. Ru particles remain three-dimensional without wetting the TiO 2 -A support within the coverage range studied (0.06–0.98 Ru nm –2 ). In contrast, at low coverages (<0.25 Ru nm –2 ), Ru wets the CeO 2 support to form flat, disordered structures. The higher wettability of CeO 2 than TiO 2 -A is attributed to oxygen vacancies in the near-surface region. The shape difference between small Ru particles or clusters on the two supports leads to drastically contrasting catalytic reactivities in polyolefin hydrogenolysis, despite similar diameters. As a result, this work highlights the implications of metal–support wetting, or cluster shape, on catalytic behaviors of small metal clusters, while establishing the foundation for future systematic studies of such a phenomenon in realistic systems, by delivering a multitechnique methodology and revealing governing fundamental principles.

42 ENGINEERING↗
DOE OSTI2021

Correspondence between neuroevolution and gradient descent

Abstract We show analytically that training a neural network by conditioned stochastic mutation or neuroevolution of its weights is equivalent, in the limit of small mutations, to gradient descent on the loss function in the presence of Gaussian white noise. Averaged over independent realizations of the learning process, neuroevolution is equivalent to gradient descent on the loss function. We use numerical simulation to show that this correspondence can be observed for finite mutations, for shallow and deep neural networks. Our results provide a connection between two families of neural-network training methods that are usually considered to be fundamentally different.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Learning in continuous action space for developing high dimensional potential energy models

Reinforcement learning (RL) approaches that combine a tree search with deep learning have found remarkable success in searching exorbitantly large, albeit discrete action spaces, as in chess, Shogi and Go. Many real-world materials discovery and design applications, however, involve multi-dimensional search problems and learning domains that have continuous action spaces. Exploring high-dimensional potential energy models of materials is an example. Traditionally, these searches are time consuming (often several years for a single bulk system) and driven by human intuition and/or expertise and more recently by global/local optimization searches that have issues with convergence and/or do not scale well with the search dimensionality. Here, in a departure from discrete action and other gradient-based approaches, we introduce a RL strategy based on decision trees that incorporates modified rewards for improved exploration, efficient sampling during playouts and a “window scaling scheme" for enhanced exploitation, to enable efficient and scalable search for continuous action space problems. Using high-dimensional artificial landscapes and control RL problems, we successfully benchmark our approach against popular global optimization schemes and state of the art policy gradient methods, respectively. We demonstrate its efficacy to parameterize potential models (physics based and high-dimensional neural networks) for 54 different elemental systems across the periodic table as well as alloys. We analyze error trends across different elements in the latent space and trace their origin to elemental structural diversity and the smoothness of the element energy surface. Broadly, our RL strategy will be applicable to many other physical science problems involving search over continuous action spaces.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Phonon thermal Hall effect in a metallic spin ice

It has become common knowledge that phonons can generate thermal Hall effect in a wide variety of materials, although the underlying mechanism is still controversial. We study longitudinal κ xx and transverse κ xy thermal conductivity in Pr 2 Ir 2 O 7 , which is a metallic analog of spin ice. Despite the presence of mobile charge carriers, we find that both κ xx and κ xy are dominated by phonons. A T/H scaling of κ xx unambiguously reveals that longitudinal heat current is substantially impeded by resonant scattering of phonons on paramagnetic spins. Upon cooling, the resonant scattering is strongly affected by a development of spin ice correlation and κ xx deviates from the scaling in an anisotropic way with respect to field directions. Strikingly, a set of the κ xx and κ xy data clearly shows that κ xy correlates with κ xx in its response to magnetic field including a success of the T/H scaling and its failure at low temperature. This remarkable correlation provides solid evidence that an indispensable role is played by spin-phonon scattering not only for hindering the longitudinal heat conduction, but also for generating the transverse response.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Improved machine learning algorithm for predicting ground state properties

Finding the ground state of a quantum many-body system is a fundamental problem in quantum physics. In this work, we give a classical machine learning (ML) algorithm for predicting ground state properties with an inductive bias encoding geometric locality. The proposed ML model can efficiently predict ground state properties of an n-qubit gapped local Hamiltonian after learning from only $\mathcal{O}$(log(n)) data about other Hamiltonians in the same quantum phase of matter. This improves substantially upon previous results that require $\mathcal{O}$(n c ) data for a large constant c. Furthermore, the training and prediction time of the proposed ML model scale as $\mathcal{O}$(n log n) in the number of qubits n. Numerical experiments on physical systems with up to 45 qubits confirm the favorable scaling in predicting ground state properties using a small training dataset.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Efficient and universal characterization of atomic structures through a topological graph order parameter

Abstract A graph-based order parameter, based on the topology of the graph itself, is introduced for the characterization of atomistic structures. The order parameter is universal to any material/chemical system and is transferable to all structural geometries. Four sets of data are used to validate both the generalizability and accuracy of the algorithm: (1) liquid lithium configurations spanning up to 300 GPa, (2) condensed phases of carbon along with nanotubes and buckyballs at ambient and high temperature, (3) a diverse set of aluminum configurations including surfaces, compressed and expanded lattices, point defects, grain boundaries, liquids, nanoparticles, all at nonzero temperatures, and (4) eleven niobium oxide crystal phases generated with ab initio molecular dynamics. We compare our proposed method to existing, state-of-the-art methods for the cases of aluminum and niobium oxide. Our order parameter uniquely classifies every configuration and outperforms all studied existing methods, opening the door for its use in a multitude of complex application spaces that can require fine structure-level characterization of atomistic graphs.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Sampling frequency thresholds for the quantum advantage of the quantum approximate optimization algorithm

We compare the performance of the Quantum Approximate Optimization Algorithm (QAOA) with state-of-the-art classical solvers Gurobi and MQLib to solve the MaxCut problem on 3-regular graphs. We identify the minimum noiseless sampling frequency and depth p required for a quantum device to outperform classical algorithms. There is potential for quantum advantage on hundreds of qubits and moderate depth with a sampling frequency of 10 kHz. We observe, however, that classical heuristic solvers are capable of producing high-quality approximate solutions in linear time complexity. In order to match this quality for large graph sizes N, a quantum device must support depth p > 11. Additionally, multi-shot QAOA is not efficient on large graphs, indicating that QAOA p ≤ 11 does not scale with N. These results limit achieving quantum advantage for QAOA MaxCut on 3-regular graphs. Other problems, such as different graphs, weighted MaxCut, and 3-SAT, may be better suited for achieving quantum advantage on near-term quantum devices.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Degree of Quantumness in Quantum Synchronization

We introduce the concept of degree of quantumness in quantum synchronization, a measure of the quantum nature of synchronization in quantum systems. Following techniques from quantum information, we propose the number of non-commuting observables that synchronize as a measure of quantumness. This fgure of merit is compatible with already existing synchronization measurements, and it captures diferent physical properties. We illustrate it in a quantum system consisting of two weakly interacting cavity-qubit systems, which are coupled via the exchange of bosonic excitations between the cavities. Moreover, we study the synchronization of the expectation values of the Pauli operators and we propose a feasible superconducting circuit setup. Finally, we discuss the degree of quantumness in the synchronization between two quantum van der Pol oscillators.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Three-State Majority-vote Model on Scale-Free Networks and the Unitary Relation for Critical Exponents

Abstract We investigate the three-state majority-vote model for opinion dynamics on scale-free and regular networks. In this model, an individual selects an opinion equal to the opinion of the majority of its neighbors with probability 1 − q , and different to it with probability q . The parameter q is called the noise parameter of the model. We build a network of interactions where z neighbors are selected by each added site in the system, a preferential attachment network with degree distribution k − λ , where λ = 3 for a large number of nodes N . In this work, z is called the growth parameter. Using finite-size scaling analysis, we obtain that the critical exponents $$\beta /\bar{\nu }$$ β / ν ¯ and $$\gamma /\bar{\nu }$$ γ / ν ¯ associated with the magnetization and the susceptibility, respectively. Using Monte Carlo simulations, we calculate the critical noise parameter q c as a function of z for the scale-free networks and obtain the phase diagram of the model. We find that the critical exponents add up to unity when using a special volumetric scaling, regardless of the dimension of the network of interactions. We verify this result by obtaining the critical noise and the critical exponents for the two and three-state majority-vote model on cubic lattice networks.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Re-evaluation of the TSL for Yttrium Hydride

Yttrium hydride (YH x ) is of interest as a high-temperature moderator material in advanced nuclear reactor systems because of its superior ability to retain hydrogen at elevated temperatures. Thermal neutron scattering laws (TSL) for hydrogen bound in yttrium hydride (H-YH 2 ) and yttrium bound in yttrium hydride (Y-YH 2 ) were previously evaluated by Naval Nuclear Laboratory using the ab initio approach and released in ENDF/B-VIII.0. In that work, density functional theory, incorporating the generalized gradient approximation (GGA) for the exchange-correlation energy, was used to simulate the face-centered cubic structure of YH 2 and calculate the interatomic Hellmann-Feynman forces for a 2×2×2 supercell containing 96 atoms. Lattice dynamics calculations using PHONON were used to determine the phonon density of states. The calculated phonon density of states for H and Y in YH 2 were then used to prepare H-YH 2 and Y-YH 2 TSL evaluations, in the incoherent approximation, using the LEAPR module of NJOY. In addition, elastic scattering was assumed to be incoherent for both H and Y. While the incoherent elastic scattering approximation is appropriate for H-YH 2 , it introduces an undesirable approximation for Y-YH 2 . In this work, we re-evaluate the TSL for Y-YH 2 using FLASSH (Full Law Analysis Scattering System Hub). Y-YH 2 is evaluated using the FLASSH generalized coherent elastic scattering capability in order to capture the Bragg peaks associated with the YH 2 crystal structure which were neglected in the prior NJOY-based evaluation due to limitations in LEAPR. An experimental approach to validate the Y-YH 2 TSL using neutron transmission measurements is discussed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗