Search NASA⌕ Search

SEARCH · Search NASA

Results for “Structural Model Decomposition”

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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Dynamic mode decomposition for extrapolating nonequilibrium Green's-function dynamics

The Hartree-Fock generalized Kadanoff-Baym ansatz (HF-GKBA) offers an approximate numerical procedure for propagating the two-time nonequilibrium Green's function (NEGF). Here, using the GW self-energy, we compare the HF-GKBA to exact results for a variety of systems with long- and short-range interactions, different two-body interaction strengths, and various nonequilibrium preparations. We find excellent agreement between the HF-GKBA and exact time evolution in models when more realistic long-range exponentially decaying interactions are considered. Furthermore, this agreement persists for long times and for intermediate to strong interaction strengths. In large systems, HF-GKBA becomes prohibitively expensive for long-time evolutions. For this reason, we look at the use of dynamical mode decomposition (DMD) to reconstruct long-time NEGF trajectories from a sample of the initial trajectory. Using no more than 16% of the total time evolution, we reconstruct the total trajectory with high fidelity. Our results show the potential for DMD to be used in conjunction with HF-GKBA to calculate long-time trajectories in large-scale systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Electron and Electron-Nuclear Dynamics of Open Quantum Subsystems

Leveraging an open-subsystem formulation of Density Functional Theory (DFT) we aim at describing periodic and molecular systems alike, including their electronic and nuclear dynamics. Subsystem DFT enables first principles simulations to approach realistic time- and length-scales, and most importantly sheds light on the dynamical behavior of complex systems. Taking subsystem DFT to the time domain allows us to inspect the electron dynamics of condensed-phase systems in real time. In liquids and interfaces, we observe all the relevant regimes proper of non-Markovian open quantum system dynamics, such as electronic energy transfer, and screening. In addition, the ab-initio modeling of system-bath interactions brought us to observe and justify the holographic time-dependent electron density theorem. Contrary to interactions between molecular (finite) systems, when molecules interact with metal or semiconductor surfaces the electron dynamics is strongly non-Markovian with dramatic repercussions to the molecule’s response to external perturbations. Metals and semiconductors typically have large polarizabilities, and even in a regime of low coupling their effect on impinging molecular species is significant–line broadening, peak shift, and intensity borrowing are observed, characterized, and explained in terms of inter-subsystem dynamical interactions and a many-body decomposition of the system’s density-density response function in a way that transcends the canons of Fermi Golden Rule.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2021

Investigation of $\gamma$-irradiated polyvinylidene fluoride and its acute toxicity

Polyvinylidene fluoride (PVDF) is an industrial polymer with many applications. The EPR spectrum of polyvinylidene fluoride (PVDF) irradiated by γ-rays at low temperature (-196°C) provides evidence for the presence of radicals formed by loss of a fluorine atom from PVDF, and polyenyl radicals. Polyenyl radicals dominate in the polymer irradiated at room temperature. Peroxide radicals formed as a result of oxidation of primary free radicals by atmospheric oxygen are observed in the EPR spectrum of a polymer irradiated in air, or when air is introduced to a polymer irradiated in a vacuum. Carbonization and oxidation of γ-irradiated PVDF were observed in the XPS spectrum. A model for polymer degradation based on correlated G3(MP2) molecular orbital theory electronic structure calculations is described in this report. The presence of HF in the products and a low yield of products with the composition CxFyHz indicate the absence of noticeable thermal destruction of the main polymer chain upon heating to 220°C of both the initial and γ-irradiated PVDF. Although any HF released during the PVDF dehydrofluorination would be toxic, a study of the acute toxicity of γ-irradiated PVDF showed that the polymer does not exhibit toxicity after a single intragastric route of administration to mice.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Near-wake dynamics of a vertical-axis turbine

Cross-flow, or vertical-axis, turbines are a promising technology for capturing kinetic energy in wind or flowing water and their inherently unsteady fluid mechanics present unique opportunities for control optimization of individual rotors or arrays. To explore the potential for beneficial interactions between turbines in an array, as well as to characterize important cycle-to-cycle variations, coherent structures in the wake of a single two-bladed cross-flow turbine are examined using planar stereo particle image velocimetry in a water channel experiment. There are three main objectives in the present work. First, the mean wake structure of this high chord-to-radius ratio rotor is described, compared with previous studies, and a simple explanation for observed wake deflection is presented. Second, the unsteady flow is then analysed via the triple decomposition, with the periodic component extracted using a combination of traditional techniques and a novel implementation of the optimized dynamic mode decomposition. The latter method is shown to outperform conditional averaging and Fourier methods, as well as uncover frequencies suggesting a transition to bluff-body shedding in the far wake. Third, vorticity and finite-time Lyapunov exponents are then employed to further analyse the oscillatory wake component. Vortex streets on both sides of the wake are identified, and their formation mechanisms and effects on the mean flow are discussed. Strong axial (vertical) flow is observed in vortical structures shed on the retreating side of the rotor where the blades travel downstream. Time-resolved tracking of these vortices is performed, which demonstrates that vortex trajectories have significant rotation-to-rotation variation within one diameter downstream. This variability suggests it would be challenging to harness or avoid such structures at greater downstream distances.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Dynamic Barriers to Crystallization of Calcium Barium Carbonates

Metastable carbonates play important roles in geochemistry, biomineralization and serve as model systems for nonclassical theories of nucleation and growth. Balcite (Ca 0.5 Ba 0.5 CO 3 ) is a remarkable carbonate phase that is isostructural with a high-temperature modification of calcite (CaCO 3 ), yet can be synthesized at ambient conditions. Here, we investigate crystallization pathways in the Ba-Ca-CO 3 -H 2 O system, with a focus on the transformation of amorphous calcium barium carbonate (ACBC) to balcite and subsequent decomposition into the equilibrium calcite (CaCO3) and witherite (BaCO 3 ) phases. Density functional theory calculations show that balcite is an unstable solid solution (Ca 1-x BaxCO 3 , R 3¯ m) in the range 0.17 < x < 0.5, but is accessible through the amorphous ACBC precursor for x ≲ 0.5, and predict its decomposition into calcite and witherite. We confirm this pathway experimentally but found demixing to proceed slowly and remain incomplete even after 9 months. Nucleation kinetics of balcite from ACBC was assessed using a microfluidic assay, where increasing barium content led to a surprising increase in the balcite nucleation rate, despite decreasing thermodynamic driving force. We attribute crystallization rates that dramatically accelerate with time to changes in interfacial structure and composition during coarsening of the amorphous precipitate. By carefully quantifying the thermodynamic and kinetic contributions in the multistep crystallization of a metastable carbonate, we produce insights that allow us to better interpret the formation and persistence of metastable minerals in natural and synthetic environments.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Management of Risk and Uncertainty Through Optimized Co-Operation of Transmission Systems and Microgrids With Responsive Loads (Final Report)

The evolution of the power system to the reliable, efficient and sustainable system of the future will involve development of both demand- and supply-side technology and operations. Ambitious national and state-level goals around the decarbonization of electricity relies on the integration of very high levels of renewable resources, most of which are variable and intermittent. The use of demand response is an ideal approach to counterbalance the intermittency of renewable generation and brings the consumer into the spotlight. Though individual consumers are interconnected at the low-voltage distribution system, these resources are typically modeled as variables at the transmission network level. Demand-side participation cannot be leveraged effectively without explicitly including the distribution system dynamics in the optimization-based wholesale market operations. This project grew from a vision for co-optimized interaction of distribution systems, or microgrids, with the high-voltage transmission system. In this framework, microgrids encompass consumers, distributed renewables and storage. The energy management system of the lower voltage system (distribution or microgrid) can also sell (buy) excess (necessary) energy from the transmission system. Until recently, very little research had been conducted on the co-optimization of these two systems due to computational limitations. However, advances in computational capabilities, and the judicious use of decomposition methods and innovative approximation methods for high-dimension dynamic programming made this goal a viable objective for this project, leading to a fundamental shift in the ability to integrate and fully utilize demand-side resources. To this end, the modeling framework developed introduces a novel co-optimization framework, to include the operations of both the transmission and distribution systems (or microgrids) in operational decision making. This framework was used to analyze renewable and distributed generation along with responsive demand and to compare the capability of co-optimized systems to perform with higher levels of variable renewables. An ideal microgrid is defined as an electric entity capable of operating in both interconnected (with the high-voltage grid) and islanded mode. As such, the microgrid should incorporate generating units (traditional units and intermittent) and if needed, exchange power with the high-voltage grid. The interplay between the microgrid and high-voltage grid motivated the development of the co-optimization approach to ensure efficient performance of the interconnected network. Results show that the use of a bi-level optimization approach is an appropriate structure, capable of co-optimizing a transmission system with multiple distribution systems and microgrids. While increasing the number of connected systems provides increasing flexibility for renewables integration this can also the economic benefits to the low-voltage subsystems with each additional system connected. Comparison of a traditional single-level decision structure with the co-optimization approach illustrates a reduction in overall system cost under co-optimization, while specific cost allocations to transmission and distribution systems are changed.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2020

Differences in substrate use linked to divergent carbon flow during litter decomposition

ABSTRACT Discovering widespread microbial processes that create variation in soil carbon (C) cycling within ecosystems may improve soil C modeling. Toward this end, we screened 206 soil communities decomposing plant litter in a common garden microcosm environment and examined features linked to divergent patterns of C flow. C flow was measured as carbon dioxide (CO2) and dissolved organic carbon (DOC) from 44-days of litter decomposition. Two large groups of microbial communities representing ‘high’ and ‘low’ DOC phenotypes from original soil and 44-day microcosm samples were down-selected for fungal and bacterial profiling. Metatranscriptomes were also sequenced from a smaller subset of communities in each group. The two groups exhibited differences in average rate of CO2 production, demonstrating that the divergent patterns of C flow arose from innate functional constraints on C metabolism, not a time-dependent artefact. To infer functional constraints, we identified features – traits at the organism, pathway or gene level – linked to the high and low DOC phenotypes using RNA-Seq approaches and machine learning approaches. Substrate use differed across the high and low DOC phenotypes. Additional features suggested that divergent patterns of C flow may be driven in part by differences in organism interactions that affect DOC abundance directly or indirectly by controlling community structure.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Effect of transport mechanism on the coarsening of bicontinuous structures: A comparison between bulk and surface diffusion

Here, coarsening of bicontinuous microstructures is observed in a variety of systems, such as nanoporous metals and mixtures that have undergone spinodal decomposition. To better understand the morphological evolution of these structures during coarsening, we compare the morphologies resulting from two different coarsening mechanisms, surface and bulk diffusion. We perform phase-field simulations of coarsening via each mechanism in a two-phase mixture at nominal volume fractions of 50%-50% and 36%-64%, and the simulated structures are characterized in terms of topology (genus density), the interfacial shape distribution, structure factor, and autocorrelations of phase and mean curvature. We observe self-similar evolution of morphology and topology and agreement with the expected power laws for dynamic scaling, in which the characteristic length scale increases over time proportionally to t 1/4 for surface diffusion and t 1/3 for bulk diffusion. While we observe the expected difference in the coarsening kinetics, we find that differences in self-similar morphology due to coarsening mechanism are relatively small, although typically they are larger at 36% volume fraction than at 50% volume fraction. In particular, we find that bicontinuous structures coarsened via surface diffusion have lower scaled genus densities than structures coarsened via bulk diffusion. We also compare the self-similar morphologies to those in literature and to two model bicontinuous structures, namely, constant-mean-curvature surfaces based on the Schoen G minimal surface and random leveled-wave structures. The average scaled mean curvatures of these model structures agree reasonably with those of the coarsened structures at both 36% and 50%, but we find substantial disagreements in the scaled genus densities and the standard deviations of mean curvature.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Three-dimensional Hydrodynamic Simulations of Convective Nuclear Burning in Massive Stars Near Iron Core Collapse

Nonspherical structure in massive stars at the point of iron core collapse can have a qualitative impact on the properties of the ensuing core-collapse supernova explosions and the multimessenger signals they produce. Strong perturbations can aid successful explosions by strengthening turbulence in the postshock region. Here we report on a set of 4π 3D hydrodynamic simulations of O- and Si-shell burning in massive star models of varied initial masses using MESA and the FLASH simulation framework. We evolve four separate 3D models for roughly the final 10 minutes prior to and including iron core collapse. We consider initial 1D MESA models with masses of 14, 20, and 25 M ⊙ to survey a range of O/Si-shell density and compositional configurations. Further, we characterize the convective shells in our 3D models and compare them to the corresponding 1D models. In general, we find that the angle-average convective speeds in our 3D simulations near collapse are three to four times larger than the convective speeds predicted by MESA at the same epoch for our chosen mixing length parameter of α MLT = 1.5. In three of our simulations, we observe significant power in the spherical harmonic decomposition of the radial velocity field at harmonic indices of ℓ = 1–3 near collapse. Our results suggest that large-scale modes are common in massive stars near collapse and should be considered a key aspect of presupernova progenitor models.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Is the Trotterized UCCSD Ansatz Chemically Well-Defined?

The variational quantum eigensolver (VQE) has emerged as one of the most enticing near-term quantum algorithms that can be used to simulate many-body systems such as molecular electronic structures. Serving as an attractive ansatz in the VQE algorithm, unitary coupled cluster (UCC) theory has seen a renewed interest in recent literature. Yet, unlike the original classical UCC theory, implementation on a quantum computer requires a finite-order Suzuki-Trotter decomposition to separate the exponentials of the large sum of Pauli operators. While previous literature has recognized the non-uniqueness of different orderings of the operators in the Trotterized form of UCC methods, the question of whether or not different orderings matter at the chemical scale has not been addressed. In this letter, we explore the effect of operator ordering on the Trotterized UCCSD ansatz, as well as the much more compact k-UpCCGSD ansatz recently proposed by Lee et al. We observe a significant, system-dependent variation in the energies of Trotterizations with different operator orderings. The energy variations occur on a chemical scale, sometimes on the order of hundreds of kcal/mol. This letter establishes the need to define not only the operators present in the ansatz, but also the order in which they appear. This is necessary for adhering to the quantum chemical notion of a ``model chemistry'', in addition to the general importance of scientific reproducibility. As a final note, we suggest a useful strategy to select out of the combinatorial number of possibilities, a single well-defined and effective ordering of the operators.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Assessment and Non-Destructive Evaluation of the Influence of Residual Solvent on a Two-Part Epoxy-Based Adhesive Using Ultrasonics

Polymers are increasingly being used in higher demanding applications due to their ability to tailor the properties of structures while allowing for a weight and cost reduction. Solvents play an important role in the manufacture of polymeric structures since they allow for a reduction in the polymer’s viscosity or assist with the dispersion of fillers into the polymer matrix. However, the incorrect removal of the solvent affects both the physical and chemical properties of polymeric materials. The presence of residual solvent can also negatively affect the curing kinetics and the final quality of polymers. Destructive testing is mainly performed to characterize the properties of these materials. However, this type of testing involves using lab-type equipment that cannot be taken in-field to perform in situ testing and requires a specific sample preparation. Here, a method is presented to non-destructively evaluate the curing process and final viscoelastic properties of polymeric materials using ultrasonics. In this study, changes in longitudinal sound speed were detected during the curing of an aerospace epoxy adhesive as a result of variations in polymer chemistry. To simulate the presence of residual solvent, samples containing different weight percentages of isopropyl alcohol were manufactured and tested using ultrasonics. Thermogravimetric analysis was used to show changes in the decomposition of the adhesive due to the presence of IPA within the polymer structure. Adding 2, 4, and 6 wt.% of IPA decreased the adhesive’s lap shear strength by 40, 58, and 71%, respectively. Ultrasonics were used to show how the solvent influenced the curing process and the final sound speed of the adhesive. Young’s modulus and Poisson’s ratio were determined using both the longitudinal and shear sound speeds of the adhesive. Using ultrasonics has the potential to non-invasively characterize the quality of polymers in both an in-field and manufacturing settings, ensuring their reliability during use in demanding applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Decomposition and Algorithmic Approaches for Solving Large-Scale Process Family Design Problems

Our most recent work expands the water desalination case study from 76 variants to 10,897 variants using the equation-oriented model built in Pyomo as part of the PARETO project. Using the discretization formulation presented in Stinchfield (2024a), rather than solving for all 10,897 variants simultaneously, we decompose the formulation into subproblems containing subsets of variants from the process family. We solve the overall problem with Progressive Hedging (PH) deployed in parallel on a distributed HPC cluster using the open-source Python package mpi-sppy (Knueven et al., 2023). This approach allowed us to solve this process family design problem to ~1.5% relative optimality gap in about 5 hours; in comparison, Gurobi reached ~50% relative optimality gap in about 6 hours (Stinchfield et al., 2024b). However, this approach still requires discretization of the common unit module design ranges; additionally, PH acts as a heuristic for MILP’s with gap-closing capabilities. Ideally, we would not have to use ML surrogates or discretization to solve this problem, instead solving the process family design problem with the equation-oriented model directly to achieve the most accurate results. However, recall that we did not consider solving the MINLP directly due to complexity and size. In this work, we aim to decompose and solve this large-scale MINLP using a Structured Nonlinear Global Optimization algorithm presented by Cao and Zavala (2019).

Stinchfield, Georgia↗
DOE OSTI2025

Beyond real: alternative unitary cluster Jastrow models for molecular electronic structure calculations on near-term quantum computers

Near-term quantum devices require wavefunction ansätze that are expressive while also of shallow circuit depth in order to both accurately and efficiently simulate molecular electronic structure. While the unitary coupled cluster ansatz (e.g., UCCSD) has become a standard, the high gate count associated with the implementation of this limits its feasibility on noisy intermediate-scale quantum (NISQ) hardware. k -Fold unitary cluster Jastrow (uCJ) ansätze mitigate this challenge by providing O( kN 2 ) circuit scaling and favorable linear depth circuit implementation. Previous work has focused on the real orbitalrotation (Re-uCJ) variant of uCJ, which allows an exact (Trotter-free) implementation. Here we extend and generalize the k -fold uCJ framework by introducing two new variants, Im-uCJ and g-uCJ, which incorporate imaginary and fully complex orbital rotation operators, respectively. Similar to Re-uCJ, both of the new variants achieve quadratic gate-count scaling. Our results focus on the simplest k = 1 model, and show that the uCJ models frequently maintain energy errors within chemical accuracy (∼1 kcal mol −1 ). Both g-uCJ and Im-uCJ are more expressive in terms of capturing electron correlation and are also more accurate than the earlier Re-uCJ ansatz. We further show that Im-uCJ and g-uCJ circuits can also be implemented exactly, without any Trotter decomposition. Numerical tests using k = 1 on H 2 , H 3 + , Be 2 , C 2 H 4 , C 2 H 6 and C 6 H 6 in various basis sets confirm the practical feasibility of these shallow Jastrow-based ansätze for applications on near-term quantum hardware.

Tkachenko, Nikolay V. [University of California, B↗
DOE OSTI2025

Autonomous direct freeform fabrication strategy for multi-axis additive manufacturing

Multi-axis additive manufacturing (M-AM) enables precise material deposition along both planar and curved layers, eliminating the need for support structures through a continuous material deposition approach. In contrast to conventional 2-dimensional planar layers constrained to a fixed building orientation, the deposition on freeform layers demands the specification of guided curves to determine material deposition directions which are no longer to be fixed to a build direction. There are challenges that arise when fabricating components with multiple “build” directions, necessitating the decomposition of geometries and the specification of guided curves for the resulting volumes. Furthermore, multi-axis systems introduce heightened challenges due to an increased degree of freedom in motion. Consequently, the potential risks of collision between the deposited geometry and the motion platform become a notable concern. This research proposes a freeform layering algorithm to address the challenge of seamless transitions between planar and curved layers in the process planning of M-AM. The algorithm computes 3D “printable” layers by leveraging topological information derived from the geometry to be built and integrates collision-free manufacturability considerations into the computational process. These accumulated volumes serve as a “substrate” and support volumes for subsequent deposition, allowing later layers to be built without the need for additional support material. In conclusion, the proposed method successfully devises a freeform layering approach suitable for intricate models that demand substantial support, thus enabling the fabrication of diverse geometries in a manner previously unachievable.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Pre-conditioned BFGS-based uncertainty quantification in elastic full-waveform inversion

SUMMARY Full-waveform inversion has become an essential technique for mapping geophysical subsurface structures. However, proper uncertainty quantification is often lacking in current applications. In theory, uncertainty quantification is related to the inverse Hessian (or the posterior covariance matrix). Even for common geophysical inverse problems its calculation is beyond the computational and storage capacities of the largest high-performance computing systems. In this study, we amend the Broyden–Fletcher–Goldfarb–Shanno (BFGS) algorithm to perform uncertainty quantification for large-scale applications. For seismic inverse problems, the limited-memory BFGS (L-BFGS) method prevails as the most efficient quasi-Newton method. We aim to augment it further to obtain an approximate inverse Hessian for uncertainty quantification in FWI. To facilitate retrieval of the inverse Hessian, we combine BFGS (essentially a full-history L-BFGS) with randomized singular value decomposition to determine a low-rank approximation of the inverse Hessian. Setting the rank number equal to the number of iterations makes this solution efficient and memory-affordable even for large-scale problems. Furthermore, based on the Gauss–Newton method, we formulate different initial, diagonal Hessian matrices as pre-conditioners for the inverse scheme and compare their performances in elastic FWI applications. We highlight our approach with the elastic Marmousi benchmark model, demonstrating the applicability of pre-conditioned BFGS for large-scale FWI and uncertainty quantification.

58 GEOSCIENCES↗
DOE OSTI2022

Raman and infrared spectra of plutonium (IV) oxalate and its thermal degradation products

For over 80 years, plutonium dioxide has been routinely produced via thermal decomposition of hydrated plutonium(IV) oxalate. Despite the longstanding utility of this process, the chemical structures of starting materials and intermediates produced during this thermal conversion remain ill-defined. To help resolve this uncertainty, we measured high-resolution Raman and infrared spectra of Pu(C 2 O 4 ) 2 ·6H 2 O that was heated to 25, 100, 220, 250, 350, and 450 °C in air. Our measurements show that Pu(C 2 O 4 ) 2 ·6H 2 O has a rich vibrational spectrum with at least 15 Raman bands between 180 cm -1 and 1900 cm -1 and 9 infrared bands between 800 cm -1 and 4000 cm -1 . As Pu(C 2 O 4 ) 2 ·6H 2 O is heated, water is liberated, and the oxalate ligand decomposes to produce plutonium oxycarbide species. When heated to 350 °C or higher, vibrational spectra are consistent with PuO 2 with some residual carbon-containing species. Full vibrational spectra, powder X-ray diffraction, and scanning electron microscopy measurements of Pu(C 2 O 4 ) 2 ·6H 2 O and its thermal degradation products are presented herein along with approximate assignments for observed spectral bands. These data can be used to validate and potentially improve existing computational models that describe the chemical structure of compounds produced during thermal degradation of plutonium (IV) oxalate. Given the utility of plutonium (IV) oxalate in synthesizing plutonium dioxide, these results are expected to provide value in the fields of nuclear fuel processing, nuclear nonproliferation, and nuclear forensics.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2022

Rank-reduced coupled-cluster. III. Tensor hypercontraction of the doubles amplitudes

In this work, we develop a quartic-scaling implementation of coupled-cluster singles and doubles (CCSD) based on low-rank tensor hypercontraction (THC) factorizations of both the electron repulsion integrals (ERIs) and the doubles amplitudes. This extends our rank-reduced (RR) coupled-cluster method to incorporate higher-order tensor factorizations. The THC factorization of the doubles amplitudes accounts for most of the gain in computational efficiency as it is sufficient, in conjunction with a Cholesky decomposition of the ERIs, to reduce the computational complexity of most contributions to the CCSD amplitude equations. Further THC factorization of the ERIs reduces the complexity of certain terms arising from nested commutators between the doubles excitation operator and the two-electron operator. We implement this new algorithm using graphical processing units and demonstrate that it enables CCSD calculations for molecules with 250 atoms and 2500 basis functions using a single computer node. Furthermore, we show that the new method computes correlation energies with comparable accuracy to the underlying RR-CCSD method.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

The Surface Chemistry of Methanol on Cu 3 Pd(111): Effects of Metal Alloying and Reaction with Hydrogen

Synchrotron-based ambient-pressure X-ray photoelectron spectroscopy (AP-XPS) was used to study the adsorption and surface chemistry of methanol on a Cu 3 Pd(111) model surface. The composition and morphological properties of the pristine Cu 3 Pd(111) substrate were analyzed using a combination of low-energy electron diffraction (LEED), scanning tunneling microscopy (STM), and low-energy ion scattering (LEIS). The results of ion scattering showed segregation of Pd toward the surface, with a Pd/Cu ratio close to 0.5, a larger value than the ratio of 0.33 expected for a Cu 3 Pd bimetallic structure. The surface of the alloy exhibited a good LEED pattern with three-fold long-range periodicities. In STM, clusters of palladium with hexagonal arrays and Pd- Pd distances of 2.7-2.8 Å were detected. Bonding to copper perturbed the electronic properties of the atoms in the Pd clusters, shifting their 4d states toward higher binding energy with respect to the Fermi level. Further, the valence band spectrum of Cu 3 Pd(111) exhibited a line shape that was very different from those displayed by Cu(111) or Pd(111). At low pressures, the adsorption of methanol on Cu 3 Pd(111) at 300 K mainly produced CH 3 O, CO and CH x species. AP-XPS showed that most Pd atoms in the surface of the Cu 3 Pd(111) alloy interacted with the decomposition products of methanol. No significant changes were observed in the core levels of copper upon the adsorption and dissociation of methanol, suggesting that the molecule mainly interacted with Pd sites of the alloy. Reaction with hydrogen led to fast removal of CH x , C and PdC x species from Cu 3 Pd(111) at moderate (< 450 K) temperatures and prevented a CH x → C transformation. If the stability of adsorbed CH 3 O is used as a descriptor for the hydrogenation of CO 2 to methanol, Cu 3 Pd should be a much better catalyst than monometallic palladium. This may be a consequence of electronic and ensemble effects in the alloy that moderate the reactivity of Pd sites.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗