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

Bond Energies of UO + and UC + : Guided Ion Beam and Quantum Chemical Studies of the Reactions of Uranium Cation with O 2 and CO

Here, guided ion beam tandem mass spectrometry was used to examine the kinetic energy dependent reactions of U + with O 2 and CO. In the reaction of U + with O 2 , the UO + product is formed in a barrierless and exothermic process with a reaction efficiency at low energies of k/k col =1.1±0.2, but increases at higher collision energies. Formation of both UO + and UC + in the reaction of U + with CO is endothermic. 0 K bond dissociation energies (BDEs) of D 0 (U + -O)=7.88±0.09 eV and D 0 (U + -C)=4.03±0.13 eV were determined by analyzing the kinetic energy dependent cross sections in the latter endothermic reactions. These values are within experimental uncertainty of previously reported experimental values. Additionally, the electronic states of UO + and UC + and the potential energy surfaces for the reactions were explored by quantum chemical calculations. The former include a full Feller-Peterson-Dixon composite approach with correlation contributions up to CCSDT(Q) for UO and UO + , yielding D 0 (U-O)=7.82 eV and D 0 (U + -O)=7.99 eV, as well as more approximate CCSD(T) calculations where a semi-empirical model was used to estimate spin-orbit energy contributions, which are generally found to improve the agreement with experiment. Both experimental BDEs are observed to be close to those of their transition metal congeners, ScL + , YL + , and GdL + (L=O and C).

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

Macroscale properties and atomic-scale mechanisms of ash removal in low-temperature hydrothermal carbonization

Biogenic ash is a significant impediment to the utilization of agricultural residues in biofuel production. Such challenge can be addressed by various treatments, as demonstrated in this study on the experimental and computational mechanisms involved in the hydrothermal treatment (HT) of wheat straw. A combination of classical (all-atom) molecular dynamics simulations of cellulose carrying silica and calcium species, along with first principles quantum chemical calculations, indicates the dissociation of inorganics from the cellulose with increased HT temperature. Here, this observation is confirmed by experimental evidence of effective ash removal by HT, showing at least 50% removal of sulfur, chlorine, potassium, and calcium, and 12.5% of silica, leading to a reduced total ash content (from 6.7% to 4.2%).Changes in structural features upon HT, such as surface cellular structure and porosity, were revealed, accompanied by an increased specific surface area (from 1.17 to 6.34 m 2 /g). Our simulations suggest that silica binds tightly to the hydrophobic face of cellulose at room temperature, but HT significantly reduces the binding free energy of association with both hydrophobic and hydrophilic surfaces. Most significantly, ash removal leads to an increased calorific value, rising from approximately 16 MJ/kg to about 19 MJ/kg, along with improved thermal behavior. The improved integration combustion index parameter S indicates that the combustion properties improve with ash removal efficiency. The proposed atomic-level mechanism for the observed removal of inorganics during mild HT underscores the potential of such treatment in producing energy-dense wheat straw, a widely available agricultural residue.

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

Re-evaluation of rate constants for the reaction N 2 H 4 (+ M) $^→_←$ NH 2 + NH 2 (+ M)

Rate constants for the dissociation/recombination reaction N 2 H 4 (+ M) $^→_←$ NH 2 + NH 2 (+ M) are determined by a combination of quantum-chemical calculations and statistical unimolecular rate theory. Between 1100 and 2500 K, limiting low-pressure rate constants for hydrazine dissociation in the bath gas Ar of k 0 = [Ar] 6.1 × 10 20 (T/1000 K) –7.3 exp (-34490 K/T) cm 3 mol –1 s –1 , limiting high-pressure rate constants of k ∞ = 7.6 × 10 16 (T/1000 K) –1.0 exp(-33600 K/T) s –1 , and center broadening factors of the falloff curves (between 1100 and 1600 K in Ar) of F cent = 0.71 exp(- T/1460 K) + 0.29 exp(- T/21 K) + exp(-13400 K/T) were calculated. Using equilibrium constants K C = 1.7 × 10 3 (T/1000 K) –1.5 exp(-33,460 K/T) mol cm –3 , between 300 and 600 K limiting low-pressure rate constants for the reverse recombination of NH 2 radicals in the bath gas N 2 of k rec,0 = [N 2 ] 4.4 × 10 20 (T/400 K) –6.9 exp (-1630 K/T) cm 6 mol –2 s –1 , limiting high-pressure rate constants of k rec,∞ = 4.4 × 10 13 (T/1000 K) 0.44 exp(-140 K/T) cm 3 mol -1 s –1 , and center broadening factors (between 300 and 600 K in N 2 ) of F cent = exp(- T/1130 K) + exp (-10340 K/T) were obtained. A comparison with experimental results for hydrazine dissociation from the literature suggests incomplete falloff extrapolations toward k ∞ and experimental problems in the determination of k 0 at temperatures above about 1600 K. Finally, implications of the present re-evaluated rate constants for the modeling of high temperature ammonia oxidation kinetics are discussed, showing an only small influence of their precise values on the overall properties of the process.

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

Computational and experimental study of different brines in temperature swing solvent extraction desalination with amine solvents

Rapid global urbanization and high-salinity wastewater disposal from industrial activities have exerted significant pressure on water resources. Over the past few years, temperature swing solvent extraction (TSSE) has been identified as a promising technique to desalinate hypersaline brines. Despite its potential, the TSSE desalination literature has been mainly based on empirical insights, and the limited molecular simulation studies have primarily focused on NaCl brines. Herein, we use molecular dynamics (MD) simulations to study the TSSE desalination of four different brines, namely, KCl, KBr, NaCl, and NaBr using diisopropylamine as the solvent. Based on both bulk and interfacial brine-diisopropylamine MD results, here we investigate the qualitative and quantitative performance of the simulations by benchmarking these results against our experimental evaluations of these same systems. MD results provide satisfactory qualitative agreement with the experimental data of water solubilization in the organic phase and amine solubilization in the aqueous phase for the KBr, KCl, and NaBr brines. Also, the molecular mechanism of solvation of ionic species by water molecules over diisopropylamine suggested by the MD simulations is in agreement with our experimental data. However, larger qualitative and quantitative deviations were observed for the NaCl brines, and this is likely due to polarization and charge transfer effects, as quantified by our quantum chemical calculations.

42 ENGINEERING↗
DOE OSTI2026

Scaling dynamics in low-salt-rejection reverse osmosis for high-salinity produced water desalination: Mechanistic modeling and membrane autopsy

Membrane scaling remains a critical barrier to the reliable operation of desalination systems, particularly for hypersaline produced water (PW) treatment. This study fills the knowledge gap of autopsy-based model validation for PW desalination by elucidating scaling mechanisms in a Low-Salt-Rejection Reverse Osmosis (LSRRO) system through the integration of pilot-scale experimentation and complementary modeling approaches. A semi-empirical modeling framework was developed and applied to a multistage pilot LSRRO system equipped with nanofiltration and RO membranes treating high-salinity PW from the Permian Basin. Water quality analysis showed that total dissolved solids decreased from ~130,000 mg/L to ~1900 mg/L in the permeate, then further reduced to ~300 mg/L by a second-pass RO. Two different thermodynamic modeling approaches were evaluated: the first extends the LSRRO framework by incorporating system complexity and scaling phenomena, whereas the second method explicitly captures concentration polarization in localized supersaturation. Both methods illustrate the tendency for carbonate and sulfate scaling throughout the stages. Membrane autopsies revealed a silica-dominated deposit matrix, localized CaSO 4 at Stage 2, and minor barite/celestite despite their prominence in model predictions. Quantum-chemical calculations indicated silica scaling can be rationalized by favorable adsorption of H 4 SiO 4 on Fe-oxide surfaces (ΔG ≈ −44 kJ/mol), providing a kinetic pathway for interfacial inorganic polymerization even when bulk equilibrium predictions are conservative. Overall, the thermodynamic scaling modeling and membrane autopsy revealed heterogeneous, localized deposits with limited impact on LSRRO performance, while quantum analysis rationalized the thermodynamically unfavorable precipitation formation under bulk equilibrium, reconciling model–autopsy discrepancies. These insights support targeted pretreatment and silica-specific antiscalants to extend membrane lifetime and increase recovery, providing a transferable framework for hypersaline water desalination systems. The combined experimental–computational approach provides new mechanistic insight into scaling in hypersaline membrane systems and establishes a transferable framework for predicting and mitigating scaling in next-generation desalination technologies.

Low-salt-rejection reverse osmosis↗
DOE OSTI2025

Consistent thermodynamic properties for alicyclic components of jet fuels: Experimental data, estimation methods, and homologous series trends

Alkylcycloalkanes represent a significant fraction of jet fuel components. An evaluation of their thermodynamic properties, enthalpies of formation in liquid and gas phases and enthalpies of vaporization, was conducted. A combination of available experimental data, up-to-date group-contribution methods, high-level quantum-chemical calculations, and homologous series trends was used to identify outliers and to recommend the most reliable values. The group-contribution approach was found to work well for the enthalpies of vaporization. Its performance for the enthalpies of formation in the liquid and gas phases was found to be substantially less effective, especially considering notable differences in this property among stereoisomers. Computationally affordable high-level ab initio results and homologous series trend analysis appeared more reliable. In conclusion, the recommended property values for 212 individual compounds and their isomeric mixtures were provided.

09 BIOMASS FUELS↗
DOE OSTI2020

Experimental and kinetic modeling study of tetralin: A naphtheno-aromatic fuel for gasoline, jet and diesel surrogates

Distillate fuels contain significant proportions of naphtheno-aromatic components and tetralin is a suitable surrogate component to represent this molecular moiety. The presence of aromatic and naphthyl rings makes kinetic modeling of tetralin very challenging. Primary radicals formed during the oxidation of tetralin can be aryl, benzylic or paraffinic in nature. Using available information on reaction paths and rate constants of naphthenes and alkyl-aromatics, in this study a kinetic model of tetralin has been developed with emphasis on low-temperature chemistry and high-pressure conditions. Due to the lack of high-level quantum chemical calculations on reaction pathways of tetralin, analogous rates from ab-initio studies on benzylic and paraffinic radicals have been adopted here. Some modifications to the reaction rate rules are incorporated to account for the unique characteristics of tetralin's molecular structure. Important reaction channels have been identified using reaction path and brute force sensitivity analyses. In order to investigate the model performance at low temperatures, new experiments are carried out in a rapid compression machine on blends of tetralin and 3-methylpentane. Blending of low-reactivity tetralin with a high-reactivity alkane allowed the investigation of tetralin ignition at very low temperatures (665 – 856 K). The kinetic model developed in the current study is found to predict the current experiments and literature data adequately. The new model will aid in high-fidelity surrogate predictions at engine-relevant conditions.

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

Mechanisms of Degradation of Toxic Nerve Agents: Quantum-chemical Insight into Interactions of Sarin and Soman with Molybdenum Dioxide

With the urgent need to provide an efficient and reliable protection for people from fatal chemical and biological weapons, our fundamental understanding of how toxins interact with filters is far from complete. The situation is further complicated by natural difficulties of performing experimental measurements with lethal toxins. Unlike experiments, computational modeling offers an attractive and safe yet reliable way of studying behavior of toxic agents on a variety of substrates at a great level of detail. Here, we report DFT-based quantum-chemical calculations of adsorption and decomposition of DMMP, sarin and soman on MoO 2 (011) surface. Our calculations show that MoO 2 strongly adsorbs toxic nerve agents and quickly decomposes them. Decomposition of DMMP on the MoO 2 (011) surface proceeds via the PO-CH 3 bond breaking and a formation of a surface methoxy group. The calculated activation barrier for this reaction is 131.5 kJ mol -1 . Unlike DMMP, decomposition of sarin and soman proceeds via the dealkylation reaction yielding propene and 3,3-dimethyl-1-butene, respectively. Decomposition of sarin requires a remarkably low energy (53.7 kJ mol- 1 ), whereas the similar reaction in soman requires 50 kJ mol -1 more energy. We also make specific predictions to guide Ambient-Pressure X-ray Photoelectron Spectroscopy (APXPS) experiments on sarin interaction with MoO 2 samples. We conclude that MoO 2 serves as an efficient substrate able of degrading toxins.

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

Atomic-Scale Structural Mapping of Active Sites in Monolayer PGM-Free Catalysts by Low-Voltage 4D-STEM

Two-dimensional (2D) materials have attracted a large amount of attention in both basic and applied fields, and scanning transmission electron microscopy (STEM) is often uniquely well-suited for characterizing the atomic-scale structure of these materials [1-4]. As a result, STEM is poised to significantly impact progress on platinum group metal (PGM)-free catalysts, which are currently under intense development to enable low-cost, commercially viable hydrogen fuel cells [5]. While recent advancements have resulted in fuel cell performance comparable to Pt catalysts by some measures [6], cell durability remains a significant challenge, limiting practical applications [7]. Catalytically active sites in PGM-free materials are proposed to be FeN4 complexes embedded in a graphene lattice (Fig. 1b) within layered or other larger materials, but this is still under debate largely due to the range of potential actives sites predicted by computational methods and lack of methods for directly validating these models [5]. Fundamental insights into the atomic structure and resulting degradation pathways of proposed active sites are therefore needed to fully understand and control cell performance and durability [6].2D materials typically make ideal samples for STEM, but those within PGM-free catalysts present additional challenges since these materials are often defect-rich, with a high density of edges, dopant atoms, etc., which significantly increase susceptibility to beam damage at standard operating voltages. This makes analysis of potential FeN4 active sites particularly challenging, since a large proportion of Fe exists at edge sites where beam-induced atomic displacements can prohibit high-resolution structural characterization [6]. Conventional dark-field imaging compounds this problem by producing less signal for a given dose and being less sensitive to light elements than dose-efficient phase contrast imaging techniques such as those enabled by four-dimensional (4D)-STEM [8-10] (Fig. 1a). Consequently, active site structural analysis is often left to methods such as low-resolution imaging combined with quantum chemical calculations [6], which hinders accurate determination of reaction and degradation mechanisms.Here, we demonstrate direct atomic-scale structural mapping of FeN4 sites by performing low-voltage 4D-STEM on a model PGM-free catalyst system with many exposed monolayer regions. To accomplish this, we pair a 30 keV aberration-corrected probe with a fast pixelated detector that has optimal performance at low beam voltages [11]. This enables us to simultaneously image light and heavy elements with high signal-to-noise by center-of-mass analysis (Fig. 1c) while minimizing beam-induced atomic displacements at sensitive sites. The monolayer nature of these materials additionally allows for experimental validation by direct comparison with multislice simulations [12] of model structures (Fig. 1d-e). This work demonstrates how low-voltage 4D-STEM will provide new insights into the atomic-scale structure and degradation mechanisms of active sites in PGM-free catalysts, facilitating the development of low-cost hydrogen fuel cells and other energy conversion technologies in the future [13].

Zachman, Michael↗
DOE OSTI2025

Map of the Zintl AM 2 Pn 2 Compounds: Influence of Chemistry on Stability and Electronic Structure

The AM 2 Pn 2 (A= Ca, Sr, Ba, Yb, Mg; M = Zn, Cd, Mg; and Pn = N, P, As, Sb, Bi) family of Zintl phases has been known as thermoelectric materials and has recently gained much attention for highly promising materials for solar absorbers in single-junction and tandem solar cells. In this paper, we will, from first principles, explore the entire family of AM 2 Pn 2 compounds in terms of their ground-state structure, thermodynamic stability, and electronic structure. We also perform photoluminescence spectroscopy on bulk powder and thin film samples to verify our results, including the first measurements of the band gaps of SrCd 2 P 2 and CaCd 2 P 2 . The AM 2 Pn 2 compounds exhibit broad stability, are mostly isostructural to CaAl 2 Si 2 (P$\overline{3}$m1), and cover a wide range of band gaps from 0 to beyond 3 eV. This could make them useful for a variety of purposes, for which we propose several candidates, such as CaZn 2 N 2 for tandem top cell solar absorbers and SrCd 2 Sb 2 and CaZn 2 Sb 2 for infrared detectors. By examining the band structures of the AM 2 Pn 2 , we find that Mg 3 Sb 2 has the most promise as a thermoelectric material due to several off-Γ valence band pockets, which are unique to it among the compositions studied here.

14 SOLAR ENERGY↗
DOE OSTI2026

A Bis(arenesulfonyl) Peroxide, an Ambient-Stable Oxidant, Is a Strong p-Dopant for Organic Semiconductors

Strong p-dopants are required to dope high-ionization energy organic semiconductors for a variety of potential applications, but strong, simple one-electron oxidants are typically sensitive to reduction by atmospheric moisture and thus challenging to store or handle. Here we show that bis(3,5-bis(trifluoromethyl)benzenesulfonyl) peroxide─a dimer formed by two highly oxidizing radicals─can function as a water-stable yet powerful oxidant, cleanly reacting with some organic semiconductors to form two radical cations and two 3,5-bis-(trifluoromethyl)benzenesulfonate anions, although in other cases sulfonylation reactions can also occur. Notably, this peroxide is capable of p-doping the high-ionization energy polymer poly[(9,9-dioctylfluorene-2,7-diyl)- alt -(benzo[2,1,3]thiadiazol-4,7-diyl)] (F8BT) to afford electrical conductivities of up to 0.03 S cm –1 , while its use with electron-rich poly(3,4-dialkoxythiophene-2,5-diyl) derivatives can afford values up to 100 S cm –1 . Quantum-chemical calculations reveal the peroxide oxidant behaves in a fashion mirroring that of relatively oxygen-stable, but highly reducing, n-dopants that have been developed based on dimers of organic radicals or organometallic sandwich compounds.

doping↗
DOE OSTI2023

How Does Bending the Uranyl Unit Influence Its Spectroscopy and Luminescence?

We report bent uranyl complexes can be formed with chloride ligands and 1,10-phenanthroline (phen) ligands bound to the equatorial and axial planes of the uranyl(VI) moiety, as revealed by the crystal structures, IR and Raman spectroscopy, and quantum-chemical calculations. With the goal of probing the influence of chloride and phenanthroline coordination enforcing the bending on the absorption and emission spectra of this complex, spin–orbit time-dependent density functional theory calculations for the bare uranyl complexes as well as for the free UO 2 Cl 2 subunit and the UO 2 Cl 2 (phen) 2 complex were performed. The emission spectra have been fully simulated by ab initio methods and compared to experimental photoluminescence spectra, recorded for the first time for UO 2 Cl 2 (phen) 2 . Notably, the bending of uranyl in UO2Cl2 and UO 2 Cl 2 (phen) 2 triggers excitations of the uranyl bending mode, yielding a denser luminescence spectrum.

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

X-ray Spectroscopy Characterization of Electronic Structure and Metal–Metal Bonding in Dicobalt Complexes

Developing multimetallic complexes with tunable metal–metal interactions has long been a target of synthetic inorganic chemistry efforts due to the unique properties that such compounds can exhibit. However, understanding relationships between metal–metal bonding and chemical properties is challenging due to system-dependent factors that influence metal–metal and metal–ligand interactions, including ligand identity, coordination geometry, and metal–metal distance. In this work, we apply X-ray absorption and emission spectroscopy and quantum chemical calculations to describe electronic structure and bonding in a series of dicobalt complexes. The compounds with silane ligands and pseudo-octahedral coordination geometry exhibit Co–Co σ and multicentered bonding character, which we characterize from both the occupied and vacant perspectives via their contributions to the Co X-ray emission and absorption spectra, respectively. In contrast, the dicobalt complexes with a pseudotetrahedral coordination environment do not exhibit Co–Co bonding due to symmetry constraints on orbital overlap. We extend these insights to the theoretical evaluation of related dicobalt complexes to explain how ligand coordination and symmetry dictate the presence or absence of a Co–Co bond. In conclusion, this work highlights how fundamental insights into electronic structure and bonding through X-ray spectroscopy uncover important factors governing metal–metal interactions and guide the rational design of multimetallic complexes with tunable metal–metal bonds.

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

Eliminating Imaginary Vibrational Frequencies in Quantum-Chemical Cluster Models of Enzymatic Active Sites

In constructing finite models of enzyme active sites for quantum-chemical calculations, atoms at the periphery of the model must be constrained to stop unphysical rearrangements during geometry relaxation. A simple fixed-atom or “coordinate-lock” approach is commonly employed but leads to undesirable artifacts in the form of small imaginary frequencies. These preclude evaluation of finite-temperature free-energy corrections, limiting thermochemical calculations to enthalpies only. Full-dimensional vibrational frequency calculations are possible by replacing the fixed-atom constraints with harmonic confining potentials. Here, we compare that approach to an alternative strategy in which fixed-atom contributions to the Hessian are simply omitted. While the latter strategy does eliminate imaginary frequencies, it tends to underestimate both the zero-point energy and the vibrational entropy while introducing artificial rigidity. Harmonic confining potentials eliminate imaginary frequencies and provide a flexible means to construct active-site models that can be used in unconstrained geometry relaxations, affording better convergence of reaction energies and barrier heights with respect to the model size, as compared to models with fixed-atom constraints.

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

NCAP: Noncanonical Amino Acid Parameterization Software for CHARMM Potentials

Noncanonical Amino Acids (NCAAs) provide numerous avenues for introduction of novel functionality to peptides and proteins. NCAAs can be incorporated through solid phase synthesis or genetic code expansion in conjugation with heterologous expression of the encoded protein modification. Due to the difficulty of synthesis, wide chemical space and lack of empirically resolved structures modeling the effects of NCAA mutation is critical for rational protein design. To evaluate the structural and functional perturbations NCAAs introduce we utilize molecular potentials that describe the forces in protein structure. Most potentials such as CHARMM are designed to model canonical residues but can be parameterized in include novel NCAAs. Here, in this work, we introduce NCAP a software package to generate CHARMM compatible parameters from quantum chemical calculation. Unlike currently available tools NCAP is designed to recognize NCAA structure and automatically bridge the gap between DFT calculations and potential parameters. For our software we discuss workflow, validation against canonical parameter sets and comparison to published NCAA-protein structures.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

TeraChem Cloud: A High-Performance Computing Service for Scalable Distributed GPU-Accelerated Electronic Structure Calculations

The encapsulation and commoditization of electronic structure arise naturally as interoperability, and the use of nontraditional compute resources (e.g., new hardware accelerators, cloud computing) remains important for the computational chemistry community. Here, we present TERACHEM CLOUD, a high-performance computing service (HPCS) that offers on-demand electronic structure calculations on both traditional HPC clusters and cloud-based hardware. The framework is designed using off-the-shelf web technologies and containerization to be extremely scalable and portable. Within the HPCS model, users can quickly develop new methods and algorithms in an interactive environment on their laptop while allowing TERACHEM CLOUD to distribute ab initio calculations across all available resources. This approach greatly increases the accessibility of hardware accelerators such as graphics processing units (GPUs) and flexibility for the development of new methods as additional electronic structure packages are integrated into the framework as alternative backends. Cost-performance analysis indicates that traditional nodes are the most cost-effective long-term solution, but commercial cloud providers offer cutting-edge hardware with competitive rates for short-term large-scale calculations. We demonstrate the power of the TERACHEM CLOUD framework by carrying out several showcase calculations, including the generation of 300,000 density functional theory energy and gradient evaluations on medium-sized organic molecules and reproducing 300 fs of nonadiabatic dynamics on the B800-B850 antenna complex in LH2, with the latter demonstration using over 50 Tesla V100 GPUs in a commercial cloud environment in 8 h for approximately $1250.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗