Search NASA⌕ Search

SEARCH · Search NASA

Results for “Chemical 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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Phosphonate‐Modified UiO‐66 Brønsted Acid Catalyst and Its Use in Dehydra‐Decyclization of 2‐Methyltetrahydrofuran to Pentadienes

Abstract Phosphorus‐modified all‐silica zeolites exhibit activity and selectivity in certain Brønsted acid catalyzed reactions for biomass conversion. In an effort to achieve similar performance with catalysts having well‐defined sites, we report the incorporation of Brønsted acidity to metal–organic frameworks with the UiO‐66 topology, achieved by attaching phosphonic acid to the 1,4‐benzenedicarboxylate ligand and using it to form UiO‐66‐PO 3 H 2 by post‐synthesis modification. Characterization reveals that UiO‐66‐PO 3 H 2 retains stability similar to UiO‐66, and exhibits weak Brønsted acidity, as demonstrated by titrations, alcohol dehydration, and dehydra‐decyclization of 2‐methyltetrahydrofuran (2‐MTHF). For the later reaction, the reported catalyst exhibits site‐time yields and selectivity approaching that of phosphoric acid on all‐silica zeolites. Using solid‐state NMR and deprotonation energy calculations, the chemical environments of P and the corresponding acidities are determined.

Dorneles de Mello, Matheus↗
DOE OSTI2020

Local Structure of Glassy Lithium Phosphorus Oxynitride Thin Films: A Combined Experimental and Ab Initio Approach

Abstract Lithium phosphorus oxynitride (LiPON) is an amorphous solid‐state lithium ion conductor displaying exemplary cyclability against lithium metal anodes. There is no definitive explanation for this stability due to the limited understanding of the structure of LiPON. Herein, we provide a structural model of RF‐sputtered LiPON. Information about the short‐range structure results from 1D and 2D solid‐state NMR experiments. These results are compared with first principles chemical shielding calculations of Li‐P‐O/N crystals and ab initio molecular dynamics‐generated amorphous LiPON models to unequivocally identify the glassy structure as primarily isolated phosphate monomers with N incorporated in both apical and as bridging sites in phosphate dimers. Structural results suggest LiPON′s stability is a result of its glassy character. Free‐standing LiPON films are produced that exhibit a high degree of flexibility, highlighting the unique mechanical properties of glassy materials.

Marple, Maxwell A. T.↗
DOE OSTI2020

Local Structure of Glassy Lithium Phosphorus Oxynitride Thin Films: A Combined Experimental and Ab Initio Approach

Lithium phosphorus oxynitride (LiPON) is an amorphous solid–state lithium ion conductor displaying exemplary cyclability against lithium metal anodes. There is no definitive explanation for this stability due to the limited understanding of the structure of LiPON. Herein, we provide a structural model of RF–sputtered LiPON. Information about the short–range structure results from 1D and 2D solid–state NMR experiments. These results are compared with first principles chemical shielding calculations of Li–P–O/N crystals and ab initio molecular dynamics–generated amorphous LiPON models to unequivocally identify the glassy structure as primarily isolated phosphate monomers with N incorporated in both apical and as bridging sites in phosphate dimers. Structural results suggest LiPON's stability is a result of its glassy character. In conclusion, free–standing LiPON films are produced that exhibit a high degree of flexibility, highlighting the unique mechanical properties of glassy materials.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Implementation of real‐time TDDFT for periodic systems in the open‐source PySCF software package

Abstract We present a new implementation of real‐time time‐dependent density functional theory (RT‐TDDFT) for calculating excited‐state dynamics of periodic systems in the open‐source Python‐based PySCF software package. Our implementation uses Gaussian basis functions in a velocity gauge formalism and can be applied to periodic surfaces, condensed‐phase, and molecular systems. As representative benchmark applications, we present optical absorption calculations of various molecular and bulk systems and a real‐time simulation of field‐induced dynamics of a (ZnO) 4 molecular cluster on a periodic graphene sheet. We present representative calculations on optical response of solids to infinitesimal external fields as well as real‐time charge‐transfer dynamics induced by strong pulsed laser fields. Due to the widespread use of the Python language, our RT‐TDDFT implementation can be easily modified and provides a new capability in the PySCF code for real‐time excited‐state calculations of chemical and material systems.

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

Quantum embedding for material chemistry based on domain separation and open subsystems

Our perspective considers two theories we recently proposed to perform quantum embedding calculations for chemical systems: domain-separated density functional theory (DSDFT) and locally coupled open subsystems (LCOS). The development includes both the fundamentals of each theory as well as potential applications, some technical aspects, and related challenges. DS-DFT is suited to study intramolecular effects, where one can apply a high level of theory (based on DFT or wave function theory) to a region of interest inside a molecule or solid, and lower level theory elsewhere, with smooth switching between the regions. LCOS, in contrast, is a fragment-based embedding which offers computational advantages to study intermolecular behavior such as electron hopping, spin-environment interaction, and charge transfer excitations. Yet, both theories can exchange roles when appropriate. In addition, these theories allow for control of computational scaling of their algorithms. We investigate paths to determine the charge-transfer operator used in LCOS, and suggest an auxiliary energy minimization that can provide a practical estimate to this operator. We also briefly discuss how to implement density fitting techniques in domain separation, and how domain separation can be used for pure wavefunction-based embedding.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2025

Computational materials assessment of the D/Li-stripping neutron source as a prototypical facility for fusion materials testing

As the US fusion materials community awaits the selection and design of a fusion prototypical neutron source (FPNS), a risk reduction exercise has been conducted to (i) provide an updated materials performance evaluation using state-of-the-art computational materials modeling, (ii) expand on legacy analysis based on pure Fe to other relevant fusion structural materials types, and (iii) ensure that materials response under FPNS operational conditions is consistent with referential fusion reactor conditions. The current paper describes the efforts undertaken to assemble a comprehensive computational methodology that includes neutronics, primary damage calculations, atomistic simulations of displacement cascades, chemical inventory evolution calculations, and a computational thermodynamic analysis of emerging phases during irradiation. Our work extends existing studies in pure Fe to reduced-activation ferritic/martensitic steels, tungsten, silicon carbide, and vanadium alloys. We focus on the single-beam deuteron/lithium-stripping neutron source behind the IFMIF-DONES concept, which we assess against ITER, two DEMO designs, and an ideal pure 14-MeV flux. Our analysis indicates that, within standard uncertainties inherent to the models employed, the DONES concept adequately captures fusion conditions in the four materials analyzed. Our work is intended as a comprehensive irradiation damage analysis of fusion-representative neutron sources, to be used for further neutron source evaluation and fusion facility operation.

Marian, Jaime [Univ. of California, Los Angeles, C↗
DOE OSTI2021

DFT Mechanism Studies: Biomimetic 1,4-NADH Chemoselective, Co-factor Regeneration with [Cp*Rh(bpy)H] + , in Tandem with the Biocatalysis Pathways of a Core Model of the (HLADH)-Zn(II) Mediated Enzyme, in the Enantioselective Reduction of Achiral Ketones to Chiral S-Alcohols

In this study, Quantum Chemical (QC) calculations, utilizing Density Functional Theory (DFT), were performed to investigate the mechanistic aspects of the chemoselective catalyzed reaction of [Cp*Rh(bpy)H] + with the biomimetic NAD + analogues, N-benzylnicotinamide triflate, 1, and β-nicotinamide ribose-5'-methyl phosphate, 2, in the conversion to their 1,4-NADH analogues, 1,4-dihydro-N-benzylnicotinamide, 4, and β-1,4-dihydronicotinamide-5'-ribose methyl phosphate, 5. This reaction was in tandem with the 1,4-NADH dependent HLADH-Zn(II)- catalyzed reduction of achiral ketones to chiral S-alcohols. The [Cp*Rh(bpy)H] + complex, and not its equilibrium tautomer, [η 4 -Cp*HRh(bpy)] + , was found to control the hydride transfer during the biomimetic NAD + /1,4-NADH conversion, through the non-covalent interactions of the biomimetic co-factors with [Cp*Rh(bpy)H] + . The thermodynamics and kinetics for the chiral reduction of the Zn(II) bound ketones, 2-pentanone and 4-phenyl-2-butanone, with co-factor, 4, catalyzed by Zn(SCH 3 ) 2 (Imidazole), a core model of the Zn(II)-based catalytic center of HLADH, was also investigated by the evaluation of two possible reaction pathways: (1) formation of a ZnH from the C4-H hydride transfer of co-factor, 4, followed by reaction of the postulated ZnH with the bound 2-pentanone or 4-phenyl-2-butanone substrate, and (2), the direct C4-H transfer to the bound achiral ketone substrates, to provide the dominant chiral alcohols, S-2-pentanol or S-4-phenyl-2-butanol. The latter pathway was found most viable, and DFT calculations also revealed an essential η 2 -coordination of the 5,6 double bond of co-factor, 4, to the HLADH-Zn(II) metal ion center, upon imidazole decomplexation, providing an asymmetric differentiation of S-η 2 -5,6-1,4-NADH-Zn(II) binding. A proposed new paradigm for the Zn(II)'s non-innocent role in the HLADH-Zn(II) biocatalysis reduction mechanism, for enantioselective hydride transfer to a Zn(II) bound ketone, providing S-alcohols.

1,4 NADH co-factors↗
DOE OSTI2023

Computational and Experimental Determination of the Properties, Structure, and Stability of Peptoid Nanosheets and Nanotubes

Peptoids (N-substituted glycines) are a group of highly controllable peptidomimetic polymers. Amphiphilic diblock peptoids have been engineered to assemble crystalline nanospheres, nanofibrils, nanosheets, and nanotubes with biochemical, biomedical, and bioengineering applications. The mechanical properties of peptoid nanoaggregates and their relationship to the emergent self-assembled morphologies have been relatively unexplored and are critical for the rational design of peptoid nanomaterials. In this work, we consider a family of amphiphilic diblock peptoids consisting of a prototypical tube-former (Nbrpm6Nc6, a NH 2 -capped hydrophobic block of six N-((4-bromophenyl)methyl)glycine residues conjugated to a polar NH 3 (CH 2 ) 5 CO tail), a prototypical sheet-former (Nbrpe6Nc6, where the hydrophobic block comprises six N-((4-bromophenyl)ethyl)glycine residues), and an intermediate sequence that forms mixed structures ((NbrpeNbrpm)3Nc6). We combine all-atom molecular dynamics simulations and atomic force microscopy to determine the mechanical properties of the self-assembled 2D crystalline nanosheets and relate these properties to the observed self-assembled morphologies. We find good agreement between our computational predictions and experimental measurements of Young’s modulus of crystalline nanosheets. A computational analysis of the bending modulus along the two axes of the planar crystalline nanosheets reveals bending to be more favorable along the axis in which the peptoids stack by interdigitation of the side chains compared to that in which they form columnar crystals with π-stacked side chains. Here, we construct molecular models of nanotubes of the Nbrpm6Nc6 tube-forming peptoid and predict a stability optimum in good agreement with experimental measurements. A theoretical model of nanotube stability suggests that this optimum is a free energy minimum corresponding to a “Goldilocks” tube radius at which capillary wave fluctuations in the tube wall are minimized.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Multireference Methods for Chemistry and Materials Science: Automated Active Spaces, Efficient Dynamic Correlation, and Extended Systems

While multiconfigurational approaches have long been relegated to expert practitioners working on a case-by-case basis, recent developments have increasingly made these methods more routine and applicable to broader sets of systems. This article outlines the state-of-the-art in multiconfigurational approaches, with an emphasis on moving from delicate hand-selected pathways through configuration space toward more robust and efficient approaches to treating a host of challenging chemical systems accurately. First, we overview recent work in automated active-space selection, which has enabled increasingly large-scale applications of multireference methods to modeling vertical excitations and reactivity. Second, we highlight the increasingly efficient methods for recovering correlation energy beyond the active space, as headlined by extensions of pair-density functional theory and its role in accurate and efficient treatment of excited-state dynamics and its utilization to train machine-learned potentials. Finally, we highlight recent efforts to treat extended systems that until recently have lied beyond the traditional limits of active-space methods, giving center stage to product-form wave functions of the localized active space family of methods that allow for the computation of multiconfigurational band structures. These recent advancements point to a broader use of multireference approaches for high-impact chemical and materials science applications.

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

Kβ X-ray Emission Spectroscopy as a Probe of Cu(I) Sites: Application to the Cu(I) Site in Preprocessed Galactose Oxidase

Cu(I) active sites in metalloproteins are involved in O 2 activation, but their O 2 reactivity is difficult to study due to the Cu(I) d 10 closed shell which precludes the use of conventional spectroscopic methods. Kβ X-ray emission spectroscopy (XES) is a promising technique for investigating Cu(I) sites as it detects photons emitted by electronic transitions from occupied orbitals. Here, we demonstrate the utility of Kβ XES in probing Cu(I) sites in model complexes and a metalloprotein. Using Cu(I)Cl, emission features from double-ionization (DI) states are identified using varying incident X-ray photon energies, and a reasonable method to correct the data to remove DI contributions is presented. Kβ XES spectra of Cu(I) model complexes, having biologically relevant N/S ligands and different coordination numbers, are compared and analyzed, with the aid of density functional theory (DFT) calculations, to evaluate the sensitivity of the spectral features to the ligand environment. While the low-energy Kβ 2,5 emission feature reflects the ionization energy of ligand np valence orbitals, the high-energy Kβ 2,5 emission feature corresponds to transitions from molecular orbitals (MOs) having mainly Cu 3d character with the intensities determined by ligand-mediated d–p mixing. A Kβ XES spectrum of the Cu(I) site in preprocessed galactose oxidase (GO pre ) supports the 1Tyr/2His structural model that was determined by our previous X-ray absorption spectroscopy and DFT study. Finally, the high-energy Kβ 2,5 emission feature in the Cu(I)-GO pre data has information about the MO containing mostly Cu 3d x 2 –y 2 character that is the frontier molecular orbital (FMO) for O 2 activation, which shows the potential of Kβ XES in probing the Cu(I) FMO associated with small-molecule activation in metalloproteins.

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

Multiple N–H and C–H Hydrogen Atom Abstractions Through Coordination-Induced Bond Weakening at Fe-Amine Complexes

We report the use of the reported Fe-phthalocyanine complex, PcFe (1; Pc = 1,4,8,11,15,18,22,25-octaethoxy-phthalocyanine), to generate PcFe–amine complexes 1-(NH 3 ) 2 , 1-(MeNH 2 ) 2 , and 1-(Me 2 NH) 2 . Treatment of 1 or 1-(NH 3 ) 2 to an excess of the stable aryloxide radical, 2,4,6-tritert-butylphenoxyl radical ( t Bu ArO•), under NH 3 resulted in catalytic H atom abstraction (HAA) and C–N coupling to generate the product 4-amino-2,4,6-tritert-butylcyclohexa-2,5-dien-1-one (2) and t Bu ArOH. Exposing 1-(NH 3 ) 2 to an excess of the trityl (CPh 3 ) variant, 2,6-di-tert-butyl-4-tritylphenoxyl radical ( Tr ArO • ), under NH 3 did not lead to catalytic ammonia oxidation as previously reported in a related Ru-porphyrin complex. However, pronounced coordination-induced bond weakening of both α N–H and β C–H in the alkylamine congeners, 1-(MeNH 2 ) 2 and 1-(Me 2 NH) 2 , led to multiple HAA events yielding the unsaturated cyanide complex, 1-(MeNH 2 )(CN), and imine complex, 1-(MeN=CH 2 ) 2 , respectively. Subsequent C–N bond formation was also observed in the latter upon addition of a coordinating ligand. Detailed computational studies support an alternating mechanism involving sequential N–H and C–H HAA to generate these unsaturated products.

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

Equatorial Electronic Structure in the Uranyl Ion: Cs 2 UO 2 Cl 4 and Cs 2 UO 2 Br 4

Electric field gradient (EFG) tensors in the equatorial plane of the linear UO 2 2+ ion have been measured by nuclear magnetic resonance (NMR) and nuclear quadrupole resonance (NQR) experiments and computed by relativistic Kohn–Sham methods with and without environment embedding for Cs 2 UO 2 Cl 4 and Cs 2 UO 2 Br 4 . This approach expands the possibilities for probing the electronic structure in uranyl complexes beyond the strongly covalent U–O bonds. The combined analyses find that one of the two largest principal EFG tensor components at the halogen sites points along the U–X bond (X = Cl, Br), and the second is parallel to the UO 2 2+ ion; in Cs 2 UO 2 Cl 4 , the components are nearly equal in magnitude, whereas in Cs 2 UO 2 Br 4 , due to short-range bromide–cesium interactions, the equatorial component is dominant for one pair of Br sites and the axial component is larger for the second pair. Additionally, the directions and relative magnitudes of the field gradient principal axes are found to be sensitive to the σ and π electron donation by the ligands and the model of the environment. Chlorine-35 NQR spectra of 235 U-depleted and 235 U-enriched Cs 2 UO 2 Cl 4 exhibited no uranium-isotope-dependent shift, but the resonance of the depleted sample displayed a 58% broader line width.

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

Evolution of Bonding and Magnetism via Changes in Valence Electron Count in CuFe 2– x Co x Ge 2

A series of solid solutions, CuFe 2–x Co x Ge 2 (x = 0, 0.2, 0.4, 0.8, and 1.0), have been synthesized by arc-melting and characterized by powder X-ray and neutron diffraction, magnetic measurements, Mössbauer spectroscopy, and electronic band structure calculations. All compounds crystallize in the CuFe 2 Ge 2 structure type, which can be considered as a three-dimensional framework built of fused MGe6 octahedra and MGe 5 trigonal bipyramids (M = Fe and Co), with channels filled by rows of Cu atoms. As the Co content (x) increases, the unit cell volume decreases in an anisotropic fashion: the b and c lattice parameters decrease while the a parameter increases. The changes in all the parameters are nearly linear, thus following Vegard’s law. CuFe 2 Ge 2 exhibits two successive antiferromagnetic (AFM) orderings, corresponding to the formation of a commensurate AFM structure, followed by an incommensurate AFM structure observed at lower temperatures. Additionally, as the Co content increases, the AFM ordering temperature (TN) gradually decreases, and only one AFM transition is observed for x ≥ 0.2. The magnetic behavior of unsubstituted CuFe 2 Ge 2 was found to be sensitive to the preparation method. The temperature-dependent zero-field 57 Fe Mössbauer spectra reveal two hyperfine split components that evolve in agreement with the two consecutive AFM orderings observed in magnetic measurements. In contrast, the field-dependent spectra obtained for fields ≥ 2 T reveal a parallel arrangement of the moments associated with the two crystallographically unique metal sites. Electronic band structure calculations and chemical bonding analysis reveal a mix of strong M–M antibonding and non-bonding states at the Fermi level, in support of the overall AFM ordering observed in zero field. The substitution of Co for Fe reduces the population of the M–M antibonding states and the overall density of states at the Fermi level, thus suppressing the T N value.

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

Computational Spectroscopy of the Cr–Cr Bond in Coordination Complexes

In this work, we report the accurate computational vibrational analysis of the Cr–Cr bond in dichromium complexes using second-order multireference complete active space methods (CASPT2), allowing direct comparison with experimental spectroscopic data both to facilitate interpreting the low-energy region of the spectra and to provide insights into the nature of the bonds themselves. Recent technological development by the authors has realized such computation for the first time. Accurate simulation of the vibrational structure of these compounds has been hampered by their notorious multiconfigurational electronic structure that yields bond distances that do not correlate with bond order. Some measured Cr–Cr vibrational stretching modes, ν(Cr 2 ), have suggested weaker bonding, even for so-called ultrashort Cr–Cr bonds, while others are in line with the bond distance. Here, we optimize geometries and compute ν(Cr 2 ) with CASPT2 for three well-characterized complexes, Cr 2 (O 2 CCH 3 ) 4 (H 2 O) 2 , Cr 2 (mhp) 4 , and Cr 2 (dmp) 4 . We obtain CASPT2 harmonic ν(Cr 2 ) modes in good agreement with experiment at 282 cm –1 for Cr 2 (mhp) 4 and 353 cm –1 for Cr 2 (dmp) 4 , compute 50 Cr and 54 Cr isotope shifts, and demonstrate that the use of the so-called IPEA shift leads to improved Cr–Cr distances. Additionally, normal mode sampling was used to estimate anharmonicity along ν(Cr 2 ), leading to an anharmonic mode of 272 cm –1 for Cr 2 (mhp) 4 and 333 cm –1 for Cr 2 (dmp) 4 .

36 MATERIALS SCIENCE↗
DOE OSTI2022

Structural Complexity and Tuned Thermoelectric Properties of a Polymorph of the Zintl Phase Ca 2 CdSb 2 with a Non-centrosymmetric Monoclinic Structure

The Zintl phase Ca 2 CdSb 2 was found to be dimorphic. Besides the orthorhombic Ca 2 CdSb 2 (-o), here we report on the synthesis, the structural characterization, and the thermoelectric transport properties of its monoclinic form, Ca 2 CdSb 2 (-m), and its Lu-doped variant Ca 2–x Lu x CdSb 2 (x ≈ 0.02). The monoclinic structure exhibits complex structural characteristics and constitutes a new structure type with the non-centrosymmetric space group Cm (Z = 30). The electrical resistivity ρ(T) measured on single crystals of both phases portrays a transition from a semiconductor to a degenerate p-type semiconductor upon doping with Lu and with an attendant change in the Hall carrier concentration nH from 7.15 × 10 18 to 2.30 × 10 19 cm –3 at 300 K. The Seebeck coefficient S(T) of both phases are comparable and indicate a hole-dominated carrier transport mechanism with magnitudes of 133 and 116 μV/K at 600 K for Ca 2 CdSb 2 (-m) and Ca 2–x Lu x CdSb 2 , respectively. The convoluted atomic bonding with an attendant large unit cell volume of ~4365 Å 3 drives a putative low thermal conductivity in these materials resulting in a power factor PF of 1.63 μW/cm K 2 and an estimated thermoelectric figure of merit zT of ~0.5 for Ca 2–x Lu x CdSb 2 at 600 K. Furthermore, differential scanning calorimetry results reveal the stability of these phases up to about 960 K, making them candidates for moderate temperature thermoelectric materials.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Four-Coordinate Fe N 2 and Imido Complexes Supported by a Hemilabile NNC Heteroscorpionate Ligand

Inspired by mechanistic proposals for N 2 reduction at the nitrogenase FeMo cofactor, we report herein a new, strongly σ-donating heteroscorpionate ligand featuring two weak-field pyrazoles and an alkyl donor. This ligand supports four-coordinate Fe(I)-N 2 , Fe(II)-Cl, and Fe(III)-imido complexes, which we have characterized using a variety of spectroscopic and computational methods. Structural and quantum mechanical analysis reveal the nature of the Fe–C bonds to be essentially invariant between the complexes, with conversion between the (formally) low-valent Fe-N 2 and high-valent Fe-imido complexes mediated by pyrazole hemilability. This presents a useful strategy for substrate reduction at such low-coordinate centers and suggests a mechanism by which FeMoco might accommodate the binding of both π-acidic and π-basic nitrogenous substrates.

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

Pu 5f Occupation in Plutonium Dioxide

Herein an actinide N 4,5 branching ratio analysis of PuO 2 and UO 2 has been performed, including measurements with a scanning transmission electron microscope at the Advanced Light Source. It is shown that the 5f occupation of the Pu in plutonium dioxide is n = 5. This is contrary to recent results from another technique.

36 MATERIALS SCIENCE↗