SEARCH · Search NASA

Results for “Forcing”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 397 records · Page 22

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Recalibration of missing low-frequency variability and trends in the North Atlantic Oscillation

Abstract Multi-decadal trends in the wintertime North Atlantic Oscillation (NAO) are under-represented by coupled general circulation models (CGCMs), consistent with a lack of autocorrelation in their NAO index series. This study proposes and tests two simple “reddening” approaches for correcting this problem in simulated indices based on simple one parameter short-term (AR; Auto-Regressive order 1) and long-term (FD; Fractional-Difference) time series filters. Using CGCMs from the Coupled Model Intercomparison Project Phase 6 (CMIP6), the FD filter successfully improves the autocorrelation structure of the NAO, and in turn the simulation of extreme trends, while the AR filter is less successful. The 1963–1993 NAO trend is the maximum 31-year trend in the historical period. Raw CGCMs underestimate the likelihood of this trend by a factor of ten but this discrepancy is corrected after reddening. CMIP6 future projections show that long-term (2024–2094) NAO ensemble mean trends systematically increase with the magnitude of radiative forcing: -2.4 to 3.5 hPa/century for low-to-high forcing after reddening (more than double the range using raw output). The related likelihood of future maximum 31year trends comparable to 1963–1993 ranges from 3 to 7% whereas none of these CMIP6 projections simulate this without reddening. Near-term projections of the next 31 years (2024–2054) are less sensitive than long term trends to the future scenario, showing weak-to-no forced trend. However, reddening increases the ensemble range by 74% (to +/-1 standard deviation/decade), which could increase/decrease regional climate change signals in the Northern Hemisphere by magnitudes that are underestimated when using raw CGCM output.

Meteorology & Atmospheric Sciences↗

DOE OSTI2024

Embedded symmetric positive semi-definite machine-learned elements for reduced-order modeling in finite-element simulations with application to threaded fasteners

Here, we present a machine-learning strategy for finite element analysis of solid mechanics wherein we replace complex portions of a computational domain with a data-driven surrogate. In the proposed strategy, we decompose a computational domain into an “outer” coarse-scale domain that we resolve using a finite element method (FEM) and an “inner” fine-scale domain. We then develop a machine-learned (ML) model for the impact of the inner domain on the outer domain. In essence, for solid mechanics, our machine-learned surrogate performs static condensation of the inner domain degrees of freedom. This is achieved by learning the map from displacements on the inner-outer domain interface boundary to forces contributed by the inner domain to the outer domain on the same interface boundary. We consider two such mappings, one that directly maps from displacements to forces without constraints, and one that maps from displacements to forces by virtue of learning a symmetric positive semi-definite (SPSD) stiffness matrix. We demonstrate, in a simplified setting, that learning an SPSD stiffness matrix results in a coarse-scale problem that is well-posed with a unique solution. We present numerical experiments on several exemplars, ranging from finite deformations of a cube to finite deformations with contact of a fastener-bushing geometry. We demonstrate that enforcing an SPSD stiffness matrix drastically improves the robustness and accuracy of FEM–ML coupled simulations, and that the resulting methods can accurately characterize out-of-sample loading configurations with significant speedups over the standard FEM simulations.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2026

Impedance Based Fixture Neutralization (FINE) Method to Replicate Field Environments in Laboratory Settings

Test article vibration test responses will vary from one test laboratory to another due to differences in shaker and fixture dynamic characteristics. This is also seen between the field and laboratory due to the differing dynamic characteristics between field assembly loading and boundary conditions compared to the laboratory configuration. This work introduces a technique called Impedance Fixture Neutralization which customizes input forces to cause consistent responses for the Device Under Test across different vibration testing conditions. The customized force neutralizes the dynamic variations between configurations for the Device Under Test. The device responses can be replicated in several situations: when using a different fixture with the same attachment points, when the test article mounting location on the fixture changes, when the force location changes, or any combination of these situations. Impedance Fixture Neutralization uses the uncoupled dynamic characteristics of the test article and excitation fixtures to customize the input thereby causing the same test article responses between two mounting configurations (either field to laboratory or two different laboratory configurations). The application of the technique is shown using an analytical model of a two-beam system and analytically using experimental FRFs from plate and frame component hardware characterization tests. In both cases a device under test is attached to two test fixtures wherein the dynamic differences in the configurations are neutralized.

Laboratory↗

DOE OSTI2025

Operando synchrotron X-ray analysis of melt pool dynamics in an Al-Sn immiscible alloy

The melt flow in an Al-50vol.% Sn immiscible alloy, produced by single-track laser melting of Al and Sn elemental powders, was studied in real time. High-speed synchrotron X-ray imaging was used to track Al and Sn liquids’ movements, and also to examine elemental distributions in the laser tracks, complimented by electron microscopy after solidification. Key aspects, including melt pool geometry, keyhole instability, and flow dynamics (flow pattern and velocity), were examined using digital image analysis. Relatively deeper melt pools formed at 400 W and 300 mm/s exhibited greater stability, with smooth surfaces, consistent outward flow, and minor vortices near the keyhole. In contrast, shallower pools produced at higher scanning speeds (>500 mm/s) demonstrated greater instability with increased surface waviness, and stronger velocity fluctuations, leading to numerous micro-vortices and increased Al-Sn heterogeneity. Velocity scale estimations, supported by experimental observations, examined the roles of vapour pressure, Marangoni effect, buoyancy, inertial, and surface tension forces in the flow. The results revealed that vapour pressure and mechanical waves dominated at high scanning speeds (shallow pools), while Marangoni forces were equally significant in deep pools at lower speeds (300 mm/s). Buoyancy was found to have minimal impact in both cases. Furthermore, the interaction between inertial and surface tension forces played a critical role in determining the degree of waviness of the pools’ surfaces. These findings offer valuable insights into melt pool dynamics during laser processing of immiscible alloys and other metallic systems using elemental powders, and provide guidance for developing high-fidelity computational fluid dynamics models.

immiscible alloys↗

DOE OSTI2026

Maximizing machine learning interatomic potential transferability for the discovery of the novel stellated octadecagon Bi18-Pt24 cage structure

Achieving true transferability remains the central challenge for Machine Learning Interatomic Potentials (ML-IAPs) in modeling complex bimetallic nanoclusters across their vast potential energy surfaces. We systematically investigate data selection strategies to optimize the Chebyshev Interaction Model for Efficient Simulation (ChIMES) potential for the Bi-Pt nanoclusters by comparing three innovative sampling methods: Principal Component Analysis (PCA)/k-means (structural diversity), t-distributedStochasticNeighborEmbedding (t-SNE)/k-means (force-space diversity), and hierarchical clustering. Quantitatively, the PCA/k-means strategy proved most effective for global accuracy, yielding the lowest force errors and achieving energy root mean square errors (RMSE) values competitive with Density Functional Theory (DFT), demonstrating excellent accuracy (19.16meV/atom). Structural validation on 34 unique DFT-optimized isomers further confirmed the potential’s high fidelity, with the best model PCA/k-means reproducing structures with an average root mean square deviation (RMSD) of 0.10 Å. However, the t-SNE methods, by maximizing diversity in the force space, demonstrated superior extrapolative power, leading to the more precise prediction of a novel stellated octadecagon Bi18⁢Pt24 cage structure, demonstrating the potential for exploring previously unseen morphologies. Our results establish a clear methodology for strategic data sampling that successfully maximizes ML-IAP transferability, providing an accurate and computationally efficient tool that accelerates the theoretical discovery of complex bimetallic architectures.

Vangheluwe, Raphaël [Université Paris-Saclay, CNRS↗

DOE OSTI2024

Analytically differentiable metrics for phase stability

Here, in this work, a long-established but sparsely documented method of obtaining semi-analytic derivatives of thermodynamic properties with respect to equilibrium conditions is briefly reviewed and rigorously derived. This procedure is then leveraged to construct general forms of derivatives of the residual driving force, a metric for measuring phase stability used in CALPHAD model optimization, with respect to overall system and individual phase compositions. Applied examples – calculating heat capacity in the Al-Fe system, thermodynamic factors in the Nb-V-W system, and residual driving force derivatives in the Ni-Ti system – demonstrate the versatility, accuracy, and extensibility of this method. Using the developed method, residual driving force gradients can be applied directly in CALPHAD model optimizers, as well as in materials design frameworks, to identify regions of phase stability with an efficient, gradient-based approach.

36 MATERIALS SCIENCE↗

DOE OSTI2025

A cohesive zone treatment for the material point method involving problems of large deformation and damage

A new algorithm is described that permits the use of cohesive zones in the material point method for problems involving large deformation and fracture. In contrast to previous cohesive zone implementations, this method does not utilize massless surface-element particles. Instead, cohesive tractions are computed using the shape function mappings from a reference grid configuration in combination with explicitly defined particle surface normals and surface positions. These normals and relative surface positions are updated each time step according to particle deformation. The tractions are converted to cohesive forces using the nodal areas and mapped back to particles using the same reference shape function mappings. These forces are then remapped by conventional particle-to-grid interpolation as external forces using the current-configuration shape-function mappings. This allows highly compliant cohesive zones to function over jump displacements larger than a grid cell. Upon damage, these interfaces can revert to conventional multi-field contact surfaces. This approach is general and readily applies to two and three dimensions as well as being compatible with damage-field gradient partitioning offering exceptional computational flexibility. The framework for this method enables other capabilities, such as improved contact precision using explicitly defined surface normals and positions, and a method to mitigate spurious material damage at weak discontinuities between stiff brittle materials and soft or compliant materials.

Cohesive zone↗

DOE OSTI2025

Comparison of interlaminar damage modeling strategies for hybrid composite/aluminum laminates subjected to low-velocity impact

Low-velocity impact of hybrid metal-composite structures was investigated experimentally and computationally. Composite laminates consisting of 2D woven glass fiber reinforced polymer (GFRP) and carbon fiber reinforced polymer (CFRP) were joined with a 6061-T6 aluminum plate using an epoxy adhesive. Two variations of the structure were studied; one consisting of all plies oriented at 0° and one consisting of all plies oriented at 45°. A drop tower was used to impact structures at a range of energies, including energies above and below the threshold at which the aluminum layer was perforated. Numerical simulations were implemented using Sierra/SM, an in-house transient dynamics finite element code developed at Sandia National Laboratories. A Hosford plasticity model was used to describe the response of the aluminum layer. A newly implemented orthotropic continuum damage mechanics (CDM) constitutive model was used to represent the composite laminate. This 3D-CDM model was compared to a cohesive zone model (2D-CDM/CZM) to investigate efficacy of aluminum perforation energy prediction, delamination prediction, and computational cost. Accuracy of each model was evaluated using the experimental results. Each showed good agreement with the tests for both the force and velocity histories, as well as the observed damage mechanisms. The 2D-CDM/CZM model was marginally more accurate in capturing both the composite and aluminum behavior — this model averaged error percentages of -11.2% and 10.8% for residual velocity and peak force, respectively. Meanwhile, the 3D-CDM model predictions yielded average error percentages of -35.5% (velocity) and 22.6% (force). However, the 3D-CDM model generally resulted in a decreased computational cost; the average run time was 14% shorter than the 2D-CDM/CZM model and 3x as many timesteps per hour were computed using the same computational resources. In conclusion, new experimental data on the impact and perforation resistance of metal-composite laminates is presented in addition to numerical predictions of the impact behavior.

Carbon fiber↗

DOE OSTI2026

Reduced-order CFD modeling of cryogenic hydrogen isotope extrusion for pellet fueling

This study presents a reduced-order model (ROM) for computational fluid dynamics (CFD) simulations of cryogenic hydrogen isotope extrusions, focusing on protium (H₂) and deuterium (D₂) piston extruders. Using a 2D axisymmetric ROM in ANSYS-Polyflow, significant computational savings were achieved (runtime reduced from 9∼24 h to 3∼5 min), with extrusion force discrepancies between the 2D ROM and 3D models being on the order of 1%. Parametric studies identified optimal cutoff shear rates in the viscosity model (0.01/s for H₂ and 0.001/s for D₂), providing recommendations for future simulations. Finally, a comprehensive comparison of ROM results with experimental data was performed across varying geometries, cryogenic materials, temperatures, extrusion lengths, and piston velocities. Predictions at low extrusion temperatures met the objective of providing quick and efficient solutions with an acceptable extrusion force error of approximately 10% or less, validating the effectiveness of the 2D ROM approach. However, at high temperatures closer to the triple point, extrusion force error grows, which necessitates developing an improved model that accounts for temperature effects, e.g. melting. Nevertheless, the findings still represent a significant improvement in efficiency of CFD modeling of cryogenic hydrogenic extrusion. The ROM framework can also be extended to tritium (T2) and screw extruders, which will ultimately provide a fast and effective tool for optimizing pellet injector design for ITER and future reactor systems.

Fan, Joy [ORNL] (ORCID:0000000229751735)↗

DOE OSTI2024

Multiparametric AFM Insights into Electron Transport Mechanisms in Biomemristors

Biomemristors have attracted significant attention due to their applications in biodegradable and flexible electronic devices, and their prospective functions in information storage, neural synapses, and neuromorphic computation. Hence, understanding biomaterials' resistive switching (RS) mechanism is essential to design novel biomemristors with enhanced performance. In the last decade, various atomic force microscopy (AFM) modes, such as Kelvin probe force microscopy (KPFM) and conductive atomic force microscopy (cAFM), in situ/operando visualize the electron transport behaviors within biomemristors at the nanoscale. This review provides a comprehensive discussion of multiparametric AFM techniques to characterize biomemristors and inspire new ideas in developing next-generation consumer electronics. We also summarize the progress of revealing the internal mechanism of biomemristors with functional AFM to visualize and understand the switching behavior. Finally, the opportunities and challenges of further utilizing multiparametric AFM to explore the electron transport mechanisms in biomemristors are presented.

Shi, Chenyang↗

DOE OSTI2025

An airfoil-based synthetic actuator disk model for wind turbine aerodynamic and structural analysis

Here, this study introduces an airfoil-based refinement technique to enhance the Actuator Disk Model (ADM) for improved wind turbine aerodynamic load prediction and structural simulation in conjunction with Large Eddy Simulations of the wind flow. While ADM offers higher computational efficiency than the more detailed but resource-intensive Actuator Line Model (ALM), it traditionally lacks the resolution needed to capture the localized blade forces accurately. To address this limitation, we introduce a refinement technique that uses airfoil-specific data and employs interpolation-based grid point refinement, achieving ALM-comparable accuracy while preserving ADM's efficiency. Unlike conventional ADM that provides only rotor-disk averaged forces, our synthetic method tracks transient aerodynamic load variations over multiple blade revolutions, allowing us to calculate the distributions of maximum and minimum loads during typical cycles. Applied to the NREL 5 MW reference turbine, our enhanced ADM accurately predicts key aerodynamic parameters (angle of attack, axial velocity, lift, drag, axial and tangential forces along the blades) as well as structural responses (blade tip deflection, maximum stress, and stress concentration). Our results show that the tip deflection ranges from 2.33m (3.69 % of blade length) to 4.28m (6.79 %), with maximum stress concentration occurring near the blade root. This research demonstrates that a refined synthetic ADM approach can serve as a computationally efficient alternative for both aerodynamic analysis and structural simulation of wind turbine blades subjected to realistic wind fields.

17 WIND ENERGY↗

DOE OSTI2026

Numerical Analysis of Offshore Wind-Farm-Induced Drag Effects on Coastal Upwelling Dynamics

Wind farms extract momentum from the atmospheric flow, generating wind-speed deficits both within the plant, and extending downstream. When located offshore, these deficits modulate air-sea coupling, potentially impacting coastal upwelling in sensitive regions. We investigate impacts of wind farms on coastal upwelling using kilometre-scale, three-way-coupled simulations with the coupled ocean-atmosphere-wave-sediment transport system for the US West Coast. Wind-farm effects are represented by a generalised turbine drag formulation, an idealised, height-dependent body force whose magnitude is systematically varied. This approach isolates the leading-order fluid-dynamical response in a realistic coastal configuration. The atmospheric adjustment exhibits an approximately linear relation between drag force and wind-speed deficit, with wakes that expand downstream and increase in magnitude as drag increases. An empirical orthogonal function analysis of sea-surface-temperature anomalies reveals the emergence of a canonical dipole pattern under strong drag forcing. Subsurface diagnostics show consistent shoaling of the mixed layer and suppressed upward velocities in areas near wind-farm region, accompanied by compensating enhancements of shoaling closer to the coast. These results identify turbine drag as a control parameter in assessing interactions between wind-farm wake and coastal upwelling and provide scaling relationships for understanding offshore wind-farm effects on the coastal circulation dynamics.

16 TIDAL AND WAVE POWER↗

DOE OSTI2024

Magneto-Stokes flow in a shallow free-surface annulus

In this study, we analyse ‘magneto-Stokes’ flow, a fundamental magnetohydrodynamic (MHD) flow that shares the cylindrical-annular geometry of the Taylor–Couette cell but uses applied electromagnetic forces to circulate a free-surface layer of electrolyte at low Reynolds numbers. The first complete, analytical solution for time-dependent magneto-Stokes flow is presented and validated with coupled laboratory and numerical experiments. Three regimes are distinguished (shallow-layer, transitional and deep-layer flow regimes), and their influence on the efficiency of microscale mixing is clarified. The solution in the shallow-layer limit belongs to a newly identified class of MHD potential flows, and thus induces mixing without the aid of axial vorticity. We show that these shallow-layer magneto-Stokes flows can still augment mixing in distinct Taylor dispersion and advection-dominated mixing regimes. The existence of enhanced mixing across all three distinguished flow regimes is predicted by asymptotic scaling laws and supported by three-dimensional numerical simulations. Mixing enhancement is initiated with the least electromagnetic forcing in channels with order-unity depth-to-gap-width ratios. If the strength of the electromagnetic forcing is not a constraint, then shallow-layer flows can still yield the shortest mixing times in the advection-dominated limit. Our robust description of momentum evolution and mixing of passive tracers makes the annular magneto-Stokes system fit for use as an MHD reference flow.

58 GEOSCIENCES↗

DOE OSTI2025

Sedimentation and shear-induced dynamics of spheroids in fluids with spatial viscosity variations

A generalized reciprocal theorem is used to relate the force and torque induced on a particle in an inertia-less fluid with small variation in viscosity to integrals involving Stokes flow fields and the spatial dependence of viscosity. These resistivity expressions are analytically evaluated using spheroidal harmonics and then used to obtain the mobility of the spheroid during sedimentation, and in linear flows, of a fluid with linear viscosity stratification. The coupling between the rotational and translational motion induced by stratification rotates the spheroid’s centerline, creating a variety of rotational and translational dynamics dependent upon the particle’s aspect ratio, κ, and the component of the stratification unit vector in the gravity direction, d g . Spheroids with 0.55 ⪅ κ ⪅ 2.0 exhibit the largest variety of settling behaviors. Interestingly, this range covers most microplastics and typical microorganisms. One of the modes include a stable orientation dependent only on κ and d g , but independent of initial orientation, thus allowing for the potential control of settling angles and sedimentation rates. In a simple shear flow, cross-streamline migration occurs due to the stratification-induced force generated on the particle. Similarly, a particle no longer stays at the stagnation point of a uniaxial extensional flow. While fully analytical results are obtained for spheroids, numerical simulations provide a source of validation. These simulations also provide additional insights into the stratification-induced force- and torque-producing mechanisms through the stratification-induced stress, which is not accessed in the reciprocal theorem-based analytical calculations.

Geophysical and Geological Flows: Stratified flows↗

DOE OSTI2026

Shear-free, inhomogeneous turbulence in a stably stratified fluid

High-resolution large eddy simulations are conducted of locally forced, shear-free turbulence in the presence of an initially sharp density interface. The simulations are reminiscent of oscillating grid turbulence experiments used to isolate the effect of turbulent diffusion and entrainment from background shear. By simulating such a flow we avoid common challenges of the experiments such as secondary-flow contamination due to sidewall effects and the inevitable interaction of the stratifying agent and forcing region. To address the latter concern, we add a heating term (potential energy sink) to the governing equations in the forcing layer, thereby preventing a heat flux through the source region. This modification sets up a continuous stratification in the mixed layer that is often assumed to be negligible in experiments. Despite this difference, we are able to make meaningful comparisons in terms of the overall entrainment rate, which varies as a power law with a turbulent Richardson number. Two exponents, −2 and −1, are measured depending on the definition of the Richardson number and entrainment rate used. The definition leading to −1 is consistent with most experiments, and we argue it is the superior choice if one is able to measure the relevant quantities. We also verify the self-similar scaling of turbulence velocity and length scales in the homogeneous fluid and propose ‘inner’ and ‘outer’ scalings for the stratified cases based on a local Froude number. The detailed scaling results are useful for turbulence model validation.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2026

Solvent Stabilization of Protic Oxonium/Ammonium Intermediates in Cation Radical Cyclization Reactions Investigated via Computational Approaches

Anodic intramolecular cyclization reactions have substantial synthetic utility for formation of cyclic carbon–carbon or carbon–heteroatom bonds. For cases of intramolecular trapping of a cation radical by a protic nucleophile, the cyclization step coincides with a substantial increase in substrate acidity and thus may exhibit particularly pronounced solvent effects. In this computational work, we employ both quantum chemical (QM) and quantum mechanics/molecular mechanics (QM/MM) methods to compute solvent effects on free energy profiles for cyclization and deprotonation reaction steps for cation radical intermediates of substrates representative for anodic intramolecular cyclizations. We find substantial solvent contribution to the thermodynamic driving force for cation radical cyclization; for example, methanol and tetrahydrofuran solvents provide ∼30–35 kJ/mol driving force to form cyclic oxonium cation radicals and ∼15–25 kJ/mol driving force to form cyclic ammonium cation radicals, compared to baseline reactions in dichloromethane solvent. Given that these solvent shifts are on par with the innate cyclization reaction thermodynamics, the choice of solvent plays a crucial role in promoting/driving the cation radical cyclization step. Methanol is particularly effective at facilitating rapid deprotonation of the cyclic cation radical intermediate, which may lead to the full electrochemical process (e.g., second electron transfer) proceeding heterogeneously at the anode.

Khan, Shahriar [ORNL] (ORCID:0000000289138430)↗

DOE OSTI2025

How Thermodynamic, Electronic, and Steric Factors Influence Mesitylcopper Oligomers

Mesitylcopper (CuMes) is a highly versatile organocopper reagent used in both organic and inorganic syntheses. It has previously been shown that CuMes exists as a tetrameric or pentameric cyclic oligomer [CuMes] n (n = 4, 5), both in solution and in the solid state. The bonding arrangement between the [CuMes] units has qualitatively been described as localized three-center twoelectron (3c-2e) bonds. However, the electronic, structural, and thermodynamic forces driving this aggregation are still not well understood. For this reason, we employed density functional theory (DFT) calculations to study mesitylcopper as a monomeric [CuMes] unit and [CuMes]n oligomers with n = 2 to n = 7. We found that there is a strong electronic driving force for aggregation caused by strong mixing between the Cu’s d orbitals and Mes’s π orbitals in oligomers larger than the dimer. This mixing is only optimized in oligomers with n ≥ 3, where the mesityl group is no longer bonded to a single copper center but instead becomes a bridging ligand. Beyond the trimer, steric and entropic factors become relevant for determining the relative stabilities of the different aggregates, with midsized oligomers (n = 4−5) having the optimal balance between the electronic Cu−C bonding character, Cu···Cu attractive forces, entropy, reduced internal ring strain, and reduced steric interactions between the mesityl groups.

Copper↗

DOE OSTI2024

A Molecular View of Methane Activation on Ni(111) through Enhanced Sampling and Machine Learning

A combination of machine learned interatomic potentials (MLIPs) and enhanced sampling simulations is used to investigate the activation of methane on a Ni(111) surface. The work entails the development and iterative refinement of MLIPs, initially trained on a dataset constructed via ab initio molecular dynamics (AIMD) simulations, supplemented by adaptive biasing forces, to enrich the sampling of catalytically relevant configurations. Our results reveal that by incorporating collective variables that capture the behavior of the reactant molecule, as well as additional frames that describe the dynamic response of the catalytic surface, it is possible to enhance considerably the accuracy of predicted energies and forces. By employing enhanced sampling schemes in the refinement of the MLIP, we systematically explore the potential energy surface, leading to a refined MLIP capable of predicting DFT-level energies and forces and replicating key geometric characteristics of the catalytic system. The resulting free energy landscapes at several temperatures provide a detailed view of the thermodynamics and dynamics of methane activation. Specifically, as methane approaches and dissociates on the catalytic surface, the process involves the dynamic interplay of CH 4 and the Ni catalyst that includes both enthalpic and entropic contributions. The progression towards the transition state involves an CH 4 moiety that is increasingly restrained in its ability to rotate or translate, while the stage following the transition state is characterized by a notable rise of the Ni atom that interacts with the cleaved C–H bond. Furthermore, this leads to an increase in the mobility of the adsorbed species, a feature that becomes more pronounced at higher temperatures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗