Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Density functional 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 1,063 records · Page 59
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Exploring Structure-Sensitive Relations for Small Species Adsorption Using Machine Learning

Accurate prediction of adsorption energies on heterogeneous catalyst surfaces is crucial to predicting reactivity and screening materials. Adsorption linear scaling relations have been developed extensively but often lack accuracy and apply to one adsorbate and a single binding site type at a time. These facts undermine their ability to predict structure sensitivity and optimal catalyst structure. Using machine learning on nearly 300 density functional theory calculations, we demonstrate that generalized coordination number scaling relations hold well for oxygen- and high-valency carbon-binding species but fail for others. Here we reveal that the valency and the electronic coupling of a species with the surface, along with the site type and its coordination environment, are critical for small species adsorption. The model simultaneously predicts the adsorption energy and preferred site and significantly outperforms linear scalings in accuracy. It can expose the structure sensitivity of chemical reactions and enable enhanced catalyst activity via engineering particle shape and facet defects. The generality of our methodology is validated by training the model with transition metal data and transferring it to predict adsorption energies on single-atom alloys.

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

First Principles Nonadiabatic Excited-State Molecular Dynamics in NWChem

Computational simulation of non-adiabatic molecular dynamics is an indispensable tool for understanding complex photoinduced processes such as internal conversion, energy transfer, charge separation, and spatial localization of excitons, to name a few. Here, we report an implementation of the fewest-switches surface hopping algorithm in the NWChem computational chemistry program. Here the surface hopping method is combined with linear-response time-dependent density functional theory calculations of adiabatic excited state potential energy surfaces. To treat quantum transitions between arbitrary electronic Born{Oppenheimer states, we have implemented both numerical and analytical differentiation schemes for derivative non-adiabatic couplings. A numerical approach for the time-derivative non-adiabatic couplings together with an analytical method for calculating non-adiabatic coupling vectors is an efficient combination for surface hopping approaches. Additionally, electronic decoherence schemes and a state reassigned unavoided crossings algorithm are also implemented to improve the accuracy of the simulated dynamics and to handle trivial unavoided crossings. We apply our code to study the ultrafast decay of photoexcited benzene, including a detailed analysis of the potential energy surface, population decay time scales, and vibrational coordinates coupled to the excitation dynamics. The development presented in this work is a baseline for future implementations of more sophisticated frameworks for simulating of non-adiabatic molecular dynamics in NWChem.

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

Ab Initio Study of Vibronic and Magnetic 5f-to-5f and Dipole-Allowed 5f-to-6d and Charge-Transfer Transitions in [UX 6 ] n – (X = Cl, Br; n = 1, 2)

The absorption spectra of the octahedral [UX 6 ] n– (X = Cl, Br; n = 1, 2) complexes in the near-infrared (NIR) and UV–vis spectral regions were studied theoretically, using a relativistic restricted active space second-order perturbation theory (RASPT2) wavefunction framework, with the spin–orbit (SO) coupling treated by state interaction, in conjunction with Kohn–Sham density functional theory calculations for determining the vibrational normal modes. The electric-dipole-allowed ligand-to-metal charge-transfer (LMCT) and 5f-to-6d transitions, and the electric-dipole-forbidden 5f-to-5f ligand field (LF) transitions, are thereby obtained within the same theoretical framework. For the 5f-to-5f LF transitions, the observed absorption intensity is mostly due to vibronic coupling with low-energy electric-dipole-allowed transitions, but in some cases, the magnetic dipole intensity of the purely electronic transition has comparable intensity to the vibronic transitions. Here, experimental LF spectra of 5f 2 open-shell systems have been reported decades back, but ab initio calculations of their vibronic intensity have not yet been reported in the literature. Although the LF spectra for the 5f 2 systems can be assigned in detail, based on the calculations, the spectra are very complex and the underlying electronic states are strongly multiconfigurational. Therefore, the usefulness of the LF spectra beyond serving as a “fingerprint” of the LF and the metal oxidation state appears to be limited.

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

Beyond Single-Reference Fixed-Node Approximation in Ab Initio Diffusion Monte Carlo Using Antisymmetrized Geminal Power Applied to Systems with Hundreds of Electrons

Diffusion Monte Carlo (DMC) is an exact technique to project out the ground state (GS) of a Hamiltonian. Since the GS is always bosonic, in Fermionic systems, the projection needs to be carried out while imposing antisymmetric constraints, which is a nondeterministic polynomial hard problem. In practice, therefore, the application of DMC on electronic structure problems is made by employing the fixed-node (FN) approximation, consisting of performing DMC with the constraint of having a fixed, predefined nodal surface. How do we get the nodal surface? The typical approach, applied in systems having up to hundreds or even thousands of electrons, is to obtain the nodal surface from a preliminary mean-field approach (typically, a density functional theory calculation) used to obtain a single Slater determinant. This is known as single reference. In this paper, we propose a new approach, applicable to systems as large as the C 60 fullerene, which improves the nodes by going beyond the single reference. In practice, we employ an implicitly multireference ansatz (antisymmetrized geminal power wave function constraint with molecular orbitals), initialized on the preliminary mean-field approach, which is relaxed by optimizing a few parameters of the wave function determining the nodal surface by minimizing the FN-DMC energy. We highlight the improvements of the proposed approach over the standard single-reference method on several examples and, where feasible, the computational gain over the standard multireference ansatz, which makes the methods applicable to large systems. We also show that physical properties relying on relative energies, such as binding energies, are affordable and reliable within the proposed scheme.

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

Efficient Mixed-Precision Matrix Factorization of the Inverse Overlap Matrix in Electronic Structure Calculations with AI-Hardware and GPUs

In recent years, a new kind of accelerated hardware has gained popularity in the artificial intelligence (AI) community which enables extremely high-performance tensor contractions in reduced precision for deep neural network calculations. In this article, we exploit Nvidia Tensor cores, a prototypical example of such AI-hardware, to develop a mixed precision approach for computing a dense matrix factorization of the inverse overlap matrix in electronic structure theory, S –1 . This factorization of S –1 , written as ZZT = S –1 , is used to transform the general matrix eigenvalue problem into a standard matrix eigenvalue problem. Here we present a mixed precision iterative refinement algorithm where Z is given recursively using matrix–matrix multiplications and can be computed with high performance on Tensor cores. To understand the performance and accuracy of Tensor cores, comparisons are made to GPU-only implementations in single and double precision. Additionally, we propose a nonparametric stopping criteria which is robust in the face of lower precision floating point operations. The algorithm is particularly useful when we have a good initial guess to Z, for example, from previous time steps in quantum-mechanical molecular dynamics simulations or from a previous iteration in a geometry optimization.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Photophysical Properties and Electronic Structure of Zinc(II) Porphyrins Bearing 0–4 meso -Phenyl Substituents: Zinc Porphine to Zinc Tetraphenylporphyrin (ZnTPP)

Six zinc(II) porphyrins bearing 0–4 meso-phenyl substituents have been examined spectroscopically and theoretically. Comparisons with previously examined free base analogues afford a deep understanding of the electronic and photophysical effects of systematic addition of phenyl groups in porphyrins containing a central zinc(II) ion versus two hydrogen atoms. Trends in the wavelengths and relative intensities of the absorption bands are generally consistent with predictions from time-dependent density functional theory calculations and simulations from Gouterman’s four-orbital model. These trends derive from a preferential effect of the meso-phenyl groups to raise the energy of the highest occupied molecular orbital. The calculations reveal additional insights, such as a progressive increase in oscillator strength in the violet-red (B-Q) absorption manifold with increasing number of phenyls. Progressive addition of 0–4 phenyl substituents to the zinc porphyrins in O 2 -free toluene engenders a reduction in the measured lifetime of the lowest singlet excited state (2.5–2.1 ns), an increase in the S 1 → S o fluorescence yield (0.022–0.030), a decrease in the yield of S 1 → T 1 intersystem crossing (0.93–0.88), and an increase in the yield of S 1 → S o internal conversion (0.048–0.090). The derived rate constants for S 1 decay reveal significant differences in the photophysical properties of the zinc chelates versus free base forms. Here, the unexpected finding of a larger rate constant for internal conversion for zinc chelates versus free bases is particularly exemplary. Collectively, the findings afford fundamental insights into the photophysical properties and electronic structure of meso-phenylporphyrins, which are widely used as benchmarks for tetrapyrrole-based architectures in solar energy and life sciences research.

14 SOLAR ENERGY↗
DOE OSTI2021

Pulse Radiolysis Investigation of Radicals Derived from Water-Soluble Cyanine Dyes: Implications for Super-resolution Microscopy

Light-induced blinking, an inherent feature of many forms of super-resolution microscopy, has been linked to transient reduction of the fluorescent cyanine dye used as an imaging agent. There is, however, only scant literature information related to one-electron reduced cyanine dyes, especially in an aqueous environment. Here, we examine a small series of cyanine dyes, possessing disparate π-conjugation lengths, under selective reducing or oxidizing conditions. The experiment allows recording of both differential absorption spectra and decay kinetics of the resultant one-electron reduced or oxidized transient species in water. Relative to the ground state, absorption transitions for the various radicals are weak and somewhat broadened but do allow correlation with the π-conjugation length. In all cases, absorption maxima lie to the blue of the main ground-state transition. Under anaerobic conditions, the transient species decay on the microsecond to millisecond time scale, with the mean lifetime depending on molecular structure, radiation dose, and dye concentration. The experimental absorption spectra recorded for the one-electron reduced radicals and the presumed dimer cation radical compare well to spectra obtained from time-dependent density functional theory calculations. Furthermore, the results allow conclusions to be drawn regarding the plausibility of the reduced species being responsible for light-induced blinking in direct stochastic optical reconstruction microscopy.

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

Electron Binding Energy Spectra of Al n Pt – Clusters—A Combined Experimental and Computational Study

Results of size-selected electron photo-detachment experiments and density functional theory calculations on anionic Al n Pt – , n = 1–7, clusters are presented and analyzed. The measured and calculated spectra of electron binding energies are, overall, in excellent accord with each other. The analysis reveals the general importance of accounting for the multiplicity of structural forms of a given-size cluster that can contribute to its measured spectrum, especially when the clusters are fluxional and/or the conditions of the experiment allow for structural transitions. Here, we show that for the systems studied here, the size-specific peculiarities of the measured spectra can be understood in terms of the combined contributions of corresponding different accessible stable equilibrium conformations, bona-fide transition-state configurations, and electronic-crossing structures that may play the role of effective barriers in electronically nonadiabatic dynamics.

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

Density Matrix Implementation of the Fermi–Löwdin Orbital Self-Interaction Correction Method

The Fermi–Löwdin orbital self-interaction correction (FLOSIC) method effectively provides a transformation from canonical orbitals to localized Fermi–Löwdin orbitals which are used to remove the self-interaction error in the Perdew–Zunger (PZ) framework. This transformation is solely determined by a set of points in space, called Fermi–Löwdin descriptors (FODs), and the occupied canonical orbitals or the density matrix. In this work, we provide a detailed workflow for the implementation of the FLOSIC method for removal of self-interaction error in DFT calculations in an orbital-by-orbital basis that takes advantage of the unitary invariant nature of the FLOSIC method. In this way, it is possible to cast the self-consistent energy minimization at fixed FODs in the same manner than standard Kohn–Sham with one additional term in the Kohn–Sham Hamiltonian that introduces the PZ self-interaction correction. Each energy minimization iteration is divided in two substeps, one for the density matrix and one for the FODs. Expressions for the effective Kohn–Sham matrix and FOD gradients are provided such that its implementation is suitable for most electronic structure codes. Here, we analyze the convergence characteristics of the algorithm and present applications for the evaluation of NMR shielding constants and real-time time-dependent DFT simulations based on the Liouville–von Neumann equation to calculate excitation energies.

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

Competition between Water–Water Hydrogen Bonds and Water–π Bonds in Pyrene–Water Cluster Anions

We present infrared spectra and density functional theory calculations of hydrated pyrene anion clusters with up to four water molecules. The experimental spectra were acquired using infrared Ar messenger photodissociation spectroscopy. Water molecules form clusters on the surface of the pyrene, forming hydrogen bonds with the pi-system. The structures of the water clusters and their interaction with the π-system are encoded in the OH stretching vibrational modes. Here, we find that the interactions between water molecules are stronger than they are between water molecules and the π-system. While all clusters show multiple conformers, three- and four-membered rings are the lowest energy structures in the larger hydrates.

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

Univariate Prediction of Hammett Parameters and Select Relative Reaction Rates Using Loewdin Atomic Charges

Loewdin charges from density functional theory calculations were used here to obtain general, univariate linear correlations for the prediction of experimental Hammett parameters and relative reaction rates. While previous studies have established that Hirshfeld and CM5 charges perform strongly as univariate predictors, the near-ubiquitous Loewdin charges have not yet been evaluated. To this end, we assess the predictive capability of Loewdin charges for three chemical systems. First, we show that Loewdin charges outperform Hirshfeld and CM5 charges for Hammett parameter prediction. Second, we see that Loewdin charges generally perform comparably to Hirshfeld charges for predicting the relative rates of olefin cleavage by photoexcited nitroarenes. The single case of poor correlation, between relative rates and the Loewdin charges on nitrogen sites, is ameliorated when considering the net charge on the NO 2 group. Third, we show that Loewdin, Hirshfeld, and CM5 charges all perform very well for generating correlations for relative reaction rates for C–H activation of 9-(4-X-phenyl)-9H-fluorene substrates by a transition metal catalyst. The equations generated throughout the study enable the prediction of Hammett parameters and relative reaction rates. Finally, these tools can accelerate synthetic and experimental studies by enabling the in silico prediction of uncharacterized chemical properties.

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

Energetics of Base–Acid Pairs for the Design of High-Temperature Fuel Cell Polymer Electrolytes

The interaction energy of base–acid plays a key role in acid retention of phosphoric acid (PA)-doped polymer electrolytes under fuel cell operating conditions. Here, we investigate the energetics of proton-accepting and hydroxide-donating organic bases using density functional theory calculations. Because of their weak basicity, proton-accepting organic bases such as benzimidazole have relatively low interaction energy with the acid in the absence of water (15.3–28.0 kcal mol –1 ). Energetics of the proton-accepting base–PA complex increases by adding water, indicating that the interactions in the base–acid complex strengthen in the presence of water. On the other hand, hydroxide-donating organic bases, such as tetramethylammonium hydroxide, have high interaction energy with PA (~110 kcal mol –1 ), which remains high in the presence of water. The chemical shifts of 31P NMR support the energetics of the base–acid complexes. Overall, this study further discusses the benefit of incorporating hydroxide-donating organic bases into the polymeric structure over proton-accepting bases as a way to increase acid retention.

25 ENERGY STORAGE↗
DOE OSTI2023

Red-Colored Circularly Polarized Luminescence from a Benzo-Fused BODIPY-BINOL Complex

Red emission with sharp bandwidth and high quantum yield are desired characteristics for organic chromophores in optoelectronic, spintronic, and biomedical applications. Here, we observe circularly polarized luminescence (CPL) with these characteristics from a benzo-fused BODIPY-BINOL complex (1). Using time-resolved optical spectroscopy, electrochemistry, and density functional theory calculations, we showed that the emissive excited state of 1 does not have a charge-transfer character, unlike the regular BODIPY counterpart (2). In conclusion, the rigidity and the lack of charge-transfer character make this class of molecules an appealing platform for CPL-active molecules in the red spectral region with ample room for improvement in the dissymmetry factor and brightness.

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

Intermolecular Proton Transfer Enabled Reactive CO 2 Capture by the Malononitrile Anion

Task-specific ionic liquids (ILs) employing carbanions represent a new class of ILs for carbon capture. The deprotonated malononitrile carbanion, [CH(CN) 2 ] - , has shown close to equimolar capacity for reactive CO 2 capture. Although the formation of the [C(CN) 2 COOH] - carboxylic acid was found to be the final product, how the hydrogen atom on the [CH(CN) 2 ] - carbanion transfers to the carboxylate group as a proton has not been fully understood. In this work, we employ density functional theory calculations with an implicit solvation model to investigate the proton transfer mechanisms in forming carboxylic acid from the reaction of the [CH(CN) 2 ] - carbanion with CO 2 . We find that the intramolecular proton-transfer pathway in [CH(CN) 2 COO] - to form [C(CN) 2 COOH] - is unlikely due to the high energy barrier of 152 kJ/mol. Instead, the intermolecular proton transfer pathway between two [CH(CN) 2 COO] - anions is more feasible to form two molecules of [C(CN) 2 COOH] - , with a significantly lower activation energy of 50 kJ/mol. Moreover, the [C(CN) 2 COOH] - dimer is further stabilized by the intermolecular hydrogen bonds of the two –COOH groups in the Z-configuration of the π-conjugated planar geometry. This insight of reactive CO 2 capture enabled by intermolecular proton transfer will be useful in designing novel carbanions and ILs for carbon capture and conversion.

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

Ab Initio Molecular Dynamics Study of the Reduction of Acetone by the Hydrated Electron

We have investigated the reaction dynamics of the reduction of acetone by the solvated electron in water using ab initio molecular dynamics simulations at 298 and 373 K. The rate constants derived from the simulations are consistent with experimental observations that the reaction has a low activation energy. Detailed analyses carried out to shed light on the mechanism of the reaction show that solvent reorganization plays a key role as a reaction coordinate, as expected from the Marcus electron transfer theory. Furthermore, constrained density functional theory calculations indicate that the electronic coupling is large, placing the reaction in the adiabatic limit. Indeed, the activation energies and rate constants of the simulations are in accordance with the predictions of adiabatic Marcus theory.

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

Real-Space Local Dynamics in 1,2,3-Triazole Using Inelastic Neutron Scattering

Enhancing proton transport in polymer electrolytes is crucial for advancing next-generation solid-state batteries, yet our understanding of proton conductivity in nonaqueous environments remains limited due to a lack of atomic-scale insights. Here, in this study, we investigated the atomic-scale dynamics of 1,2,3-triazole, a small molecule capable of dynamic hydrogen bonding, as a model system for proton hopping in nonaqueous environments. Using the real-space correlation function determined by the double Fourier transformation of inelastic neutron scattering spectra, we identified that the self-motion of protons and intermolecular dynamics occur on comparable time scales. Furthermore, we observed that the activation energy associated with the intermolecular dynamics matches the energy barrier for molecular rotations determined through Density Functional Theory calculations. These findings underscore the importance of controlling molecular dynamics at the atomic scale to control proton transport. Additionally, we demonstrated that intermolecular dynamics in systems involving protons can be studied using inelastic neutron scattering even without deuteration, thereby providing a broader avenue for studying atomic-scale dynamics in soft matter systems.

Shinohara, Yuya [Oak Ridge National Laboratory (OR↗
DOE OSTI2020

Predicting Elastic Properties of Materials from Electronic Charge Density Using 3D Deep Convolutional Neural Networks

Materials representation plays a key role in machine learning-based prediction of materials properties and new materials discovery. Currently both graph and three-dimensional (3D) voxel representation methods are based on the heterogeneous elements of the crystal structures. Here, we propose to use electronic charge density (ECD) as a generic unified 3D descriptor for materials property prediction with the advantage of possessing close relation with the physical and chemical properties of materials. We developed an ECD-based 3D convolutional neural networks (CNNs) for predicting the elastic properties of materials, in which CNNs can learn effective hierarchical features with multiple convolving and pooling operations. Extensive benchmark experiments over 2170 $Fm\bar3m$ face-centered-cubic materials show that our ECD-based CNNs can achieve good performance for elasticity prediction. Especially, our CNN models based on the fusion of elemental Materials-Agnostic Platform for Informatics and Exploration features and ECD descriptors achieved the best fivefold cross-validation performance. More importantly, we showed that our ECD-based CNN models can achieve significantly better extrapolation performance when evaluated over nonredundant data sets, where there are few neighbor-training samples around test samples. As an additional validation, we evaluated the predictive performance of our models on 329 materials of space group $Fm\bar3m$ by comparing to density functional theory calculated values, which shows a better prediction power of our model for bulk modulus than shear modulus. Because of the unified representation power of ECD, it is expected that our ECD-based CNN approach can also be applied to predict other physical and chemical properties of crystalline materials.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Nanoscale Mg–B via Surfactant Ball Milling of MgB 2 : Morphology, Composition, and Improved Hydrogen Storage Properties

Metal borides have attracted the attention of researchers due to their useful physical properties and unique ability to form high hydrogen-capacity metal borohydrides. We demonstrate improved hydrogen storage properties of a nanoscale Mg–B material made by surfactant ball milling MgB 2 in a mixture of heptane, oleic acid, and oleylamine. Transmission electron microscopy data show that Mg–B nanoplatelets are produced with sizes ranging from 5 to 50 nm, which agglomerate upon ethanol washing to produce an agglomerated nanoscale Mg–B material of micron-sized particles with some surfactant still remaining. X-ray diffraction measurements reveal a two-component material where 32% of the solid is a strained crystalline solid maintaining the hexagonal structure with the remainder being amorphous. Fourier transform infrared shows that the oleate binds in a “bridge-bonding” fashion preferentially to magnesium rather than boron, which is confirmed by density functional theory calculations. The Mg–B nanoscale material is deficient in boron relative to bulk MgB 2 with a Mg–B ratio of ~1:0.75. The nanoscale MgB 0.75 material has a disrupted B–B ring network as indicated by X-ray absorption measurements. Hydrogenation experiments at 700 bar and 280 °C show that it partially hydrogenates at temperatures 100 °C below the threshold for bulk MgB 2 hydrogenation. In addition, upon heating to 200 °C, the H–H bond-breaking ability increases ~10-fold according to hydrogen–deuterium exchange experiments due to desorption of oleate at the surface. This behavior would make the nanoscale Mg–B material useful as an additive where rapid H–H bond breaking is needed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗