Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “DFT 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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Discovery of a Ferromagnetic Nickel Chalcogenide Nanocluster Ni 3 S 3 H(PEt 3 ) 5

Atomically precise ligated nanoclusters (NC) are promising cluster-based materials with novel molecular architectures and tunable magnetic properties. Herein, the synthesis and characterization of a nickel sulfide NC Ni 3 S 3 H(PEt 3 ) 5 (PEt 3 = triethylphosphine) with distinct magnetic properties are reported. Magnetization measurements reveal its magnetic moment of 1.5 µ B in the solid phase, consistent with the existence of one unpaired electron predicted by density functional theory (DFT) calculations. Additionally, experimental measurements indicate the presence of ferromagnetic ordering within each Ni 3 S 3 H(PEt 3 ) 5 NC and strong coercivity at temperatures below 20 K. Ion mobility-mass spectrometry is employed in conjunction with DFT calculations and collision cross-section simulations to investigate the structure of the isolated Ni 3 S 3 H(PEt 3 ) 5 . Theoretical studies show that [Ni 3 S 3 H(PEt 3 ) 5 ] + has a planar Ni 3 S 3 core where three Ni atoms are arranged in a triangle with three bridging S atoms residing in the same plane. This structure is preserved in both solution and solid phases, which is confirmed by spectroscopic studies of Ni 3 S 3 H(PEt 3 ) 5 . Additionally, DFT calculations indicate that all spins at the Ni sites are aligned parallel, confirming the presence of ferromagnetic coupling. Overall, this study provides key insights into the structure and magnetic properties of Ni 3 S 3 H(PEt 3 ) 5 , which will facilitate the design of new NC-based magnetic materials.

Nickel sulfide nanocluster↗
DOE OSTI2023

Electronic structure and energetics of a heterodimeric BChl g ′/Chl a ′ special pair generated by exposure of Heliomicrobium modesticaldum to dioxygen

Heliobacteria are anoxygenic phototrophs that have a Type I homodimeric reaction center containing bacteriochlorophyll g (BChl g). Previous experimental studies have shown that in the presence of light and dioxygen, BChl g is converted into 8 1 -OH-chlorophyll a F (hereafter Chl a F ), with an accompanying loss of light-driven charge separation. These studies suggest that the reaction center only loses the ability to transfer electrons once both BChl g′ molecules of the P 800 special pair have been converted to Chl a F ′. Here, the present work confirms that the partially converted BChl g′/Chl a F ′ special pair remains functional in samples exposed to dioxygen by demonstrating its presence using hyperfine couplings obtained from Q-band 1 H ENDOR, 2D 14 N HYSCORE and DFT methods. The DFT calculations of the BChl g′/BChl g′ homodimeric primary donor, which are based on the recently published X-ray crystal structure, predict that the unpaired electron spin is equally delocalized over both BChl g′ molecules and provide an excellent match to the experimental hyperfine couplings of the anaerobic samples. Exposure to dioxygen leads to substantial changes in the hyperfine interactions, indicative of greater localization of the unpaired electron spin. The measured hyperfine couplings are reproduced in the DFT calculations by replacing one of the BChl g′ molecules of the primary donor with a Chl a F ′ molecule. The calculations reveal that the spin density becomes localized on BChl g′ in the heterodimeric primary donor. Time-dependent DFT calculations demonstrate that conversion of either or both of the accessory BChl g molecules and/or one of the BChl g′ molecules of P800 to Chl a F ′ results in minor effects on the energy of the charge-separated states. In contrast, if both of the BChl g′ molecules of P 800 are converted a large increase in the energy of the charge-separated state occurs. This suggests that the reaction center remains functional when only one half of the dimer is converted, however, conversion of both halves of the P 800 dimer leads to loss of function

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

Dispersion and stability mechanism of Pt nanoparticles on transition-metal oxides

The heterogeneous catalysts of Pt/transition-metal oxides are typically synthesized through calcination at 500 °C, and Pt nanoparticles are uniformly and highly dispersed when hydrogen peroxide (H 2 O 2 ) is applied before calcination. The influence of H 2 O 2 on the dispersion and the stability of Pt nanoparticles on titania-incorporated fumed silica (Pt/Ti–FS) supports was examined using X-ray absorption fine structure (XAFS) measurements at the Pt L 3 and Ti K edges as well as density functional theory (DFT) calculations. The local structural and chemical properties around Pt and Ti atoms of Pt/Ti–FS with and without H 2 O 2 treatment were monitored using in-situ XAFS during heating from room temperature to 500 °C. XAFS revealed that the Pt nanoparticles of H 2 O 2 -Pt/Ti–FS are highly stable and that the Ti atoms of H 2 O 2 -Pt/Ti–FS support form into a distorted-anatase TiO 2 . DFT calculations showed that Pt atoms bond more stably to oxidized–TiO 2 surfaces than they do to bare- and reduced–TiO 2 surfaces. XAFS measurements and DFT calculations clarified that the presence of extra oxygen atoms due to the H 2 O 2 treatment plays a critical role in the strong bonding of Pt atoms to TiO 2 surfaces.

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

Accelerating self-consistent field iterations in Kohn-Sham density functional theory using a low-rank approximation of the dielectric matrix

We present an efficient preconditioning technique for accelerating the fixed-point iteration in real-space Kohn-Sham density functional theory (DFT) calculations. The preconditioner uses a low-rank approximation of the dielectric matrix (LRDM) based on Gâteaux derivatives of the residual of fixed-point iteration along appropriately chosen direction functions. We develop a computationally efficient method to evaluate these Gâteaux derivatives in conjunction with the Chebyshev filtered subspace iteration procedure, an approach widely used in large-scale Kohn-Sham DFT calculations. Further, we propose a variant of LRDM preconditioner based on adaptive accumulation of low-rank approximations from previous self-consistent field iterations, and also extend the LRDM preconditioner to spin-polarized Kohn-Sham DFT calculations. We demonstrate the robustness and efficiency of the LRDM preconditioner against other widely used preconditioners on a range of benchmark systems with sizes ranging from ~100 to 1100 atoms (~500–20,000 electrons). The benchmark systems include various combinations of metal-insulating-semiconducting heterogeneous material systems, nanoparticles with localized d orbitals near the Fermi energy, nanofilm with metal dopants, and magnetic systems. In all benchmark systems, the LRDM preconditioner converges robustly within 20–30 iterations. In contrast, other widely used preconditioners show slow convergence in many cases, as well as divergence of the fixed-point iteration in some cases. Lastly, we demonstrate the computational efficiency afforded by the LRDM method, with up to 3.4-fold reduction in computational cost for the total ground-state calculation compared to other preconditioners.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Insights on carbon dioxide adsorption in a flexible 1-D coordination polymer from in situ X-ray scattering and density functional theory

A combination of in situ small-angle X-ray scattering (SAXS) microstructure characterization over scales ranging from 1 nm to 10 µm and powder X-ray diffraction (XRD) structure characterization under various gas/pressure/temperature conditions with density functional theory (DFT) calculations provides new insights for the CO2 sorption behavior of a one-dimensional porous coordination polymer: catena-bis­(di­benzoyl­methanato)(4,4′-bi­pyridyl)­nickel(II), denoted NiDBM-Bpy. The NiDBM-Bpy chains are held together by van der Waals forces, but the structure of guest-free NiDBM-Bpy is unsolved due to a lack of suitable crystals and high-quality powder XRD patterns. Nevertheless, SAXS and powder XRD can follow microstructural and structural changes as a function of gas pressure, composition and temperature. Both mixed-gas flow and static supercritical CO2 regimes are explored experimentally. DFT calculations are used to model the structural variation associated with XRD changes and hysteresis in the sorption isotherms. XRD and DFT calculations suggest that an orthorhombic Fddd structure with two CO2 per Ni emerges following a transition from the structure with lower molar volume and symmetry that exists without CO2 present.

Allen, Andrew J.↗
DOE OSTI2021

Identifying Order and Disorder in Double Four-Membered Rings via Raman Spectroscopy during Crystallization of LTA Zeolite

Fluoride (F – ) has been essential for the synthesis of low-defect siliceous zeolites. It has been hypothesized that F – balances the positive charges from organic structure-directing agents (OSDAs) and stabilizes key building units during zeolite crystallization such as the double four-membered ring (D4R). However, due to the lack of characterization techniques for investigating medium-range structures, including rings and cages formed during zeolite crystallization, the roles of F – in stabilizing building units and maintaining local charge balance during zeolite assembly are not yet fully understood. Here, the crystallization of siliceous Linde type A (LTA) zeolite in the presence of F – was investigated using Raman spectroscopy and periodic density functional theory (DFT) calculations. We have discovered that the F – -filled double four-membered ring (F – /D4R) and the empty D4R exhibit rather distinct Raman features. Both the F – /D4R and empty D4R are formed in the LTA zeolite synthesized in the presence of F – using 1,2-dimethyl-3-(4-methylbenzyl) imidazolium as the OSDA. Observed Raman bands of the F – /D4R and empty D4R, along with predictive DFT calculations on LTA supercells, reveal an ordered distribution of these two D4R units in the final as-made LTA zeolite. The discovery of these distinct Raman signatures of F – /D4R and empty D4R units opens an interesting new window for studying defects in the D4R during zeolite formation. In particular, we have observed variation in Raman intensities of F – /D4R and empty D4R bands during LTA crystallization; periodic DFT calculations indicate that the observed Raman behavior is consistent with empty D4R units containing one or two Si vacancies surrounded by Q 3 Si—consistent also with solid-state nuclear magnetic resonance measurement. These defects appear to heal during further crystallization, leading to the formation of defect-free LTA zeolite crystals. Altogether, our results provide deeper understanding on the roles of F – in charge balancing and stabilizing intact D4R units during zeolite formation.

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

Fermi surface topology and magnetotransport properties of superconducting Pd 3 Bi 2 Se 2

Pd 3 Bi 2 Se 2 is a rare realization of a superconducting metal with a non-zero topological invariant. Here, in this study, we report the growth of high-quality single crystals of layered Pd 3 Bi 2 Se 2 with a superconducting transition at T c ≈ 0.80 K and upper critical fields of ~10 mT and ~5 mT for the in plane and out-of-plane directions, respectively. Our density functional theory (DFT) calculations reveal three pairs of doubly degenerate bands crossing the Fermi level all displaying clear three dimensional dispersion consistent with the overall low electronic anisotropy (<2). The multiband electronic nature of Pd 3 Bi 2 Se 2 is evident in magneto-transport measurements, yielding a sign changing Hall resistivity at low temperatures. The magnetoresistance is non-saturating and follows Kohler’s scaling rule. We interpret the magneto-transport data in terms of open orbits that are revealed in the DFT calculated Fermi surface. de Haas-van Alphen (dHvA) oscillation measurements using torque magnetometry on single crystals yield four frequencies for out-of-plane fields: F α = (150 ± 26) T, F β = (293 ± 10) T, F γ = (375 ± 20) T and F η = (1017 ± 12) T, with the low frequency dominating the spectrum. Through the measurement of angular dependent dHvA oscillations and DFT calculations we identify the F α frequency with an approximately ellipsoidal electron pocket centered on the L 2 point of the Brillouin zone. Lifshitz-Kosevich analysis of the dHvA oscillations reveals a small cyclotron effective mass m* = (0.11 ± 0.02)m 0 and a nontrivial Berry phase for the dominant orbit. The presence of nontrivial topology in a bulk superconductor positions Pd 3 Bi 2 Se 2 as a potential candidate for exploring topological superconductivity.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Rapid Design of Cost-Effective Refractory High Entropy Alloys Strengthened by Precipitation

Funded by the ARPA-E ULTIMATE program, this project aims to develop cost-effective refractory high entropy alloys (RHEAs) for gas turbine blades applications at 1300 degree Celsius and above. The alloy design focuses on carbides precipitation strengthening in a strong and ductile solid solution BCC matrix. The goal is to achieve balanced mechanical properties including room temperature ductility & fracture toughness and high temperature strength and creep resistance, while maintaining comparable oxidation resistance and densities to Ni-base superalloys. To accelerate alloy design, multiscale modeling including CALPHAD and first-principles density functional theory (DFT) calculations and machine learning are performed. Specifically, high throughput phase diagram calculations and screening are carried out using CALPHAD; DFT calculations are performed to predict intrinsic ductility, grain boundary strengthening potential, coefficient of thermal expansion (CTE) and temperature-dependent elastic constants. Temperature dependent yield stress, creep strain and CTE are also predicted using machine learning. Down-selected alloys are synthesized in small buttons of about 250g using arc melting for rapid evaluation on microstructure and mechanical properties before employing plasma arc melting and additive manufacturing for producing large ingots. Preliminary computational and experimental results of this project will be presented and future direction in RHEA development will be discussed.

high temperature materials↗
DOE OSTI2021

CO and H 2 adsorption on Au-Ni bimetallic surfaces: a combined experimental and DFT theoretical study

Au-Ni bimetallic thin films were grown on refractory metal substrates. CO and H 2 adsorption on Au-Ni bimetallic surfaces have been studied by a combination of in situ polarization modulation infrared reflection absorption spectroscopy (PM-IRRAS), temperature-programmed desorption (TPD), and density functional theory (DFT) calculations. It is found CO desorption peak shifts from 413 K on pure Ni surfaces to 293 K on the isolated Ni atoms formed by alloying with Au atoms. The sharp decrease of CO desorption temperature on Au-Ni surfaces with increasing Au coverage is caused by the change of the favored CO adsorption sites from bridge/hollow sites on pure Ni surfaces to Ni top sites on Au-Ni bimetallic surfaces. In situ PM-IRRAS shows two CO adsorption bands on Au sites at 2119 cm −1 and 2103 cm −1 on Au-Ni surfaces at 80 K, which are due to CO bound on under-coordinated Au atoms and electron negatively charged Au sites modified with nearby Ni atoms, respectively. Even with the Au-Ni surface temperature at as low as 100 K, CO adsorption induced Ni surface segregation has been observed by in situ PM-IRRAS. Furthermore, DFT calculation results discover the adsorption energy of CO on Ni top sites continues to decrease with increasing Au coverage due to the geometric ensemble effect and the lowered d-band center after Ni alloying with Au. H 2 desorption temperature decreases from 363 K on pure Ni thin films to 302 K with increasing Au coverage to 0.6 ML. A new H 2 peak appears at around 170 K on the Au-Ni surfaces with Au coverages between 0.6 ML and 0.9 ML. This new H 2 TPD peak is assigned to H 2 desorption from the totally isolated Ni sites. With Au coverage above 1.5 ML, there is no any H 2 desorption detected. Finally, the combined surface science studies and DFT calculations provide new insights into the surface structure-activity correlation of Ni-base bimetallic surface alloys.

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

Using density functional theory to construct multiphase equations of state: tin as an example

We present a general methodology for using density functional theory (DFT) calculations as a basis for constructing multiphase equations of state (EOS). Focusing on tin as an example, we discuss the full process of generating a tabular EOS, starting from first-principles calculations and arriving at a full EOS in agreement with existing experimental data. We begin by describing DFT calculations for the five solid phases and liquid phase of tin, including cold curve, phonon, and DFT-based molecular dynamics calculations. We then discuss methods for incorporating DFT results into materials models used in OpenSesame, the program used to create a full tabular multiphase EOS. Next, we outline a general strategy for adjusting model parameters in OpenSesame to fit to a variety of experimental data, including isobaric data, isothermal data, shock data, solid-solid phase boundary measurements, and measurements of the melt curve. We end with a discussion of the advantages of using DFT data as a foundation for EOS generation and discuss potential strategies for automated EOS generation based on the fitting process highlighted in this work.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Hydrodechlorination of 1,2-Dichloroethane on Platinum Catalysts: Insights from Reaction Kinetics Experiments, Density Functional Theory, and Microkinetic Modeling

Catalytic hydrodechlorination is a promising strategy for treating industrial 1,2-dichloroethane wastes, for which Pt and Pt-based alloy catalysts are widely used. Here, we performed a detailed mechanistic study for 1,2-dichloroethane hydrodechlorination on Pt using a synergistic approach combining density functional theory (DFT) calculations, reaction kinetics experiments, and microkinetic modeling. Using planewave DFT calculations, we evaluated the reaction energy and activation energy barrier of each elementary step involved in the reaction network on Pt(111). The calculated energetics were then incorporated into a comprehensive mean-field microkinetic model accounting for a total of 65 elementary steps. The model-predicted reaction rates were compared with the results from our reaction kinetics experiments on SiO2-supported Pt catalysts. Our results indicated that the hydrodechlorination of 1,2-dichloroethane on Pt(111) starts with a H-removal step; then, it proceeds through a sequence of alternating dechlorination and dehydrogenation steps until vinylidene (CH 2 C*) is formed; finally, CH 2 C* is hydrogenated to the final product, ethane, sequentially via vinyl (CH 2 CH*), ethylene, and ethyl (CH 3 CH 2 *) intermediates. After model parameter adjustments, we achieved good agreement between our theoretical model and experimental results; the adjustments to the calculated parameters are consistent with the typically anticipated coverage effects. Furthermore, our study offers valuable mechanistic insights, which are useful for improving catalysts for this chemistry.

1,2-dichloroethane↗
DOE OSTI2023

Revisiting two thiophosphate compounds constituting d 0 transition metal HfP 2 S 6 and d 10 transition metal α-Ag 4 P 2 S 6 as multifunctional materials for combining second harmonic generation response and photocurrent response

Two acentric thiophosphate compounds, HfP 2 S 6 and α-Ag 4 P 2 S 6 , are revisited and studied as infrared nonlinear optical materials. HfP 2 S 6 and α-Ag 4 P 2 S 6 were structurally characterized without any property measurements. Here, in this work, HfP 2 S 6 and α-Ag 4 P 2 S 6 were synthesized via high temperature salt flux reactions. Low-temperature polymorph acentric α-Ag 4 P 2 S 6 was purified and grown as mm-sized crystals with the aid of AgBr flux. The AgBr flux was revealed to play an important role in stabilizing the acentric α-Ag 4 P 2 S 6 . The acentric α-Ag 4 P 2 S 6 transforms to centrosymmetric β-Ag 4 P 2 S 6 at 850(5) K, which is revealed by differential scanning calorimetry (DSC) analysis and powder X-ray diffraction experiments. HfP 2 S 6 and α-Ag 4 P 2 S 6 are discovered by UV-vis spectrum measurements as indirect bandgap semiconductors with bandgaps of 2.2(1) eV and 2.5(1) eV, respectively, which is supported by DFT calculations and TB-LMTO-ASA calculations. The bonding pictures of α-Ag 4 P 2 S 6 were studied by crystal orbital Hamilton population calculations (COHP) coupled with electron localization function (ELF) analysis. DFT calculations predict that HfP 2 S 6 and α-Ag 4 P 2 S 6 would exhibit different optical performances regardless of being constructed from identical [P 2 S 6 ] motifs. HfP 2 S 6 exhibits a low second harmonic generation (SHG) response, ~0.21 × AGS (for the sample of particle size of 25 μm). α-Ag 4 P 2 S 6 possesses moderate SHG response, ~0.61 × AGS (for the sample of particle size of 225 μm) coupled with a high laser damage threshold (LDT) of ~3.2 × AGS. Characteristics of high ambient stability, moderate bandgap and SHG response, type-I phase-matching capability, and high LDT together with the easy growth of large crystals make α-Ag 4 P 2 S 6 attractive for future infrared nonlinear optical applications. Photocurrent measurements found that α-Ag 4 P 2 S 6 and β-Ag 4 P 2 S 6 have high photocurrent response, 165 nA cm -2 and 135 nA cm -2 , respectively. α-Ag 4 P 2 S 6 is a new multifunctional material of the ternary Ag–P–S system, which combines nonlinear optical (NLO) properties and photocurrent response.

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

Impact of solvation on the GW quasiparticle spectra of molecules

First-principles calculations for electrochemistry require accurate treatment of both electronic structure and solvation. The perturbative GW approximation starting from density functional theory (DFT) calculations accurately models materials systems with varying dimensionality. Continuum solvation models enable efficient treatment of solvation effects in DFT calculations, but their applications with beyond-DFT electronic structure methods such as GW have been limited. Here, we introduce the frequency-dependent liquid polarizability from a nonlocal continuum solvation model in the screened Coulomb interaction of full-frequency GW calculations with a solvated DFT starting point. In this study, we show that the liquid screening contributions substantially reduce the HOMO–LUMO gap of molecules by 3–5 eV, while solvent effects on the DFT starting point negligibly impact the GW gap. The resulting framework facilitates the simultaneous electronic and solvation accuracy needed for first-principles electrochemistry.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Ligand Many-Body Expansion as a General Approach for Accelerating Transition Metal Complex Discovery

Methods that accelerate the evaluation of molecular properties are essential for chemical discovery. While some degree of ligand additivity has been established for transition metal complexes, it is underutilized in asymmetric complexes, such as the square pyramidal coordination geometries highly relevant to catalysis. To develop predictive methods beyond simple additivity, we apply a many-body expansion to octahedral and square pyramidal complexes and introduce a correction based on adjacent ligands (i.e., the cis interaction model). We first test the cis interaction model on adiabatic spin-splitting energies of octahedral Fe(II) complexes, predicting DFT-calculated values of unseen binary complexes to within an average of 1.4 kcal/mol. Uncertainty analysis reveals the optimal basis, comprising the homoleptic and mer symmetric complexes. We next show that the cis model (i.e., the cis interaction model solved for the optimal basis) infers both DFT- and CCSD(T)-calculated model catalytic reaction energies to within 1 kcal/mol on average. The cis model predicts low-symmetry complexes with reaction energies outside the range of binary complex reaction energies. We observe that trans interactions are unnecessary for most monodentate systems but can be important for some combinations of ligands, such as complexes containing a mixture of bidentate and monodentate ligands. Lastly, we demonstrate that the cis model may be combined with Δ-learning to predict CCSD(T) reaction energies from exhaustively calculated DFT reaction energies and the same fraction of CCSD(T) reaction energies needed for the cis model, achieving around 30% of the error from using the CCSD(T) reaction energies in the cis model alone.

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

Analysis of Intermediates and Products from the Dehydrogenation of Mg(BH 4 ) 2

The thermodynamic properties of key compounds, Mg(B 3 H 8 ) 2 , MgB 2 H 6 , MgB 10 H 10 , Mg(B 11 H 14 ) 2 , Mg 3 (B 3 H 6 ) 2 , and MgB 12 H 12 , proposed to be formed in the release of hydrogen from magnesium borohydride Mg(BH 4 ) 2 and uptake of hydrogen by MgB 2 , have been investigated using solid–state density functional theory (DFT) calculations. Further, more accurate treatment of cell–size effects to the entropies were also investigated, in order to improve the accuracy of the thermodynamic properties of complex borohydrides. We find that the zero–point energy corrections can lower the electronic energies of reaction by 20–30 kJ/(mol H 2 ) for these intermediates, while adding the thermal and entropy contributions results in the total decrease up to ~50 kJ/(mol H 2 ). Although our treatment lowers the calculated formation energy of Mg(B 3 H 8 ) 2 , it is still too high to explain the experimental observation of B 3 H 8 – . We discuss possible reasons for this disparity and propose that the formation of B 3 H 8 – and H – in a disordered amorphous phase has a large energy difference compared to the phase–separated Mg(B 3 H 8 ) 2 and MgH 2 considered in calculations. Comparison of the experimental and NMR chemical shifts calculated within a DFT approach for known species Mg(BH 4 ) 2 , Mg(B 3 H 8 ) 2 , Mg(B 11 H 14 ) 2 , MgB 10 H 10 and MgB 12 H 12 provides validation for predicting the chemical shifts of the other compounds which are yet to be confirmed experimentally. These include MgB 2 H 6 and the proposed tri–anion species Mg 3 (B 3 H 6 ) 2 that both have favorable thermodynamics for reversible hydrogen storage in Mg(BH 4 ) 2 without the formation of MgH 2 as a co–product which could phase–separate and inhibit rehydrogenation.

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

Simulations and Experiments Reveal Effect of Nanopores on Helium Diffusion in Quartz

The diffusion properties of noble gases in minerals are widely used to reconstruct the thermal histories of rocks. In this work, we combine density functional theory (DFT) calculations with laboratory experiments to investigate controls on helium diffusion in quartz. DFT calculations for perfect α-quartz predict substantially lower activation energies and frequency factors for helium diffusion than observed in laboratory experiments, especially in the [001] direction. These results imply that no helium could be retained in quartz at Earth surface temperatures, which conflicts with observations of partial cosmogenic 3 He retention over geologic time scales. Here, we implement a model of helium diffusion in α-quartz modulated by nanopore defects that disrupt energetically favorable diffusion pathways. In this model, we find that laboratory-determined diffusivities can be most closely reproduced when a helium atom encounters ~70 nanopore sites per million interstitial sites. The results of our model indicate that diffusion of helium in natural quartz, like other noble gases in other minerals, can be significantly modulated by extended defects.

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

The Structural Characterisation and DFT-Aided Interpretation of Vibrational Spectra for Cyclo(l-Cys-d-Cys) Cyclic Dipeptide in a Solid State

Cyclic dipeptides with two intramolecular peptide bonds forming a six-membered 2,5-diketopiperazine ring are gaining significant attention due to their biological and chemical properties. Small changes in the local geometry of such molecules (from cis to trans) can lead to significant structural differences. This work presents the results of a study of cyclo(l-Cys-d-Cys), a dipeptide comprising two cysteine molecules in opposite chiral configurations, with the functional groups situated at both sides of the diketopiperazine ring. X-ray diffraction (XRD) experiment revealed that the molecule crystallises in the $P$-1 space group, which includes the centre of inversion. The IR and Raman vibrational spectra of the molecule were acquired and interpreted in terms of the potential energy distribution (PED) according to the results of density functional theory (DFT) calculations. The DFT-assisted analysis of energy frameworks for the hydrogen bond network within molecular crystals was performed to support the interpretation of X-ray structural data. The optimisation of the computational model based on three-molecule geometry sections from the crystallographic structure, selected to appropriately reflect the intermolecular interactions responsible for the formation of 1D molecular tapes in cyclo(l-Cys-d-Cys) crystal, allowed for better correspondence between theoretical and experimental vibrational spectra. This work can be considered the first complete structural characterisation of cyclo(l-Cys-d-Cys), complemented via vibrational spectroscopy results with full band assignment aided with the use of the DFT method.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Toward an Understanding of Linear Scaling Relations through Energy Decomposition Analysis

The discovery of linear scaling relations has fundamentally changed the field of heterogeneous catalysis. The scaling relations have been rationalized based on the d-band theory, specifically a separation of sp and d electron contributions to adsorption energies. Within the framework of energy decomposition analysis, a full understanding of such a separation would require one to further break down the adsorption energy into distinct energy components such as electrostatics, polarization, charge transfer, and van der Waals interactions, and to examine the sp and d contributions to each of them. As a step in this direction, we analyzed the interaction energy between CH x (x = 1–4) adsorbates and fcc(100) transition metal surfaces (M = Cu, Ag, Au, Rh, and Pt), with the surfaces represented both as slabs in plane-wave density functional theory (pw-DFT) calculations and as atomic clusters in atomic-orbital basis density functional theory (ao-DFT) calculations. Through an absolutely localized molecular orbital (ALMO) based energy decomposition analysis of the ao-DFT adsorption energy, each of the interaction energy components (electrostatics, polarization, van der Waals, and charge transfer) was found to follow its own scaling relations, with an intricate interplay among these energy components yielding the overall scaling relations for the total adsorption energies. Using the recently introduced ALMO-based polarization and charge-transfer analysis schemes, we further dissected polarization into metal surface and adsorbate contributions, and charge transfer into metal → adsorbate and adsorbate → metal contributions. The contributions from the sp and d electrons of the metal to these terms were further quantified, and the dominant role of the metal d electrons was reaffirmed. These results shed light on how CHx adsorbates interact with metal surfaces and further reveal the physical origin of the scaling relations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗