Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “binding energy difference”

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

ORNL_AISD-Ex: Quantum chemical prediction of UV/Vis absorption spectra for over 10 million organic molecules

We performed calculations of electronic excitation energies and associated oscillator strengths based on the time-dependent density-functional tight-binding (TD-DFTB) method [1]. The SMILES (Simplified molecular-input line-entry system) strings of the molecules from the AISD HOMO-LUMO database [2] were converted to a 3D atomistic structure and stored in a PDB file after preliminary geometry optimization using the Merck Molecular Force Field (MMFF94) in RDKit [3,4]. The primary information stored in the PDB file archive consists of Cartesian coordinates for each atom of the molecule in their 3D location in space, along with summary information about the structure, sequence, and experiment. We then performed molecular geometry optimization using the density-functional tight-binding (DFTB) method [5] in the electronic ground state, followed by single-point excited states calculations, as described below. We note that, since RDKit employs a random choice for the generation of molecular conformers, the molecular geometries obtained in this dataset could be different from the ones that were generated when the AISD HOMO-LUMO dataset was generated. The computed excitation energies and associated oscillator strengths can be converted to predict UV/Vis absorption spectra, where excitation energies correspond to absorption peak positions, and oscillator strengths are a good measure of the probability of absorption of visible or UV light in transitions between electronic ground and excited states. The conversion of SMILES strings to 3D Cartesian coordinates of fully DFTB-optimized molecules was successful for 10,502,904 out of 10,502,917 molecules. For these molecules, both geometry optimizations and excited states calculations were successful. The DFTB calculations did not complete for 13 molecules of the original AISD HOMO-LUMO dataset. We still provide information about the geometry of these molecules. The molecules are diverse for chemical compositions (which span 5 non-hydrogen elements: oxygen, carbon, nitrogen, fluorine, sulfur) and molecular size (the smallest molecule contains 5 non-hydrogen atoms, and the largest molecule contains 71 non-hydrogen atoms). The DFTB method [5] is an approximation to density functional theory (DFT), utilizing a minimal basis set in conjunction with a two-center approximation to the electronic Hamiltonian and overlap matrix elements. The DFTB total energy is the sum of an electronic and a repulsive energy contribution, and their calculation requires optimized electronic parameters and diatomic repulsive potential energy functions. All DFTB calculations were performed using the DFTB+ code [6] (version 21.2) and the wrapper for DFTB+ in the Atomic Simulation Environment (ASE) (version 3.22.1) [7], which performed an internal conversion of Cartesian coordinates from PDB to the .gen file format. For the geometry optimizations on the electronic ground state potential energy surface of the molecules, we have chosen the third-order DFTB (DFTB3) method [5c] and employed the matching 3ob set of electronic parameters and repulsive potentials [8]. The empirical γ-damping for hydrogen bond correction, and Grimme's D3 empirical dispersion correction with Becke-Johnson damping (D3(BJ)) [9] dispersion correction was included to improve the description of non-covalent interactions. For excited states single-point energy calculations, we employed the TD-DFTB method in conjunction with the DFTB2 method [5b] and the matching mio [5b,10] and halorg [11] parameter sets. We opted to request the simultaneous calculation of 50 excited states for singlet transition to investigate sufficient number of excited states, based on linear response theory using the Casida equation [Ref: T. A. Niehaus, S. Suhai, F. Della Sala, P Lugli, M. Elstner, G. Seifert, and Th. Frauenheim. Tight-binding approach to time-dependent density-functional response theory. Phys. Rev. B, 63:085108, 2001] and the ARPACK diagonalizer [R. B. Lehoucq, D. C. Sorensen, and C. Yang. Arpack users guide: Solution of large-scale eigenvalue problems by implicitly restarted arnoldi methods, 1997. 46, 51]. The dataset contains 1001 tar.gz files. Tar files are named as “ornl_aisd_ex_1.tar.gz†through “ornl_aisd_ex_1000.tar.gzâ€. Additionally, the 13 failed molecules are in “ornl_aisd_ex_unprocessed.tar.gzâ€. Except for the tar files listed below, each tar file contains 10,500 molecules. Tar files numbered 34, 121, 128, 352, 360, 429, 495, 509, 518, 627, 676, 668, and 862 contain 10,499 molecules each. The last tar file numbered 1000 contains 13,417 molecules. The total size of the uncompressed dataset is over 283 Gigabytes. The code for calculating the electronic excitation energies and statistical analysis of the dataset is provided at the following GitLab repository: https://github.com/ORNL/Analysis-of-Large-Scale-Molecular-Datasets-with-Python Calculating the UV spectrum of a molecule requires performing 3 main operations: 1. Converting the smiles string representation of a molecule into a geometric structure where each atom is assigned XYZ coordinates. The geometric structure is written to the file smiles.pdb. 2. Using smiles.pdb to compute the relaxed geometry of the molecule, which corresponds with the position of the atoms at the position of equilibrium at the ground state. This generates the files band.out, detailed.out, and geo_end.gen. 3. Using geo_end.gen to calculate the UV spectrum of the molecule which is written into the file EXC.DAT. Every molecule in the dataset has its own directory. The files contained in each molecule directory are as follows: 1. geo_end.gen 2. detailed.out 3. band.out 4. EXC.DAT 5. smiles.pdb REFERENCES [1] Niehaus, T. A.; Suhai, S.; Della Salla, F.; Lugli, P.; Elstner, M.; Seifert, G.; Frauenheim, Th. Tight-binding approach to time-dependent density-functional response theory. Phys. Rev. B, 2001, 63, 085108/1-9. [2] Blanchard, A.; Gounley, J.; Metha, K.; Yoo, P.; Irle, S. AISD HOMO-LUMO. DOI: 10.13139/ORNLNCCS/1869409 [3] RDKit: Cheminformatics and Machine Learning Software. 2013, [http://www.rdkit.org] [4] Tosco, P.; Stiefl, N. and Landrum, G. Bringing the MMFF force field to the RDKit: implementation and validation. J Cheminform. 2014, 6, 1–4. [5] a) Porezag, D.; Frauenheim, T.; Kohler, T.; Seifert, G.; Kaschner, Construction of tight-binding-like potentials on the basis of density-functional theory: Application to carbon, R. Phys. Rev. B 1995, 51, 12947-12957; b) Elstner, M.; Porezag, D.; Jungnickel, G.; Elsner, J.; Haugk, M.; Frauenheim, Th.; Suhai, S.; Seifert, G.; Phys. Rev. B 1998, 58, 7260-7268; c) Gaus, M.; Cui, Q.; Elstner, M. DFTB3: Extension of the Self-Consistent-Charge Density-Functional Tight-Binding Method (SCC-DFTB), J. Chem. Theory Comput. 2011, 7, 931-948; d) Cui, Q.; Elstner, M. Density functional tight binding: values of semi-empirical methods in an ab initio era, Phys. Chem. Chem. Phys. 2014, 16, 14368-14377. [6] Hourahine, B. et al. DFTB+, a software package for efficient approximate density functional theory based atomistic simulations, J. Chem. Phys. 2020, 152, 124101/1-19. [7] Larsen, A. H. et al. The atomic simulation environment—a Python library for working with atoms. J. Phys.: Cond. Matter 2017, 29, 273002. [8] Kubillus, M.; Kubar, T.; Gaus, M.; Rezac, J.; Elstner, M. Parameterization of the DFTB3 Method for Br, Ca, Cl, F, I, K, and Na in Organic and Biological Systems, J. Chem. Theory Comput. 2015, 11, 332-342. [9] Brandenburg, J. G.; Grimme, S. Accurate Modeling of Organic Molecular Crystals by Dispersion-Corrected Density Functional Tight Binding (DFTB), J. Phys. Chem. Lett. 2014, 5, 1785−1789. [10] a) Niehaus, T. A.; Elstner, M.; Frauenheim, Th.; Suhai, S. Application of an approximate density-functional method to sulfur containing compounds. J. Mol. Struct.: THEOCHEM 2001, 541, 185-94; b) Elstner, M.; Hobza, P.; Frauenheim, Th.; Suhai, S.; Kaxiras, E. Hydrogen bonding and stacking interactions of nucleic acid base pairs: A density-functional-theory based treatment. J. Chem. Phys. 2001, 114, 5149-55. [11] Kubar, T.; Bodrog, Z.; Gaus, M.; Köhler, C.; Aradi, B.; Frauenheim, Th.; Elstner, M. Parametrization of the SCC-DFTB Method for Halogens. J. Chem. Theory Comput. 2013, 9, 2939-49.

36 MATERIALS SCIENCE↗
DOE OSTI2025

First-Principles Study of Molecular Hydrogen Activation by Defects in Boron Nitride

Here, we used density functional theory simulations in combination with ab initio thermodynamics to determine the H 2 partial pressure (p H2 )-dependent energetics associated with H 2 activation and recovery at various defect sites in hexagonal boron nitride (h-BN). We found that some defects are very reactive with hydrogen, thereby definitely trapping hydrogen in defective h-BN. However, depending on hydrogen partial pressure, less reactive defect sites can be populated. Because of the lower binding capability of these sites, they would allow hydrogen to be recycled and recovered. For small defect sizes, we found that hydrogen preferentially binds to nitrogen sites by forming N–H bonds, and if no N sites are available then boron sites would be the next to bind hydrogen. Hydrogen dissociation via frustrated Lewis pair is found to be more favorable than forming only N–H bonds but only if the defect size is large enough to accommodate steric effects. For specific conditions such as T = 400 K, p H2 = 1 bar, and only considering one molecular H 2 per defect, three defects, namely, the N monovacancy, 3V(1B2N), and hexagonal 6V(3B3N) could play a role in both the activation and recycling of H 2 as they would be reacting enough to allow a favorable splitting of H 2 while not binding too strongly to allow its recovery. More broadly, a range of p H2 and hydrogen loading conditions were investigated for different types of defects and the finding suggests that p H2 could be used to fine-tune the Gibbs free energy of hydrogenation, thereby allowing several types of defects at different hydrogen loading contents to play a role in the activation/recovery process of H 2 in defective h-BN.

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

Identifying Support Effects in Au-Catalyzed CO Oxidation

Some catalytic oxide supports are more equal than others, with numerous variable properties ranging from crystal symmetry to surface chemistry and electronic structure. As a consequence, it is often very difficult to determine which of these act as the driver of performance changes observed in catalysis. In this work, we hold many of these variable properties constant with structurally similar LnScO 3 (Ln = La, Sm, and Nd) nanoparticle supports with cuboidal shapes and a common Sc-rich surface termination. Using CO oxidation over supported Au nanoparticles as a probe reaction, we observe higher activation energy and a slower rate using NdScO 3 as the support material. This change is found to correlate to the strength of CO 2 binding to the support surface, identified by temperature-programmed desorption measurements. Furthermore, the change is due to differences in the 4f electrons of the lanthanide cations, the cations’ Lewis acidity, and the inductive effect they impose.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2000

Chemisorption and Diffusion of H on a Graphene Sheet and Single-Wall Carbon Nanotubes

Recent experiments on hydrogen storage in single wall nanotubes and nanotube bundles have reported large fractional weight of stored molecular hydrogen which are not in agreement with theoretical estimates based of simulation of hydrogen storage by physisorption mechanisms. Hydrogen storage in catalytically doped nanotube bundles indicate that atomic H might undergo chemisorption changing the basic nature of the storage mechanism under investigation by many groups. Using a generalized tight-binding molecular dynamics (GTBMD) method for reactive C-H dynamics, we investigate chemisorption and diffusion of atomic H on graphene sheet and C nanotubes. Effective potential energy surfaces (EPS) for chemisorption and diffusion are calculated for graphene sheet and nanotubes of different curvatures. Analysis of the activation barriers and quantum rate constants, computed via wave-packet dynamics method, will be discussed in this presentation.

Srivastava, Deepak↗
DOE OSTI2021

Correlated Response of Protein Side-Chain Fluctuations and Conformational Entropy to Ligand Binding

The heterogeneous fast side-chain dynamics of proteins plays crucial roles in molecular recognition and binding. Site-specific NMR experiments quantify these motions by measuring the model-free order parameter (O axis 2 ) on a scale of 0 (most flexible) to 1 (least flexible) for each methyl-containing residue of proteins. In this work, we have examined ligand-induced variations in the fast side-chain dynamics and conformational entropy of calmodulin (CaM) using five different CaM–peptide complexes. O axis 2 of CaM in the ligand-free (O axis,U 2 ) and ligand-bound (O axis,B 2 ) states are calculated from molecular dynamics trajectories and conformational energy surfaces obtained using the adaptive biasing force (ABF) method. ΔO axis 2 = O axis,B 2 – O axis,U 2 follows a Gaussian-like unimodal distribution whose second moment is a potential indicator of the binding affinity of these complexes. The probability for the binding-induced O axis,U 2 → O axis,B 2 transition decreases with increasing magnitude of ΔO axis 2 , indicating that large flexibility changes are improbable for side chains of CaM after ligand binding. A linear correlation established between ΔOaxis2 and the conformational entropy change of the protein makes possible the determination of the conformational entropy of binding of protein–ligand complexes. The results not only underscore the functional importance of fast side-chain fluctuations but also highlight key motional and thermodynamic correlates of protein–ligand binding.

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

Coordination-Induced Weakening of N–H Bonds Driven by Bimetallic Cooperativity in Zr/Co Compounds

The bond dissociation free energy (BDFE) of the element-hydrogen bonds of protic substrates have been found to decrease upon metal coordination. Herein, an early/late heterobimetallic complex is used to examine the impact on the BDFE N−H when the substrate binding site and the redox-active site are two different metals that are spatially separated. A tris- (phosphinoamide) framework is used to link a d 0 ZrIV center with an accessible substrate binding site to a coordinatively saturated redox-active Co center, which serves as an appended electron reservoir. A series of aniline, amido, and imido Zr/Co model compounds were synthesized starting from the Zr IV /Co −I aniline adduct PhH 2 N−Zr(MesNP i Pr 2 ) 3 CoCN t Bu (2). 2,4,6-tristert- butylphenoxyl radical ( t Bu 3 ArO • ) was used to abstract one or two H atoms and produce the amido and imido complexes PhHNZr( MesNP i Pr 2 ) 3 CoCN t Bu (3) and PhN≡Zr(MesNP i Pr 2 ) 3 CoCNtBu (4), respectively. Using open-circuit potential measurements, the BDFE N−H within 2 and 3 were determined to be 37 kcal/mol (2) and 55 kcal/mol (3). Cyclic voltammetry measurements were conducted to determine the Co I/0 and Co 0/−I redox potentials. The pK a s were then estimated using the Bordwell equation to provide further insight into the thermochemical aspects of the observed proton coupled electron transfer (PCET) reactions.

Bond dissociation free energy↗
DOE OSTI2024

How Does HF-DFT Achieve Chemical Accuracy for Water Clusters?

Bolstered by recent calculations of exact functional-driven errors (FEs) and density-driven errors (DEs) of semilocal density functionals in the water dimer binding energy , we investigate approximate FEs and DEs in neutral water clusters containing up to 20 monomers, charged water clusters, and alkali- and halide-water clusters. Our proxy for the exact density is r 2 SCAN 50, a 50% global hybrid of exact exchange with r 2 SCAN, which may be less correct than r 2 SCAN for the compact water monomer but importantly more correct for long-range electron transfers in the noncompact water clusters. We show that SCAN makes substantially larger FEs for neutral water clusters than r 2 SCAN, while both make essentially the same DEs. Unlike the case for barrier heights, these FEs are small in a relative sense and become large in an absolute sense only due to an increase in cluster size. SCAN@HF, short for SCAN evaluated on the Hartree–Fock (HF) density, produces a cancellation of errors that makes it chemically accurate for predicting the absolute binding energies of water clusters. Likewise, adding a long-range dispersion correction to r 2 SCAN@HF, as in the composite method HF-r 2 SCAN-DC4, makes its FE more negative than in r 2 SCAN@HF, permitting a near-perfect cancellation of FE and DE. r 2 SCAN by itself (and even more so, r 2 SCAN evaluated on the r 2 SCAN 50 density), is almost perfect for the energy differences between water hexamers, and thus probably also for liquid water away from the boiling point. Thus, the accuracy of composite methods like SCAN@HF and HF-r 2 SCAN-DC4 is not due to the HF density being closer to the exact density, but to a compensation of errors from its greater degree of localization. We also give an argument for the approximate reliability of this unconventional error cancellation for diverse molecular properties. Lastly, we confirm this unconventional error cancellation for the SCAN description of the water trimer via Kohn–Sham inversion of the CCSD(T) density.

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

Dynamics-Based Peptide–MHC Binding Optimization by a Convolutional Variational Autoencoder: A Use-Case Model for CASTELO

An unsolved challenge in the development of antigen-specific immunotherapies is determining the optimal antigens to target. Comprehension of antigen–major histocompatibility complex (MHC) binding is paramount toward achieving this goal. Here, we apply CASTELO, a combined machine learning-molecular dynamics (ML-MD) approach, to identify per-residue antigen binding contributions and then design novel antigens of increased MHC-II binding affinity for a type 1 diabetes-implicated system. We build upon a small-molecule lead optimization algorithm by training a convolutional variational autoencoder (CVAE) on MD trajectories of 48 different systems across four antigens and four HLA serotypes. We develop several new machine learning metrics including a structure-based anchor residue classification model as well as cluster comparison scores. ML-MD predictions agree well with experimental binding results and free energy perturbation-predicted binding affinities. Moreover, ML-MD metrics are independent of traditional MD stability metrics such as contact area and root-mean-square fluctuations (RMSF), which do not reflect binding affinity data. Finally, our work supports the role of structure-based deep learning techniques in antigen-specific immunotherapy design.

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

Monte Carlo MP2-F12 for Noncovalent Interactions: The C 60 Dimer

A scalable stochastic algorithm is presented that can evaluate explicitly correlated (F12) second-order many-body perturbation (MP2) energies of weak, noncovalent, intermolecular interactions. It first transforms the formulas of the MP2 and F12 energy differences into a short sum of high-dimensional integrals of Green’s functions in real space and imaginary time. Furthermore, these integrals are then evaluated by the Monte Carlo method augmented by parallel execution, redundant-walker convergence acceleration, direct-sampling autocorrelation elimination, and control-variate error reduction. By sharing electron-pair walkers across the supermolecule and its subsystems spanned by the joint basis set, the statistical uncertainty is reduced by one to 2 orders of magnitude in the MP2 binding energy corrected for the basis-set incompleteness and superposition errors. The method predicts the MP2-F12/aug-cc-pVDZ binding energy of 19.1 ± 4.0 kcal mol –1 for the C 60 dimer at the center distance of 9.748 Å.

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

Combined multiplet theory and experiment for the Fe 2p and 3p XPS of FeO and Fe 2 O 3

The Al K alpha, 1486.6eV, based XPS of Fe 2p and Fe 3p for Fe(III) in Fe 2 O 3 and Fe(II) in FeO are compared with theoretical predictions based on ab initio wavefunctions that accurately treat the final, core-hole, multiplets. The principal objectives of this comparison are to understand the multiplet structure and to evaluate the use of both the 2p and 3p spectra in determining oxidation state. In order to properly interpret the features of these spectra and to use the XPS to provide atomistic insights as well as atomic composition, it is necessary to understand the origin of the multiplet energies and intensities. The theoretical treatment takes into account the ligand field and spin-orbit splittings, the covalent mixing of ligand and Fe 3d orbitals, and the angular momentum coupling of the open shell electrons. These effects lead to the distribution of XPS intensity into a large number of final, ionic, states that are only partly resolved with energies spread over a wide range of binding energies, BEs. For this reason, it is necessary to record the Fe 2p and 3p XPS spectra over a wide energy range which includes all the multiplets in the theoretical treatment as well as additional shake satellites. We also evaluate the effects of differing assumptions concerning the extrinsic background subtraction, to be sure our experimental spectrum may be fairly compared to the theory. In this work, we conclude that the Fe 3p XPS provides an additional means for distinguishing Fe(III) and Fe(II) oxidation states beyond just using the Fe 2p spectrum. In particular, with the use of the Fe 3p XPS, the depth of the material probed is about 1.5 times greater than for the Fe 2p XPS. In addition, a new type of atomic many-body effect that involves excitations into orbitals that have Fe f,$\ell$=3, symmetry has been shown to be important for the Fe 3p XPS.

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

Reduction of vanadium diffusivity within copper grain boundaries due to enhanced binding

We investigate permeation of vanadium (V) through copper (Cu) layers in V/Cu/niobium (Nb) and molybdenum (Mo)/Cu/V tri-layer composites and compare it with the permeation behavior of Nb. Despite their low solubility in single-crystal Cu, both V and Nb permeate readily due to short-circuit transport through Cu grain boundaries (GBs). However, Nb permeates preferentially to V. Atomistic simulations suggest that this difference is due to lower V diffusivity at general GBs in Cu, compared to Nb. Here, the lower V diffusivity, in turn, arises from stronger V binding to the boundary, which manifests as a greater equilibrium GB concentration as well as a higher activation energy for GB diffusion.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Understanding 2p core-level excitons of late transition metals by analysis of mixed-valence copper in a metal–organic framework

The L 2,3 -edge X-ray absorption spectra of late transition metals such as Cu, Ag, and Au exhibit absorption onsets lower in energy for higher oxidation states, which is at odds with the measured spectra of earlier transition metals. Time-dependent density functional theory calculations for Cu 2+ /Cu + reveal a larger 2p core-exciton binding energy for Cu 2+ , overshadowing shifts in single-particle excitation energies with respect to Cu + . Finally, we explore this phenomenon in a Cu + metal–organic framework with ~12% Cu 2+ defects and find that corrections with self-consistent excited-state total energy differences provide accurate XAS peak alignment.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1990

Vacancy formation and extraction energies in semiconductor compounds and alloys

A model for calculating the extraction energies and vacancy-formation energies in diamond-cubic and zinc-blende semiconductors is developed on the basis of Harrison's (1980, 1983) tight-binding theory. The extraction energies provide a reference from which other final states of the removed atoms can be calculated. The results of calculations show that, in a given compound, the calculated extraction energies are larger for the anion than for the cation, with the difference between the cation and the anion being larger in the II-VI than in the III-V compounds. This is in agreement with experimental results.

Berding, M. A.↗
NASA NTRS1986

Comparative isosteric ion adsorption for minerals

A comparative isosteric ion adsorption study for minerals (kaolinite, rutile, and quartz) was performed in aqueous solutions of CaCl2, LaCl3, and Th(NO3)4 in the presence of the neutral salt NaCl. It was observed that the concentration of Ca(2+) ions required to produce a standard reduction in the electrophoretic mobility of mineral particles was always appreciably greater than the concentration required for the Th(4+) ions. The effectiveness of adsorption of the cations differed from particle to particle and showed that ion adsorption on a mineral surface depends, among other things, on the nature of the mineral surface and on the particular adsorbed cation. The number of cation binding sites on mineral surfaces and the electrochemical free energies of cation adsorption were calculated. It was found that the adsorption energy of La(3+) and Th(4+) ions on rutile, kaolinite, and quartz was greater than that of Ca(2+) on these minerals.

Omenyi, Samuel N.↗
DOE OSTI2024

Chemical changes from N-doped graphene and Metal-Organic Frameworks to N-G/MOF composites for improved electrocatalytic activity

Integrating N doped graphene (N G) with Metal Organic Frameworks (MOFs) enhances catalytic activity for oxygen reduction reaction (ORR), often exceeding both the performances of their precursors and, in some cases, even PGM based catalysts. However, the factors driving this improved catalytic activity in N G/MOF composites remain unexplored, particularly from the perspective of the chemical changes. To investigate the chemical changes in N G and MOF upon their integration and the implications of these changes on ORR catalytic activity, an N G/MOF was synthesized from N G with a MOF (ZIF 8) following a mechanochemical wet ball milling process. The N G, ZIF 8, and N G/MOF samples were examined for changes in elemental composition, chemical state of carbon, different nitrogen and carbon bonds, and other chemical interactions. In the N G/MOF catalyst, compared to its N G and ZIF 8 precursors, the relative oxygen content increased, indicating the formation of additional oxygen containing groups. The C 1s peak shifted to a lower binding energy in N G/MOF, suggesting changes in the overall chemical or oxidation state of the carbon atoms. Besides, the increase in pyridinic N functional groups in N G/MOF points to the formation of additional active sites. Furthermore, the formation of C Zn bonds in N G/MOF suggests the probable emergence of single atom Zn sites, while the increase in C=O bonds points to the formation of carboxyl or carbonyl groups. These chemical changes could be linked to the enhanced electrocatalytic activity of the N G/MOF composite for ORR. In conclusion, this study may also be beneficial for other research focused on developing composite catalysts involving various N G and MOF materials.

25 ENERGY STORAGE↗
DOE OSTI2020

Accuracy of XAS theory for unraveling structural changes of adsorbates: CO on Ni(100)

Studying surface reactions using ultrafast optical pump and X-ray probe experiments relies on accurate calculations of X-ray spectra of adsorbates for correct identification of the spectral signatures and their dynamical evolution. We show that experimental XAS can be well reproduced for the different binding sites in a static prototype system CO/Ni(100) at a standard DFT (GGA) level of theory, using a plane-wave basis and pseudopotentials. This validates its utility in analyzing ultrafast X-ray probe experiments. The accuracy of computed relative core level binding energies is about 0.2 eV, representing a lower limit for which spectral features can be resolved with this method. We also show that the commonly used Z +1 approximation gives very good core binding energy shifts overall. However, we find a discrepancy for CO adsorbed in the hollow site, which we assign to the significantly stronger hybridization in hollow bonding than in on-top.

74 ATOMIC AND MOLECULAR PHYSICS↗
NASA NTRS2002

An Accurate Quartic Force Field, Fundamental Frequencies, and Binding Energy for the High Energy Density Material T(d)N4

The CCSD(T) method has been used to compute a highly accurate quartic force field and fundamental frequencies for all N-14 and N-15 isotopomers of the high energy density material T(sub d)N(sub 4). The computed fundamental frequencies show beyond doubt that the bands observed in a matrix isolation experiment by Radziszewski and coworkers are not due to different isotopomers of T(sub d)N(sub 4). The most sophisticated thermochemical calculations to date yield a N(sub 4) -> 2N(sub 2) heat of reaction of 182.22 +/- 0.5 kcal/mol at 0 K (180.64 +/- 0.5 at 298 K). It is hoped that the data reported herein will aid in the ultimate detection of T(sub d)N(sub 4).

Lee, Timothy J.↗
DOE OSTI2021

Photoluminescence Dynamics Defined by Exciton Trapping Potential of Coupled Defect States in DNA-Functionalized Carbon Nanotubes

Chemical reactions between semiconducting single-wall carbon nanotubes (SWCNTs) and single-stranded DNA (ssDNA) achieve spatially patterned covalent functionalization sites and create coupled fluorescent quantum defects on the nanotube surface, tailoring SWCNT photophysics for applications such as single-photon emitters in quantum information technologies. The evaluation of relaxation dynamics of photoluminescence (PL) from those coupled quantum defects is essential for understanding the nanotube electronic structure and beneficial to the design of quantum light emitters. Here, we measured the PL decay for ssDNA-functionalized SWCNTs as a function of the guanine content of the ssDNA oligo that dictates the red-shifting of their PL emission peaks relative to the band-edge exciton. We then correlate the observed dependence of PL decay dynamics on energy red-shifts to the exciton potential energy landscape, which is modeled using first-principles approaches based upon the morphology of ssDNA-altered SWCNTs obtained by atomic force microscopy (AFM) imaging. In this work, our simulations illustrate that the multiple guanine defects introduced within a single ssDNA strand strongly interact to create a deep exciton trapping well, acting as a single hybrid trap. The emission decay from the distinctive trapping potential landscape is found to be biexponential for ssDNA-modified SWCNTs. We attributed the fast time component of the biexponential PL decay to the redistribution of exciton population among the lowest energy bright states and a manifold of dark states emerging from the coupling of multiple guanine defects. The long lifetime component in the biexponential decay, on the other hand, is attributed to the redistribution of exciton population among different exciton trapping sites that arise from the binding of multiple ssDNA strands along the nanotube axis. AFM measurements indicate that those trapping sites are separated on average by ~8 nm along the nanotube axis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗