Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “quantum mechanical calculations”

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 595 records · Page 33

Enhanced spectral purity of WSe 2 quantum emitters via conformal organic adlayers

Quantum emitters in solid-state materials are typically embedded in the bulk of their hosts, making their electronic transitions inaccessible to surface modification. In contrast, two-dimensional materials, with their all-surface nature, offer a platform for tuning quantum emitters via chemical functionalization. Because of its semiconducting properties that enable electrical addressability, monolayer WSe 2 is a promising candidate for quantum emission, although the complex interplay between point defects and the localized strain needed to activate quantum emission leads to poor spectral purity. Here, we demonstrate that functionalizing monolayer WSe 2 with conformal adlayers of 3,4,9,10-perylenetetracarboxylic dianhydride (PTCDA) improves quantum emission spectral purity. Optical spectroscopy reveals that PTCDA functionalization lowers defect activation energies by 10 meV and induces a 30 nm redshift in quantum emission wavelength, while preserving the bright and dark exciton energies of monolayer WSe 2 . First-principles calculations corroborate these findings, thus providing molecular-level insight into the underlying mechanism of enhanced spectral purity.

Ananth, Riddhi [Northwestern Univ., Evanston, IL (↗

Fission Product Evaluation (Quarterly Report)

The LLNL team has developed a theoretical framework to compute the fission fragment charge and mass distributions of odd-mass nuclei, which can be formed in neutron-induced reactions on even-even nuclei. For example, n+ 238 U leads to the compound nucleus 239 U which has 147 neutrons. Odd-mass and odd-odd nuclei are spin-polarized: the spin of the nucleus is not equal to 0 as in even-even nuclei (a consequence of pairing correlations) but is in fact set by a single particle. In practice, that single particle can occupy different quantum-mechanical configurations, each with different spin. This affects the deformation properties of the nucleus and, as a result, the distribution of fission fragments. These distributions require two inputs for the calculation: a potential energy surface (= the total energy of the nucleus as a function of the deformation of the nuclear shape, often parametrized with 2 degrees of freedom) and the collective inertia tensor associated with each point of the PES. In an odd-mass or odd-odd nucleus, each quantity depends on the configuration of the single particle. We computed of the order of 20 to 30 different such configurations to extract the PES for spins K=1/2, 3/2, 5/2, 7/2, and 9/2 of 239 U. We then computed the primary yields for each such PES, as shown in Fig. 1 below. These yields were weighed with the spin distribution associated with the reaction 238 U+n and the resulting distribution used as input to FREYA to compute the independent yields shown in Fig. 2. These yields are then used as inputs to the FIER code to extract the cumulative yields and compare with LLNL-LANL measurements at TUNL.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Impact of solvation on the GW quasiparticle spectra of molecules

First-principles calculations for electrochemistry require accurate treatment of both electronic structure and solvation. The perturbative GW approximation starting from density functional theory (DFT) calculations accurately models materials systems with varying dimensionality. Continuum solvation models enable efficient treatment of solvation effects in DFT calculations, but their applications with beyond-DFT electronic structure methods such as GW have been limited. Here, we introduce the frequency-dependent liquid polarizability from a nonlocal continuum solvation model in the screened Coulomb interaction of full-frequency GW calculations with a solvated DFT starting point. In this study, we show that the liquid screening contributions substantially reduce the HOMO–LUMO gap of molecules by 3–5 eV, while solvent effects on the DFT starting point negligibly impact the GW gap. The resulting framework facilitates the simultaneous electronic and solvation accuracy needed for first-principles electrochemistry.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Vibrational Entropy of Crystalline Solids from Covariance of Atomic Displacements

The vibrational entropy of a solid at finite temperature is investigated from the perspective of information theory. Ab initio molecular dynamics (AIMD) simulations generate ensembles of atomic configurations at finite temperature from which we obtain the N-body distribution of atomic displacements, ρN. We calculate the information-theoretic entropy from the expectation value of lnρN. At a first level of approximation, treating individual atomic displacements independently, our method may be applied using Debye–Waller B-factors, allowing diffraction experiments to obtain an upper bound on the thermodynamic entropy. At the next level of approximation we correct the overestimation through inclusion of displacement covariances. We apply this approach to elemental body-centered cubic sodium and face-centered cubic aluminum, showing good agreement with experimental values above the Debye temperatures of the metals. Below the Debye temperatures, we extract an effective vibrational density of states from eigenvalues of the covariance matrix, and then evaluate the entropy quantum mechanically, again yielding good agreement with experiment down to low temperatures. Our method readily generalizes to complex solids, as we demonstrate for a high entropy alloy. Further, our method applies in cases where the quasiharmonic approximation fails, as we demonstrate by calculating the HCP/BCC transition in Ti.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Probing Oxidation-Driven Amorphized Surfaces in a Ta(110) Film for Superconducting Qubit

Recent advances in superconducting qubit technology have led to significant progress in quantum computing, but the challenge of achieving a long coherence time remains. Despite the excellent lifetime performance that tantalum (Ta) based qubits have demonstrated to date, the majority of superconducting qubit systems, including Ta-based qubits, are generally believed to have uncontrolled surface oxidation as the primary source of the two-level system loss in two-dimensional transmon qubits. Therefore, atomic-scale insight into the surface oxidation process is needed to make progress toward a practical quantum processor. Here, in this study, the surface oxidation mechanism of native Ta films and its potential impact on the lifetime of superconducting qubits were investigated using advanced scanning transmission electron microscopy (STEM) techniques combined with density functional theory calculations. The results suggest an atomistic model of the oxidized Ta(110) surface, showing that oxygen atoms tend to penetrate the Ta surface and accumulate between the two outermost Ta atomic planes; oxygen accumulation at the level exceeding a 1:1 O/Ta ratio drives disordering and, eventually, the formation of an amorphous Ta 2 O 5 phase. In addition, we discuss how the formation of a noninsulating ordered TaO 1–δ (δ < 0.1) suboxide layer could further contribute to the losses of superconducting qubits. Subsurface oxidation leads to charge redistribution and electric polarization, potentially causing quasiparticle loss and decreased current-carrying capacity, thus affecting superconducting qubit coherence. The findings enhance the comprehension of the realistic factors that might influence the performance of superconducting qubits, thus providing valuable guidance for the development of future quantum computing hardware.

36 MATERIALS SCIENCE↗

Prediction of nonlinear optical properties of organic materials. General theoretical considerations

The prediction of nonlinear optical properties of organic materials is geared to assist materials scientists in the selection of good candidate molecules. A brief summary of the quantum mechanical methods used for estimating hyperpolarizabilities will be presented. The advantages and limitations of each technique will be discussed. Particular attention will be given to the finite-field method for calculating first and second order hyperpolarizabilities, since this method is better suited for large molecules. Corrections for dynamic fields and bulk effects will be discussed in detail, focusing on solvent effects, conformational isomerization, core effects, dispersion, and hydrogen bonding. Several results will be compared with data obtained from third-harmonic-generation (THG) and dc-induced second harmonic generation (EFISH) measurements. These comparisons will demonstrate the qualitative ability of the method to predict the relative strengths of hyperpolarizabilities of a class of compounds. The future application of molecular mechanics, as well as other techniques, in the study of bulk properties and solid state defects will be addressed. The relationship between large values for nonlinear optical properties and large conjugation lengths is well known, and is particularly important for third-order processes. For this reason, the materials with the largest observed nonresonant third-order properties are conjugated polymers. An example of this type of polymer is polydiacetylene. One of the problems in dealing with polydiacetylene is that substituents which may enhance its nonlinear properties may ultimately prevent it from polymerizing. A model which attempts to predict the likelihood of solid-state polymerization is considered, along with the implications of the assumptions that are used. Calculations of the third-order optical properties and their relationship to first-order properties and energy gaps will be discussed. The relationship between monomeric and polymeric third-order optical properties will also be considered.

Cardelino, B.↗

Spectrum of end of the world branes in holographic BCFTs

We study overlaps between two regularized boundary states in conformal field theories. Regularized boundary states are dual to end of the world branes in an AdS black hole via the AdS/BCFT. Thus they can be regarded as microstates of a single sided black hole. Owing to the open-closed duality, such an overlap between two different regularized boundary states is exponentially suppressed as $\langle$ψ a |ψ b $\rangle$~e -O (h$^{(min)}_{ab}$) , where h$^{(min)}_{ab}$ is the lowest energy of open strings which connect two different boundaries a and b. Our gravity dual analysis leads to h$^{(min)}_{ab}$ = c/24 for a pure AdS 3 gravity. This shows that a holographic boundary state is a random vector among all left-right symmetric states, whose number is given by a square root of the number of all black hole microstates. We also perform a similar computation in higher dimensions, and find that h$^{(min)}_{ab}$ depends on the tensions of the branes. In our analysis of holographic boundary states, the off diagonal elements of the inner products can be computed directly from on-shell gravity actions, as opposed to earlier calculations of inner products of microstates in two dimensional gravity.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Real-space Green’s function approach to photoelectron diffraction

We discuss the real-space Green’s function (RSGF) approach for core-level photoelectron diffraction (PD) and its comparison to x-ray absorption spectra (XAS). We focus especially on the separable Green’s function formalism for efficient calculations of high-order multiple scattering. Computational details such as scattering potentials, self-energy effects, Debye–Waller factors, and inelastic losses are discussed. Finally, we consider prospects for improving ab initio PD calculations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Importance learning estimator for the site-averaged turnover frequency of a disordered solid catalyst

For disordered catalysts such as atomically dispersed “single-atom” metals on amorphous silica, the active sites inherit different properties from their quenched-disordered local environments. The observed kinetics are site-averages, typically dominated by a small fraction of highly active sites. Standard sampling methods require expensive ab initio calculations at an intractable number of sites to converge on the siteaveraged kinetics. We present a new method that efficiently estimates the site-averaged turnover frequency (TOF). The new estimator uses the same importance learning algorithm [Vandervelden et al., React. Chem. Eng. 5, 77 (2020)] that we previously used to compute the siteaveraged activation energy. We demonstrate the method by computing the site-averaged TOF for a simple disordered lattice model of an amorphous catalyst. The results show that with the importance learning algorithm, the site-averaged TOF and activation energy can now be obtained concurrently with orders of magnitude reduction in required ab initio calculations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Iterative subspace algorithms for finite-temperature solution of Dyson equation

One-particle Green’s functions obtained from the self-consistent solution of the Dyson equation can be employed in the evaluation of spectroscopic and thermodynamic properties for both molecules and solids. However, typical acceleration techniques used in the traditional quantum chemistry self-consistent algorithms cannot be easily deployed for the Green’s function methods because of a non-convex grand potential functional and a non-idempotent density matrix. Moreover, the optimization problem can become more challenging due to the inclusion of correlation effects, changing chemical potential, and fluctuations of the number of particles. In this paper, we study acceleration techniques to target the self-consistent solution of the Dyson equation directly. We use the direct inversion in the iterative subspace (DIIS), the least-squared commutator in the iterative subspace (LCIIS), and the Krylov space accelerated inexact Newton method (KAIN). We observe that the definition of the residual has a significant impact on the convergence of the iterative procedure. Based on the Dyson equation, we generalize the concept of the commutator residual used in DIIS and LCIIS and compare it with the difference residual used in DIIS and KAIN. The commutator residuals outperform the difference residuals for all considered molecular and solid systems within both GW and GF2. For a number of bond-breaking problems, we found that an easily obtained high-temperature solution with effectively suppressed correlations is a very effective starting point for reaching convergence of the problematic low-temperature solutions through a sequential reduction of temperature during calculations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Mass Property Calculator

A mass property calculator has been developed to compute the moment of inertia properties of an assemblage of parts that make up a system. The calculator can take input from spreadsheets or Creo mass property files or it can be interfaced with Phoenix Integration Model Center. The input must include the centroidal moments of inertia of each part with respect to its local coordinates, the location of the centroid of each part in the system coordinates and the Euler angles needed to rotate from the part coordinates to the system coordinates. The output includes the system total mass, centroid and mass moment of inertia properties. The input/output capabilities allow the calculator to interface with external optimizers. In addition to describing the calculator, this document serves as its user's manual. The up-to-date version of the calculator can be found in the Git repository https://cee-gitlab.sandia.gov/cj?ete/mass-properties-calculator.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

The synthesis and characterization of an iron(VII) nitrido complex

Abstract Complexes of iron in high oxidation states are captivating research subjects due to their pivotal role as active intermediates in numerous catalytic processes. Structural and spectroscopic studies of well-defined model complexes often provide evidence of these intermediates. In addition to the fundamental molecular and electronic structure insights gained by these complexes, their reactivity also affects our understanding of catalytic reaction mechanisms for small molecule and bond-activation chemistry. Here, we report the synthesis, structural and spectroscopic characterization of a stable, octahedral Fe(VI) nitrido complex and an authenticated, unique Fe(VII) species, prepared by one-electron oxidation. The super-oxidized Fe(VII) nitride rearranges to an Fe(V) imide through an intramolecular amination mechanism and ligand exchange, which is characterized spectroscopically and computationally. This enables combined reactivity and stability studies on a single molecular system of a rare high-valent complex redox pair. Quantum chemical calculations complement the spectroscopic parameters and provide evidence for a diamagnetic ( S = 0) d 2 Fe(VI) and a genuine S = 1/2, d 1 Fe(VII) configuration of these super-oxidized nitrido complexes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Low-energy pathways lead to self-healing defects in CsPbBr 3

Self-regulation of free charge carriers in perovskites via Schottky defect formation has been posited as the origin of the well-known defect-tolerance of metal halide perovskite materials. Understanding the mechanisms of self-regulation promises to lead to the fabrication of better performing solar cell materials with higher efficiencies. We investigate many such mechanisms here for CsPbBr 3 , a popular representative of a more commercially viable all-inorganic metal halide perovskite. We investigate different atomic-level mechanisms and pathways of the diffusion and recombination of neutral and charged interstitials and vacancies (Schottky pairs) in CsPbBr 3 . We use nudged elastic band calculations and ab initio-derived pseudopotentials within quantum ESPRESSO to determine energies of formation and migration and hence the activation energies for these defects. While halide vacancies are known to exhibit low formation energies, the migration of interstitials is less studied. Our calculations uncover interstitial defect pathways capable of producing an activation energy at, or below, the single experimental value of 0.53 eV observed for the slow, temperature-dependent recovery of light-induced conductivity in bulk CsPbBr 3 . Our work reveals the existence of a low-energy diffusion pathway involving a concerted “domino effect” of interstitials, with the net result that interstitials can diffuse more readily over longer distances than expected. This observation suggests that defect self-healing can be promoted if the “domino effect” strategy can be engaged.

14 SOLAR ENERGY↗

Ab initio investigation of a hypersonic double cone experiment

This article presents a direct molecular simulation (DMS) of a reactive Mach 8.2 oxygen flow over a double cone geometry. The free stream conditions and article configuration generate a flow with thermal and chemical nonequilibrium, which are common attributes of hypersonic flight. This scenario was first studied experimentally at Calspan University of Buffalo Research Center’s test facility. DMS is a particle method that uses quantum mechanically derived interaction potentials to simulate molecular collisions within a flow field. Since these interaction potentials are the only modeling inputs used in the simulation, all flow features can solely be attributed to the ab initio potential energy surfaces. Hence, providing a comparison of a hypersonic ground test and numerical data anchored to quantum mechanics. The fundamental nature of DMS is leveraged to investigate molecular level mechanisms prevalent in the flow, and comparisons with lower fidelity simulations are presented to highlight the role of these first principles calculations as benchmark solutions.

Science & Technology - Other Topics↗

Mechanism and dynamics of fatty acid photodecarboxylase

Photoenzymes are rare biocatalysts driven by absorption of a photon at each catalytic cycle; they inspire development of artificial photoenzymes with valuable activities. Fatty acid photodecarboxylase (FAP) is a natural photoenzyme that has potential applications in the bio-based production of hydrocarbons, yet its mechanism is far from fully understood. RATIONALE To elucidate the mechanism of FAP, we studied the wild-type (WT) enzyme from Chlorella variabilis (CvFAP) and variants with altered active-site residues using a wealth of techniques, including static and time-resolved crystallography and spectroscopy, as well as biochemical and computational approaches. RESULTS A 1.8-Å-resolution CvFAP x-ray crystal structure revealed a dense hydrogen-bonding network positioning the fatty acid carboxyl group in the vicinity of the flavin adenine dinucleotide (FAD) cofactor. Structures solved from free electron laser and low-dose synchrotron x-ray crystal data further highlighted an unusual bent shape of the oxidized flavin chromophore, and showed that the bending angle (14°) did not change upon photon absorption (step 1) or throughout the photocycle. Calculations showed that bending substantially affected the energy levels of the flavin. Structural and spectroscopic analysis of WT and mutant proteins targeting two conserved active-site residues, R451 and C432, demonstrated that both residues were crucial for proper positioning of the substrate and water molecules and for oxidation of the fatty acid carboxylate by 1 FAD* (~300 ps in WT FAP) to form FAD ∙– (step 2). Time-resolved infrared spectroscopy demonstrated that decarboxylation occured quasi-instantaneously upon this forward electron transfer, consistent with barrierless bond cleavage predicted by quantum chemistry calculations and with snapshots obtained by time-resolved crystallography. Transient absorption spectroscopy in H 2 O and D 2 O buffers indicated that back electron transfer from FAD ∙– was coupled to and limited by transfer of an exchangeable proton or hydrogen atom (step 3). Unexpectedly, concomitant with FAD ∙– reoxidation (to a red-shifted form FAD RS ) in 100 ns, most of the CO 2 product was converted, most likely into bicarbonate (as inferred from FTIR spectra of the cryotrapped FAD RS intermediate). Calculations indicated that this catalytic transformation involved an active-site water molecule. Cryo-Fourier transform infrared spectroscopy studies suggested that bicarbonate formation (step 4) was preceded by deprotonation of an arginine residue (step 3). At room temperature, the remaining CO 2 left the protein in 1.5 μs (step 4'). The observation of residual electron density close to C432 in electron density maps derived from time-resolved and cryocrystallography data suggests that this residue may play a role in stabilizing CO 2 and/or bicarbonate. Three routes for alkane formation were identified by quantum chemistry calculations; the one shown in the figure is favored by the ensemble of experimental data. CONCLUSION Finally, we provide a detailed and comprehensive characterization of light-driven hydrocarbon formation by FAP, which uses a remarkably complex mechanism including unique catalytic steps. We anticipate that our results will help to expand the green chemistry toolkit.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗