Search NASA⌕ Search

SEARCH · Search NASA

Results for “Density functional theory simulations”

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

Multiple incommensurate magnetic states in the kagome antiferromagnet Na 2 Mn 3 Cl 8

The kagome lattice can host exotic magnetic phases arising from frustrated and competing magnetic interactions. However, relatively few insulating kagome materials exhibit incommensurate magnetic ordering. Here, we present a study of the magnetic structures and interactions of antiferromagnetic Na 2 Mn 3 Cl 8 with an undistorted Mn 2+ kagome network. Using neutron-diffraction and bulk magnetic measurements, we show that Na 2 Mn 3 Cl 8 hosts two different incommensurate magnetic states, which develop at T N1 = 1.6 K and T N2 = 0.6 K. Magnetic Rietveld refinements indicate magnetic propagation vectors of the form q=(q x , q y , $\frac{3}{2}$), and our neutron-diffraction data can be well described by cycloidal magnetic structures. By optimizing exchange parameters against magnetic diffuse-scattering data, we show that the spin Hamiltonian contains ferromagnetic nearest-neighbor and antiferromagnetic third-neighbor Heisenberg interactions, with a significant contribution from long-ranged dipolar coupling. This experimentally determined interaction model is compared with density-functional-theory simulations. Using classical Monte Carlo simulations, we show that these competing interactions explain the experimental observation of multiple incommensurate magnetic phases and may stabilize multi-q states. Here, our results expand the known range of magnetic behavior on the kagome lattice.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Anomalous stability of non-van der Waals bonded B4C nanosheets through surface reconstruction

Boron carbide (B4C) has been well studied both theoretically and experimentally in its bulk form due to its exceptional hardness and use as a high-temperature thermoelectric. However, the properties of its two-dimensional nanosheets are not well established. In this paper, using van der Waals-corrected density-functional theory simulations, we show that bulk B4C can be cleaved along different directions to form B4C nanosheets with low formation energies. We find that there is minimal dependence of formation energies on cleavage planes and surface terminations, even though the bulk is not van der Waals layered. This anomalous stability of B4C nanosheets is found to be a result of surface reconstructions that are unique to B-rich systems. While the density of states of the bulk B4C indicate that it is a semiconductor, the B4C nanosheets are found to be predominantly metallic. We attribute this metallic behavior to a redistribution of charges on the surface bonds of the films. The Seebeck coefficients of the B4C films remain comparable to those of the bulk and are nearly constant as a function of temperature. Our results provide guidance for experimental synthesis efforts and future application of B4C nanosheets in nanoelectronic and thermoelectric applications.

Physics↗
DOE OSTI2020

Amorphous Ag 2-x Cu x S quantum dots: “all-in-one” theranostic nanomedicines for near-infrared fluorescence/photoacoustics dual-modal-imaging-guided photothermal therapy

Integrating dual-modal imaging and photothermal effect within a single nanosystem is an enormous challenge for current nanomedicine application. In this study, non-toxic amorphous Ag 2-x Cu x S quantum dots (QDs) possessing near-infrared fluorescence (NIRF), enhanced photothermal and photoacoustic (PA) performance, as “one-in-all” nanomedicines, are designed and obtained. Cu doping in amorphous Ag 2 S QDs causes red-shift of NIRF to the lowest tissue absorption wavelength in so called “first biowindow” (~820 nm) and improves photothermal conversion efficiency to 44.0% simultaneously. Density functional theory simulation results suggest the fluorescence red-shift and enhanced photothermal/photoacoustics properties can be attributed to the generation of intragap states introduced by Cu doping. The amorphous Ag 2-x Cu x S QDs exhibit high long-term biocompatibility. Under a single laser irradiation with relatively low power density, an effective tumor ablation induced by high photothermal effect of amorphous Ag 2-x Cu x S QDs is achieved in vivo. Simultaneously, the PTT process can be guided by high quality NIRF/PA dual-modal-imaging.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Design of enhanced porous organic cage solubility in Type 2 porous liquids

Decreasing cost of technologies for direct air capture of carbon can be achieved through the design of new materials with high CO 2 selectivity that can be incorporated into existing industrial processes. An emerging class of materials for these applications are porous liquids (PLs). PLs are mixtures of porous hosts and solvents with intrinsic porosity due to steric exclusion of solvent from inside the porous host. It is currently unknown how solvent -porous host interactions affect porous host solubility in the bulk solvent. Here, in this study, density functional theory simulations were used to investigate interactions between nine solvents and a CC13 porous organic cage (POC). Calculations identified that solvent molecules were the most stable when placed either inside the CC13 POC or in the pore window compared to interfacial binding sites. Structural changes to the CC13 POC correlated with reported experimental solubilities, including expansion of the CC13 POC with solvent molecule infiltration and expansion or contraction of the pore window. Based on these results, new PL design guidelines should include compositions with (1) high concentrations of POCs with flexible cage structures that can expand when solvated and (2) solvent molecule-POC combinations that contract the pore window during solvent molecule-host binding.

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

CH 4 Activation over Perovskite Catalysts: True Density and Reactivity of Active Sites

The high thermal stability of perovskites has drawn attention toward their applications for catalytic CH 4 activation and conversion, typically occurring at high temperatures. The reaction rates of perovskite catalysts for CH 4 combustion, however, trail behind those of noble metal catalysts. Ways to optimize the performance of perovskite catalysts are destined to trial-and-error approaches unless their complex reconstructed surfaces are correlated with fundamental kinetic parameters. Discerning the intrinsic activity of surface catalytic sites and the density of those sites is crucial to rationally envision complex metal oxides with enhanced catalytic performance. Here, the present work presents a detailed kinetic analysis of catalytic CH 4 combustion over a set of seven perovskites (SrTiO 3 , SrZrO 3 , SrFeO 3 , LaFeO 3 , LaInO 3 , LaCoO 3 , LaMnO 3 ) with various surface terminations. Steady-state isotopic transient kinetic analysis was employed to measure turnover frequency (TOF) and density of surface intermediates (N) under operando conditions. Top surface characterization elucidated performance-structure relationships between near-monolayer surface composition and intrinsic reactivity of the catalysts. By using a chemical etching procedure to expose Fe-sites at the top surface of LaFeO 3 (LaFeO 3 , HNO 3 ), its TOF was increased 4-fold, compared with the unmodified sample, although N on the surface of LaFeO 3 , HNO 3 decreased. Density functional theory simulations corroborated that surface Fe-termination and La-Fe termination offer lower energetic barriers for CH 4 activation when compared with La-termination. In general, surface reconstruction is shown as a tool to tune TOF and N to improve reaction rates. This work fills a gap in current kinetic studies of perovskites through a careful assessment and discussion of the density and intrinsic reactivity of active sites for methane combustion over well-characterized reconstructed perovskite surfaces.

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

Integrating Machine Learning Potential and X-ray Absorption Spectroscopy for Predicting the Chemical Speciation of Disordered Carbon Nitrides

Precise determination of atomic structural information in functional materials holds transformative potential and broad implications for emerging technologies. Spectroscopic techniques, such as X-ray absorption near-edge structure (XANES), have been widely used for material characterization; however, extracting chemical information from experimental probes remains a significant challenge, particularly for disordered materials. We present an integrated approach that combines atomic simulations, data-driven techniques, and experimental measurements to investigate chemical speciation of amorphous carbon nitride systems as a case study. Here, we discuss the development of machine learning potentials that can efficiently explore the vast configuration space of amorphous carbon nitrides. By employing statistical methods, this structural database enables the elucidation of the most representative local structures and how they evolve with chemical compositions and density. Density functional theory simulations are used to establish a correlation between the local structure and spectroscopic signatures, which then serve as the basis for interpreting and extracting chemical content from experimental data. Although our framework is specifically demonstrated for XANES and carbon nitrides, the approach described herein is readily adaptable as applied to other experimental characterization probes and materials classes.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Improved Kelbg Potentials for Z > 1 and Application to Carbon Plasmas

In this work, we present a general form for the electron‐ion diffractive potential derived from the quantum pair density matrix and fit to the improved Kelbg potential for atomic numbers up to $Z = 54$. We apply classical molecular dynamics using the improved Kelbg potential for carbon with various forms of the Pauli potential to compute internal energies and pressures for hot, dense plasma conditions. Our results are compared to an equation of state model based on path integral Monte Carlo and density functional theory simulations to examine the extent to which the improved Kelbg potential reproduces the internal energy and pressure of carbon plasmas. The regions of validity for carbon agree generally with those derived previously for hydrogen once pressure ionization effects are incorporated. Based on our carbon results and previously published hydrogen studies, we discuss the general applicability and limitations of these potentials for equation of state studies in warm dense matter and high energy density plasmas.

general physics↗
DOE OSTI2023

Ab Initio Simulations of Tritium Diffusion in Al-Rich Iron Aluminide Coating Phases

Surface coatings of steels used in extreme conditions and corrosive environments generally aim to provide protection and increased durability. In the case of tritium-producing burnable absorber rods (TPBARs) used in nuclear reactors, a 316 stainless steel has been coated with Al. Scanning transmission electron microscopy (STEM) characterization of the coating found three Al-rich (>60 atom % Al) iron aluminide alloys identified as hexagonal FeNiAl 5 , monoclinic Fe 4 Al 13 , and orthorhombic Fe 2 Al 5 . Density functional theory simulations using nudged elastic band have been performed to investigate the diffusion of interstitial tritium in each Al-rich iron aluminide phase. While FeNiAl 5 and Fe 4 Al 13 can be viewed as the stacking of two layers, the structural peculiarity of Fe 2 Al 5 is that channels of variable Al vacancy content are present along the c-axis. Therefore, three stoichiometries for Fe 2 Al x phase, namely, Fe 2 Al 4 , Fe 2 Al 5 , and Fe 2 Al 6 , have been considered to evaluate the impact of Al vacancy concentration on tritium diffusion behavior. Altogether, we found that at 600 K, tritium diffusion decreases from a faster rate in the channels of Fe 2 Al x phases (D T ≤ 10 –11 m 2 ·s –1 ) to Fe 4 Al 13 (D T ≈ 10–12 m2·s–1), and finally in FeNiAl5 (D T ≈ 10 –13 m 2 ·s –1 ). Here, we also find that interstitial tritium generally diffuses faster in Fe–Al coating phases than in the tritium breeding material γ-LiAlO 2 (D T ≈ 10 –14 m 2 ·s –1 ) but slightly slower than in 316 stainless steel (D T ≈ 10 –10 m 2 ·s –1 ).

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

Complex Oxides under Simulated Electric Field: Determinants of Defect Polarization in AB O 3 Perovskites

Abstract Polarization of ionic and electronic defects in response to high electric fields plays an essential role in determining properties of materials in applications such as memristive devices. However, isolating the polarization response of individual defects has been challenging for both models and measurements. Here the authors quantify the nonlinear dielectric response of neutral oxygen vacancies, comprised of strongly localized electrons at an oxygen vacancy site, in perovskite oxides of the form AB O 3 . Their approach implements a computationally efficient local Hubbard U correction in density functional theory simulations. These calculations indicate that the electric dipole moment of this defect is correlated positively with the lattice volume, which they varied by elastic strain and by A‐site cation species. In addition, the dipole of the neutral oxygen vacancy under electric field increases with increasing reducibility of the B‐site cation. The predicted relationship among point defect polarization, mechanical strain, and transition metal chemistry provides insights for the properties of memristive materials and devices under high electric fields.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Ferroelectric Fractals: Switching Mechanism of Wurtzite AlN

The advent of wurtzite ferroelectrics is enabling new ferroelectric devices for computer memory that have the potential to bypass the von Neumann bottleneck due to their robust polarization and silicon compatibility. However, the atomistic switching mechanism of wurtzites is still undetermined due to the limitations of density functional theory simulation size and experimental temporal and spatial resolution. Thus, physics-informed materials engineering to reduce coercive field and breakdown in these devices has been limited. In this work, the atomistic mechanism of domain wall migration and domain growth in aluminum nitride-based wurtzites is uncovered using molecular dynamics and Monte Carlo simulations. We reveal the anomalous switching mechanism of fast 1D single columns of atoms propagating from a slow-moving 2D fractallike domain wall. We find that the critical nucleus is a single aluminum ion that breaks its bond with one nitrogen and bonds to another nitrogen; this creates a cascade that flips atoms directly only in the same column, due to the extreme locality (sharpness) of the domain walls in wurtzites. We further show how the fractallike shape of the domain wall in the 2D plane breaks assumptions in the Kolmogorov, Avrami, and Ishibashi (KAI) model and leads to the anomalously fast switching in wurtzite structured ferroelectrics.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Atomic cluster expansion potential for large scale simulations of hydrocarbons under shock compression

We present an Atomic Cluster Expansion (ACE) machine learned potential developed for high-fidelity atomistic simulations of hydrocarbons, targeting pressures and temperatures near and above supercritical fluid regimes for molecular fluids. A diverse set of stoichiometries were covered in training, including 1:0 (pure carbon), 1:4 (methane), and 1:1 (benzene), and rich bonding environments sampled at supercritical temperatures, hydrogen rich, reactive mixtures where metastable stoichiometries arise, including 1:2 (ethylene) and 1:3 (ethane). A high-fidelity training database was constructed by performing large-scale quantum molecular dynamic simulations [density functional theory (DFT) MD] of diamond, graphite, methane, and benzene. A novel approach to selecting structures from DFT MD is also presented, which allows for the rapid selection of unique DFT MD frames from complex trajectories. Comparisons to DFT and experimental data demonstrate that the presented ACE potential accurately reproduces isotherms, carbon melting curves, radial distribution functions, and shock Hugoniots for carbon and hydrocarbon systems for pressures up to 100 GPa and temperatures up to 6000 K for hydrocarbon systems and up to 9000 K for pure carbon systems. This work delivers a potential that can be used for accurate, large-scale simulations of shocked hydrocarbons and demonstrates a methodology for fitting and validating machine learning interatomic potentials to complex molecular environments, which can be applied to energetic materials in future works.

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

Predicted thermophysical properties of UN, PuN, and (U,Pu)N

Molecular dynamics and density functional theory simulations are used to predict the lattice and electronic contributions of thermophysical properties for UN, PuN, and mixed (U,Pu)N systems. The properties predicted include the lattice parameter, linear thermal expansion, enthalpy, and specific heat capacity, as a function of temperature. The simulation predictions for high temperature specific heat capacity are compared against experimental measurements to understand the behavior, and why differences in the experimental measurements are observed. The influence of adding U vacancies, N interstitials, and Pu to UN is also examined. For this, a new PuN potential parameter set is developed and used with the Kocevski UN potential, enabling the dynamics of mixed (U,Pu)N systems to be studied. How defects impact the thermophysical properties is important for understanding fuel behavior under different reactor conditions, and these mechanistic predictions can be used to support fuel performance codes where data is scarce.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Bayesian blacksmithing: discovering thermomechanical properties and deformation mechanisms in high-entropy refractory alloys

Finding alloys with specific design properties is challenging due to the large number of possible compositions and the complex interactions between elements. This study introduces a multi-objective Bayesian optimization approach guiding molecular dynamics simulations for discovering high-performance refractory alloys with both targeted intrinsic static thermomechanical properties and also deformation mechanisms occurring during dynamic loading. The objective functions are aiming for excellent thermomechanical stability via a high bulk modulus, a low thermal expansion, a high heat capacity, and for a resilient deformation mechanism maximizing the retention of the BCC phase after shock loading. Contrasting two optimization procedures, we show that the Pareto-optimal solutions are confined to a small performance space when the property objectives display a cooperative relationship. Conversely, the Pareto front is much broader in the performance space when these properties have antagonistic relationships. Density functional theory simulations validate these findings and unveil underlying atomic-bond changes driving property improvements.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Electronically induced defect creation at semiconductor/oxide interface revealed by time-dependent density functional theory

Carrier induced defect creation at the semiconductor-oxide interface has been known as the origin of electronic device degradation for a long time, but how exactly the interface lattice can be damaged by carriers (especially low-energy ones) remains unclear. Here we carry out real-time time-dependent density functional theory simulations on concrete Si/SiO 2 interfaces to study the interaction between excited electrons and interface bonds. We show that the normal interface Si-H bonds are generally resistant to electrons due to the delocalized nature and high energy level of the Si-H antibonding states, and due to the high-energy barrier to break the Si-H bond. However, if an additional hydrogen atom exists by attaching to a nearby oxygen atom (forming a “Si-H···H-O” complex), the Si-H bond will be greatly weakened, including the reduction of energy barrier for bond breaking, and the lowering of the antibonding state energy level which favors electron injection. Together with the multiple vibrational excitation process, the corresponding Si-H bond can be broken much more easily. Thus we propose that the Si-H···H-O complex will be the center for defect creation and device degradation. Finally, we also explain why such a center might be relatively easy to form during the hydrogen annealing process.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Cu Based Dilute Alloys for Tuning the C 2+ Selectivity of Electrochemical CO 2 Reduction

Electrochemical CO 2 reduction is a promising technology for replacing fossil fuel feedstocks in the chemical industry but further improvements in catalyst selectivity need to be made. So far, only copper-based catalysts have shown efficient conversion of CO 2 into the desired multi-carbon (C 2+ ) products. This work explores Cu-based dilute alloys to systematically tune the energy landscape of CO 2 electrolysis toward C 2+ products. Selection of the dilute alloy components is guided by grand canonical density functional theory simulations using the calculated binding energies of the reaction intermediates CO*, CHO*, and OCCO* dimer as descriptors for the selectivity toward C 2+ products. Here, a physical vapor deposition catalyst testing platform is employed to isolate the effect of alloy composition on the C 2+ /C 1 product branching ratio without interference from catalyst morphology or catalyst integration. Six dilute alloy catalysts are prepared and tested with respect to their C 2+ /C 1 product ratio using different electrolyzer environments including selected tests in a 100-cm 2 electrolyzer. Consistent with theory, CuAl, CuB, CuGa and especially CuSc show increased selectivity toward C 2+ products by making CO dimerization energetically more favorable on the dominant Cu facets, demonstrating the power of using the dilute alloy approach to tune the selectivity of CO 2 electrolysis.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Antiferromagnetic Inverse Spinel Oxide LiCoVO 4 with Spin-Polarized Channels for Water Oxidation

Exploring highly efficient catalysts for the oxygen evolution reaction (OER) is essential for water electrolysis. Cost-effective transition-metal oxides with reasonable activity are raising attention. Recently, OER reactants' and products' differing spin configurations have been thought to cause slow reaction kinetics. Catalysts with magnetically polarized channels could selectively remove electrons with opposite magnetic moment and conserve overall spin during OER, enhancing triplet state oxygen molecule evolution. Herein, antiferromagnetic inverse spinel oxide LiCoVO 4 is found to contain d 7 Co 2+ ions that can be stabilized under active octahedral sites, possessing high spin states S = 3/2 (t 2g 5 e g 2 ). With high spin configuration, each Co 2+ ion has an ideal magnetic moment of 3 µ B , allowing the edge-shared Co 2+ octahedra in spinel to be magnetically polarized. Density functional theory simulation results show that the layered antiferromagnetic LiCoVO 4 studied contains magnetically polarized channels. The average magnetic moment (µ ave ) per transition-metal atom in the spin conduction channel is around 2.66 µ B . Such channels are able to enhance the selective removal of spin-oriented electrons from the reactants during the OER, which facilitates the accumulation of appropriate magnetic moments for triplet oxygen molecule evolution. Finally, the LiCoVO 4 reported has been identified as an oxide catalyst with excellent OER activity.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Additive‐Free, Low‐Temperature Crystallization of Stable α‐FAPbI 3 Perovskite

Abstract Formamidinium lead triiodide (FAPbI 3 ) is attractive for photovoltaic devices due to its optimal bandgap at around 1.45 eV and improved thermal stability compared with methylammonium‐based perovskites. Crystallization of phase‐pure α‐FAPbI 3 conventionally requires high‐temperature thermal annealing at 150 °C whilst the obtained α‐FAPbI 3 is metastable at room temperature. Here, aerosol‐assisted crystallization (AAC) is reported, which converts yellow δ‐FAPbI 3 into black α‐FAPbI 3 at only 100 °C using precursor solutions containing only lead iodide and formamidinium iodide with no chemical additives. The obtained α‐FAPbI 3 exhibits remarkably enhanced stability compared to the 150 °C annealed counterparts, in combination with improvements in film crystallinity and photoluminescence yield. Using X‐ray diffraction, X‐ray scattering, and density functional theory simulation, it is identified that relaxation of residual tensile strains, achieved through the lower annealing temperature and post‐crystallization crystal growth during AAC, is the key factor that facilitates the formation of phase‐stable α‐FAPbI 3 . This overcomes the strain‐induced lattice expansion that is known to cause the metastability of α‐FAPbI 3 . Accordingly, pure FAPbI 3 p–i–n solar cells are reported, facilitated by the low‐temperature (≤100 °C) AAC processing, which demonstrates increases of both power conversion efficiency and operational stability compared to devices fabricated using 150 °C annealed films.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Achieving Electrode Smoothing by Controlling the Nucleation Phase of Metal Deposition Through Polymer‐Substrate Binding

Polymer additives [like polyethylene oxide (PEO)] are widely used for smooth electrode deposition in aqueous zinc and many other battery systems. However, the precise mechanism by which they regulate morphology and suppress dendrite formation remains unclear. In this study, the knowledge gap is addressed by using in situ electrochemical atomic force microscopy to directly observe the interfacial evolution during Zn electrodeposition and polymer adsorption on Cu substrates in the presence of varying concentrations of ZnSO 4 and PEO. Contrary to previous literature assumptions, which emphasize the binding to the growing Zn crystal surfaces or Zn 2+ ions, the results demonstrate that PEO smooths Zn films by promoting nucleation of (002)-oriented Zn platelets through interactions with the Cu substrate. Density functional theory simulations support this finding by showing that PEO adsorption on Cu modifies the interfacial energy of Zn/Cu/electrolyte interfaces, favoring the stabilization of Zn (002) on the Cu substrate, as well as confines Zn electrodeposition to a narrow near-surface region. In conclusion, these findings elucidate a novel design principle for electrode smoothing, emphasizing the importance of substrate selection paired with polymer additives that exhibit an attractive interaction with the substrate but minimal interaction with growing crystals, offering a mechanistic perspective for improved battery performance.

Electrodeposition↗