SEARCH · Search NASA

Results for “Describing function”

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 307 records · Page 17

All sources NASA NTRS DOE OSTI

DOE OSTI2019

Fermi-Löwdin orbital self-interaction correction using the strongly constrained and appropriately normed meta-GGA functional

Despite the success of density functional approximations (DFAs) in describing the electronic properties of many-electron systems, the most widely used approximations suffer from self-interaction errors (SIE) that limit their predictive power. Herein we describe the effects of removing SIE from the strongly constrained and appropriately normed (SCAN) meta-generalized gradient approximation (GGA) using the Fermi-Lowdin Orbital Self-Interaction Correction (FLOSIC) method. FLOSIC is a size-extensive implementation of the Perdew-Zunger self-interaction correction (PZ-SIC) formalism. We find that FLOSIC-SCAN calculations require careful treatment of numerical details and describe an integration grid that yields reliable accuracy with this approach. We investigate the performance of FLOSIC-SCAN for predicting a wide array of properties and find that it provides better results than FLOSIC-LDA and FLOSIC-PBE in nearly all cases. Furthermore, it gives better predictions than SCAN for orbital energies and dissociation energies where self-interaction effects are known to be important, but total energies and atomization energies are made worse. For these properties, we also investigate the use of the self-consistent FLOSIC-SCAN density in the SCAN functional and find that this DFA@FLOSIC-DFA approach yields improved results compared to pure, self-consistent SCAN calculations. Thus FLOSIC-SCAN provides improved results over the parent SCAN functional in cases where SIEs are dominant, and even when they are not, if the SCAN@FLOSIC-SCAN method is used.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Effective description of sub-maximal chaos: stringy effects for SYK scrambling

It has been proposed that the exponential decay and subsequent power law saturation of out-of-time-order correlation functions can be universally described by collective ‘scramblon’ modes. We develop this idea from a path integral perspective in several examples, thereby establishing a general formalism. After reformulating previous work on the Schwarzian theory and identity conformal blocks in two-dimensional CFTs relevant for systems in the infinite coupling limit with maximal quantum Lyapunov exponent, we focus on theories with sub-maximal chaos: we study the large-q limit of the SYK quantum dot and chain, both of which are amenable to analytical treatment at finite coupling. In both cases we identify the relevant scramblon modes, derive their effective action, and find bilocal vertex functions, thus constructing an effective description of chaos. The final results can be matched in detail to stringy corrections to the gravitational eikonal S-matrix in holographic CFTs, including a stringy Regge trajectory, bulk to boundary propagators, and multi-string effects that are unexplored holographically.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2022

J-PLUS: Spectral evolution of white dwarfs by PDF analysis

We estimated the spectral evolution of white dwarfs with effective temperature using the Javalambre Photometric Local Universe Survey (J-PLUS) second data release (DR2), which provides 12 photometric optical passbands over 2176 deg 2 . We analyzed 5926 white dwarfs with r ≤ 19.5 mag in common between a white dwarf catalog defined from Gaia EDR3 and J-PLUS DR2. We performed a Bayesian analysis by comparing the observed J-PLUS photometry with theoretical models of hydrogen- and helium-dominated atmospheres. We estimated the probability distribution functions for effective temperature (T eff ), surface gravity, parallax, and composition; and the probability of having a H-dominated atmosphere (p H ) for each source. We applied a prior in parallax, using Gaia EDR3 measurements as a reference, and derived a self-consistent prior for the atmospheric composition as a function of T eff . We described the fraction of white dwarfs with a He-dominated atmosphere (f He ) with a linear function of the effective temperature at 5000 < T eff < 30 000 K. We find f He = 0.24 ± 0.01 at T eff = 10 000 K, a change rate along the cooling sequence of 0.14 ± 0.02 per 10 kK, and a minimum He-dominated fraction of 0.08 ± 0.02 at the high-temperature end. We tested the obtained p H by comparison with spectroscopic classifications, finding that it is reliable. We estimated the mass distribution for the 351 sources with distance d < 100 pc, mass M > 0.45 M ⊙ , and T eff > 6000 K. The result for H-dominated white dwarfs agrees with previous studies, with a dominant M = 0.59 M ⊙ peak and the presence of an excess at M ~ 0.8 M ⊙ . This high-mass excess is absent in the He-dominated distribution, which presents a single peak. The J-PLUS optical data provide a reliable statistical classification of white dwarfs into H- and He-dominated atmospheres. We find a 21 ± 3% increase in the fraction of He-dominated white dwarfs from T eff = 20 000 K to T eff = 5000 K.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

Heterosynaptic plasticity in biomembrane memristors controlled by pH

Abstract In biology, heterosynaptic plasticity maintains homeostasis in synaptic inputs during associative learning and memory, and initiates long-term changes in synaptic strengths that nonspecifically modulate different synapse types. In bioinspired neuromorphic circuits, heterosynaptic plasticity may be used to extend the functionality of two-terminal, biomimetic memristors. In this article, we explore how changes in the pH of droplet interface bilayer aqueous solutions modulate the memristive responses of a lipid bilayer membrane in the pH range 4.97–7.40. Surprisingly, we did not find conclusive evidence for pH-dependent shifts in the voltage thresholds ( V* ) needed for alamethicin ion channel formation in the membrane. However, we did observe a clear modulation in the dynamics of pore formation with pH in time-dependent, pulsed voltage experiments. Moreover, at the same voltage, lowering the pH resulted in higher steady-state currents because of increased numbers of conductive peptide ion channels in the membrane. This was due to increased partitioning of alamethicin monomers into the membrane at pH 4.97, which is below the pKa (~5.3–5.7) of carboxylate groups on the glutamate residues of the peptide, making the monomers more hydrophobic. Neutralization of the negative charges on these residues, under acidic conditions, increased the concentration of peptide monomers in the membrane, shifting the equilibrium concentrations of peptide aggregate assemblies in the membrane to favor greater numbers of larger, increasingly more conductive pores. It also increased the relaxation time constants for pore formation and decay, and enhanced short-term facilitation and depression of the switching characteristics of the device. Modulating these thresholds globally and independently of alamethicin concentration and applied voltage will enable the assembly of neuromorphic computational circuitry with enhanced functionality. Impact statement We describe how to use pH as a modulatory “interneuron” that changes the voltage-dependent memristance of alamethicin ion channels in lipid bilayers by changing the structure and dynamical properties of the bilayer. Having the ability to independently control the threshold levels for pore conduction from voltage or ion channel concentration enables additional levels of programmability in a neuromorphic system. In this article, we note that barriers to conduction from membrane-bound ion channels can be lowered by reducing solution pH, resulting in higher currents, and enhanced short-term learning behavior in the form of paired-pulse facilitation. Tuning threshold values with environmental variables, such as pH, provide additional training and learning algorithms that can be used to elicit complex functionality within spiking neural networks. Graphical abstract

36 MATERIALS SCIENCE↗

DOE OSTI2020

Automation of Active Space Selection for Multireference Methods via Machine Learning on Chemical Bond Dissociation

Predicting and understanding the chemical bond is one of the major challenges of computational quantum chemistry. Kohn–Sham density functional theory (KS-DFT) is the most common method, but approximate density functionals may not be able to describe systems where multiple electronic configurations are equally important. Multiconfigurational wave functions, on the other hand, can provide a detailed understanding of the electronic structures and chemical bonds of such systems. In the complete active space self-consistent field (CASSCF) method, one performs a full configuration interaction calculation in an active space consisting of active electrons and active orbitals. However, CASSCF and its variants require the selection of these active spaces. This choice is not black box; it requires significant experience and testing by the user, and thus active space methods are not considered particularly user-friendly and are employed only by a minority of quantum chemists. Our goal is to popularize these methods by making it easier to make good active space choices. Here, we present a machine learning protocol that performs an automated selection of active spaces for chemical bond dissociation calculations of main group diatomic molecules. The protocol shows high prediction performance for a given target system as long as a properly correlated system is chosen for training. Good active spaces are correctly predicted with a considerably better success rate than random guess (larger than 80% precision for most systems studied). Our automated machine learning protocol shows that a “black-box” mode is possible for facilitating and accelerating the large-scale calculations on multireference systems where single-reference methods such as KS-DFT cannot be applied.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Feasibility of a Novel Density Functional Method Outside the Kohn-Sham Framework for Modeling Global Potential Energy Surfaces of Molecular Chemical Reactions (Final Technical Report)

The project aimed towards the construction of computational methods capable of modeling the global potential energy surface (PES) of small molecules, molecular ions, and radicals — including the parts of the PES which correspond to chemical reactions, and the reaction paths connecting reactants, intermediates, and reaction products. This goal may seem humble at first glance, but for chemical systems with more than about six atoms in total, at the time the project was proposed, there were no established theoretical methods capable of simulating such systems reliably, not even for small molecules in the gas phase which are electronically benign. This restriction severely hampers our ability to model and control chemical processes under harsh conditions. With the goal of constructing a method capable of modeling such chemical systems, we proposed to pursue a novel approach towards a Multi-Configuration (MC) DFT outside the traditional frameworks of Kohn-Sham theory and other methods of coupling wave function theory with DFT. Rather than being a complete active space (CAS) method, the proposed DOCI-DFT would employ a special restricted form of the active space wave function, called Doubly-Occupied CI (DOCI)—this wave function form is sufficient to describe not only heterolytic, but also homolytic bond dissociation processes at the zeroth order (i.e., as active space wave function). There are no other standard mean field methods which can do so. In the original proposal, we also outlined strong formal and practical arguments speaking for this method. The project proposal was accepted by the Department of Energy and provided two years of seed funding for one graduate student, as well as two years of two weeks PI summer-salary for the PI. However, ultimately the project could not be effectively pursued due to severe interferences outside the context scientific problems and was cancelled. Apart from two weeks of PI salary (and associated fringe & overhead costs issues by the performing institution), no costs were charged against the award, with the entire rest of funding returned to the Department of Energy. This report reiterates the primary background information regarding the project, its original goals, and outlines the preliminary work performed during the two weeks of DOE funding charged.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Driving mode analysis—How uncertain functional inputs propagate to an output

Abstract Driving mode analysis elucidates how correlated features of uncertain functional inputs jointly propagate to produce uncertainty in the output of a computation. Uncertain input functions are decomposed into three terms: the mean functions, a zero‐mean driving mode, and zero‐mean residual. The random driving mode varies along a single direction, having fixed functional shape and random scale. It is uncorrelated with the residual, and under linear error propagation, it produces an output variance equal to that of the full input uncertainty. Finally, the driving mode best represents how input uncertainties propagate to the output because it minimizes expected squared Mahalanobis distance amongst competitors. These characteristics recommend interpretation of the driving mode as the single‐degree‐of‐freedom component of input uncertainty that drives output uncertainty. We derive the functional driving mode, show its superiority to other seemingly sensible definitions, and demonstrate the utility of driving mode analysis in an application. The application is the simulation of neutron transport in criticality experiments. The uncertain input functions are nuclear data that describe how Pu reacts to bombardment by neutrons. Visualization of the driving mode helps scientists understand what aspects of correlated functional uncertainty have effects that either reinforce or cancel one another in propagating to the output of the simulation.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

Bootstrapping a stress-tensor form factor through eight loops

We bootstrap the three-point form factor of the chiral stress-tensor multiplet in planar $\mathcal{N}$ = 4 supersymmetric Yang-Mills theory at six, seven, and eight loops, using boundary data from the form factor operator product expansion. This may represent the highest perturbative order to which multi-variate quantities in a unitary four-dimensional quantum field theory have been computed. In computing this form factor, we observe and employ new restrictions on pairs and triples of adjacent letters in the symbol. We provide details about the function space required to describe the form factor through eight loops. Plotting the results on various lines provides striking numerical evidence for a finite radius of convergence of perturbation theory. By the principle of maximal transcendentality, our results are expected to give the highest weight part of the gg → Hg and H → ggg amplitudes in the heavy-top limit of QCD through eight loops. These results were also recently used to discover a new antipodal duality between this form factor and a six-point amplitude in the same theory.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2023

Measurement of the Higgs boson inclusive and differential fiducial production cross sections in the diphoton decay channel with pp collisions at $\sqrt{s}$ = 13 TeV

The measurements of the inclusive and differential fiducial cross sections of the Higgs boson decaying to a pair of photons are presented. The analysis is performed using proton-proton collisions data recorded with the CMS detector at the LHC at a centre-of-mass energy of 13 TeV and corresponding to an integrated luminosity of 137 fb -1 . The inclusive fiducial cross section is measured to be σ fid = 73.4$^{+5.4}_{-5.3}$ (stat)$a^{+2.4}_{-2.2 }$(syst) fb, in agreement with the standard model expectation of 75.4 ± 4.1 fb. The measurements are also performed in fiducial regions targeting different production modes and as function of several observables describing the diphoton system, the number of additional jets present in the event, and other kinematic observables. Two double differential measurements are performed. No significant deviations from the standard model expectations are observed.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2024

Phase Selection During Solidification and Solid-State Phase Transformations in an Al-10Ce-8Mn (wt pct) Alloy

In multicomponent Al-Ce alloys, and especially after additive manufacturing (AM), complex and metastable solidification microstructures are frequently observed. Here, in this research, the relationship between solidification conditions and phase selection are explored for an Al-10Ce-8Mn (wt pct) alloy using a systematic study of laser melting conditions. Three solidification modes were observed: primary Al 10 Mn 2 Ce; primary Al 20 Mn 2 Ce; and eutectic FCC Al + Al 20 Mn 2 Ce. These solidification modes were correlated to specific liquid-solid interface velocities using a simple thermal model, showing the change in primary solidification phase for low (< 6.8 × 10 −4 m/s), moderate (between 8.2 × 10 −4 and 5.9 × 10 −2 m/s) and high solidification velocities (> 6.2 x 10 −2 m/s) for the above three solidification microstructures, respectively. These results were rationalized by using interface response function (IRF) theory to describe the solidification undercooling for the possible primary intermetallic phases. The implication of the local phase selection from differing solidification conditions is summarized by a comparison of hardness which demonstrates the potential variance of Vickers hardness from 101 to 242 (VHV) by changing the laser velocity from 1 to 83 mm/s. Interestingly, on heat treatment at 400°C, the decomposition pathways of the solidification microstructure and hardness were also found to be different, thereby opening multiple pathways for spatial microstructure and property control within AM components.

Sisco, Kevin↗

DOE OSTI2022

Concurrent Precipitation of Nb(C,N) and Metastable M 23 C 6 in Alloy 347H at 700°C and 750°C: Computer Simulations and Comparison to Experiment

Here, we present our results for the concurrent precipitation of metastable M 23 C 6 , Nb(C,N) secondary precipitates, and the Nb(C,N) primary crystals in 347H austenitic stainless steel. For precipitation modeling, we have accounted for the elastic contribution to interfacial energy, and for the Fe-spin-polarization for NbC/Fe and M 23 C 6 /Fe interfacial energy values: for NbC/Fe ~ 0.63 J/m 2 . For M 23 C 6 precipitates, an error function was used to describe the interfacial energy growth with particle size. In precipitation simulations, the average size of the primary Nb(C,N) particles remained ~ 1 μm at 700°C and ~ 0.3 μm at 750°C. The M 23 C 6 precipitates at 750°C dissolved after 120 h (our simulations) compared to 300 h (experiments). The Nb(C,N)/Fe interfacial energy was not affected by the nitrogen additions. With these modifications, reasonable agreement with the available experimental data was obtained, which allows using them in the development of the 2nd-phase particle-informed creep theory.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Electrochemical grand potential-based phase-field simulation of electric field-assisted sintering

Here, an electrochemical grand potential functional was proposed to describe the sintering of an ionic ceramic green body. The resultant phase-field description enables simulation of the consolidation of an arbitrary number of granular particles and their interactions with the surrounding void phase. The model includes the effects of charged vacancies and the associated interactions between internal and applied electric fields. Defect segregation to grain boundaries is also accounted for, as well as enhanced interfacial defect mobilities. The model was parameterized for Y 2 O 3 . Simulations of two-particle systems showed that the applied electric field had an increasingly important impact on neck growth as particle size increased. A sudden rapid increase in temperature occurred for larger field strengths, which has been reported to be correlated to the onset of a flash event in flash sintering. Simulations of many particles showed that internal heat generation by Joule heating was localized at particle–particle contacts (grain boundaries), even though their conductivities were lower than nearby internal particle-void interfaces. A percolative path for ionic charge across the green body and the ceramic sintered solid was thus defined, accelerating the Joule heating process as the porosity of the green body is removed.

36 MATERIALS SCIENCE↗

DOE OSTI2019

A solvable problem in statistical mechanics: The dipole-type Hamiltonian mean field model

The present study documents a type of mean field approximation inspired by the dipole interaction model, which is analytically solved in the canonical and microcanonical ensembles. The current calculations were derived from the Hamiltonian mean field model developments. After describing the canonical partition function, the free and internal energies, magnetization, and specific heat are derived and graphically depicted. In the microcanonical ensemble, the entropy is calculated as well as other thermodynamic properties. The system shows a second-order phase transition emphasizing that both methods coincide, which is only valid only in equilibrium. In addition, the current model represents a nonsymmetric Hamiltonian mean field model that shows a phase transition.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2019

Bootstrap-based confidence interval estimation for thermal security region of bulk power grid

Thermal security region (TSR) is a powerful tool for monitoring and controlling the thermal security of bulk power grid with high penetration of renewable energy. One of the major challenges of TSR in practical application is to obtain the exact analytical expression of the TSR boundary (TSRB) because TSRB is a nonlinear hypersurface described by an implicit function. For this reason, the TSRB is usually approximated by hyperplane instead of getting the exact analytical expression. Traditionally, the coefficients of the hyperplane are estimated with a set of original TSRB points using least square estimation (LSE). However, LSE is a point estimation method, which is incapable of evaluating the quality of the estimation. In addition, the accuracy and computational efficiency of the hyperplane approximation are influenced by the number of TSRB points significantly. Here, a bootstrap based confidence interval estimation is proposed to estimate not only the coefficients of TSRB approximation hyperplane, but also the standard deviations and confidence intervals of the coefficients for evaluating the quality and reliability of the approximation results. First, empirical distribution functions (EDFs) of the coefficients of TSRB approximation hyperplane are approximated from a set of original TSRB points by using residual resampling bootstrap method. Then, the EDFs of the coefficients are employed to estimate the coefficients of the approximation hyperplane. Meanwhile, the standard deviations and confidence intervals of the estimated coefficients of the hyperplanes are also calculated for evaluating the quality and reliability of the approximation. The proposed approach is tested on the China Southern Power Grid (CSG). Results of simulations validate the accuracy and efficiency of the proposed approach in approximating the TSRB.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2022

Characterization of Multianode Photomultiplier Tubes for use in the CLAS12 RICH detector

In this work, we present results of the detailed study of several hundred Hamamatsu H12700 Multianode Photomultiplier Tubes (MaPMTs), characterizing their response to the Cherenkov light photons in the second Ring Imaging Cherenkov detector, a part of the CLAS12 upgrade at Jefferson Lab. The total number of pixels studied was 25536. The single photoelectron spectra were measured for each pixel at different high voltages and light intensities of the laser test setup. Using the same dedicated front-end electronics as in the first RICH detector, the setup allowed us to characterize each pixel’s properties such as gain, quantum efficiency, signal crosstalk between neighboring pixels, and determine the signal threshold values to optimize their efficiency to detect Cherenkov photons. A recently published state-of-the-art mathematical model, describing photon detector response functions measured in low light conditions, was extended to include the description of the crosstalk contributions to the spectra. The database of extracted parameters will be used for the final selection of the MaPMTs, their arrangement in the new RICH detector, and the optimization of the operational settings of the front-end electronics. The results show that the characteristics of the H12700 MaPMTs satisfy our requirements for the position-sensitive single photoelectron detectors.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2022

Heat transport and temperature boundary-layer profiles in closed turbulent Rayleigh–Bénard convection with slippery conducting surfaces

Here, we report direct numerical simulations (DNS) of the Nusselt number Nu , the vertical profiles of mean temperature Θ(z) and temperature variance Ω(z) across the thermal boundary layer (BL) in closed turbulent Rayleigh–Bénard convection (RBC) with slippery conducting surfaces ( z is the vertical distance from the bottom surface). The DNS study was conducted in three RBC samples: a three-dimensional cuboid with length L=H and width W=H/4 ( H is the sample height), and two-dimensional rectangles with aspect ratios Γ≡L/H=1 and 10. The slip length b for top and bottom plates varied from 0 to ∞ . The Rayleigh numbers Ra were in the range 10 6 $&LessSlantEqual;$Ra$&LessSlantEqual;$10 10 and the Prandtl number Pr was fixed at 4.3 . As b increases, the normalised Nu/Nu 0 ( Nu 0 is the global heat transport for b =0 ) from the three samples for different Ra and Γ can be well described by the same function Nu/Nu 0 = N 0 tanh( b /λ 0 )+1 , with N 0 =0.8±0.03 . Here λ 0 ≡ L /(2 Nu 0 ) is the thermal boundary layer thickness for b =0 . Considering the BL fluctuations for Pr >1, one can derive solutions of temperature profiles Θ(z) and Ω(z) near the thermal BL for b $&GreaterSlantEqual;$ 0 . When b =0 , the solutions are equivalent to those reported by Shishkina et al. ( Phys. Rev. Lett. , vol. 114, 2015, 114302) and Wang et al. ( Phys. Rev. Fluids , vol. 1, 2016, 082301(R)), respectively, for no-slip plates. For b >0, the derived solutions are in excellent agreement with our DNS data for slippery plates.

42 ENGINEERING↗

DOE OSTI2024

Static Subspace Approximation for Random Phase Approximation Correlation Energies: Implementation and Performance

Developing theoretical understanding of complex reactions and processes at interfaces requires using methods that go beyond semilocal density functional theory to accurately describe the interactions between solvent, reactants and substrates. Methods based on many-body perturbation theory, such as the random phase approximation (RPA), have previously been limited due to their computational complexity. However, this is now a surmountable barrier due to the advances in computational power available, in particular through modern GPU-based supercomputers. In this work, we describe the implementation of RPA calculations within BerkeleyGW and show its favorable computational performance on large complex systems relevant for catalysis and electrochemistry applications. Our implementation builds off of the static subspace approximation which, by employing a compressed representation of the frequency dependent polarizability, enables the evaluation of the RPA correlation energy with significant acceleration and systematically controllable accuracy. We find that the computational cost of calculating the RPA correlation energy scales only linearly with system size for systems containing up to 50 thousand bands, and is expected to scale quadratically thereafter. We also show excellent strong scaling results across several supercomputers, demonstrating the performance and portability of this implementation.

algorithmic development↗

DOE OSTI2020

In Situ Studies of Methanol Decomposition Over Cu(111) and Cu2O/Cu(111): Effects of Reactant Pressure, Surface Morphology, and Hot Spots of Active Sites

The dissociative adsorption of methanol was investigated on Cu(111) and ultrathin Cu 2 O films. We employed synchrotron-based Ambient Pressure X-ray Photoelectron Spectroscopy (AP-XPS) and Scanning Tunneling Microscopy (STM) to study the dynamics of gas–solid interactions, and calculations based on Density Functional Theory (DFT) were used to examine the reaction path. C 1s XPS spectra revealed that methanol underwent dissociative adsorption on plain Cu(111) to form methoxy (CH 3 O), formaldehyde (H 2 CO), and formate (HCOO) at a pressure range of 0.5–10 mTorr, with these species remaining on the surface after evacuation. This was accompanied by the appearance of a low coverage (~0.05 ML) of O ads in the O 1s which can be considered a highly active site for methanol activation. The high activity is apparent by a coverage of 0.8 ML of methoxy at room temperature. STM was unable to image these species at room temperature as they were highly mobile on metallic copper. In contrast, for CH 3 OH on Cu 2 O/Cu(111), STM showed clear hot spots for reaction and a complex array of adsorption structures. On the oxide substrate, there was decomposition of methanol to H 2 CO, CH 3 O, HCOO, and hydrocarbon species (CH x ) due to the subsequent interactions of methanol with lattice oxygen. Cu(111) remained entirely saturated with decomposition products under 10 mTorr of methanol (θ ≈ 0.97 ML), whereas the Cu 2 O overlayer was saturated at a much lower coverage (θ ≈ 0.30 ML). STM revealed rows and step edges of Cu 2 O decorated with decomposition products and metallic Cu islands ~5 nm in size. The difference in activity between Cu(111) and Cu 2 O/Cu(111) is attributed to the significant amount of O present on the oxide surface. Additionally, Density Functional theory (DFT) calculations described the XPS measurements well, showing a likely methanol dissociation to *CH 3 O and therefore a surface reduction. More importantly, the DFT results revealed that it was the chemisorbed oxygen on Cu 2 O/Cu(111) which oxidized the dissociated *CH 3 O to *HCOO and eventually CO 2 , while the reaction only involving upper oxygen on the Cu 2 O hexagonal ring led to the formation of H 2 CO.

36 MATERIALS SCIENCE↗