Search NASA⌕ Search

SEARCH · Search NASA

Results for “density functional theory (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 217 records · Page 12
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 OSTI2024

INSPIRED: Inelastic neutron scattering prediction for instantaneous results and experimental design

Inelastic neutron scattering (INS) has unique advantages in probing how atoms vibrate and how the vibrations propagate and interact. Such dynamic information is crucial in understanding various material properties, from heat capacity, thermal conductivity, phase transitions, and chemical reactions to more exotic quantum behavior. The analysis and interpretation of the INS spectra often start from a model structure of the sample, followed by a series of calculations to obtain the simulated spectra to compare with experiments. The conventional way to perform such calculations usually requires significant time, computing resources, and specialized expertise. Here, we present a new program named INSPIRED (Inelastic Neutron Scattering Prediction for Instantaneous Results and Experimental Design), which enables users to perform rapid INS simulations in several different ways on their personal computers in just a few clicks, with the crystal structure as the only input file. Specifically, the users can choose a pre-trained symmetry-aware neural network (coupled with an autoencoder) to predict the phonon density of states (DOS), 1D S(E) and 2D S(|Q|,E) spectra for any given structure. One can also choose an existing density functional theory (DFT) calculation from a database (containing over 12,000 crystals), and quickly obtain the simulated INS spectra for single crystals and powders. It is also possible to use pre-trained universal machine learning force fields to relax a given crystal structure, calculate the phonon dispersion and DOS, and, subsequently, the INS spectra. All these functions are implemented with a PyQt graphic user interface. Finally, we expect these new tools will benefit broad user communities and significantly improve the efficiency of experiment design, execution, and data analysis for INS.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

TaTi_BCC_SolidSolution_128atoms_VASP6

We performed density functional theory (DFT) calculations for body-centered-cubic (BCC) structures with 128 lattices sites of solid solution binary alloys tantalum-titanium (Ta-Ti). The electronic structures of alloys have been calculated using Vienna Ab initio Simulation Package (VASP). Within this package the DFT approach is used to reduce many-body Schrodinger equation to set of single particle Kohn-Sham (KS) equations. The generalized electronic exchange-correlation functional is described by generalized gradient approximation with the Perdew-Burke-Ernzerhof parametrization. The electron-ion interactions is described by pseudopotentials developed within the plane-wave basis projector augmented-wave (PAW) approach. These pseudopotentials are available at the VASP portal (http://cms.mpi.univie.ac.at/vasp/). Our calculations have been run with the pseudopotentials treating s and p semi-core states as valence in case for the elements Hf and Nb. The electronic densities and potentials are expanded over plane-waves with energy cutoff of 350 eV. 2x2x2 k-mesh and normal precision were used. The alloys were modeled by supercell containing 128 randomly distributed atoms. At initial step the atoms occupy perfect bcc lattice cites. This initial structure was optimized until energy changes less than 1e-6 eV, while forces acting on atoms don't exceed 1e-2 eV/angstrom. The electron-ion interaction is described by PAW pseudopotentials. The calculations have been collected by sampling chemical compositions across the entire compositional range. The chemical compositions have been sampled by progressively changing the number of atoms per constituent by 4. For each chemical composition of binaries and ternaries, the first-principle calculations have been run for 100 randomized arrangements of the constituents on the BCC lattice sites. We collected data for a total of 3,100 randomized atomic structures over 31 chemical compositions. The calculations have been collected on Air Force HPC11 cluster using the VASP 6.5.1. Additional methodology and file structure information is available in the dataset README.txt file.

36 MATERIALS SCIENCE↗
DOE OSTI2024

VZr_BCC_SolidSolution_128atoms_VASP6

We performed density functional theory (DFT) calculations for body-centered-cubic (BCC) structures with 128 lattices sites of solid solution binary alloys vanadium-zirconium (V-Zr). The electronic structures of alloys have been calculated using Vienna Ab initio Simulation Package (VASP). Within this package the DFT approach is used to reduce many-body Schrodinger equation to set of single particle Kohn-Sham (KS) equations. The generalized electronic exchange-correlation functional is described by generalized gradient approximation with the Perdew-Burke-Ernzerhof parametrization. The electron-ion interactions is described by pseudopotentials developed within the plane-wave basis projector augmented-wave (PAW) approach \cite{PAW}. These pseudopotentials are available at the VASP portal (http://cms.mpi.univie.ac.at/vasp/). Our calculations have been run with the pseudopotentials treating s and p semi-core states as valence in case for the elements V and Zr. The electronic densities and potentials are expanded over plane-waves with energy cutoff of 350 eV. 2x2x2 k-mesh and normal precision were used. The alloys were modeled by supercell containing 128 randomly distributed atoms. At initial step the atoms occupy perfect bcc lattice cites. This initial structure was optimized until energy changes less than 1e-6 eV, while forces acting on atoms don't exceed 1e-2 eV/angstrom. The electron-ion interaction is described by PAW pseudopotentials. The calculations have been collected by sampling chemical compositions across the entire compositional range. The chemical compositions have been sampled by progressively changing the number of atoms per constituent by 4. For each chemical composition of binaries and ternaries, the first-principle calculations have been run for 100 randomized arrangements of the constituents on the BCC lattice sites. We collected data for a total of 3,100 randomized atomic structures over 31 chemical compositions.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Theoretical Investigation of the Hydrogenation of Cyclohexene Catalyzed by Supported Single-Atom Sites on Redox Noninnocent LiMn 2 O 4 and Li 2 Mn 2 O 4 Surfaces

Here, the tuning of catalyst activity via stereoelectronic modulation of the active-site structure remains a grand challenge in heterogeneous catalysis. In homogeneous catalysis, the redox noninnocent ligands can be introduced to organometallic fragments to donate electrons to the metal center and fine-tune the catalytic activity. Analogously, lithium-ion battery materials, such as lithium manganese oxide (LMO), lithium titanium oxides (LTO), etc., can serve as redox noninnocent catalyst support to modulate the electronic structures of the active site via lithiation, hence tuning the catalytic activity of supported active sites. Experimentally, the Ni single-atom site was supported on LiMn 2 O 4 via oxidative grafting and exhibited no catalytic activity toward the hydrogenation of cyclohexene. After introducing additional Li into the catalyst support forming Li 2 Mn 2 O 4 , the Ni single-atom site becomes active, with the catalytic rates increasing as a function of the lithiation for Li/Mn ratios >0.9. In this paper, density functional theory (DFT) calculations are performed to study the Ni site structure via X-ray absorption near edge structure (XANES) simulations and investigate the electronic properties of Ni single-atom site before and after the addition of intercalated lithium in the LMO spinel structure. Furthermore, the study of the reaction mechanism is also carried out to understand the thermodynamically and kinetically favored pathways. XANES simulation suggests that the Ni single-atom site is likely to stay in the Li channel of the spinel support structure and form an octahedral structure. After Li intercalation, the Ni site becomes less positively charged, indicating the partial reduction of the Ni site. The simulated reaction energy profile over the LiMn 2 O 4 support exhibits high-energy barriers (1.07 eV) for the hydrogenation of cyclohexene; however, the Li 2 Mn 2 O 4 support is able to better stabilize low-coordinated ion sites and improve ion mobility, leading to lower overall energy barriers (0.64 eV). The reduced and low-coordinated ion site, thus, can better stabilize the reaction intermediates and promote hydrogenation reaction. Similarly, other transition metal ions (Fe, Co, and Cu) are also considered over the LiMn 2 O 4 and Li 2 Mn 2 O 4 catalyst supports for hydrogenation reaction.

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

Mechanistic Investigations of Gas-Phase Catalytic Hydrogenation in Metal–Organic Frameworks: Cooperative Activity of the Metal and Linker Sites in Cu x Rh 3– x (BTC) 2

We report the Cu x Rh 3–x (BTC) 2 catalyst (abbreviated CuRhBTC, BTC 3– = benzene tricarboxylate) provides excellent dispersion of active metal sites coupled with well-defined, robust structures for propylene hydrogenation reactions. This material therefore serves as a unique prototype for understanding catalytic activity in metal organic frameworks (MOFs). The mechanism of gas-phase hydrogenation at the bimetallic metal nodes of a MOF has been investigated in detail for the first time using in situ spectroscopy and diffraction experiments combined with density functional theory (DFT) calculations. The reaction occurs via a cooperative process in which the metal and linker sites play complementary roles; specifically, H 2 is dissociated at a Rh 2+ site with a missing Rh–O bond, while protonation of the decoordinated carboxylate linker stabilizes the active sites and promotes H 2 dissociation. In situ X-ray diffraction experiments show that the crystalline structure of the MOF is retained under reaction conditions at 20–100 °C. In situ Raman spectroscopy and diffuse reflectance infrared Fourier transform spectroscopy (DRIFTS) experiments demonstrate that propylene adsorbs at both Rh 2+ and Cu 2+ sites via π bonding. Cu 2+ is catalytically inactive, but at Rh 2+ sites, a propyl intermediate is observed when H 2 is introduced into the propylene feed. Furthermore, the appearance of the O–H stretch of COOH at ~3690 cm –1 in the DRIFT spectra is characteristic of defects consisting of missing Rh–O bonds. These experimental results are in general agreement with a reaction mechanism proposed by DFT, in which the decoordinated carboxylate linker is protonated, and the active Rh 2+ site remains available for readsorption of reactants in the subsequent catalytic cycle.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Ba 1−x Sr x FeO 3−δ as an improved oxygen storage material for chemical looping air separation: a computational and experimental study

Chemical looping air separation (CLAS) is a promising technology to generate oxygen-rich gas streams to enable efficient carbon dioxide capture during fossil fuel combustion or gasification. CLAS relies on the capture and release of oxygen from the atmosphere using the redox properties of an oxygen-selective solid oxide carrier. This study investigates the redox characteristics of Ba 1−x Sr x FeO 3−δ (0.0 ≤ x ≤ 0.417, 0.0 ≤ δ ≤ 0.5) using a combination of density functional theory (DFT) calculations and experimental verification using X-ray diffraction, thermogravimetric analysis, and oxygen-temperature-programmed desorption. The DFT computed energies of the Ba 1−x Sr x FeO 3−δ perovskites reveal a composition-dependent transition from hexagonal to cubic phases as the Sr-concentration or oxygen vacancy concentration increases. Oxygen vacancy formation energies of the cubic perovskites are found to be lower than those of their hexagonal counterparts. A low oxygen diffusion barrier of ∼1 eV combined with the thermodynamic preference of Ba 1−x Sr x FeO 3−δ compositions that form in a cubic phase suggests them as promising candidates for oxygen storage applications. The experimental results corroborate this finding by identifying Ba 0.75 Sr 0.25 FeO 3−δ in the cubic phase as an optimal composition offering low-temperature oxygen storage capacities comparable to that of the state-of-the-art Sr 0.75 Ca 0.25 FeO 3−δ perovskite oxygen storage material at 325 °C and 350 °C.

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

Ab initio thermodynamics of Ni and Co incorporation in Mg hydroxide, carbonate, and hydroxycarbonate minerals

Ni and Co are critical elements needed for modern technologies, and a better understanding of the ability of Mg-based minerals to incorporate these elements would benefit strategy development for Ni and Co recovery from mafic and ultramafic deposits. Here, in this work, we performed density functional theory (DFT) calculations of Ni and Co incorporation in six potential products of the carbonation of mafic and ultramafic silicates: brucite (Mg(OH) 2 ), magnesite (MgCO 3 ), nesquehonite (MgCO 3 ⸱3H 2 O), lansfordite (MgCO 3 ⸱5H 2 O), artinite (Mg 2 CO 3 (OH) 2 ⸱3H 2 O), and hydromagnesite (Mg 5 (CO 3 ) 4 (OH) 2 ⸱4H 2 O). The DFT results were used in an ab initio thermodynamics framework to explore the pH 2 O–pCO 2 conditions at which the Mg-based minerals were predicted to be thermodynamically stable and to quantify the Gibbs free energy of Ni and Co substitution at Mg sites. Among the six Mg-based minerals, brucite and magnesite were predicted to have the lowest Ni and Co substitution free energy. An analysis of the effect of temperature indicated that, at low temperature (<100 K), brucite more readily accommodated Ni and Co, while, at higher temperature (>335 K), magnesite more favorably incorporated Ni and Co. Between 100 K and 335 K, Ni was predicted to preferentially substitute for Mg in brucite and Co for Mg in magnesite, thus leading to a driving force for separating Ni and Co in conditions where brucite and magnesite both form. Insights gained in this work could therefore help select experimental conditions that either promote or inhibit incorporation of these critical elements into Mg-based mineral phases.

Critical elements↗
DOE OSTI2024

High-entropy Li-rich layered oxide cathode for Li-ion batteries

High-entropy oxides (HEOs) are emerging as promising cathode materials for Li-ion batteries (LIBs) due to their stable solid-state phase and compositional flexibility. Herein, we investigate the structural and electrochemical properties of a novel non-equimolar high-entropy cathode material, termed high-entropy Li-rich layered oxide (HE-LLO, Li 1.15 Na 0.05 Ni 0.19 Mn 0.56 Fe 0.02 Mg 0.02 Al 0.02 O 1.97 F 0.03 ), in comparison to a pristine Li-rich layered oxide (PR-LLO, Li1.2Ni0.2Mn0.6O2). The incorporation of multiple cations (Na + , Al 3+ , Mg 2+ , Fe 3+ ) and anion (F - ) into HE-LLO introduces compositional diversity, enhancing structural stability through the entropy stabilization effect. Theoretical calculations confirm a significantly higher configurational entropy in HE-LLO compared to PRLLO, supporting its high-entropy nature. Electrochemical evaluations demonstrate that HE-LLO exhibits considerable capacity retention, preserving 76.8 % of its discharge capacity at 0.5C after 200 cycles, compared to only 36.2 % for PR-LLO. Even under high-temperature conditions, HE-LLO outperformed PR-LLO, maintaining 76.1 % of its discharge capacity after 100 cycles at 5C, while PR-LLO retained only 12.4 %. These enhancements are attributed to the improved phase reversibility and higher Li + ion diffusion coefficients of HE-LLO, validated by ex-situ characterizations using a synchrotron X-ray technique, along with density functional theory (DFT) calculations. In conclusion, these findings highlight the promise of non-equimolar HEOs as a novel design strategy for highperformance cathode materials.

25 ENERGY STORAGE↗
DOE OSTI2023

Utilizing metal cation reactions with carbonyl sulfide to remove isobaric interferences in tandem inductively coupled plasma mass spectrometry analyses

Carbonyl sulfide (OCS) was used as a reaction gas to investigate gas phase metal (M + ) ion-molecule reactions using the Agilent 8900 inductively coupled plasma tandem mass spectrometer (ICP-MS/MS) to yield insight on how this gas may be used to remove isobaric interferences in analytical measurements. The experimental work was paired with density functional theory (DFT) calculations of the reaction enthalpy to predict whether M + will react with OCS. A multi-element standard containing 46 elements ranging from 9 to 208 u was analyzed in the presence and absence of OCS. When a reaction was observed, the dominant product was the sulfide (MS + ). Oxide products were also observed for many M + but formation was less efficient with OCS than previously observed with other reaction gases. This is likely due to the weaker OC-S bond that makes MS + formation more favorable. Increasing the flow rate from 0.1 to 0.2 mL/min (corresponding to a change in reaction gas pressure from 0.35 to 0.53 Pa (2.6 to 4.0 mTorr)) generally resulted in greater MS + production, including the secondary product MS 2 + for a few cations. The early lanthanide series ions (La + , Ce + , Pr + and Nd + ) produced greater quantities of MO + at the higher pressure, although MS + products were still the dominant product. The DFT-predicted reaction enthalpies were consistent with the observed sulfide formation, with an accuracy >90%; however, model predictions were less accurate for the minor and higher order products (< 77% for MO + ). Finally, the work presented here continues a systematic study of ion-molecule reactions in ICP-MS/MS to understand and develop new and novel ways to analyze complex mixtures with minimal pre-analysis treatment.

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

Structural anisotropy in Sb thin films

Sb thin films have attracted wide interest due to their tunable band structure, topological phases, high electron mobility, and thermoelectric properties. We successfully grow epitaxial Sb thin films on a closely lattice-matched GaSb(001) surface by molecular beam epitaxy. We find a novel anisotropic directional dependence on their structural, morphological, and electronic properties. The origin of the anisotropic features is elucidated using first-principles density functional theory (DFT) calculations. The growth regime of crystalline and amorphous Sb thin films was determined by mapping the surface reconstruction phase diagram of the GaSb(001) surface under Sb 2 flux, with confirmation of structural characterizations. Crystalline Sb thin films show a rhombohedral crystal structure along the rhombohedral (211) surface orientation parallel to the cubic (001) surface orientation of the GaSb substrate. At this coherent interface, Sb atoms are aligned with the GaSb lattice along the [1̄10] crystallographic direction but are not aligned well along the [110] crystallographic direction, which results in anisotropic features in reflection of high-energy electron diffraction patterns, misfit dislocation formation, surface morphology, and transport properties. Our DFT calculations show that the preferential orientation of the rhombohedral Sb (211) plane may originate from the GaSb surface, where Sb atoms align with the Ga and Sb atoms on the reconstructed surface. The formation energy calculations confirm the stability of the experimentally observed structures. Our results provide optimal film growth conditions for further studies of novel properties of Bi 1-x Sb x thin films with similar lattice parameters and an identical crystal structure, as well as functional heterostructures of them with III–V semiconductor layers along the (001) surface orientation, supported by a theoretical understanding of the anisotropic film orientation.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Adsorbate-induced adatom formation on Au-Cu bimetallic alloys and its possible consequences for CO 2 electroreduction

The adsorbate-induced formation of sub-nanometer clusters on transition-metal single crystals observed in previous high-pressure microscopic studies hinted at the in-situ formation of unique active sites even on large nanoparticle catalysts. We propose that the adatom formation energy can be used as an energetic descriptor for the initial step toward the adsorbate-induced metal-cluster formation process. This descriptor can be efficiently computed using density functional theory (DFT) calculations and applied for screening and identification of metal catalysts where this phenomenon may play an important role in generating active sites in-situ. As a proof of concept, here, we construct an adatom formation energy database for three Au x Cu y alloys (x:y = 3:1, 1:1, or 1:3) and eighteen adsorbates (H, C, N, O, F, S, Cl, Br, I, CH x , NH x (x = 1 – 3), CO, NO, and OH) commonly involved in catalytic reactions. The energetics of adatom formation were examined in all cases where the (111) terrace, (211) step-edge, and (874) kink were the sources of the adatom. We demonstrate that the presence of an adsorbate could alter not only the energetics for adatom formation but also the elemental nature of the preferred adatom being formed. Using our database, we identified promising systems which favor adsorbate-induced adatom formation under near-ambient conditions. Specifically, CO-induced adatom formation on all three Au-Cu alloy surfaces could occur under CO 2 electroreduction (CO 2 RR) conditions. This phenomenon offers a qualitative explanation for the experimentally observed CO 2 RR activity on Au-Cu alloy catalysts. As a result, our methodology offers an easily expandable and efficient approach for large-scale catalyst screening with regards to adatom/cluster formation under reaction conditions and provides insight into the possible nature of active sites on alloy catalysts from a novel perspective.

Active site↗
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 OSTI2023

Lattice dynamics and thermodynamics for δ -plutonium from density functional theory

Here, we present results from density functional theory (DFT) calculations of the lattice dynamics (phonons) and thermodynamics for δ-phase plutonium. The fully relativistic electronic structure is calculated assuming a three-dimensional noncollinear magnetic structure in conjunction with DFT and the general gradient approximation for the electron exchange and correlation interactions. The electronic-structure model is further enhanced by addressing strong orbital-orbital coupling via the conventional orbital-polarization (OP) scheme as has been successfully done for plutonium. The temperature dependence of the phonons is calculated within the self-consistent ab initio lattice dynamics approach. The obtained phonons compare very well with measurements although a modest overestimation of the transverse L-point [ξξξ] phonon is acknowledged. Calculated thermal vibration amplitudes and the associated Debye-Waller temperatures are close to experiments. Lattice, electronic, and magnetic contributions to the heat capacity are predicted and consistent to a few percent with that deduced from experimental data. Good agreement is only achieved when a magnetic contribution to the specific heat is recognized. The parameter-free DFT+OP electronic model is thus capable of predicting phonon properties and thermodynamic behavior of δ-phase plutonium rather accurately.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Ambient pressure x-ray photoelectron spectroscopy study of water formation and adsorption under two-dimensional silica and aluminosilicate layers on Pd(111)

Ambient pressure x-ray photoelectron spectroscopy (AP-XPS) supported by density functional theory (DFT) calculations was used to characterize the interaction of water with two-dimensional (2D) silica and aluminosilicate bilayers on Pd(111). Starting with oxygen adsorbed at the SiO2/Pd interface, exposure to water caused the SiO 2 -derived XPS peaks to shift to higher binding energy and the removal of an O 1s feature associated with interfacial adsorbed oxygen. These observations were attributed to the formation of a mixed water–hydroxyl interface, which eliminates the interfacial dipolar layer, and its associated electrostatic potential, created by adsorbed oxygen. Interfacial oxygen also reacted with H 2 to produce adsorbed water which also caused an upward binding energy shift of the SiO 2 peaks. Spectra recorded under 0.5 Torr water revealed additional water adsorption and a further shift of the overlayer peaks to higher binding energy. Incorporating Al into the 2D material caused the bilayer peaks to shift to lower binding energy which could be explained by electron donation from the metal to the bilayer. Although the stronger interaction between the bilayer and Pd substrate should restrict interfacial adsorption and reaction, similar trends were observed for water and hydrogen exposure to interfacial adsorbed oxygen. Less water adsorption was observed at the aluminosilicate interface which is a consequence of Al strengthening the bond to the metal substrate. The results reveal how the sensitivity of XPS to interfacial dipoles can be exploited to distinguish reactions taking place in confined spaces under 2D layers and how tuning the composition of the 2D layer can impact such reactions.

Density functional theory↗
DOE OSTI2024

Field-induced magnetization of δ -plutonium from density functional theory

Here, we present results from density functional theory (DFT) calculations of magnetization, induced by an external magnetic field, for δ-phase plutonium. The fully relativistic electronic structure accounts for Zeeman splitting effects through a Hamiltonian that couples the magnetic field to both spin and orbital magnetic moments. The electronic-structure model is further improved by an extension to DFT in terms of an orbital-orbital coupling via the conventional orbital-polarization method, as has routinely been done for plutonium. The response to the applied magnetic field is shown to be weak in the DFT model, with induced magnetic moments of the order of in magnetic fields up to 30 T. These results are in accord with recent assessments from x-ray magnetic circular dichroism measurements in magnetic field on δ-plutonium. Somewhat surprisingly, a model assuming δ-plutonium to be absent magnetic moments shows stronger response to an external static magnetic field than models allowing for antiferromagnetism or magnetic disorder.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Gaussian approximation potentials for accurate thermal properties of two-dimensional materials

Two-dimensional materials (2DMs) continue to attract a lot of attention, particularly for their extreme flexibility and superior thermal properties. Molecular dynamics simulations are among the most powerful methods for computing these properties, but their reliability depends on the accuracy of interatomic interactions. While first principles approaches provide the most accurate description of interatomic forces, they are computationally expensive. In contrast, classical force fields are computationally efficient, but have limited accuracy in interatomic force description. Machine learning interatomic potentials, such as Gaussian Approximation Potentials, trained on density functional theory (DFT) calculations offer a compromise by providing both accurate estimation and computational efficiency. Here, in this work, we present a systematic procedure to develop Gaussian approximation potentials for selected 2DMs, graphene, buckled silicene, and h-XN (X = B, Al, and Ga, as binary compounds) structures. We validate our approach through calculations that require various levels of accuracy in interatomic interactions. The calculated phonon dispersion curves and lattice thermal conductivity, obtained through harmonic and anharmonic force constants (including fourth order) are in excellent agreement with DFT results. HIPHIVE calculations, in which the generated GAP potentials were used to compute higher-order force constants instead of DFT, demonstrated the first-principles level accuracy of the potentials for interatomic force description. Molecular dynamics simulations based on phonon density of states calculations, which agree closely with DFT-based calculations, also show the success of the generated potentials in high-temperature simulations.

2D Materials↗
DOE OSTI2025

Predicting interface structure using the minima hopping method

Here, we adapt the minima hopping method (MHM) to the problem of interfacial structure prediction and apply it to study a canonical problem, the tilt grain boundaries in SrTiO 3 . Our method employs a hybrid approach by first exploring the potential energy surface (PES) of different grain boundary samplings with an empirical force field, among which the fifteen candidates with lower energies are then refined using ab initio density functional theory (DFT) calculations. During the exploratory stage, we bias the search using a local order parameter to primarily sample various reconstructions in the vicinity of the interface, while preserving the crystallinity of the bulk regions. We further enhance the search by incorporating initial structures with rigid body displacements to account for translational variations between bulk phases, enabling the MHM to effectively generate both stoichiometric and nonstoichiometric SrTiO 3 Σ⁢3(111)[110] and Σ⁢3(112)[110] grain boundaries. From an algorithmic standpoint, MHM outperforms earlier studies based on genetic algorithms (GA) by identifying more stable interfacial structures of several SrTiO 3 grain boundaries. The performance of the present implementation of the MHM approach is primarily limited by exploring an approximate description of the PES with a rather simple Buckingham potential. This limitation leads to variations in performance when compared to approaches utilizing more advanced surrogate PES models, such as direct DFT-PES sampling or GA with the embedded atom method (EAM). Despite the present limitations, the MHM approach is able to yield interfacial structures with comparable or lower interfacial energies in specific cases, such as Σ⁢3(111)[110] Γ=1, ±0.5 and Σ⁢3(112)[110] Γ= ±1, −2, underscoring the robustness of the MHM approach even with a simple approximation of the DFT PES. The MHM interfacial structure prediction method thus offers an efficient approach to understanding the grain boundaries and heterointerfaces at the atomic scale, providing an important prerequisite for effective materials design.

density functional theory↗