Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Classical mechanics”

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 847 records · Page 47
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

SIERRA Code Coupling Module: Arpeggio User Manual - Version 5.28

The SNL Sierra Mechanics code suite is designed to enable simulation of complex multiphysics scenarios. The code suite is composed of several specialized applications which can operate either in standalone mode or coupled with each other. Arpeggio is a supported utility that enables loose coupling of the various Sierra Mechanics applications by providing access to Framework services that facilitate the coupling. More importantly Arpeggio orchestrates the execution of applications that participate in the coupling. This document describes the various components of Arpeggio and their operability. The intent of the document is to provide a fast path for analysts interested in coupled applications via simple examples of its usage.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

A theoretical study of the upper bound of surface elevation variance in the Phillips initial stage during wind-wave generation

The resonance mechanism in the initial of wind-wave generation proposed by Phillips is a foundation of wind-wave generation theory, but a precise theoretical quantification of wave energy growth in this initial stage has not been obtained yet after more than six decades of research. In this study, we aim to address this knowledge gap by developing an analytical approach based on a novel complex analysis method to theoretically investigate the temporal evolution of the wave energy in the Phillips initial stage. We quantitatively derive and analyse the growth behaviour of the surface wave energy and obtain an analytical solution for its upper bound. Our result highlights the crucial effects of surface tension. Because the phase velocity of gravity–capillary waves has a minimal value at a critical wavenumber, gravity–capillary waves and gravity waves (which neglect surface tension) exhibit distinct resonance curve properties and wave energy growth behaviours. For gravity waves, the resonance curve extends indefinitely; for gravity–capillary waves, it either forms a finite-length curve or does not exist, depending on the wind speed. The leading-order term of the upper-bound solution of the energy of gravity waves increases linearly over time, while for gravity–capillary waves, the term increases linearly over time under strong wind conditions but remains finite under weak wind conditions. This theoretical study provides an analytical framework for the generation of wind-waves in the Phillips initial stage, which may inspire further theoretical, numerical and experimental research.

58 GEOSCIENCES↗
DOE OSTI2021

Toward witnessing molecular exciton entanglement from spectroscopy

Entanglement is a defining feature of quantum mechanics that can be a resource in engineered and natural systems, but measuring entanglement in experiment remains elusive especially for large chemical systems. Most practical approaches require determining and measuring a suitable entanglement witness which provides some level of information about the entanglement structure of the probed state. A fundamental quantity of quantum metrology is the quantum Fisher information (QFI), which is a rigorous witness of multipartite entanglement that can be evaluated from linear response functions for certain states. Here, in this paper, we explore measuring the QFI of molecular exciton states of the first-excitation subspace from spectroscopy. In particular, we utilize the fact that the linear response of a pure state subject to a weak electric field over all possible driving frequencies encodes the variance of the collective dipole moment in the probed state, which is a valid measure for QFI. The systems that are investigated include the molecular dimer, N-site linear aggregate with nearest-neighbor coupling, and N-site circular aggregate, all modeled as a collection of interacting qubits. Our theoretical analysis shows that the variance of the collective dipole moment in the brightest dipole-allowed eigenstate is the maximum QFI. The optical response of a thermally equilibrated state in the first-excitation subspace is also a valid QFI. Theoretical predictions of the measured QFI for realistic linear dye aggregates as a function of temperature and energetic disorder due to static variations of the host matrix show that two- to three-partite entanglement is realizable. This paper lays some groundwork and inspires measurement of multipartite entanglement of molecular excitons with ultrafast pump-probe experiments.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Kinetics of TATB and LX-17 Decomposition by Thermal Analysis

The kinetics of TATB decomposition have been investigated by thermal analysis (thermogravimetric analysis, TGA, and differential scanning calorimetry, DSC) at both LLNL and the University of Alabama Birmingham. This report summarizes progress in understanding from those results within the context of identifying the global reaction mechanism, i.e., sequential, autocatalytic, competitive, or some mixture thereof. A particularly important issue is how pressure affects both the reaction rate and enthalpy. The effect of pressure helps identify dominant mechanisms and construct and constrain a global mechanism that will work over a wide range of conditions. This memorandum is simply a snapshot of current progress and not a definitive interpretation. I have written it to organize work of others and my interpretations in a systematic manner with the intention of stimulating further advances in understanding. SDT data comes from Jen Montgomery at LLNL, and high-pressure DSC comes from Victoria Stanford and Sergey Vyazovkin at University of Alabama at Birmingham. Major observations and conclusions, some new and some not, are (1) outgassing below 150 °C prior to decomposition releases adsorbed water vapor, (2) outgassing between 150 and 250 °C may come primarily from intercalated water (3) initial thermal decomposition above 250 °C or so involves dehydration of TATB to create furazans, (4) quickly heating to higher temperatures in open pans with flowing gas leads primarily to TATB evaporation between 270 and 390 °C, (5) using a pinhole crucible at atmospheric pressure largely eliminates TATB sublimation, (6) at typical thermal analysis heat rates, the main stage of decomposition is described by a single mass loss process, and decomposition of residual material is spread over a wide range of thermal severity, (7) at typical thermal analysis heat rates and constant temperatures above 300 °C, heat release occurs in two overlapping peaks, with the first one a little faster than mass loss and the second on a little slower than mass loss, (8) total heat release increases significantly with heating rate and substantially with increased pressure, and (9) isothermal experiments at 340 °C indicate that the first exotherm increases by more with LX-17 than TATB, which suggests a possible chemical role of Kel-F decomposition.

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

Symposium: Quantum Theory, Quantum Materials, Quantum Computing (Final Report, DE-SC0022551)

Quantum mechanics is fundamental to understanding, predicting, and controlling many physical processes in a vast range of scientific and technological areas. In molecules and materials, the quantum mechanics of Coulombic many-body systems determines far more than the ground state and bulk properties. Quantum mechanics governs magnetism, electron and spin transport, chemical reactions and catalysis, superconductivity, photosynthesis, etc. Each plays a vital role in modern society. As digital technologies push both temporal and length-scale limits, fundamental quantum principles become the current focus for a new class of technologies: quantum computing, quantum sensing, quantum encryption, etc. This oncoming quantum information sciences wave is manifest recently in frontier foci of all six Department of Energy Office of Science programs. The Sanibel Symposium offers a proven, influential platform for facilitating and advancing cross-fertilization among theoretical and computational chemists and physicists, computer scientists, and materials engineers. The Symposia have a long, successful record of driving progress on theory and computation of materials, nanostructures, and their molecular constituents. The Symposia thus are well-suited for advancing the quantum information science thrust in the Office of Basic Energy Sciences (BES). With that QIS thrust as the focus, we propose five sessions of the 61st Sanibel Symposium on the theme Quantum Theory, Quantum Materials, Quantum Computing. From a longer perspective, this thematic program for 2022 is the second year of a five year plan of program segments addressed at challenges in quantum information sciences and quantum materials. The rest of the Symposium program will be of considerable relevance, with, for example, sessions on quantum molecular and material dynamics, and new directions in coupled cluster theory. New ideas for methods and applications can be expected to emerge from the stimulus of the formal presentations and informal discussions. The proposed award is targeted at students, post-docs, and junior scientists. Conscious effort, characteristic of the Sanibel Symposia, to stimulate the participation of women and under-represented minorities will continue.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Pseudogap and Fermi arc induced by Fermi surface nesting in a centrosymmetric skyrmion magnet

Skyrmions in noncentrosymmetric materials are believed to occur due to the Dzyaloshinskii-Moriya interaction. By contrast, the skyrmion formation mechanism in centrosymmetric materials remains elusive. Here, in this work, we reveal the intrinsic electronic structure of the centrosymmetric GdRu 2 Si 2 by selectively measuring magnetic domains using angle-resolved photoemission spectroscopy (ARPES). We found robust Fermi surface (FS) nesting, consistent with the magnetic modulation q -vector detected by the previous resonant x-ray scattering measurements. The pseudogap opens at the nested FS portions, which vary for different magnetic domains. The anomalous pseudogap disconnects the FS to generate Fermi arcs with twofold symmetry. These results indicate that the Ruderman-Kittel-Kasuya-Yosida (RKKY) interaction plays a decisive role in generating the screw spin modulation responsible for the skyrmion formation in GdRu 2 Si 2 . Furthermore, we demonstrate the flexible nature of magnetism in GdRu 2 Si 2 by manipulating magnetic domains with magnetic field and temperature cyclings, providing potential future applications for data storage and processing devices.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Injection mechanisms in a III -nitride light-emitting diode as seen by self-emissive electron microscopy

Here, we report on the investigation of an electrically biased high efficiency green III-nitride light-emitting diode (LED) by electron emission microscopy (EEM) using a low-energy electron microscope (LEEM). The surface of the LED was activated to negative electron affinity via deposition of a submonolayer of Cs. With the illumination column of the LEEM turned off, upon electrical injection of the LED, we directly image the hot electrons generated by eeh Auger-Meitner nonradiative processes that diffuse through the top p-Ga N layer and emit out the surface of the biased LED. By determining the source of emitted electrons using complementary electron emission spectroscopy measurements, EEM allows us to effectively map the carrier density within the LED. Using EEM, we observed nonelectron emitting regions with a density of approximately 3 × 10 8 cm -2 , identified as V-shaped defects (V-defects). This is confirmed through the corresponding dark spots of panchromatic cathodoluminescence measurements of the same sample and by plan-view transmission electron microscopy. The absence of electron emission at the sidewall of the V-defects can be attributed to several factors, including reduced carrier density in the sidewall quantum wells due to carriers traveling fast through the semipolar sidewalls before being injected into the planar quantum wells, the reduced population of hot electrons surviving diffusion through the thicker p-GaN filling in the V-defect before emission onto vacuum, and a smaller Auger-Meitner coefficient for the low In content semipolar sidewall quantum wells. The stronger electron emission observed at the ridges of most V-defects compared to the planar quantum well regions indicates larger local injected carrier densities, confirming that V-defect sidewalls allow for strong lateral carrier injection when compared to the weaker vertical injection away from the V-defect as evidenced by the weaker electron emission intensity away from the V-defects.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Towards satellite tests combining general relativity and quantum mechanics through quantum optical interferometry: progress on the deep space quantum link

The Deep Space Quantum Link (DSQL) is a space-mission concept that aims to explore the interplay between general relativity and quantum mechanics using quantum optical interferometry. This mission concept was formally presented to the United States National Academy of Science Decadal Survey as a research campaign for Fundamental Physics in 2022. Since then, advances have been made in the space-based quantum optical technologies required to conduct a DSQL-type mission. In addition, other research efforts have defined alternative measurement concepts to explore the same scientific questions motivating the DSQL mission. This paper serves as an update to the community on the status of the DSQL mission concept and related research and technology development efforts.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Enabling Large-Scale Condensed-Phase Hybrid Density Functional Theory Based Ab Initio Molecular Dynamics. 1. Theory, Algorithm, and Performance

By including a fraction of exact exchange (EXX), hybrid functionals reduce the self-interaction error in semilocal density functional theory (DFT) and thereby furnish a more accurate and reliable description of the underlying electronic structure in systems throughout biology, chemistry, physics, and materials science. However, the high computational cost associated with the evaluation of all required EXX quantities has limited the applicability of hybrid DFT in the treatment of large molecules and complex condensed-phase materials. To overcome this limitation, we describe a linear-scaling approach that utilizes a local representation of the occupied orbitals (e.g., maximally localized Wannier functions (MLWFs)) to exploit the sparsity in the real-space evaluation of the quantum mechanical exchange interaction in finite-gap systems. In this work, we present a detailed description of the theoretical and algorithmic advances required to perform MLWF-based ab initio molecular dynamics (AIMD) simulations of large-scale condensed-phase systems of interest at the hybrid DFT level. We focus our theoretical discussion on the integration of this approach into the framework of Car–Parrinello AIMD, and highlight the central role played by the MLWF-product potential (i.e., the solution of Poisson’s equation for each corresponding MLWF-product density) in the evaluation of the EXX energy and wave function forces. We then provide a comprehensive description of the exx algorithm implemented in the open-source Quantum ESPRESSO program, which employs a hybrid MPI/OpenMP parallelization scheme to efficiently utilize the high-performance computing (HPC) resources available on current- and next-generation supercomputer architectures. Furthermore, this is followed by a critical assessment of the accuracy and parallel performance (e.g., strong and weak scaling) of this approach when AIMD simulations of liquid water are performed in the canonical (NVT) ensemble. With access to HPC resources, we demonstrate that exx enables hybrid DFT-based AIMD simulations of condensed-phase systems containing 500–1000 atoms (e.g., (H₂O)₂₅₆) with a wall time cost that is comparable to that of semilocal DFT. In doing so, exx takes us one step closer to routinely performing AIMD simulations of complex and large-scale condensed-phase systems for sufficiently long time scales at the hybrid DFT level of theory.

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

Impact of Laboratory-Accelerated Aging Methods to Study Alkali–Silica Reaction and Reinforcement Corrosion on the Properties of Concrete

This study focuses on two separate investigations of the main aging mechanisms: alkali–silica reactivity (ASR) and the corrosion of reinforcing steel (rebar) concrete, both of which may result in a premature failure to meet the serviceability or strength requirements of a concrete structure. However, these processes occur very slowly, spanning decades. The impact of direct chemical additives to fresh concrete to accelerate ASR and the corrosion of reinforcing steel on the fresh and hardened properties of the ensuing material are investigated to inform the potential use of chemicals in large-scale studies. The deterioration of reinforced concrete (RC) is determined by means of expansion, cracking, bulk diffusivity and surface resistivity measurements, and compressive, split tensile and flexural strength tests. The results indicate that the addition of sodium hydroxide and calcium chloride can effectively accelerate the crack formation and propagation in concrete due to ASR and the corrosion of rebar, respectively. The ASR-induced cracks maintained a constant crack width from 0.05 mm to 0.1 mm over the measurement period regardless of the intensity of aging acceleration. Adding 4% chloride by weight of cement for accelerating rebar corrosion resulted in an average crack that was 82% larger than in the case of ASR accelerated with the addition of sodium hydroxide. The addition of alkali resulted in an increase in early-age (7-day) strength. At a total alkali loading of 2.98 kg/m3, 3.84 kg/m3 and 5.57 kg/m3, the 28-day compressive strength of concrete decreased by 3%, 10% and 24%, respectively. Similarly, a higher early-age strength and a lower later-age strength was observed for the concrete in the presence of corrosive calcium chloride. The results from this research are expected to inform future studies on the long-term performance of RC structures under accelerated ASR and corrosion.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Single-Step Conversion of Ethanol to n-Butene over Ag-ZrO2/SiO2 Catalysts

Ethanol is a promising platform molecule for production of a variety of fuels and chemicals. Of particular interest is producing middle distillate fuels (i.e., jet and diesel blendstock) from renewable ethanol feedstock. State-of-the-art alcohol-to-jet technology requires multiple process steps based on catalytic dehydration of ethanol to form ethylene, followed by sometimes a multi-step oligomerization, and then hydrotreatment and distillation. Here we report on a new catalytic route in which ethanol is directly converted to n-butene (1- and 2-butene mixtures) over Ag-ZrO2/SBA-16, thus offering the potential for a reduction in the number of required processing steps versus conventional alcohol-to-jet technology. This catalyst system provides the balanced metal and Lewis acid sites required to selectively facilitate a cascading sequence of reactions that includes dehydrogenation, aldol condensation, Meerwein–Ponndorf–Verley reduction, dehydration, and hydrogenation. High conversion and selectivity toward either n-butene or 1,3-butadiene is achieved by tuning the hydrogen feed partial pressure and other process/catalyst parameters. With sufficient hydrogen partial pressure 1,3-butadiene is completely and selectively hydrogenated to form n-butene. The reaction mechanism was elucidated through operando-nuclear magnetic resonance investigations coupled with reactivity measurements. Combined experimental-computational investigation reveals how changes in silver and zirconium composition and the silver oxidation state affects reactivity under controlled hydrogen partial pressures and after prolonged run times. Finally, catalyst effectiveness also was demonstrated when using wet ethanol feed, thus highlighting process flexibility in terms of feedstock purity requirements. This work was financially supported by the U.S. Department of Energy (DOE), Office of Energy Efficiency and Renewable Energy, Bioenergy Technologies Office, and was performed at the Pacific Northwest National Laboratory (PNNL) under Contract No. DE-AC05-76RL01830 and the National Renewable Energy Laboratory under Contract No. DE-AC36- 08GO28308. Part of the work conducted by S. A. Akhade was performed under the auspices of the U.S. DOE at Lawrence Livermore National Laboratory under Contract No. DE-AC52-07NA27344. This work was partly supported through the PNNL-WSU Distinguished Graduate Research Program for ADW. NMR and XPS experiments were performed using EMSL (grid.436923.9), a DOE Office of Science User Facility sponsored by the Office of Biological and Environmental Research

Dagle, Vanessa↗
DOE OSTI2021

Broadband 2DES detection of vibrational coherence in the S x state of canthaxanthin

The nonadiabatic mechanism that mediates nonradiative decay of the bright S 2 state to the dark S 1 state of carotenoids involves population of a bridging intermediate state, S x , in several examples. The nature of S x remains to be determined definitively, but it has been recently suggested that S x corresponds to conformationally distorted molecules evolving along out-of-plane coordinates of the isoprenoid backbone near a low barrier between planar and distorted conformations on the S 2 potential surface. Here, the electronic and vibrational dynamics accompanying the formation of S x in toluene solutions of the ketocarotenoid canthaxanthin (CAN) are characterized with broadband two-dimensional electronic spectroscopy (2DES) with 7.8 fs excitation pulses and detection of the linear polarization components of the third-order nonlinear optical signal. A stimulated-emission cross peak in the 2DES spectrum accompanies the formation of S x in <20 fs following excitation of the main absorption band. S x is prepared instantaneously, however, with excitation of hot-band transitions associated with distorted conformations of CAN’s isoprenoid backbone in the low frequency onset of the main absorption band. Vibrational coherence oscillation maps and modulated anisotropy transients show that S x undergoes displacements from the Franck–Condon S 2 state along out-of-plane coordinates as it passes to the S 1 state. The results are consistent with the conclusion that CAN’s carbonyl-substituted β-ionone rings impart an intramolecular charge-transfer character that frictionally slows the passage from S x to S 1 compared to carotenoids lacking carbonyl substitution. Despite the longer lifetime, the S 1 state of CAN is formed with retention of vibrational coherence after passing through a conical intersection seam with the S x state.

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

Projector-Free Capped-Fragment Scheme within Density Functional Embedding Theory for Covalent and Ionic Compounds

Quantum-mechanics-(QM)-based simulations now routinely aid in understanding and even discovering new chemistries involving molecules and materials exhibiting desired functionalities. Ab initio correlated wavefunction (CW) theories systematically improve QM methods, with many exhibiting high accuracy. However, execution of CW methods requires expensive computations that typically scale poorly with system size. Divide-and-conquer approaches partition large systems into smaller fragments; a lower level of theory treats fragment interactions while a preferred higher level of theory describes the important fragment. These methods offer ways to incorporate CWs into chemical simulations of large systems, e.g., biomolecules, surfaces, large inorganic clusters, bulk crystals, etc. Here we propose a partitioning protocol that utilizes capping atoms to saturate severed covalent bonds at fragment interfaces and density functional embedding theory (DFET) to describe fragment interactions. The capping groups in each fragment provide an ad-hoc potential that approximates the effects of the environment. An embedding potential optimized via DFET then serves as an augmentation of the capping group to simulate the effects of the environment. We concurrently use an auxiliary fragment (a separate system comprised of only the combined capping groups) to account for, and thereby correct, the electron density contributions of all the capping groups added to all of the fragments. This method depends only on the capped-subsystem and auxiliary-fragment electron densities, forgoing, as with the original DFET developed for metallic systems, orbital-based projector approaches that determine a nonlocal action of the embedding potential onto the fragment electrons. By using an auxiliary fragment, the method maintains a purely electron-density-dependent embedding potential, substantially lessening the cost and leading to simpler implementation. In this work, we demonstrate the utility of our capped-DFET and ensuing capped embedded CW method in two contrasting systems, namely, an organic molecule and an ionic metal oxide cluster.

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

Losses of CO and CO 2 upon collision-activated dissociation of substituted 2-methoxyphenoxides after methyl radical loss

Some deprotonated, substituted 2-methoxyphenols fragment by methyl radical loss followed by competing losses of CO and CO 2 upon collisional activation in linear quadrupole ion trap mass spectrometers. These reactions were examined experimentally and computationally in order to determine their mechanisms. For deprotonated vanillin, the CO loss was found to involve ring contraction, with the highest free energy barrier of 58.0 kcal mol -1 (for the syn rotamer; here, syn refers to the relationship between the aldehyde oxygen atom and the oxygen atom of the methoxy group) for the entire process. The atoms lost in this fragmentation are the oxygen atom that was bound to the first eliminated methyl group and the aromatic carbon atom bound to this oxygen. Examination of carbon-13 labeled vanillin supports these assignments. Examination of several model compounds revealed that this reaction requires the presence of an electron-withdrawing substituent in the para-position relative to the phenol moiety. In contrast, the CO 2 loss from deprotonated vanillin occurs via ring opening followed by re-cyclization and then ring contraction, leading to the loss of CO 2 in a process wherein the highest free energy barrier is 84.5 kcal mol -1 (for the syn isomer). The atoms lost in this fragmentation are the carbon and oxygen atoms from the phenoxide group and the oxygen atom that was bound to the first eliminated methyl group, which is supported by examination of carbon-13 labeled vanillin. Despite the higher total free energy requirement, the CO 2 loss is competitive with CO loss, possibly due to favorable entropy and low activation energy (free energy barrier of 48.9 kcal mol -1 ) for the first reaction step (the analogous value for CO loss is 58.0 kcal mol -1 ). The extent of CO 2 loss is strongly affected by substituents – it either is not observed or is very slow for the other compounds studied here. For example, it is substantially less favorable than CO loss (highest free energy barrier 60.9 kcal mol -1 ) for deprotonated acetovanillone for which the first reaction step for CO 2 loss has a substantially greater free energy barrier (60.6 kcal mol-1) than for vanillin.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Illuminating the Photophysics of Oxygen Atom Transfer in Rare Earth Metal-Organic Complexes

Industrial separations of rare earth metals are reliant on selectivity based on ionic radii, which are time, energy, and resource consuming.1 The unique and discontinuous manifold of energy states associated with the f-electron configurations of the rare earth ions2 offer an opportunity to overcome some of these limitations by incorporating light into the separation process. Here, we use transient absorption spectroscopy to understand the differences in mechanism of light-initiated oxygen atom transfer (OAT) in yttrium (Y) and dysprosium (Dy) rare earth complexes. These metal-organic complexes are functionalized with hexfluoroacetyl-acetonate (hfac) and 4-methylmorpholine-N-oxide (NMMO). After light absorption in the presence of triphenyl phosphine (TPP), the NMMO in the Y-NMMO complex is replaced by the TPPO OAT reaction product. Because the Dy f-electron manifold is accessible to the ligand excited states, while there are no accessible states associated with f-electrons for Y, there are significant differences in reactivity. Transient absorption spectroscopy (TAS) was used to monitor changes after 340 nm photoexcitation of Dy/Y-NMMO or the product complex Dy/Y-TPPO. TAS is a pump-probe technique where a white light probe pulse is delayed relative to a visible pump pulse, which photoexcites a sample. Differences between the ground and excited states are monitored over time. With TAS, we access time scales from 100 fs-400 us to determine the photophysics of the complexes of interest. Photoexcitation of Y-NMMO (Fig. 1A) creates a broad photoinduced absorption (PIA) that evolves into a sharp positive feature over 4.2 ps centered at 450 nm. A secondary feature centered at 550 nm grows in with a 200 ps time constant and remains through the 5 ns time window of our experiment. Control experiments for Y-TPPO reveal that the feature at 450 nm is due to internal dynamics of the hfac ligand, while the 550 nm shoulder is due to the presence of the NMMO. Comparing kinetics at 550 nm for Y-NMMO and Dy-NMMO (Fig. 1B) shows that while Y exhibits a distinct growth in the NMMO feature, the Dy decays almost completely within 5 ns. Control experiments with Dy-TPPO prove that the ligand excited states centered on the hfac ligand are deactivated via energy transfer to the Dy manifold within 200 ps. This deactivation leads to significant differences in OAT reactivity between Y and Dy.

CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS↗
DOE OSTI2020

Benchmarking an Embedded Adaptive Sampling Configuration Interaction Method for Surface Reactions: H 2 Desorption from and CH 4 Dissociation on Cu(111)

Embedded (emb-) correlated wavefunction (CW) theory enables accurate assessments of both ground- and excited-state reaction mechanisms involved in heterogeneous catalysis. Embedded multireference second-order perturbation theory (emb-MRPT2) based on reference wavefunctions generated via embedded complete active space self-consistent field (emb-CASSCF) theory is currently state-of-the-art. However, the factorial scaling of CASSCF limits the size of active space and the complexity of systems that can be studied. In this work, we assess the efficacy of an alternative CW method, adaptive sampling configuration interaction (ASCI)–which enables large active spaces to be used–for studying surface reactions. We couple ASCI with density functional embedding theory (DFET) and benchmark its performance for two reactions: H 2 desorption from and CH 4 dissociation on the Cu(111) surface. Unlike embedded complete active space second-order perturbation theory (emb-CASPT2) that accurately reproduces a measured H 2 desorption barrier, embedded ASCI, using a very large active space (though one that still comprises a small portion of the full set of orbitals) fails to do so. Adding an extra correlation term from embedded Møller–Plesset second-order perturbation theory (emb-MP2) improves the desorption barrier and endothermicity predictions. Thus, the inaccuracy of embedded ASCI comes from the missing dynamic correlation from the many other electrons and orbitals not included in the active space. For CH 4 dissociation, again embedded ASCI overestimates the dissociation barrier compared to emb-CASPT2 predictions. Adding dynamic correlation from emb-MP2 helps correct the barrier. However, this composite approach suffers from double counting of correlation within embedded ASCI followed by emb-MP2 calculations. We therefore conclude that the state-of-the-art emb-MRPT2 based on reference wavefunctions generated via emb-CASSCF remains the method of choice for studying surface reactions. emb-ASCI is useful when large active spaces beyond the limit of emb-CASSCF are essential, such as to study complex surface reactions with significant multiconfigurational character (static correlation) but weak dynamic correlation.

08 HYDROGEN↗
DOE OSTI2020

ChIMES: A Machine-Learned Interatomic Model Targeting Improved Description of Condensed Phase Chemistry in Energetic Materials

In this report we detail completion of a Physics and Engineering Model Level Two Milestone targeting improved reactive interatomic potentials (IAPs) for energetic materials (EM) through machine learning. The specific goals of this milestone were to develop, validate, and document a new reactive molecular dynamics method for EM, based on machine learning by (1) generating databases of first-principles-derived forces, stresses, and energies for HN3 and 3,4-bis(3-nitrofurazan- 4-yl)furoxan (DNTF) (2) generate atomistic force fields from these databases via ML, and (3) benchmark model performance against first principles calculations. These goals were achieved by (1) further developing a machine learned reactive IAP and generation approach (i.e. the Chebyshev Interaction Model for Efficient Simulation or “ChIMES”), for which resulting IAPs can approach the predictive power of quantum-mechanical approaches at a fraction of the computational expense, and (2) applying the ChIMES framework to develop models for HN3 and DNTF. We find that for simple energetic materials like HN3, high accuracy ChIMES models can be obtained through application of a fitting approach that does not use active machine learning. We demonstrate the suitability of ChIMES models for simulations involving EM by using the HN3 model in multiscale shock technique simulations to predict the HN3 Chapman-Jouguet detonation state and investigate chemical evolution out to 1 ns following shock compression. This model is then used in larger direct shock (DS) simulations for a preliminary investigation of how bubbles (i.e. voids) influence material response under shock compression. We find that more complex EM (i.e. DNTF) necessitate a more sophisticated fitting approach, and develop a new active learning method and python tool to meet this challenge. We demonstrate that this fitting approach yields ChIMES models that out-perform commonly used standard reactive IAPs as well as semi-empirical quantum methods, and discuss the systematic improvability of these actively learned ChIMES models. We also describe challenges related to model development for EM such as DNTF, for which few experimental or previous simulation data are available (e.g. which could otherwise inform generation of training data). To overcome this issue, we establish a semi-empirical quantum ChIMES capability which can be used to efficiently map out relevant thermodynamic and configurational space, and generate ChIMES-IAP training data in a multiscale manner. We also show that these semi-empirical quantum ChIMES models can be used to generate predictions for the shock Hugoniot (the Hugoniot is the locus of thermodynamic states found in a shocked material) equation of state, investigate related thermochemistry, and explore carbon condensation following shock compression. This work represents a substantial advance in our atomistic modeling capability for EM that will provide much needed information on the chemistry of detonation for continued development of continuum models based on the Cheetah thermochemical code.

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

Metal-to-Ligand Charge-Transfer Spectrum of a Ru-Bipyridine-Sensitized TiO 2 Cluster from Embedded Multiconfigurational Excited-State Theory

Understanding optical properties of the dye molecule in dye-sensitized solar cells (DSSCs) from first-principles quantum mechanics can contribute to improving the efficiency of such devices. While density functional theory (DFT) and time-dependent DFT have been pivotal in simulating optoelectronic properties of photoanodes used in DSSCs at the atomic scale, questions remain regarding DFT’s adequacy and accuracy to furnish critical information needed to understand the various excited-state processes involved. In this work, we simulate the absorption spectra of a dye-sensitized solar cell analogue, comprised of a Ru-bipyridine (Ru-bpy) dye molecule and a small TiO2 cluster via DFT and via an accurate embedded correlated wavefunction (CW) theory. We generated CW spectra for the adsorbed Ru-bpy dye via a recently introduced capped density functional embedding theory or capped-DFET (to generate the embedding potential that accounts for the interaction of the molecule and the TiO 2 cluster). We then combined capped-DFET with the accurate but expensive multiconfigurational complete active space second-order perturbation theory (CASPT2)—embedded CASPT2. Because the CW theory is conducted on only a portion of the total system in the presence of an embedding potential that describes that portion’s interaction with its environment, we efficiently obtain CW-quality predictions that reflect local properties of the entire system. Specifically, for example, with capped-DFET and embedded CW theory, we can simulate accurately a plethora of metal-to-ligand charge-transfer excited properties at a manageable computational cost. Here, we predict detailed electronic spectra within the visible region, featuring the lowest three singlet and triplet excited states, along with predictions of the singlets’ lifetimes. We illustrated these results using a Jablonski diagram that show the relative energy position of the singlet and longer-lived triplet excited states and analyzed and proposed relaxation paths for the excited state corresponding to the most intense but short-lived absorption (interconversion, intersystem crossing, fluorescence, and phosphorescence) that may lead to longer-lived excited states necessary for efficient charge separation required to generate current in solar cells.

14 SOLAR ENERGY↗