Search NASASearch

SEARCH · Search NASA

Results for “Elastic Search”

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 55 records · Page 3

Models, In situ, and Remote Sensing of Aerosols (MIRA)

There is a natural partitioning of scientific interest amongst three specialties of aerosol re-search: modeling, in situ measurements, and remote sensing. The broader aerosol community benefits when these groups interact, and this strengthens overall scientific impact on climate and air quality research and predictions. The new MIRA working group establishes a forum for identifying collaborations and improving discussions amongst specialties and across regional boundaries. One area of keen interest is uniting satellite and ground-based lidar groups with other aerosol disciplines. Elastic backscatter lidars depend upon a priori knowledge of aerosol properties to convert measured lidar profiles into aerosol extinction profiles. Acquiring additional insight on aerosol properties and transport is highly valuable to these lidar groups to improve the data quality and aid in their scientific interpretation. Another area of interest is to facilitate the incorporation of measurements into global aerosol models. Modelers need aerosol optical look-up tables that enable quick conversions of hydrated (and dry) polydisperse size distributions into aerosol optical properties (extinction, scattering, etc.), but many of the existing tables are based upon outdated measurements. Thus, MIRA is building the Tables of Aerosol Optics (TAO), which is a community collection of aerosol optical calculations. This expands on the historical efforts of Shettle and Fenn, d’Almeida, GADS, OPAC, etc., except that TAO seeks continual input from the community. Thus, as aerosol measurement and computational techniques advance, so does the TAO collective. Since TAO is a community collective, it is expected that users will optionally upload computations for aerosol type as well as computations for all of the traditional aerosol species(like amm sulfate, amm nitrate, organics, etc.). The near-term purpose of the MIRA working group is to: *Characterize regional aerosol lidar ratios to support improvements of aerosol extinction profiles and understanding of aerosol typing. *Create TAO --a community cooperative of aerosol optical tables. *Facilitate international communications between aerosol measurement and modeling groups. *Encourage the use of regional knowledge to develop and improve remote sensing techniques for current and future backscatter lidars located in space. *Enable and foster communication between the scientists who run global aerosol models and scientists who analyze space-based lidar data. We will discuss how a collaborative aerosol working group will be organized on this topic. Those who are interested in MIRA can sign up for the MIRA email listserver at https://forms.gle/qdbCnngzNJc5YJi57.

Greg Schuster

Models, In situ, and Remote sensing of Aerosols (MIRA): Formation of an International Working Group

There is a natural partitioning of scientific interest amongst three specialties of aerosol re- search: modeling, in situ measurements, and remote sensing. The broader aerosol community benefits when these groups interact, and this strengthens overall scientific impact on climate and air quality research and predictions. The new MIRA working group establishes a forum for identifying collaborations and improving discussions amongst specialties and across regional boundaries. One area of keen interest is uniting satellite and ground-based lidar groups with other aerosol disciplines. Elastic backscatter lidars depend upon a priori knowledge of aerosol properties to convert measured lidar profiles into aerosol extinction profiles. Acquiring additional insight on aerosol properties and transport is highly valuable to these lidar groups to improve the data quality and aid in their scientific interpretation. Another area of interest is to facilitate the incorporation of measurements into global aerosol models. Modelers need aerosol optical look-up tables that enable quick conversions of hydrated (and dry) polydisperse size distributions into aerosol optical properties (extinction, scattering, etc.), but many of the existing tables are based upon outdated measurements. Thus, MIRA is building the Tables of Aerosol Optics (TAO), which is a community collection of aerosol optical calculations. This expands on the historical efforts of Shettle and Fenn, d’Almeida, GADS, OPAC, etc., except that TAO seeks continual input from the community. Thus, as aerosol measurement and computational techniques advance, so does the TAO collective. Since TAO is a community collective, it is expected that users will optionally upload computations for aerosol type as well as computations for all of the traditional aerosol species (like amm sulfate, amm nitrate, organics, etc.). The near-term purpose of the MIRA working group is to: * Characterize regional aerosol lidar ratios to support improvements of aerosol extinction profiles and understanding of aerosol typing. * Create TAO -- a community cooperative of aerosol optical tables. * Facilitate international communications between aerosol measurement and modeling groups. * Encourage the use of regional knowledge to develop and improve remote sensing techniques for current and future backscatter lidars located in space. * Enable and foster communication between the scientists who run global aerosol models and scientists who analyze space-based lidar data. We will discuss how a collaborative aerosol working group will be organized on this topic. Send email to calipso_v5alr-join@lists.nasa.gov with the word `subscribe' in the subject line to join MIRA.

Gregory Schuster

MHOST: An efficient finite element program for inelastic analysis of solids and structures

An efficient finite element program for 3-D inelastic analysis of gas turbine hot section components was constructed and validated. A novel mixed iterative solution strategy is derived from the augmented Hu-Washizu variational principle in order to nodally interpolate coordinates, displacements, deformation, strains, stresses and material properties. A series of increasingly sophisticated material models incorporated in MHOST include elasticity, secant plasticity, infinitesimal and finite deformation plasticity, creep and unified viscoplastic constitutive model proposed by Walker. A library of high performance elements is built into this computer program utilizing the concepts of selective reduced integrations and independent strain interpolations. A family of efficient solution algorithms is implemented in MHOST for linear and nonlinear equation solution including the classical Newton-Raphson, modified, quasi and secant Newton methods with optional line search and the conjugate gradient method.

Nakazawa, S.

Hierarchical Reinforcement Learning of a Short-Range Bond-Order Potential for Silica: Analytic Embedding of Coordination with Classical Efficiency

Reinforcement learning (RL) has recently emerged as a data-efficient strategy to parametrize short-range interatomic potentials. Building on our past RL optimization of pairwise silica models, we extend the framework to a bond-order (Tersoff-type) potential that provides an analytic embedding of local coordination through a three-body term. A hierarchical RL workflow combining continuous-action Monte Carlo Tree Search and property-based rewards efficiently explores the 26-dimensional parameter space, sequentially optimizing lattice parameters, densities, angles, and cohesive energies of 21 silica polymorphs. The resulting models, Q-Tersoff and ML-Tersoff, reproduce the energetic ordering of low-energy phases and capture the angular correlations and amorphous structure factors of silica with improved fidelity over pairwise force fields, while remaining orders of magnitude faster than high-dimensional machine-learned potentials. Both models underperform for elastic constants and high-energy frameworks, delineating the limits of the current analytic form. The approach establishes a general and interpretable route to angle-aware, short-range potentials that bridge physics-based and machine-learned descriptions of silicate materials.

36 MATERIALS SCIENCE

Estimating lithospheric properties at Atla Regio, Venus

Magellan spehrical harmonic gravity and topography models are used to estimate lithospheric properties at Alta Regio, Venus, a proposed hotspot with dynamic support from mantle plume(s). Global spherical harmonic and local representations of the gravity field share common properties in the Atla region interms of their spectral behavior over a wavelength band from approximately 2100 to approximately 700 km. The estimated free-air admittance spectrum displays a rather featureless long-wavelength portion followed by a sharp rise at wavelengths shorter than about 1000 km. This sharp rise requires significant flexural support of short-wavelength structures. The Bouguer coherence also displays a sharp drop in this wavelength band, indicating a finite flexural rigidity of the lithosphere. A simple model for lithospheric loading from above and below is introduced (D. W. Forsyth, 1985) with four parameters: f, the ratio of bottom loading to top loading; z(sub m), crustal thickness; z(sub l) depth to bottom loading source; and T(sub e) elastic lithosphere thickness. A dual-mode compensation model is introduced in which the shorter wavelengths (lambda approximately less than 1000 km) might be explained best by a predominance of top loading by the large shield volcanoes Maat Mons, Ozza Mons, and Sapas Mons, and the longer wavelengths (lambda approximately greater than 1500 km) might be explained best by a deep depth of compensation, possibly representing bottom loading by a dynamic source. A Monte Carlo inversion technique is introduced to thoroughly search out the four-space of the model parameters and to examine parameter correlation in the solutions. Venus either is a considerabe deficient in heat sources relative to Earth, or the thermal lithosphere is overthickened in response to an earlier episode of significant heat loss from the planet.

Phillips, Roger J.

Possible flexural signatures around Olympus and Ascraeus Montes, Mars

The effective elastic thickness (h) of the lithosphere provides a measure of the thermal and mechanical state of a planet's shallow interior. An estimate of h in the vicinity of a feature that constitutes a load on a planetary surface can be determined from the flexural response of the lithosphere to the load. This approach has been applied to Mars by calculating radial stresses associated with lithospheric flexure associated with surface loads, and comparing the results to the positions of circumferential graben surrounding the major Martian shield volcanoes and mascon basins. However, many prominent surface loads on Mars, most notably the Olympus Mons volcano, do not exhibit flexural graben. In these instances application of the above method can provide only a lower limit of effective elastic thickness. An alternative method of determining h is to calculate the vertical displacements associated with the flexural loading and to compare the amplitude and shape of the flexural profile to observed topography. This method has not been applied to flexural problems on Mars because of the poor resolution of Martian topographic data. However, previous analyses have shown that the lithosphere around major volcanic shields should exhibit vertical deflections of order 1 km over horizontal baselines of order 100 km. We were thus motivated to search for the presence of flexural troughs in the existing Mars topography data.

Zuber, M. T.

Measuring short-range correlations and quasi-elastic cross sections in A(e,e’) at x > 1 and modest Q 2

We present results from the Jefferson Lab E08-014 experiment, investigating short-range correlations (SRC) through measurements of inclusive quasi-elastic scattering from 2 H, 3 He, 4 He, 12 C, 40 Ca, and 48 Ca. The kinematics were selected to isolate scattering from SRCs, yielding a plateau in the A/ 2 H cross-section ratios due to the universal two-body structure of the 2N-SRCs in light and heavy nuclei. We observe approximate plateaus in the A/ 2 H ratios and provide the first extractions of the A/ 2 H ratio for 40 Ca and 48 Ca. We also examine the A/ 3 He ratio, aiming to identify three-nucleon SRCs (3N-SRCs). Following the approach for isolating 2N-SRCs, searching for 3N-SRC dominance involved measuring the A/ 3 He cross section ratio at modest-to-large Q 2 values and looking for a plateau ratios for x ≳ 2.5. This was not observed in the data, and in fact increasing Q 2 values moved the data further away from the predicted plateau. We show here that, when analyzed in terms of the struck nucleon’s light-cone momentum, the data exhibited the expected trend, progressively approaching the predicted 3N-SRC plateau. These observations suggest that future measurements at higher energies may facilitate a definitive isolation and identification of 3N-SRCs.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS

Algebraic Multigrid with Filtering: An Efficient Preconditioner for Interior Point Methods in Large-Scale Contact Mechanics Optimization

Large-scale contact mechanics simulations are crucial in many engineering fields such as structural design and manufacturing. In the frictionless case, contact can be modeled by minimizing an energy functional; however, these problems are often nonlinear, nonconvex, and increasingly difficult to solve as mesh resolution increases. In this work, we employ a Newton-based interior-point (IP) filter line-search method, an effective approach for large-scale constrained optimization. While this method converges rapidly, each iteration requires solving a large saddle-point linear system that becomes ill-conditioned as the optimization process converges, largely due to IP treatment of the contact constraints. Such ill-conditioning can hinder solver scalability and increase iteration counts with mesh refinement. Here, to address this, we introduce a novel preconditioner, algebraic multigrid with filtering (AMGF), tailored to the Schur complement of the saddle-point system. Building on the classical AMG solver, commonly used for elasticity, we augment it with a specialized subspace correction that filters near null space components introduced by contact interface constraints. Through theoretical analysis and numerical experiments on a range of linear and nonlinear contact problems, we demonstrate that the proposed solver achieves mesh independent convergence and maintains robustness against the ill-conditioning that notoriously plagues IP methods. These results indicate that AMGF makes contact mechanics simulations more tractable and broadens the applicability of Newton-based IP methods in challenging engineering scenarios. More broadly, AMGF is well suited for problems, optimization or otherwise, where solver performance is limited by a low-dimensional subspace, such as those arising from localized constraints, interface conditions, or model heterogeneities. This makes the method widely applicable beyond contact mechanics and constrained optimization.

Mathematics and Computing

An Integrated Data Analytics Platform

An Integrated Science Data Analytics Platform is an environment that enables the confluence of resources for scientific investigation. It harmonizes data, tools and computational resources which subsequently enable the research community to focus on the investigation rather than spending time on security, data preparation, management, etc. OceanWorks is a NASA technology integration project to establish a cloud-based Integrated Ocean Science Data Analytics Platform at NASA’s Physical Oceanography Distributed Active Archive Center (PO.DAAC) for big ocean science. It focuses on advancement and maturity by bringing together several NASA open-source, big data projects for parallel analytics, anomaly detection, in-situ to satellite data matchup, quality-screened data subsetting, search relevancy, and data discovery. Our communities are relying on data distributed through data centers such as the PO.DAAC, COAPS, NCAR, and many others to conduct their research. In typical investigations, scientists would engage in: search for data, evaluate the relevance of that data, download it, and then apply algorithms to identify trends. Such workflow cannot scale if the research involves a massive amount of data or multi-variate measurements. NASA’s Surface Water and Ocean Topography (SWOT) mission is expected to produce massive amount of observational data during its 3-year nominal mission. Collections like SWOT challenges all existing Earth Science data archival, distribution and analysis paradigms. In this paper, we will discuss how OceanWorks enhances the analysis of physical ocean data where the computation is done on an elastic cloud platform next to the archive to deliver fast, web-accessible services for working with oceanographic measurements.

Yang, Chaowei

Unleashing the power of EFT in neutrino-nucleus scattering

Neutrino physics is advancing into a precision era with the construction of new experiments, particularly in the few GeV energy range. Within this energy range, neutrinos exhibit diverse interactions with nucleons and nuclei. This study delves in particular into neutrino-nucleus quasi-elastic cross sections, taking into account both standard and, for the first time, non-standard interactions, all within the framework of effective field theory (EFT). The main uncertainties in these cross sections stem from uncertainties in the nucleon-level form factors, and from the approximations necessary to solve the nuclear many-body problem. We explore how these uncertainties influence the potential of neutrino experiments to probe new physics introduced by left-handed, right-handed, scalar, pseudoscalar, and tensor interactions. For some of these interactions the cross section is enhanced, making long-baseline experiments an excellent place to search for them. Our results, including tabulated cross sections for all interaction types and all neutrino flavors, can serve as the foundation for such searches.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS

Ultrasonic/Sonic Driller/Corer (USDC) as a Subsurface Sampler and Sensors Platform for Planetary Exploration Applications

The search for existing or past life in the Universe is one of the most important objectives of NASA's mission. For this purpose, effective instruments that can sample and conduct in-situ astrobiology analysis are being developed. In support of this objective, a series of novel mechanisms that are driven by an Ultrasonic/Sonic actuator have been developed to probe and sample rocks, ice and soil. This mechanism is driven by an ultrasonic piezoelectric actuator that impacts a bit at sonic frequencies through the use of an intermediate free-mass. Ultrasonic/Sonic Driller/Corer (USDC) devices were made that can produce both core and powdered cuttings, operate as a sounder to emit elastic waves and serve as a platform for sensors. For planetary exploration, this mechanism has the important advantage of requiring low axial force, virtually no torque, and can be duty cycled for operation at low average power. The advantage of requiring low axial load allows overcoming a major limitation of planetary sampling in low gravity environments or when operating from lightweight robots and rovers. The ability to operate at duty cycling with low average power produces a minimum temperature rise allowing for control of the sample integrity and preventing damage to potential biological markers in the acquired sample. The development of the USDC is being pursued on various fronts ranging from analytical modeling to mechanisms improvements while considering a wide range of potential applications. While developing the analytical capability to predict and optimize its performance, efforts are made to enhance its capability to drill at higher power and high speed. Taking advantage of the fact that the bit does not require rotation, sensors (e.g., thermocouple and fiberoptics) were integrated into the bit to examine the borehole during drilling. The sounding effect of the drill was used to emit elastic waves in order to evaluate the surface characteristics of rocks. Since the USDC is driven by piezoelectric actuation mechanism it can designed to operate at extreme temperature environments from very cold as on Titan and Europa to very hot as on Venus. In this paper, a review of the latest development and applications of the USDC will be given.

sampling

Application of machine learning to discover new intermetallic catalysts for the hydrogen evolution and the oxygen reduction reactions

The adsorption energies for hydrogen, oxygen, and hydroxyl were calculated by means of density functional theory on the lowest energy surface of 24 pure metals and 332 binary intermetallic compounds with stoichiometries AB, A 2 B, and A 3 B taking into account the effect of biaxial elastic strains. This information was used to train two random forest regression models, one for the hydrogen adsorption and another for the oxygen and hydroxyl adsorption, based on 9 descriptors that characterized the geometrical and chemical features of the adsorption site as well as the applied strain. All the descriptors for each compound in the models could be obtained from physico-chemical databases. The random forest models were used to predict the adsorption energy for hydrogen, oxygen, and hydroxyl of ≈2700 binary intermetallic compounds with stoichiometries AB, A 2 B, and A 3 B made of metallic elements, excluding those that were environmentally hazardous, radioactive, or toxic. This information was used to search for potential good catalysts for the HER and ORR from the criteria that their adsorption energy for H and O/OH, respectively, should be close to that of Pt. Further, this investigation shows that the suitably trained machine learning models can predict adsorption energies with an accuracy not far away from density functional theory calculations with minimum computational cost from descriptors that are readily available in physico-chemical databases for any compound. Moreover, the strategy presented in this paper can be easily extended to other compounds and catalytic reactions, and is expected to foster the use of ML methods in catalysis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH

Novel cross section ratios as possible signals of saturation in ultraperipheral collisions

We propose new cross section ratios in ultraperipheral A + A and p + A collisions (UPCs) as a potential experimental signal of saturation physics. We consider the ratio R 1 of elastic vector meson photoproduction cross section to the inclusive hadron or jet photoproduction cross section. The ratio can be measured in the γ + A and γ + p collisions taking place in the UPCs. We label the ratios R 1 ( A ) and R 1 ( p ) , respectively. Constructing the double ratio R UPC ( A ) = R 1 ( A ) / R 1 ( p ) , and performing a small- x calculation both in the quasiclassical approximation and by including small- x evolution, we observe that R UPC ( A ) exhibits a markedly different dependence on the nuclear atomic number inside and outside the saturation region. This result indicates that R UPC ( A ) and R 1 ( A ) measurements in UPCs may help in the experimental searches for the evidence of saturation physics. Published by the American Physical Society 2024

73 NUCLEAR PHYSICS AND RADIATION PHYSICS

Improved Joint FHC+RHC Sterile Analysis with Nu-on-E Sample

The NOvA experiment measures neutrino oscillations in a nearly pure muon neutrino beam. The neutral current (NC) and the muon neutrino charged current (CC) interaction disappearance are studied to search for active-to-sterile neutrino oscillations. NOvA consists of two functionally identical detectors placed 14 mrad off-axis of Fermilab's NuMI beamline. The Near Detector (ND) is located 100 m underground on-site at Fermilab, while the Far Detector (FD) sits on the surface 810 km away at Ash River, Minnesota. The NuMI beam is capable of running in forward or reverse horn current modes, creating a predominantly neutrino or antineutrino beam, respectively. For this analysis, a simultaneous near and far detector joint fit approach was utilized. Systematic uncertainties constrain the ND while the statistical uncertainties constrain the FD. For the upcoming analysis, we plan on including a dedicated neutrino-on-electron elastic scattering sample. This interaction is advantageous as it is subject to minimal cross-section uncertainties, which will help reduce degeneracies between other uncertainty categories. This thesis presents the sterile neutrino oscillation sensitivity results using 27.2e+20 POT in FHC mode, 12.6e+20 in RHC mode.

Haejun, Stella Oh [Cincinnati U.]

NOvA Sterile Neutrino Search with Additional Systematic Constraints

The NOvA experiment measures neutrino oscillations in a nearly pure muon neutrino beam. The neutral current (NC) and the muon neutrino charged current (CC) interaction disappearance are studied to search for active-to-sterile neutrino oscillations. NOvA consists of two functionally identical detectors placed 14 mrad off-axis of the Fermi National Accelerator Laboratory's NuMI beamline. The Near Detector (ND) is located 100 meters underground on-site at Fermilab, while the Far Detector (FD) is situated on the surface 810 kilometers away at Ash River, Minnesota. The NuMI beam is capable of operating in either forward or reverse horn current modes, producing predominantly neutrino or antineutrino beams, respectively. For this analysis, a simultaneous near and far detector joint fit approach was utilized. Systematic uncertainties constrain the ND while the statistical uncertainties constrain the FD. In addition to the dual baseline approach, the neutrino-on-electron elastic scattering sample was included. This interaction is advantageous as it is subject to minimal cross-section uncertainties, which will help reduce degeneracies between other uncertainty categories. This thesis presents the sterile neutrino oscillation sensitivity results using 27.2e+20 POT in FHC mode, 12.6e+20 in RHC mode, with the nu-on-e sample of 13.6e+20 POT in FHC.

Oh, Haejun [Cincinnati U.]

A comprehensive review on valorization of chestnut processing wastes into bio‐based composites and bioplastics

Abstract This review examines the characterization and utilization of chestnut processing wastes (35%) in the production of bioplastics and biocomposites. In this review, a Web of Science search without any publishing year restriction on the biochemical compositions of all the components of Castanea sativa . The obtaining of bioplastics and biocomposites based on C. sativa was reviewed. First, it highlights the biochemical composition and antioxidant properties of chestnut fruit, shell, burrs, leaves, flowers, and wood focusing on the most important compounds, such as phenolic acids, flavonoids, carbohydrates, Klason lignin, cellulose, and glucan, which can enhance the properties of these materials. Then the review covers using several chestnut extracts and fillers in bioplastics production through solvent casting technique. The mechanical, structural, bioactive properties, and moisture content were optimized through the composition and production. The color, UV absorption, antioxidant, and antimicrobial activity were also discussed. Biocomposites reinforced with chestnut burs, shells, or wood flour increased the intended properties. The enhancements in tensile strength, elastic modulus, and the effects of a pre‐treatment were evaluated. Additionally, it discusses material recovery, recycling, and reuse, particularly how it affects the biodegradability of composites incorporating chestnut waste residues. Highlights Chestnut fruit, shells, and burrs are rich in starch, lignin, and cellulose. The waste of chestnut processing can be used in bioplastics and biocomposites. Chestnut‐based films and biocomposites exhibit promising mechanical properties. The antimicrobial activity, making films, and composites proper for food packaging. New techniques boost performance, offering alternatives to conventional plastics.

Silva, Simão B. [REQUIMTE/LAQV, ISEP, Polytechnic

Nuclear-level effective theory of 𝜇 → 𝑒 conversion: Inelastic process

Mu2e and COMET will search for electrons produced via the neutrinoless conversion of stopped muons bound in 1s atomic orbits of 27 Al, improving existing limits on charged lepton flavor violation (CLFV) by roughly four orders of magnitude. Conventionally, 𝜇 → 𝑒 conversion experiments are optimized to detect electrons originating from transitions where the nucleus remains in the ground state, thereby maximizing the energy of the outgoing electron. Clearly, detection of a positive signal in forthcoming experiments would stimulate additional work—including subsequent conversion experiments using complementary nuclear targets—to further constrain the new physics responsible for CLFV. Here we argue that additional information can be extracted without the need for additional experiments, by considering inelastic conversion in 27 Al. Transitions to low-lying nuclear excited states can modify the near-endpoint spectrum of conversion electrons, with the ratio of the elastic and inelastic responses being sensitive to the underlying CLFV operator. We extend the nuclear effective theory of 𝜇 → 𝑒 conversion to the inelastic case, which adds five new response functions to the six that arise for the elastic process. We evaluate these nuclear response functions in 27 Al and calculate the resulting conversion-electron signal, taking into account the resolution anticipated in Mu2e/COMET. We find that 27 Al is an excellent target choice from the perspective of the new information that can be obtained from inelastic 𝜇 → 𝑒 conversion.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS

Direct Deoxygenation of Phenol over Fe-Based Bimetallic Surfaces Using On-the-Fly Surrogate Models

We present an accelerated nudged elastic band (NEB) study of phenol direct deoxygenation (DDO) on Fe-based bimetallic surfaces using a recently developed Gaussian process regression (GPR) calculator. Our test calculations demonstrate that the GPR calculator achieves up to 3 times speedup compared to conventional density functional theory calculations while maintaining high accuracy, with energy barrier errors below 0.015 eV. Using GPR-NEB, we systematically examine the DDO mechanism on pure Fe(110) and surfaces modified with Co and Ni in both top and subsurface layers. Our results show that subsurface Co and Ni substitutions preserve favorable thermodynamics and kinetics for both C–O bond cleavage and C–H bond formation, comparable to those on the pure Fe(110) surface. In contrast, top-layer substitutions generally increase the C–O bond cleavage barrier, render the step endothermic, and result in significantly higher reverse reaction rates, making DDO unfavorable on these surfaces. This work demonstrates the effectiveness of GRR-accelerated transition state searches for complex surface reactions and provides insights into rational design of bimetallic catalysts for selective deoxygenation.

Aromatic compounds