Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “calculator”

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

Accurate calculation of many-body energies in water clusters using a classical geometry-dependent induction model

We incorporate geometry-dependent distributed multipole and polarizability surfaces into an induction model that is used to describe the 3- and 4-body terms of the interaction between water molecules. The expansion is carried out up to hexadecapole with the multipoles distributed on the atom sites. Dipole-dipole, dipole-quadrupole, and quadrupole-quadrupole distributed polarizabilities are used to represent the response of the multipoles to an electric field. We compare the model against two large databases consisting of 43,844 3-body terms and 3,603 4-body terms obtained from high level ab initio calculations previously used to fit the MB-pol and q-AQUA interaction potentials. The classical induction model with no adjustable parameters reproduces the ab-initio 3- and 4-body terms contained in these teo Databases with a Root-Mean-Square-Error (RMSE) of 0.104/0.058 and a Mean-Absolute-Error (MAE) of 0.054/0.026 kcal/mol, respectively, results that are on a par with those obtained 1 by fitting the same data using tens of thousands of Permutationally Invariant Polynomials (PIPs). This demonstrates the accuracy of this physically motivated model in describing the 3- and 4-body terms in the interactions between water molecules with no adjustable parameters. The triple-dipole-dispersion energy was included in the 3-body energy and was found to be small but not quite negligible. The model represents a practical, efficient and transferable approach for obtaining accurate non-additive interactions for multi-component systems without the need of performing tens of thousands of high level electronic structure calculations and fitting them with tens of thousands of PIPs.

molecular analysis water, aqueous↗
DOE OSTI2024

Energetic and Spectroscopic Properties of Astrophysically Relevant MgC 4 H Radicals Using High-Level Ab Initio Calculations

Considering the importance of magnesium-bearing hydrocarbon molecules (MgC n H; n = 2, 4, and 6) in the carbon-rich circumstellar envelopes (e.g., IRC+10216), a total of 28 constitutional isomers of MgC 4 H have been theoretically investigated using density functional theory (DFT) and coupled-cluster methods. The zero-point vibrational energy corrected relative energies at the ROCCSD(T)/cc-pCVTZ level of theory reveal that the linear isomer, 1-magnesapent-2,4-diyn-1-yl (1, 2 Σ + ), is the global minimum geometry on the MgC 4 H potential energy surface. The latter has been detected both in the laboratory and in the evolved carbon star, IRC+10216. The calculated spectroscopic data for 1 match well with the experimental observations (error ~ 0.78%) which validates our theoretical methodology. Plausible isomerization processes happening among different isomers are examined using DFT and coupled-cluster methods. CASPT2 calculations have been performed for a few isomers exhibiting multireference characteristics. The second most stable isomer, 1-ethynyl-1λ 3 -magnesacycloprop-2-ene-2,3-diyl (2, 2 A 1 , μ = 2.54 D), is 146 kJ mol –1 higher in energy than 1 and possibly the next promising candidate to be detected in the laboratory or in the interstellar medium in future.

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

Exploring Anomalous Photoelectron Angular Distributions in the Photoelectron Spectra of Gd 3 O 3 – : Study of Gd 3 O 2 – and Gd 3 O 3 – Using Photoelectron Spectroscopy and Density Functional Theory Calculations

Anion photoelectron (PE) spectra of lanthanide oxide clusters obtained previously have exhibited anomalous photoelectron angular distributions which were attributed to strong PE–valence electron (PEVE) interactions. Here, to further explore this effect, we have obtained the PE spectra of Gd 3 O 2 – and Gd 3 O 3 – , two clusters that have similarly complex electronic structures but contrasting symmetries. The spectra exhibit manifolds of detachment transitions at similar binding energies in a 0.5 eV window of energy. The electron affinity of Gd 3 O 2 is measured to be 1.29 ± 0.05 eV, and that of Gd 3 O 3 is 1.31 ± 0.05 eV. As seen in previous studies on lanthanide oxide cluster anions in lower than conventional oxidation states, transitions in spectra obtained lower photon energies are more congested than those obtained with higher photon energy, a signature of strong PEVE interactions. While the detachment transitions have predominantly parallel photoelectron angular distributions (PAD), the PAD varies across the manifold of transitions in the PE spectrum of Gd 3 O 3 – in a way that suggests four different subgroups of transitions. Results of calculations on Gd 3 O 2 – suggest kite or V-shape structures with antiferromagnetic coupling between one of the 4f 7 subshells with the two others. Calculations on Gd 3 O 3 – more definitively point to ring structures with a nearly isoenergetic ferromagnetically coupled high spin (24-tet) state and a dectet state in which one of the 4f 7 subshells is antiferromagnetically coupled with the other two. Taking these results as qualitative, we propose that strong mixing between the unperturbed states predicted computationally leads to overlapping transitions with different PADs.

anions↗
DOE OSTI2023

Phonons in Copper Diphosphide (CuP 2 ): Raman Spectroscopy and Lattice Dynamics Calculations

Copper diphosphide (CuP 2 ) is an emerging binary semiconductor with promising properties for energy conversion and storage applications. While functionality and possible applications of CuP 2 have been studied, there is a curious gap in the investigation of its vibrational properties. In this work, we provide a reference Raman spectrum of CuP 2 , with a complete analysis of all Raman active modes from both experimental and theoretical perspectives. Raman measurements have been performed on polycrystalline CuP 2 thin films with close to stoichiometric composition. Detailed deconvolution of the Raman spectrum with Lorentzian curves has allowed identification of all theoretically predicted Raman active modes (9A g and 9B g ), including their positions and symmetry assignment. Furthermore, calculations of the phonon density of states (PDOS), as well as the phonon dispersions, provide a microscopic understanding of the experimentally observed phonon lines, in addition to the assignment to the specific lattice eigenmodes. We further provide the theoretically predicted positions of the infrared (IR) active modes, along with the simulated IR spectrum from density functional theory (DFT). Overall good agreement is found between the experimental and DFT-calculated Raman spectra of CuP 2 , providing a reference platform for future investigations on this material.

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

Structural Stability of Graphene-Supported Pt Layers: Diffusion Monte Carlo and Density Functional Theory Calculations

In this study, we used a combination of diffusion Monte Carlo and density functional theory calculations to investigate the stability and interlayer binding of various layered structures of Pt atoms adsorbed on graphene. Our findings show that vertically buckled Pt monolayer and bilayer with (111)-packing order are more energetically favorable than the corresponding buckled (100) or flat (100)-packing structures. This can be attributed to the significant lattice mismatch (>10%) between pristine graphene and a free-standing (100)-packing Pt layer. Additionally, our calculations reveal that among the (100)-packing Pt layers, an incommensurate structure with a Pt/C atomic ratio less than 1/2 may be more stable than the commensurate structures registered at the bridge sites, which aligns with recent experimental findings of incommensurate (100)-packing Pt layers on graphene. The interlayer binding between the Pt layer and graphene is found to be primarily driven by van der Waals interaction, except for the AA-stacked buckled-(100) Pt bilayer where the bottom Pt atoms show chemisorption to the graphene surface. In conclusion, this research offers a comprehensive examination of the stability and interlayer binding of metallic Pt layers, providing valuable insights for potential applications as a next-generation catalyst.

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

Understanding Pore Filling Processes and Adsorption/Desorption Hysteresis in Nanoporous Metal–Organic Frameworks: Insights from Grand Canonical Monte Carlo Simulations and Free Energy Calculations

Grand canonical Monte Carlo (GCMC) simulations were used to investigate pore filling and hysteresis in nanoporous metal-organic frameworks (MOFs). Adsorption and desorption isotherms were calculated for argon at 87 K in 1866 MOFs from the CoRE MOF database and for short n-alkanes in selected MOFs, keeping the adsorbent structure rigid. Analysis of the molecular configurations showed two different mechanisms and origins of hysteresis: one involving a transition of the adsorbate arrangement in the pores similar to a gas-to-liquid transition associated with a large change in the loading and one more similar to a liquid-to-solid transition associated with a relatively small change in the loading. Our GCMC simulations in MOFs with diverse pore topologies indicate exceptions to an empirical relationship for the minimum diameter of a cylindical pore required for hysteresis as a function of the adsorbate diameter and reduced temperature. The simulations reveal some structures where isotherms exhibit two steps in the adsorption branch and only one step in the desorption branch. Hysteresis loops with a different number of adsorption and desorption steps are not common. Here, to better understand why hysteresis is observed in the GCMC simulations, the concept of the transition probability for observing a step in the adsorption isotherm at a given pressure in a GCMC simulation is introduced. We used two different methods to calculate the transition probabilities and find that these yield comparable results. Furthermore, the transition probability provides a measure for the length of GCMC simulations to yield reliable results.

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

Predicting Device Parameters for Dye-Sensitized Solar Cells from Electronic Structure Calculations to Reproduce Experiment.

Given that improvements to the power-conversion efficiency (PCE) of dye-sensitized solar cells (DSSCs) have slowed in recent years, a means to accurately predict the device parameters yielded by trial dyes in silico, without having to synthesize them, would be extremely valuable to speed up the design process. Currently, the best-performing methods of calculating device parameters rely on a set of experimentally determined kinetic coefficients. In practice, it is very difficult to measure these kinetic parameters accurately, limiting the overall accuracy of such predictive methods. This work proposes a model to obtain key parameters such as J(SC), V-OC, and PCE using only the results from density functional theory (DFT) and time-dependent DFT calculations, noting that rates of electron-transfer steps are ultimately linked to the electronic structure of the dye center dot center dot center dot TiO2 working electrode. Six organic DSSC dyes from dissimilar chemical classes (L0, L1, L2, WS-2, WS-92, and C281) were chosen to demonstrate the power of this approach. Their a priori known experimentally determined device performance metrics served to validate our predictions. The greatest absolute error in our predicted PCE values was 0.36% relative to the experiment, while the greatest fractional error was 0.042. This indicates that the proposed model offers a dramatic improvement on previous predictive methods for DSSC device parameters, both in accuracy and in consistency. Moreover, such a predictive model has great potential to be applied to other photovoltaic applications, further enabling the design of novel, highly efficient photoactive materials.

density functional theory↗
DOE OSTI2022

Mechanistic Study of 1,2-Dichloroethane Hydrodechlorination on Cu-Rich Pt–Cu Alloys: Combining Reaction Kinetics Experiments with DFT Calculations and Microkinetic Modeling

Cu-rich Pt–Cu bimetallic catalysts are among the most promising candidates for actively catalyzing the hydrodechlorination of 1,2-dichloroethane (1,2-DCA) toward ethylene production. Combining reaction kinetics experiments with density functional theory (DFT) calculations and mean-field microkinetic modeling, we present a systematic mechanistic study for 1,2-DCA hydrodechlorination on Cu-rich Pt–Cu alloy catalysts. Our DFT (PBE+(TS+SCS)) results suggest that increasing Cu content in the Pt–Cu alloy destabilizes C 2 -species adsorption while stabilizing the binding of atomic chlorine. The reaction energetics of all the elementary steps in the 1,2-DCA reaction network were calculated on a Pt 1 Cu 3 (111) model surface. The DFT results were then used to construct a microkinetic model, and the model-predicted reaction rates were compared with our reaction kinetics experimental results on a Cu-rich SiO 2 -supported Pt–Cu alloy catalyst through a parameter estimation procedure. Both the reaction kinetics experiments and the microkinetic model after parameter adjustments yielded 100% selectivity to ethylene. The microkinetic model pointed to a reaction pathway involving two sequential chlorine-removal steps on the Pt–Cu alloy catalyst, a mechanism distinct from the one previously identified on pure Pt/SiO 2 catalysts, which involved an initial hydrogen-removal step. Adjustments to the DFT-derived parameters indicate the possible formation of chlorine-induced Cu-enriched surface sites during 1,2-DCA hydrodechlorination conditions, sites that are more active than those encountered in the bulk Pt 1 Cu 3 (111) alloy surface. Furthermore, our study offers valuable initial insights on the 1,2-DCA hydrodechlorination reaction mechanism and the nature of the active sites on PtCu bimetallic catalysts.

1,2-Dichloroethane↗
DOE OSTI2025

Mg-Ion Conduction in Antiperovskite Solid Electrolytes Revealed by 25 Mg Ultrahigh Field NMR and First-Principles Calculations

Magnesium-ion batteries hold the potential to outperform the energy density of lithium-ion batteries, given the divalent charge carried by each Mg 2+ cation, but remain in an early stage of development. Here, in this study, 25 Mg solid-state nuclear magnetic resonance (ssNMR) is used to gain insight into the local structure and Mg-ion dynamics of candidate Mg-ion solid electrolytes, the antiperovskites Mg 3 SbN and Mg 3 AsN. Using the highest available magnetic field (35.2 T) for high-resolution solid-state NMR, the largest 25 Mg quadrupole coupling constants (C Q ) yet measured of up to 22 MHz are reported and corroborated by first-principles calculations. Predicted C Q values are shown to correlate with the antiperovskite’s tolerance factor; thus, 25 Mg NMR linewidths can report on lattice distortions and phase stability of these antiperovskites. Variable-temperature 25 Mg NMR spectra demonstrate changes at elevated temperatures, ascribed to Mg-ion motional effects. 25 Mg T 1 relaxometry measurements at ultrahigh field reveal a lower activation energy for the more distorted Mg 3 AsN phase, matching computational predictions of a lower energy barrier for Mg 2+ ion migration and suggesting that additional scrutiny of antiperovskites as Mg-ion conductors is warranted. Given the inherent challenges of 25 Mg NMR, this work demonstrates the benefits of combining ultrahigh field NMR spectroscopy, advanced pulse sequences, modern signal processing, and first-principles calculations to facilitate NMR of quadrupolar nuclei as a tool to probe the local structure and ion dynamics in beyond-Li battery materials.

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

fluxfinder: An R Package for Reproducible Calculation and Initial Processing of Greenhouse Gas Fluxes From Static Chamber Measurements

Fluxes of greenhouse gases are a critical component of the earth's natural climate, but anthropogenic emissions have created an imbalance and resulted in global climate change. Quantifying the emission of these gases is vital to our understanding of their sources and sinks, both natural and anthropogenic. The static chamber method, in which a system of interest is enclosed, and gas concentrations are measured over time, is widely used to estimate fluxes of greenhouse gases. With the development of instruments such as infrared gas analyzers (IRGAs) supporting high-frequency concentration data, there is a growing need for open-source workflows to calculate fluxes. Here we present fluxfinder, an R package designed to support reproducible calculations and processing of greenhouse gas fluxes measured with the static chamber method. The package includes raw data file parsing from widely used IRGAs, metadata matching, unit conversion, flux estimations, and initial quality assurance/quality control (QA/QC). Diagnostic graphical plots provide a transparent way to differentiate between measurement issues and nonlinear behavior. The package is also designed to be easily integrated with the gasfluxes package for further fitting of nonlinear concentration-time models, allowing alternative or additional flux QA/QC. The fluxfinder package offers a flexible workflow that is easily adaptable to promote open and reproducible greenhouse gas flux estimations.

Wilson, Stephanie J.↗
DOE OSTI2021

Efficient calculation of carrier scattering rates from first principles

Abstract The electronic transport behaviour of materials determines their suitability for technological applications. We develop a computationally efficient method for calculating carrier scattering rates of solid-state semiconductors and insulators from first principles inputs. The present method extends existing polar and non-polar electron-phonon coupling, ionized impurity, and piezoelectric scattering mechanisms formulated for isotropic band structures to support highly anisotropic materials. We test the formalism by calculating the electronic transport properties of 23 semiconductors, including the large 48 atom CH 3 NH 3 PbI 3 hybrid perovskite, and comparing the results against experimental measurements and more detailed scattering simulations. The Spearman rank coefficient of mobility against experiment ( r s = 0.93) improves significantly on results obtained using a constant relaxation time approximation ( r s = 0.52). We find our approach offers similar accuracy to state-of-the art methods at approximately 1/500th the computational cost, thus enabling its use in high-throughput computational workflows for the accurate screening of carrier mobilities, lifetimes, and thermoelectric power.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Automated calculation and convergence of defect transport tensors

Defect diffusion is a key process in materials science and catalysis, but as migration mechanisms are often too complex to enumerate a priori, calculation of transport tensors typically have no measure of convergence and require significant end-user intervention. These two bottlenecks prevent high-throughput implementations essential to propagate model-form uncertainty from interatomic interactions to predictive simulations. In order to address these issues, we extend a massively parallel accelerated sampling scheme, autonomously controlled by Bayesian estimators of statewide sampling completeness, to build atomistic kinetic Monte Carlo models on a state-space irreducible under exchange and space group symmetries. Focusing on isolated defects, we derive analytic expressions for drift and diffusion coefficients, providing a convergence metric by calculating the Kullback–Leibler divergence across the ensemble of diffusion processes consistent with the sampling uncertainty. The autonomy and efficacy of the method is demonstrated on surface trimers in tungsten and Hexa-interstitials in magnesium oxide, both of which exhibit complex, correlated migration mechanisms.

36 MATERIALS SCIENCE↗
DOE OSTI2023

pyGWBSE: a high throughput workflow package for GW-BSE calculations

Abstract We develop an open-source python workflow package, py GWBSE to perform automated first-principles calculations within the GW-BSE (Bethe-Salpeter) framework. GW-BSE is a many body perturbation theory based approach to explore the quasiparticle (QP) and excitonic properties of materials. GW approximation accurately predicts bandgaps of materials by overcoming the bandgap underestimation issue of the more widely used density functional theory (DFT). BSE formalism produces absorption spectra directly comparable with experimental observations. py GWBSE package achieves complete automation of the entire multi-step GW-BSE computation, including the convergence tests of several parameters that are crucial for the accuracy of these calculations. py GWBSE is integrated with Wannier90 , to generate QP bandstructures, interpolated using the maximally-localized wannier functions. py GWBSE also enables the automated creation of databases of metadata and data, including QP and excitonic properties, which can be extremely useful for future material discovery studies in the field of ultra-wide bandgap semiconductors, electronics, photovoltaics, and photocatalysis.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Probing radical–molecule interactions with a second generation energy decomposition analysis of DFT calculations using absolutely localized molecular orbitals

Intermolecular interactions between radicals and closed-shell molecules are ubiquitous in chemical processes, ranging from the benchtop to the atmosphere and extraterrestrial space. While energy decomposition analysis (EDA) schemes for closed-shell molecules can be generalized for studying radical–molecule interactions, they face challenges arising from the unique characteristics of the electronic structure of open-shell species. In this work, we introduce additional steps that are necessary for the proper treatment of radical–molecule interactions to our previously developed unrestricted Absolutely Localized Molecular Orbital (uALMO)-EDA based on density functional theory calculations. A “polarize-then-depolarize” (PtD) scheme is used to remove arbitrariness in the definition of the frozen wavefunction, rendering the ALMO-EDA results independent of the orientation of the unpaired electron obtained from isolated fragment calculations. The contribution of radical rehybridization to polarization energies is evaluated. It is also valuable to monitor the wavefunction stability of each intermediate state, as well as their associated spin density profiles, to ensure the EDA results correspond to a desired electronic state. The examples suggest that this second generation uALMO-EDA is a useful tool for furthering one's understanding of both energetic and property changes associated with radical–molecule interactions.

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

First-principles calculations for understanding microstructures and mechanical properties of co-sputtered Al alloys

Recent experimental studies show that co-sputtering solutes with Al, together, can refine columnar grain size around few tens of nanometers and promote the formation and enhance the stability of planar defects such as stacking faults (SFs) and grain boundaries (GBs) in Al alloys. These crystal defects and fine columnar grains result in high strength, enhanced strain hardening and thermal stability of Al alloys. Using first-principles density-functional theory (DFT) calculations, we studied the role of eleven solutes in tailoring kinetics and energetics of adatoms and clusters on Al {111} surface, stable and unstable stacking fault energies, and kinetic energy barriers for the migration of defects. The calculations show that most solutes can effectively refine columnar grain size by decreasing the diffusivity of adatoms and surface clusters. These solutes do not necessarily decrease the stacking fault energy of Al alloys, but reduce the formation energy of faulted surface clusters and increase the energy barriers for the recovery of faulted surface clusters. Correspondingly, the formation of SFs is kinetically promoted during sputtering. Furthermore, solutes are segregated into the core of Shockley partial dislocations and play a pinning effect on SFs, SF arrays and twin boundaries, enhancing the thermal stability of these crystal defects. These findings provide insights into the design of high-strength Al alloys for high-temperature applications.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Informing air–carbon ablation modeling with theoretical calculations of atomic oxygen and nitrogen interacting with carbon surfaces

To understand the gas–surface chemistry above the thermal protection system of a hypersonic vehicle, it is necessary to map out the kinetics of key elementary reaction steps. In this work, extensive periodic density functional theory (DFT) calculations are performed to elucidate the interaction of atomic oxygen and nitrogen with both the basal plane and edge sites of highly oriented pyrolytic graphite (HOPG). Reaction energies and barriers are determined for adsorption, desorption, diffusion, recombination, and several reactions. These DFT results are compared with the most recent finite-rate model for air–carbon ablation. Our DFT results corroborated some of the parameters used in the model but suggest that further refinement may be necessary for others. The calculations reported here will help to establish a predictive kinetic model for the complex reaction network present under hypersonic flight conditions.

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

Relativistic coupled-cluster calculations of RaOH pertinent to spectroscopic detection and laser cooling

Here a relativistic coupled-cluster study of the low-lying electronic states in the radium monohydroxide molecule (RaOH), a radioactive polyatomic molecule of interest to laser cooling and to the search of new physics beyond the Standard Model, is reported. The level positions of the A 2 Π 1/2 and C 2 Σ states have been computed with an accuracy of around 200 cm –1 to facilitate spectroscopic observation of RaOH using laser induced fluorescence spectroscopy, thereby exploiting the systematic convergence of electron-correlation and basis-set effects in relativistic coupled-cluster calculations. The energy level for the B 2 Δ 3/2 state has also been calculated accurately to conclude that the B 2 Δ 3/2 state lies above the A 2 Π 1/2 state. This confirms X 2 Σ ↔ A 2 Π 1/2 as a promising optical cycling transition for laser cooling RaOH.

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

Refining the thermochemical properties of CF, SiF, and their cations by combining photoelectron spectroscopy, quantum chemical calculations, and the Active Thermochemical Tables approach

Fluorinated species have a pivotal role in semiconductor material chemistry and some of them have been detected beyond the Earth's atmosphere. Achieving good energy accuracy on fluorinated species using quantum chemical calculations has long been a challenge. In addition, obtaining direct experimental thermochemical quantities has also proved difficult. Here, we report the threshold photoelectron and photoion yield spectra of SiF and CF radicals generated with a fluorine reactor. The spectra were analysed with the support of ab initio calculations, resulting in new experimental values for the adiabatic ionisation energies of both CF (9.128 ± 0.006 eV) and SiF (7.379 ± 0.009 eV). Using these values, the underlying thermochemical network of Active Thermochemical Tables was updated, providing further refined enthalpies of formation and dissociation energies of CF, SiF, and their cationic counterparts.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗