Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dependability modeling”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 1,243 records · Page 69

Determinants of cognitive and brain resilience to tau pathology: a longitudinal analysis

Mechanisms of resilience against tau pathology in individuals across the Alzheimer’s disease spectrum are insufficiently understood. Longitudinal data are necessary to reveal which factors relate to preserved cognition (i.e. cognitive resilience) and brain structure (i.e. brain resilience) despite abundant tau pathology, and to clarify whether these associations are cross-sectional or longitudinal. We used a longitudinal study design to investigate the role of several demographic, biological and brain structural factors in yielding cognitive and brain resilience to tau pathology as measured with PET. In this multicentre study, we included 366 amyloid-β-positive individuals with mild cognitive impairment or Alzheimer’s disease dementia with baseline 18 F-flortaucipir-PET and longitudinal cognitive assessments. A subset (n = 200) additionally underwent longitudinal structural MRI. We used linear mixed-effects models with global cognition and cortical thickness as dependent variables to investigate determinants of cognitive resilience and brain resilience, respectively. Models assessed whether age, sex, years of education, APOE-ε4 status, intracranial volume (and cortical thickness for cognitive resilience models) modified the association of tau pathology with cognitive decline or cortical thinning. We found that the association between higher baseline tau-PET levels (quantified in a temporal meta-region of interest) and rate of cognitive decline (measured with repeated Mini-Mental State Examination) was adversely modified by older age (Stβ interaction = -0.062, P = 0.032), higher education level (Stβ interaction = -0.072, P = 0.011) and higher intracranial volume (Stβ interaction = -0.07, P = 0.016). Younger age, higher education and greater cortical thickness were associated with better cognitive performance at baseline. Greater cortical thickness was furthermore associated with slower cognitive decline independent of tau burden. Higher education also modified the negative impact of tau-PET on cortical thinning, while older age was associated with higher baseline cortical thickness and slower rate of cortical thinning independent of tau. Our analyses revealed no (cross-sectional or longitudinal) associations for sex and APOE-ε4 status on cognition and cortical thickness. In this longitudinal study of clinically impaired individuals with underlying Alzheimer’s disease neuropathological changes, we identified education as the most robust determinant of both cognitive and brain resilience against tau pathology. The observed interaction with tau burden on cognitive decline suggests that education may be protective against cognitive decline and brain atrophy at lower levels of tau pathology, with a potential depletion of resilience resources with advancing pathology. Finally, we did not find major contributions of sex to brain nor cognitive resilience, suggesting that previous links between sex and resilience might be mainly driven by cross-sectional differences.

59 BASIC BIOLOGICAL SCIENCES↗

Time-dependent orbital-free density functional theory: Background and Pauli kernel approximations

Time-dependent orbital-free density functional theory (DFT) is an efficient method for calculating the dynamic properties of large-scale quantum systems due to the low computational cost compared to standard time-dependent DFT. In this work, we formalize this method by mapping the real system of interacting fermions onto a fictitious system of noninteracting bosons. The dynamic Pauli potential and associated kernel emerge as key ingredients of time-dependent orbital-free DFT. Using the uniform electron gas as a model system, we derive an approximate frequency-dependent Pauli kernel. Pilot calculations suggest that space nonlocality is a key feature for this kernel. Nonlocal terms arise already in the second-order expansion with respect to unitless frequency and reciprocal space variable ($\frac{ω}{qk_F}$ and $\frac{q}{2k_F}$, respectively). Given the encouraging performance of the proposed kernel, we expect it will lead to more accurate orbital-free DFT simulations of nanoscale systems out of equilibrium. Additionally, the proposed path to formulate nonadiabatic Pauli kernels presents several avenues for further improvements which can be exploited in future work to improve the results.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

TRINIDI: Time-of-Flight Resonance Imaging With Neutrons for Isotopic Density Inference

Accurate reconstruction of 2D and 3D isotope densities is a desired capability with great potential impact in applications such as evaluation and development of next-generation nuclear fuels. Neutron time-of-flight (TOF) resonance imaging offers a potential approach by exploiting the characteristic neutron absorption spectra of each isotope. However, it is a major challenge to compute quantitatively accurate images due to a variety of confounding effects such as severe Poisson noise, background scatter, beam non-uniformity, absorption non-linearity, and extended source pulse duration. We present the TRINIDI algorithm which is based on a two-step process in which we first estimate the neutron flux and background counts, and then reconstruct the areal densities of each isotope and pixel. Both components are based on the inversion of a forward model that accounts for the highly non-linear absorption, energy-dependent emission profile, and Poisson noise, while also modeling the substantial spatio-temporal variation of the background and flux. Further, to do this, we formulate the non-linear inverse problem as two optimization problems that are solved in sequence. We demonstrate on both synthetic and measured data that TRINIDI can reconstruct quantitatively accurate 2D views of isotopic areal density that can then be reconstructed into quantitatively accurate 3D volumes of isotopic volumetric density.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Effect of Thermal Expansion Coefficient, Viscosity and Melting Range in Simulation of PCM Embedded Heat Exchangers With and Without Fins

Phase change material heat exchangers (PCMHX) have an important role in integrating renewable energy systems. PCMHX can offer high storage density at various temperatures of interest and can be used for grid load shifting purposes. Numerical models enable engineers to estimate PCMHX performance for different design parameters and operating conditions. Modeling phase change phenomena is challenging due to the complex time-dependent nature of the process. The accuracy of models is highly sensitive to PCM thermo-physical properties. Thermal expansion coefficient (β), viscosity (μ) and melting temperature range (MR) of a PCM are important properties, especially when natural convection is not negligible. In PCMHX modeling, using less than accurate values for these properties can have significant impact on the simulation outcomes. These properties and discussions thereof are not readily available in the literature. This paper presents a brief review of the literature and a numerical study investigating the model sensitivity to the above-mentioned properties for PCMHX with and without fins. CFD is used to evaluate the charging (melting) phenomena. The study quantifies the impact of uncertainty in these properties on the melting rate and temperature profile. Constant wall temperature was considered as heat source with no heat loss to ambient. Here, the results show that β and μ has significant effect on the melting rate and evolution of the melting front. For a non-finned domain when comparing results for different published values of β and μ, the deviation in melting time can be up to 12.9% and 57.6% respectively. For high wall temperatures, change in melting range did not impact melting time. But when the wall temperature is reduced, up to 9.8% deviation in melting time is observed.

CFD↗

Evaluation of Model Summertime Boundary Layer Cloud Development over Complex Terrain in New York State

Weather forecasting over complex terrain with diverse land cover is challenging. Utilizing the high-resolution observations from New York State Mesonet (NYSM), we are able to evaluate the surface processes of the Weather Research Forecast (WRF) Model in a detailed, scale-dependent manner. Here, in the study, possible impacts of land–atmosphere interaction on surface meteorology and boundary layer cloud development are investigated with different model resolutions, land surface models (LSMs), and planetary boundary layer (PBL) physical parameterizations. The High-Resolution Rapid Refresh, version 3 (HRRR), forecasting model is used as a reference for the sensitivity evaluation. Results show that over complex terrain, the high-resolution simulations (1 km × 60 vertical levels) generally perform better compared to low-resolution (3 km × 50 levels) in both surface meteorology and cloud fields. LSMs play a more important role in surface meteorology compared to PBL schemes. The NoahMP land surface model exhibits daytime warmer and drier biases compared to the Rapid Update Cycle (RUC) due to better prediction of the Bowen ratio in RUC. The PBL schemes would affect the convective strength in the boundary layer. The Shin–Hong (SH) scale-aware scheme tends to produce the strongest convective strength in the PBL, while the ACM2 PBL scheme rarely resolved convection even at 1-km resolution. By considering the radiation effect of subgrid-scale (SGS) clouds, the Mellor–Yamada–Nakanishi–Niino eddy diffusivity mass flux (MYNN-EDMF) predicted the highest cloud coverage and lowest surface solar radiation bias. The configuration of SGS clouds in MYNN-EDMF would not only significantly reduce shortwave radiation bias, but also affect the convection behaviors through land surface–cloud–radiation interaction.

54 ENVIRONMENTAL SCIENCES↗

Investigation of Scale-Dependent Groundwater/Surface-water Exchange in Rivers by Gradient Self-Potential Logging: Numerical Modeling and Field Experiments

Exchanges of groundwater and surface-water are fundamental to a wide range of water-supply and water-quality management issues but challenging to map beyond the reach scale. Waterborne gradient self-potential (SP) measurements are directly sensitive to water flow through riverbed sediments and can be used to infer exchange locations, direction (gain versus loss), scale, and relative changes, but to date applications to river corridor hydrology are limited. Numerical modeling and field experiments were therefore performed herein, each emphasizing waterborne gradient SP logging for identifying and locating focused vertical groundwater discharge (surface-water gain) and recharge (surface-water loss) in a river. Two and three-dimensional numerical models were constructed to simulate the polarities, appearances, and peak amplitudes of streaming-potential and electric-field anomalies on a riverbed and in the surface-water that were attributable to steady-state vertical fluxes of groundwater through high-permeability conduits in the riverbed. Effects of varied hydraulic length-scale of exchange and surface-water depth were tested through numerical modeling. Modeling results aided in data acquisition and interpretation for three repeated field experiments performed along a 1.5–2.0 km reach of the Quashnet River in Cape Cod, Massachusetts, where focused, meter-scale groundwater discharges occur at discrete locations within otherwise ubiquitous and more diffuse groundwater upwelling conditions. Strong gradient SP anomalies were repeatedly measured in the Quashnet River at previously confirmed locations of focused groundwater discharge, showing the efficacy of waterborne gradient SP logging in identifying and characterizing groundwater/surface water exchange dynamics at multiple river network scales.

Engineering↗

Modeling Clustered DNA Damage by Ionizing Radiation Using Multinomial Damage Probabilities and Energy Imparted Spectra

Simple and complex clustered DNA damage represent the critical initial damage caused by radiation. In this paper, a multinomial probability model of clustered damage is developed with probabilities dependent on the energy imparted to DNA and surrounding water molecules. The model consists of four probabilities: (A) direct damage of sugar-phosphate moieties leading to SSB, (B) OH− radical formation with subsequent SSB and BD formation, (C) direct damage to DNA bases, and (D) energy imparted to histone proteins and other molecules in a volume not leading to SSB or BD. These probabilities are augmented by introducing probabilities for the relative location of SSB using a ≤10 bp criteria for a double-strand break (DSB) and for the possible success of a radical attack that leads to SSB or BD. Model predictions for electrons, 4He, and 12C ions are compared to the experimental data and show good agreement. Thus, the developed model allows an accurate and rapid computational method to predict simple and complex clustered DNA damage as a function of radiation quality and to explore the resulting challenges to DNA repair.

Biochemistry & Molecular Biology↗

Multiscale Models for Fibril Formation: Rare Events Methods, Microkinetic Models, and Population Balances

Amyloid fibrils are thought to grow by a two-step dock-lock mechanism. However, previous simulations of fibril formation (i) overlook the bi-molecular nature of the docking step and obtain rates with first-order units, or (ii) superimpose the docked and locked states when computing the potential of mean force for association and thereby muddle the docking and locking steps. Here, we developed a simple microkinetic model with separate locking and docking steps and with the appropriate concentration dependences for each step. We constructed a simple model comprised of chiral dumbbells that retains qualitative aspects of fibril formation. We used rare events methods to predict separate docking and locking rate constants for the model. The rate constants were embedded in the microkinetic model, with the microkinetic model embedded in a population balance model for “bottom-up” multiscale fibril growth rate predictions. These were compared to “top-down” results using simulation data with the same model and multiscale framework to obtain maximum likelihood estimates of the separate lock and dock rate constants. We used the same procedures to extract separate docking and locking rate constants from experimental fibril growth data. Our multiscale strategy, embedding rate theories, and kinetic models in conservation laws should help to extract docking and locking rate constants from experimental data or long molecular simulations with correct units and without compromising the molecular description.

Shayesteh Zadeh, Armin (ORCID:0000000299896710)↗

Multiphysics Time-Integration for Turbulent Combustion at the Exascale

Turbulent reacting flow systems are often modeled with coupled time-dependent partial differential equations (PDEs). Solving such equations can easily tax the world's largest supercomputers. One pragmatic strategy for attacking such problems is to split the PDEs into components that can more easily be solved in isolation. This generic operator-splitting strategy leads to a set of ordinary differential equations (ODEs) that need to be solved as part of an "outer-loop" time-stepping approach. In many combustion applications, the ODEs to be solved can be very stiff, exhibiting timescales that span many orders of magnitude. The SUNDIALS library provides a plethora of robust time integration algorithms for solving these ODEs on exascale-capable computing hardware, yet for many complex applications (such multicomponent fuels or emissions predictions), the chemical models remain too complex to solve using reasonable resources. The Quasi-Steady State Approximation (QSSA) can be an effective tool for reducing the size and stiffness of the simulations. In this talk, I will discuss the use of the SUDIALS library of ODE solvers together with automatic code generation tools to solve complex turbulent reacting flow problems using QSSA models.

chemistry↗

Physics-Informed Neural Network Modeling of Li-Ion Batteries: Preprint

Li-ion batteries (LIB) are a promising solution to enable storage of intermittent energy sources due to their high energy density. However, LIBs are known to significantly degrade after about 1000 charge-discharge cycles. LIBs degrade following different degradation modes and at a rate that depends on the operating conditions (external temperature, load). To plan the installation of batteries, appropriate understanding and prediction capabilities of their lifecycle is needed. High-fidelity numerical models of LIBs such as the pseudo-two-dimensional (P2D) model have been shown to accurately represent the charge-discharge-cycle of an LIB given ac- curate choice of the physical parameters. Given the large number of P2D parameters, adjusting them using forward runs is intractable. This work describes the development of a physics-informed neural network (PINN) as a surrogate substitute of the P2D model that captures parameter dependence. The PINN is advantageous as it needs little to no data, and can naturally encode the dependence of every model parameter. Here, a specific training procedure is adopted to efficiently cover parameter space, handle model stiffness and enforce boundary conditions. The trained PINN is validated against numerical solutions of the P2D model, and its applicability to battery degradation modeling is discussed.

battery degradation↗

Many-body factorization and position–momentum equivalence of nuclear short-range correlations

While mean-field approximations, such as the nuclear shell model, provide a good description of many bulk nuclear properties, they fail to capture the important effects of nucleon–nucleon correlations such as the short-distance and high-momentum components of the nuclear many-body wave function1. Here, we study these components using the effective pair-based generalized contact formalism and ab initio quantum Monte Carlo calculations of nuclei from deuteron to 40 Ca. We observe a universal factorization of the many-body nuclear wave function at short distance into a strongly interacting pair and a weakly interacting residual system. The residual system distribution is consistent with that of an uncorrelated system, showing that short-distance correlation effects are predominantly embedded in two-body correlations. Spin- and isospin-dependent ‘nuclear contact terms’ are extracted in both coordinate and momentum space for different realistic nuclear potentials. The contact coefficient ratio between two different nuclei shows very little dependence on the nuclear interaction model. These findings thus allow extending the application of mean-field approximations to short-range correlated pair formation by showing that the relative abundance of short-range pairs in the nucleus is a long-range (that is, mean field) quantity that is insensitive to the short-distance nature of the nuclear force.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Active and Transfer Learning of High-Dimensional Neural Network Potentials for Transition Metals

Classical molecular dynamics (MD) simulations represent a very popular and powerful tool for materials modeling and design. The predictive power of MD hinges on the ability of the interatomic potential to capture the underlying physics and chemistry. There have been decades of seminal work on developing interatomic potentials, albeit with a focus predominantly on capturing the properties of bulk materials. Such physics-based models, while extensively deployed for predicting the dynamics and properties of nanoscale systems over the past two decades, tend to perform poorly in predicting nanoscale potential energy surfaces (PESs) when compared to high-fidelity first-principles calculations. These limitations stem from the lack of flexibility in such models, which rely on a predefined functional form. Machine learning (ML) models and approaches have emerged as a viable alternative to capture the diverse size-dependent cluster geometries, nanoscale dynamics, and the complex nanoscale PESs, without sacrificing the bulk properties. Here, in this study, we introduce an ML workflow that combines transfer and active learning strategies to develop high-dimensional neural networks (NNs) for capturing the cluster and bulk properties for several different transition metals with applications in catalysis, microelectronics, and energy storage, to name a few. Our NN first learns the bulk PES from the high-quality physics-based models in literature and subsequently augments this learning via retraining with a higher-fidelity first-principles training data set to concurrently capture both the nanoscale and bulk PES. Our workflow departs from status-quo in its ability to learn from a sparsely sampled data set that nonetheless covers a diverse range of cluster configurations from near-equilibrium to highly nonequilibrium as well as learning strategies that iteratively improve the fingerprinting depending on model fidelity. All the developed models are rigorously tested against an extensive first-principles data set of energies and forces of cluster configurations as well as several properties of bulk configurations for 10 different transition metals. Our approach is material agnostic and provides a methodology to transfer and build upon the learnings from decades of seminal work in molecular simulations on to a new generation of ML-trained potentials to accelerate materials discovery and design.

36 MATERIALS SCIENCE↗

Residual elastic strain evolution due to thermal cycling of a ceramic-metal composite (WC-Cu) via high energy X-ray diffraction and analytical modeling

Residual stress, when superimposed with in-service loading, can significantly reduce the lifetime and performance of a component. Ceramic-metal composites are susceptible to residual stresses due to the thermal expansion mismatch of the ceramic and metallic phases. The WC-Cu composite explored in the present study provides a promising combination of thermal conductivity and strength properties, while exhibiting counterintuitive improvements in strength and ductility after thermal cycling. Further, this work quantifies the evolution of the residual elastic strains as a result of processing and cyclic thermal loading in a co-continuous WC-Cu composite through experimental high energy X-ray diffraction and kinetics-based modeling. Both analyses indicate that processing-induced residual tensile stress in the copper phase is relieved upon subsequent thermal cycling, with kinetics modeling revealing the cyclic-dependent nature of the active power-law creep mechanisms. The results indicate that, through stress relaxation, this material system maintains structural stability during thermal cycling. The illustrated kinetics of relaxation can inform general material processors and designers of ceramic-metal composites to minimize detrimental residual stress and improve performance of these material systems.

36 MATERIALS SCIENCE↗

Statistical tools for a better optical model

Background: Modern statistical tools provide the ability to compare the information content of observables and provide a path to explore which experiments would be most useful to give insight into and constrain theoretical models. Purpose: Here we study three such tools in the context of nuclear reactions with the goal of constraining the optical potential. Method: The three statistical tools examined are (i) the principal component analysis, (ii) the sensitivity analysis based on derivatives, and (iii) the Bayesian evidence. We first apply these tools to a toy-model case, comparing the form of the imaginary part of the optical potential. Then we consider two different reaction observables, elastic angular distributions and polarization data for reactions on 48 Ca and 208 Pb at two different beam energies. Results: For the toy-model case, we find significant discrimination power in the sensitivities and the Bayesian evidence, showing clearly that the volume imaginary term is more useful to describe scattering at higher energies. When comparing between elastic cross sections and polarization data using realistic optical models, sensitivity studies indicate that both observables are roughly equally sensitive but the variability of the optical model parameters is strongly angle dependent. The Bayesian evidence shows some variability between the two observables, but the Bayes factor obtained is not sufficient to discriminate between angular distributions and polarization. Conclusions: From the cases considered, we conclude that, in general, elastic scattering angular distributions have similar impact in constraining the optical potential parameters compared with the polarization data. The angular ranges for the optimum experimental constraints can vary significantly with the observable considered.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Shifts in hydroclimatology of US megaregions in response to climate change

Most of the population and economic growth in the United States occurs in megaregions as the clustered metropolitan areas, whereas climate change may amplify negative impacts on water and natural resources. This study assesses shifts in regional hydroclimatology of fourteen US megaregions in response to climate change over the 21st century. Hydroclimatic projections were simulated using the Variable Infiltration Capacity (VIC) model driven by three downscaled climate models from the Multivariate Adaptive Constructed Analogs (MACA) dataset to cover driest to wettest future conditions in the conterminous United States (CONUS). Shifts in the regional hydroclimatolgy and basin characteristics of US megaregions were represented as a combination of changes in the aridity and evaporative indices using the Budyko framework and Fu's equation. Changes in the climate types of US megaregions were estimated using the Fine Gaussian Support Vector Machine (SVM) method. The results indicate that Los Angeles, San Diego, and San Francisco are more likely to experience less arid conditions with some shifts from Continental to Temperate climate type while the hydroclimatology of Houston may become drier with some shifts from Temperate to Continental climate type. Additionally, water yield is likely to decrease in Seattle. Change in the hydroclimatology of Denver and Phoenix highly depends on the selected climate model. However, the basin characteristics of Phoenix have the highest sensitivity to climate change. Overall, the hydroclimatic conditions of Los Angeles, San Diego, Phoenix, Denver, and Houston have the highest sensitivity to climate change. Understanding of future shifts in hydroclimatology of megaregions can help decision-makers to attenuate negative consequences by implementing appropriate adaptation strategies, particularly in the water-scare megaregions.

54 ENVIRONMENTAL SCIENCES↗

Updated neutrino mass constraints from galaxy clustering and CMB lensing-galaxy cross-correlation measurements

In this work, we revisit cosmological constraints on the sum of the neutrino masses Σm v from a combination of full-shape BOSS galaxy clustering [P(k)] data and measurements of the cross-correlation between Planck Cosmic Microwave Background (CMB) lensing convergence and BOSS galaxy overdensity maps [$C_{l}^{kg}$], using a simple but theoretically motivated model for the scale-dependent galaxy bias in auto- and cross-correlation measurements. We improve upon earlier related work in several respects, particularly through a more accurate treatment of the correlation and covariance between P(k) and $C_{l}^{kg}$ measurements. When combining these measurements with Planck CMB data, we find a 95% confidence level upper limit of Σm v sum of neutrino masses < 0.14 eV, while slightly weaker limits are obtained when including small-scale ACTPol CMB data, in agreement with our expectations. We confirm earlier findings that (once combined with CMB data) the full-shape information content is comparable to the geometrical information content in the reconstructed BAO peaks given the precision of current galaxy clustering data, discuss the physical significance of our inferred bias and shot noise parameters, and perform a number of robustness tests on our underlying model. While the inclusion of $C_{l}^{kg}$ measurements does not currently appear to lead to substantial improvements in the resulting Σm v of neutrino mass constraints, we expect the converse to be true for near-future galaxy clustering measurements, whose shape information content will eventually supersede the geometrical one.

79 ASTRONOMY AND ASTROPHYSICS↗

Measurements of dielectron production in Au + Au collisions at s N N = 27 , 39, and 62.4 GeV from the STAR experiment

We report systematic measurements of dielectron (e + e – ) invariant-mass M ee spectra at midrapidity in Au + Au collisions at s N N = 27 = 27, 39, and 62.4 GeV taken with the STAR detector at the Relativistic Heavy Ion Collider. For all energies studied, a significant excess yield of dielectrons is observed in the low-mass region ( 0.40 < M ee < 0.75 MeV / c 2 ) compared to hadronic cocktail simulations at freeze-out. Models that include an in-medium broadening of the ρ -meson spectral function consistently describe the observed excess. In addition, we report acceptance-corrected dielectron-excess spectra for Au + Au collisions at midrapidity ( |y ee | < 1) in the 0–80% centrality bin for each collision energy. The integrated excess yields for 0.4 < M ee < 0.75 GeV / c 2 , normalized by the charged particle multiplicity at midrapidity, are compared with previously published measurements for Au + Au at s N N = 27 = 19.6 and 200 GeV. Models that include an in-medium broadening of the ρ -meson spectral function consistently describe the observed excess. Furthermore, the normalized excess yields in the low-mass region show no significant collision energy dependence. The data, however, are consistent with model calculations that demonstrate a modest energy dependence.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Interaction of Anisole on Alumina-Supported Ni and Mo Oxide Hydrodeoxygenation Catalysts

The conversion of biomass to transportation fuels and value-added chemicals is a promising method to reduce reliance on fossil fuels. Mo-based catalysts have been shown to be highly active in the hydrodeoxygenation of biomass-derived phenolic compounds. The catalyst active phase, surface species, and the effect of adding additional metals is not comprehensively understood. Here we compare the temperature-dependent adsorption behavior of the model compound anisole 2 on an alumina-supported mixed nickel molybdenum oxide catalyst with two reference catalysts of molybdenum oxide and nickel oxide. Raman spectroscopy showed that the catalysts contain significant amounts of molybdates and molybdoaluminates, in addition to NiMoO 4 in the nickel molybdenum catalyst and MoO 3 in the molybdenum-only catalyst. Using transmission infrared spectroscopy under a controlled environment, we find that anisole chemisorbed largely through the oxygen in the methoxy group to form surface-bound phenoxy and methoxy species on all the catalysts. Ambient pressure X-ray photoelectrons spectroscopy measurements of the catalysts in anisole vapor showed reduced Mo atoms are the binding sites. The surface interaction and removal temperature of these species varied with the metal composition. The MoO x component dominated the adsorption behavior in both MoO x and NiMoO x catalysts. The formation of new aromatics, including methylated rings, depended on Ni composition. Upon adding hydrogen to induce the hydrodeoxygenation of anisole, undesirable polynuclear aromatic species were quickly formed on the Mo-based catalysts. In conclusion, these results suggest that the molybdenum oxide controls the adsorption and reactivity of the surface species, with a cooperative effect by Ni.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗