Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “density functional theory (DFT)”

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 181 records · Page 10

Experimental and DFT Calculated IR Spectra of Guests in Zeolites: Acyclic Olefins and Host-Guest Interactions

We performed experimental and periodic density functional theory (DFT) IR spectroscopy to investigate the adsorption of acyclic olefins over both acidic and non-acidic zeolites. Two conjugated polyenes, 2,4-dimethyl-1,3-pentadiene (I) and 2,6-dimethyl-2,4,6-octatriene (II) were studied to probe organic intermediates that can form during methanol conversion, and that can lead to deactivating species known collectively as “coke.” We computed vibrational spectra using zeolite-adsorbed and gas-phase models for both neutral and protonated forms of I and II, and compared these DFT results to diffuse reflectance IR Fourier transform (DRIFT) spectra of zeolite-guest systems. Our experimental and computational results are precise enough to pinpoint that the gauche s-cis conformation of species I predominates during adsorption over de-aluminated zeolite beta. Computed zeolite-adsorbed spectra of the protonated species I and II best represent the DRIFT spectra obtained after the adsorption of the olefins on HMOR at 20 ?, with computed bands at 1543 and 1562 cm–1 for molecules I+ and II+, respectively, attributed to the allylic stretching mode, ?(C=C–C+). These computed band frequencies are within 6 cm–1 of experimental data and confirm that the interaction between neutral acyclic olefins and acidic zeolites leads to protonation of the olefin. Comparison of computed spectra of the protonated species in the gas phase to those in the zeolite indicates that the electrostatic interaction between alkenyl and alkadienyl cations and negative zeolite framework does not significantly impact the position of the allylic stretching bands. These results highlight that computed spectroscopy and thermodynamics coupled with experimental spectra can be used for elucidating complex mixtures in zeolites, and certain spectral features of adsorbed olefins can be accurately modeled by gas-phase calculations. ACKNOWLEDGMENTS BM and SMA were supported by the National Science Foundation under Award # CBET-1512442. Acknowledgment is made to the Donors of the American Chemical Society Petroleum Research Fund for support (or partial support) of this research (#589978-ND5, EDH and FCJ). MDB was supported by the US Department of Energy (DOE), Office of Science, Office of Basic Energy Sciences (BES), Division of Materials Science and Engineering. CJM was supported by the DOE, Office of Science, BES, Division of Chemical Sciences, Geosciences, and Biosciences. Pacific Northwest National Laboratory (PNNL) is a multiprogram national laboratory operated by Battelle for the U.S. DOE. Computing resources were generously provided through PNNL’s institutional computing (PIC).

Manookian, Babgen↗

Conventional High-Temperature Superconductivity in Metallic, Covalently Bonded, Binary-Guest C–B Clathrates

Inspired by the synthesis of XB 3 C 3 (X = Sr, La) compounds in the bipartite sodalite clathrate structure, density functional theory (DFT) calculations are performed on members of this family containing up to two different metal atoms. A DFT chemical pressure analysis on systems with X = Mg, Ca, Sr, Ba reveals that the size of the metal cation, which can be tuned to stabilize the B–C framework, is key for their ambient-pressure dynamic stability. High-throughput density functional theory calculations on 105 Pm$\bar{3}$ symmetry XYB 6 C 6 binary-guest compounds (where X, Y are electropositive metal atoms) find 22 that are dynamically stable at 1 atm, expanding the number of potentially synthesizable phases by 19 (18 metals and 1 insulator). Here, the density of states at the Fermi level and superconducting critical temperature, T c , can be tuned by changing the average oxidation state of the metal atoms, with T c being highest for an average valence of +1.5. KPbB 6 C 6 , with an ambient-pressure Eliashberg T c of 88 K, is predicted to possess the highest T c among the studied Pm$\bar{3}$n XB 3 C 3 or Pm$\bar{3}$XYB 6 C 6 phases, and calculations suggest it may be synthesized using high-pressure high-temperature techniques and then quenched to ambient conditions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Predictive Model of Charge Mobilities in Organic Semiconductor Small Molecules with Force-Matched Potentials

Charge mobility of crystalline organic semiconductors (OSC) is limited by local dynamic disorder. Recently, the charge mobility for several high mobility OSCs, including TIPS-pentacene, were accurately predicted from a density functional theory (DFT) simulation constrained by the crystal structure and the inelastic neutron scattering spectrum, which provide direct measures of the structure and the dynamic disorder in the length scale and energy range of interest. However, the computational expense required for calculating all of the atomic and molecular forces is prohibitive. In this work, we demonstrate the use of density functional tight binding (DFTB), a semiempirical quantum mechanical method that is 2 to 3 orders of magnitude more efficient than DFT. We show that force matching a many-body interaction potential to DFT derived forces yields highly accurate DFTB models capable of reproducing the low-frequency intricacies of experimental inelastic neutron scattering (INS) spectra and accurately predicting charge mobility. We subsequently predicted charge mobilities from our DFTB model of a number of previously unstudied structural analogues to TIPS-pentacene using dynamic disorder from DFTB and transient localization theory. The approach we establish here could provide a truly rapid simulation pathway for accurate materials properties prediction, in our vision applied to new OSCs with tailored properties.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Bayesian automated weighting of aggregated DFT, MD, and experimental data for candidate thermodynamic models of aluminum with uncertainty quantification

Atomic-scale modeling methods such as density functional theory (DFT) and molecular dynamics (MD) can predict the thermodynamic properties of materials at a lower cost than experimental measurements. However, their regular usage in thermodynamic model construction is hampered by the lack of quantitative agreement with experimental measurements and the lack of uncertainty estimates on the data. To make regular usage of this atomistic simulation data, it is important to assess whether the atomistic simulation datasets, by themselves or in combination with experimental measurements, result in the same physics-informed models best supported by experimental measurements alone. Here, models of aluminum thermodynamic properties are discussed using three data sources: atomistic calculations (DFT and MD), experiments, and a combination of atomistic calculations and experiments. The study shows that, after ensuring self-consistency in predicting key invariant points, both experimental measurements and atomistic calculations can significantly contribute to an optimal model.

36 MATERIALS SCIENCE↗

Interaction of Methane with Calcium Ferrite in the Chemical Looping Partial Oxidation Application: Experimental and DFT Study

In this study, experimental studies using mass spectroscopy, thermogravimetric analysis, in situ X-ray diffraction (XRD), and theoretical density functional theory calculations were carried out to investigate the mechanism of methane partial oxidation over a CaFe 2 O 4 oxygen carrier. The XRD analysis revealed the formation of the calcium ferrite phases during methane reduction occurred in the sequence CaFe 2 O 4 →CaFe 3 O 5 + Ca 2 Fe 2 O 5 →FeO + Fe + CaO→Fe + CaO. Each of these phases was modeled using periodic density functional theory to study the activity of these phases for the reactivity with CH 4 , CO, and H 2 along with the adsorbed surface species to help elucidate the effect that phase changes during reduction have on the chemical looping products CO and H 2 . Mass spectroscopic analysis revealed the formation of H 2 O and CO 2 occurred prior to the formation of CO and H 2 , and these products may be due to the removal of lattice oxygen bound in the discrete primary phase, which contributes to a small portion (~16%) of the total transferable oxygen in the solid. The density functional theory (DFT) results suggest that the gaseous CO and H 2 are not thermodynamically favorable to react with the Ca 2 Fe 2 O 5 where the release of CO and H 2 led to an increase in selectivity for these products. The DFT results suggest that some of the adsorbed atomic hydrogen combined to form adsorbed H 2 and then gaseous H 2 on the CaFe 2 O 4 surface, while the carbon deposition and some adsorbed atomic hydrogen react to form the CaFe 3 O 5 and Ca 2 Fe 2 O 5 phases simultaneously. DFT simulations suggest that both the gaseous H 2 and adsorbed atomic H are less thermodynamically favorable to drive the deep reduction of Ca 2 Fe 2 O 5 to CaO and Fe0. The deep reduction of the calcium ferrite therefore favors the generation of syngas products (CO and H 2 ) through mechanisms of methane dehydrogenation and selective oxidation of carbon to CO. The significant difference in mechanisms for partial oxidation of methane to syngas with the CaFe 2 O 4 carrier can be correlated to the differences in the oxide phases present during the reduction process, concentration of oxygen species, and bonding in the lattice structure.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Real-space formulation of the stress tensor for O (N) density functional theory: Application to high temperature calculations

We present an accurate and efficient real-space formulation of the Hellmann–Feynman stress tensor for O (N) Kohn–Sham density functional theory (DFT). While applicable at any temperature, the formulation is most efficient at high temperature where the Fermi–Dirac distribution becomes smoother and the density matrix becomes correspondingly more localized. We first rewrite the orbital-dependent stress tensor for real-space DFT in terms of the density matrix, thereby making it amenable to O (N) methods. We then describe its evaluation within the O (N) infinite-cell Clenshaw–Curtis Spectral Quadrature (SQ) method, a technique that is applicable to metallic and insulating systems, is highly parallelizable, becomes increasingly efficient with increasing temperature, and provides results corresponding to the infinite crystal without the need of Brillouin zone integration. We demonstrate systematic convergence of the resulting formulation with respect to SQ parameters to exact diagonalization results and show convergence with respect to mesh size to the established plane wave results. Here, we employ the new formulation to compute the viscosity of hydrogen at 10 6 K from Kohn–Sham quantum molecular dynamics, where we find agreement with previous more approximate orbital-free density functional methods.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Modeling Multi-Step Organic Reactions: Can Density Functional Theory Deliver Misleading Chemistry?

Many organic reactions are characterized by a complex mechanism with a variety of transition states and intermediates of different chemical natures. Their correct and accurate theoretical characterization critically depends on the accuracy of the computational method used. In this work, we study a complex ambimodal cycloaddition with five transition states, two intermediates, and three products, and we ask whether density functional theory (DFT) can provide a correct description of this type of complex and multifaceted reaction. Our work fills a gap in that most systematic benchmarks of DFT for chemical reactions have considered much simpler reactions. Our results show that many density functionals not only lead to seriously large errors but also differ from one another in predicting whether the reaction is ambimodal. Only a few of the available functionals provide a balanced description of the complex and multifaceted reactions. The parameters varied in the tested functionals are the ingredients, the treatment of medium-range and nonlocal correlation energy, and the inclusion of Hartree–Fock exchange. Furthermore, these results show a clear need for more benchmarks on the mechanisms of large molecules in complex reactions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Charge-density based evaluation and prediction of stacking fault energies in Ni alloys from DFT and machine learning

A combination of high strength and high ductility has been observed in multi-principal element alloys due to twin formation attributed to low stacking fault energy (SFE). In the pursuit of low SFE alloys, a key bottleneck is the lack of understanding of the composition–SFE correlations that would guide tailoring SFE via alloy composition. Using density functional theory (DFT), we show that dopant radius, which have been postulated as a key descriptor for SFE in dilute alloys, does not fully explain SFE trends across different host metals. Instead, charge density is a much more central descriptor. It allows us to (1) explain contrasting SFE trends in Ni and Cu host metals due to various dopants in dilute concentrations, (2) explain the large SFE variations observed in the literature even within a given alloy composition due to the nearest neighbor environments in “model” concentrated alloys, and (3) develop a machine learning model that can be used to predict SFEs in multi-elemental alloys. Further, this model opens a possibility to use charge density as a descriptor for predicting SFE in alloys.

36 MATERIALS SCIENCE↗

OrbNet Denali: A machine learning potential for biological and organic chemistry with semi-empirical cost and DFT accuracy

We present OrbNet Denali, a machine learning model for an electronic structure that is designed as a drop-in replacement for groundstate density functional theory (DFT) energy calculations. The model is a message-passing graph neural network that uses symmetryadapted atomic orbital features from a low-cost quantum calculation to predict the energy of a molecule. OrbNet Denali is trained on a vast dataset of 2.3 × 10 6 DFT calculations on molecules and geometries. This dataset covers the most common elements in biochemistry and organic chemistry (H, Li, B, C, N, O, F, Na, Mg, Si, P, S, Cl, K, Ca, Br, and I) and charged molecules. OrbNet Denali is demonstrated on several well-established benchmark datasets, and we find that it provides accuracy that is on par with modern DFT methods while offering a speedup of up to three orders of magnitude. For the GMTKN55 benchmark set, OrbNet Denali achieves WTMAD-1 and WTMAD-2 scores of 7.19 and 9.84, on par with modern DFT functionals. For several GMTKN55 subsets, which contain chemical problems that are not present in the training set, OrbNet Denali produces a mean absolute error comparable to those of DFT methods. For the Hutchison conformer benchmark set, OrbNet Denali has a median correlation coefficient of R 2 = 0.90 compared to the reference DLPNO-CCSD(T) calculation and R 2 = 0.97 compared to the method used to generate the training data (ωB97X-D3/def2- TZVP), exceeding the performance of any other method with a similar cost. Similarly, the model reaches chemical accuracy for noncovalent interactions in the S66x10 dataset. For torsional profiles, OrbNet Denali reproduces the torsion profiles of ωB97X-D3/def2-TZVP with an average mean absolute error of 0.12 kcal/mol for the potential energy surfaces of the diverse fragments in the TorsionNet500 dataset.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Metal-to-Ligand Charge-Transfer Spectrum of a Ru-Bipyridine-Sensitized TiO 2 Cluster from Embedded Multiconfigurational Excited-State Theory

Understanding optical properties of the dye molecule in dye-sensitized solar cells (DSSCs) from first-principles quantum mechanics can contribute to improving the efficiency of such devices. While density functional theory (DFT) and time-dependent DFT have been pivotal in simulating optoelectronic properties of photoanodes used in DSSCs at the atomic scale, questions remain regarding DFT’s adequacy and accuracy to furnish critical information needed to understand the various excited-state processes involved. In this work, we simulate the absorption spectra of a dye-sensitized solar cell analogue, comprised of a Ru-bipyridine (Ru-bpy) dye molecule and a small TiO2 cluster via DFT and via an accurate embedded correlated wavefunction (CW) theory. We generated CW spectra for the adsorbed Ru-bpy dye via a recently introduced capped density functional embedding theory or capped-DFET (to generate the embedding potential that accounts for the interaction of the molecule and the TiO 2 cluster). We then combined capped-DFET with the accurate but expensive multiconfigurational complete active space second-order perturbation theory (CASPT2)—embedded CASPT2. Because the CW theory is conducted on only a portion of the total system in the presence of an embedding potential that describes that portion’s interaction with its environment, we efficiently obtain CW-quality predictions that reflect local properties of the entire system. Specifically, for example, with capped-DFET and embedded CW theory, we can simulate accurately a plethora of metal-to-ligand charge-transfer excited properties at a manageable computational cost. Here, we predict detailed electronic spectra within the visible region, featuring the lowest three singlet and triplet excited states, along with predictions of the singlets’ lifetimes. We illustrated these results using a Jablonski diagram that show the relative energy position of the singlet and longer-lived triplet excited states and analyzed and proposed relaxation paths for the excited state corresponding to the most intense but short-lived absorption (interconversion, intersystem crossing, fluorescence, and phosphorescence) that may lead to longer-lived excited states necessary for efficient charge separation required to generate current in solar cells.

14 SOLAR ENERGY↗

Solid-State Hydroxide Ion Conductivity in Silver(I) Oxide, Ag 2 O

Silver (I) oxide, Ag 2 O, precipitated as microcrystals by combining aqueous silver(I) nitrate and KOH solutions, was found to be a solid-state hydroxide ion conductor with ionic conductivity on the order of 10 -3 S/cm. The proton chemical shifts at 4.87 and -7.35 ppm measured by solid-state 1H-NMR experiments are attributed to water molecules and in-lattice OH - coordinated to silver, respectively. The lack of spinning sidebands around the 4.87 ppm peak indicates rapid reorientation on the NMR timescale, suggesting that the water molecules are adsorbed to the surface of the Ag 2 O crystals. Pulsed field gradient measurements gave similar diffusion coefficients (2 x 10 -7 cm 2 /s at 298K) for all three proton environments indicating chemical exchange between sites on the millisecond timescale. The activation energy for OH - diffusion measured by NMR (0.18 eV) was comparable to that obtained by conductivity measurements and density functional theory (DFT) electronic structure calculations. Furthermore, the calculated Pourbaix diagram of Ag 2 O is consistent with the slightly lower sample density observed in He pycnometry and thermogravimetric measurements.

25 ENERGY STORAGE↗

Thorium Monosilicide, ThSi: An Experimental and Theoretical Study

The present theoretical and experimental combination study investigates the ThSi molecule in detail. Computationally, we utilized high-level multireference and coupled-cluster levels of theory conjoined with large correlation consistent basis sets to study a series of electronic and spin–orbit states of ThSi. Here, we report potential energy curves (PECs), electron configurations at equilibrium distances, spectroscopic constants, energetics, and spin–orbit coupling effects for 16 electronic states of ThSi. The studied 16 electronic states are arranged tightly within 0.9 eV, highlighting the complexity of the electronic spectrum of ThSi. The ground electronic state of ThSi is a single-reference 1 1 Σ + state that derives from the 1σ 2 2σ 2 1π 4 electronic configuration. The Ω = 0 + spin–orbit ground state of ThSi is composed of 1 1 Σ + (47%) and 13Π (44%) electronic states. Our measured bond energy (D0) of ThSi, obtained using resonant two-photon ionization (R2PI) spectroscopy is 3.146(4) eV, where the assigned error limit is given in parentheses in units of the last quoted digits. The computed D0 of ThSi (Ω = 0 + ) at the CBS-C-CCSD(T)-δT(Q)-δDK-δSO level (3.181 eV) is in good agreement with the experimental value. Our derived enthalpy of formation for ThSi, Δ f H 0K o (ThSi(g)), is 971.8(6.0) kJ/mol. Finally, we have performed density functional theory (DFT) calculations for ThSi(1 1 Σ + ) using 16 exchange correlation functionals that span multiple rungs of “Jacob’s ladder” of density functional approximation (DFA) to assess the DFT errors on D 0 , r e , and ω e of ThSi with respect to experimental and ab initio coupled-cluster values.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Spectroscopic Insight on Neodymium Solvation in Lithium Borohydride-Supported Electrolyte

Borohydride-based electrolytes have recently emerged as promising media for the electrodeposition of electropositive metals, including rare earth (RE) elements. While the presence of supporting alkali metal cations and RE counteranions provides essential electrochemical conductivity for achieving fast metal electrodeposition, interactions between the host ligand and solvated neodymium (Nd) complexes remain unclear. This study provides insights into the coordination structure of a concentrated and directly solvated Nd salt in a lithium borohydride-supported electrolyte. Our spectroscopic results indicate that the RE coordination environment is significantly influenced by the solvation mechanism, which can vary between metathesis and complexation pathways, primarily dictated by stoichiometric factors. Under dilute conditions, nearly complete metathesis of anions leads to a high coordination number for the host ligand (borohydride), consistent with the previously reported solvated Nd speciation in chlorine-free electrolytes. In contrast, concentrated dissolution of the Nd salt in the supported electrolyte is dominated by a complexation pathway featuring a Li-ion-paired complex with a low coordination number of the host ligand. Density functional theory (DFT) calculations indicated that the observed blue shift in the borohydride vibration was the result of an increase in electron density drawn into the terminal B–H interbond region from the hydride as the coordination changed from Li to Nd. In conjunction with DFT results, vibrational analyses allowed correlation of the experimental shifts associated with changes in Nd ligation and coordination spheres, further consolidating the prevalence of highly chloride-coordinated species under concentrated conditions. In conclusion, the outcomes of this work illuminate the distinctive and heterogeneous coordination structures that the electroactive RE species can adopt at high concentrations in lithium borohydride-supported electrolytes, as a key step to comprehend the reported metal electrodeposition performance in these media.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Toward Quantum Chemical Free Energy Simulations of Platinum Nanoparticles on Titania Support

Platinum nanoparticles (Pt-NPs) supported on titania surfaces are costly but indispensable heterogeneous catalysts because of their highly effective and selective catalytic properties. Therefore, it is vital to understand their physicochemical processes during catalysis to optimize their use and to further develop better catalysts. However, simulating these dynamic processes is challenging due to the need for a reliable quantum chemical method to describe chemical bond breaking and bond formation during the processes but, at the same time, fast enough to sample a large number of configurations required to compute the corresponding free energy surfaces. Density functional theory (DFT) is often used to explore Pt-NPs; nonetheless, it is usually limited to some minimum-energy reaction pathways on static potential energy surfaces because of its high computational cost. In this work, we report a combination of the density functional tight binding (DFTB) method as a fast but reliable approximation to DFT, the steered molecular dynamics (SMD) technique, and the Jarzynski equality to construct free energy surfaces of the temperature-dependent diffusion and growth of platinum particles on a titania surface. In particular, we present the parametrization for Pt-X (X = Pt, Ti, or O) interactions in the framework of the second-order DFTB method, using a previous parametrization for titania as a basis. The optimized parameter set was used to simulate the surface diffusion of a single platinum atom (Pt 1 ) and the growth of Pt 6 from Pt 5 and Pt 1 on the rutile (110) surface at three different temperatures (T = 400, 600, 800 K). The free energy profile was constructed by using over a hundred SMD trajectories for each process. We found that increasing the temperature has a minimal effect on the formation free energy; nevertheless, it significantly reduces the free energy barrier of Pt atom migration on the TiO 2 surface and the transition state (TS) of its deposition. In a concluding remark, the methodology opens the pathway to quantum chemical free energy simulations of Pt-NPs’ temperature-dependent growth and other transformation processes on the titania support.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Theory-guided design of duplex-phase multi-principal-element alloys

Density-functional theory (DFT) is used to identify phase-equilibria in multi-principal-element and high-entropy alloys (MPEAs/HEAs), including duplex-phase and eutectic microstructures. Here, a combination of composition-dependent formation energy and electronic-structure-based ordering parameters were used to identify a transition from FCC to BCC favoring mixtures, and these predictions experimentally validated in the Al-Co-Cr-Cu-Fe-Ni system. A sharp crossover in lattice structure and dual-phase stability as a function of composition were predicted via DFT and validated experimentally. The impact of solidification kinetics and thermodynamic stability was explored experimentally using a range of techniques, from slow (castings) to rapid (laser remelting), which showed a decoupling of phase fraction from thermal history, i.e., phase fraction was found to be solidification rate-independent, enabling tuning of a multi-modal cell and grain size ranging from nanoscale through macroscale. Strength and ductility tradeoffs for select processing parameters were investigated via uniaxial tension and small-punch testing on specimens manufactured via powder-based additive manufacturing (directed-energy deposition). This work establishes a pathway for design and optimization of next-generation multiphase superalloys via tailoring of structural and chemical ordering in concentrated solid solutions.

36 MATERIALS SCIENCE↗

Universal machine learning framework for defect predictions in zinc blende semiconductors

Our article introduces a universal predictive framework for point defect formation energies and charge transition levels in a wide chemical space of zinc blende semiconductors and possible impurity atoms selected from across the periodic table. This framework was developed by leveraging high-throughput quantum mechanical simulations benchmarked using some experimental data from the literature, as well as machine learning (ML)-based regressions techniques that map unique materials descriptors to computed defect properties and yield optimized and generalizable models. Furthermore, the power and utility of these models is revealed through quick predictions for thousands of new defects and screening of low-energy impurities, which may tune the equilibrium conductivity in the semiconductor. This work presents, to our knowledge, the largest density functional theory (DFT) dataset of defect properties in semiconductors and the largest DFT+ML-based screening of point defects in semiconductors to date.

36 MATERIALS SCIENCE↗