Search NASA⌕ Search

SEARCH · Search NASA

Results for “optimized energy density functionals”

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 667 records · Page 37
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

A self-consistent first-principle based approach to model carrier mobility in organic materials

Transport through thin organic amorphous films, utilized in OLEDs and OPVs, has been a challenge to model by using ab-initio methods. Charge carrier mobility depends strongly on the disorder strength and reorganization energy, both of which are significantly affected by the details in environment of each molecule. Here we present a multi-scale approach to describe carrier mobility in which the materials morphology is generated using DEPOSIT, a Monte Carlo based atomistic simulation approach, or, alternatively by molecular dynamics calculations performed with GROMACS. From this morphology we extract the material specific hopping rates, as well as the on-site energies using a fully self-consistent embedding approach to compute the electronic structure parameters, which are then used in an analytic expression for the carrier mobility. We apply this strategy to compute the carrier mobility for a set of widely studied molecules and obtain good agreement between experiment and theory varying over several orders of magnitude in the mobility without any freely adjustable parameters. The work focuses on the quantum mechanical step of the multi-scale workflow, explains the concept along with the recently published workflow optimization, which combines density functional with semi-empirical tight binding approaches. This is followed by discussion on the analytic formula and its agreement with established percolation fits as well as kinetic Monte Carlo numerical approaches. Finally, we skatch an unified multi-disciplinary approach that integrates materials science simulation and high performance computing, developed within EU project MMM@HPC.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Temperature-dependent dielectric and energy-storage properties of Pb(Zr,Sn,Ti)O{sub 3} antiferroelectric bulk ceramics

The dielectric and energy-storage properties of Pb{sub 0.99}Nb{sub 0.02}[(Zr{sub 0.60}Sn{sub 0.40}){sub 0.95}Ti{sub 0.05}]{sub 0.98}O{sub 3} (PNZST) bulk ceramics near the antiferroelectric (AFE)-ferroelectric (FE) phase boundary are investigated as a function of temperature. Three characteristic temperatures T{sub 0}, T{sub C}, T{sub 2} are obtained from the dielectric temperature spectrum. At different temperature regions (below T{sub 0}, between T{sub 0} and T{sub C}, and above T{sub C}), three types of hysteresis loops are observed as square double loop, slim loop and linear loop, respectively. The switching fields and recoverable energy density all first increase and then decrease with increasing temperature, and reach their peak values at ∼T{sub 0}. These results provide a convenient method to optimize the working temperature of antiferroelectric electronic devices through testing the temperature dependent dielectric properties of antiferroelectric ceramics.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2016

Charge transfer mechanism in titanium-doped microporous silica for photocatalytic water-splitting applications

Solar energy conversion into chemical form is possible using artificial means. One example of a highly-efficient fuel is solar energy used to split water into oxygen and hydrogen. Efficient photocatalytic water-splitting remains an open challenge for researchers across the globe. Despite significant progress, several aspects of the reaction, including the charge transfer mechanism, are not fully clear. Density functional theory combined with density matrix equations of motion were used to identify and characterize the charge transfer mechanism involved in the dissociation of water. A simulated porous silica substrate, using periodic boundary conditions, with Ti 4+ ions embedded on the inner pore wall was found to contain electron and hole trap states that could facilitate a chemical reaction. A trap state was located within the silica substrate that lengthened relaxation time, which may favor a chemical reaction. A chemical reaction would have to occur within the window of photoexcitation; therefore, the existence of a trapping state may encourage a chemical reaction. Furthermore, this provides evidence that the silica substrate plays an integral part in the electron/hole dynamics of the system, leading to the conclusion that both components (photoactive materials and support) of heterogeneous catalytic systems are important in optimization of catalytic efficiency.

14 SOLAR ENERGY↗
DOE OSTI2025

Lewis Acid Site Engineering in Chromite Spinels Orchestrated Surface Reconstruction and Surpasses RuO 2 in Oxygen Evolution

Atomic-scale engineering of chromite spinels featuring redox-active tetrahedral A-sites and strong Cr–O covalency offers a promising route to superior platinum-group-metal-free oxygen evolution reaction (OER) catalysts. However, comprehensive studies addressing how cation substitution influences surface chemistry and governs OER activity and durability in chromite spinels remain limited. Here, in this work, a systematic investigation of the multicationic chromite series Ni x Fe y Cr 3−x−y O 4 is presented, identifying composition-dependent Lewis acidity as a descriptor of superior OER performance. It is further demonstrated that tuning surface acidity directly controls dynamic reconstruction processes and lattice-oxygen participation during spinel-based electrocatalysis. Following activation, the optimized Ni 0.8 Fe 0.3 Cr 1.9 O 4 catalyst delivers a current density of 10 mA cm −2 at an overpotential of 235 mV, surpassing RuO 2 , with excellent long-term stability. Integrating microscopic and spectroscopic analysis with operando impedance spectroscopy, it shows that activation generates an oxyhydroxide overlayer and reveals a previously unrecognized link between surface Lewis acidity and the growth kinetics and activity of these shells. Density functional theory calculations indicate that Fe incorporation at octahedral sites raises the O 2p-band center and lowers oxygen-vacancy formation energy, promoting lattice-oxygen activation and triggering reconstruction, yielding enhanced OER. This work integrates cation-driven surface-acidity modulation, acidity-governed reconstruction, and OER activity enhancement into a unified predictive framework for designing earth-abundant spinel-based catalysts.

operando impedance spectroscopy↗
DOE OSTI2014

Predicted boron-carbide compounds: A first-principles study

By using developed particle swarm optimization algorithm on crystal structural prediction, we have explored the possible crystal structures of B-C system. Their structures, stability, elastic properties, electronic structure, and chemical bonding have been investigated by first-principles calculations with density functional theory. The results show that all the predicted structures are mechanically and dynamically stable. An analysis of calculated enthalpy with pressure indicates that increasing of boron content will increase the stability of boron carbides under low pressure. Moreover, the boron carbides with rich carbon content become more stable under high pressure. The negative formation energy of predicted B{sub 5}C indicates its high stability. The density of states of B{sub 5}C show that it is p-type semiconducting. The calculated theoretical Vickers hardnesses of B-C exceed 40 GPa except B{sub 4}C, BC, and BC{sub 4}, indicating they are potential superhard materials. An analysis of Debye temperature and electronic localization function provides further understanding chemical and physical properties of boron carbide.

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

Application of Autonomous Smart Inverter Volt-VAR Function for Voltage Reduction Energy Savings and Power Quality in Electric Distribution Systems: Preprint

This paper evaluated the impact of smart inverter Volt-VAR function on voltage reduction energy saving and power quality in electric power distribution systems. A methodology to implement the voltage reduction optimization was developed by controlling the substation LTC and capacitor banks, and having smart inverters participate through their autonomous Volt-VAR control. In addition, a power quality scoring methodology was proposed and utilized to quantify the effect on power distribution system power quality. All of these methodologies were applied to a utility distribution system model to evaluate the voltage reduction energy saving and power quality under various PV penetrations and smart inverter densities.

14 SOLAR ENERGY↗
DOE OSTI2016

SU-G-BRC-15: The Potential Clinical Significance of Dose Mapping Error for Intra- Fraction Dose Mapping for Lung Cancer Patients

Purpose: To evaluate the dose-mapping error (DME) inherent to conventional dose-mapping algorithms as a function of dose-matrix resolution. Methods: As DME has been reported to be greatest where dose-gradients overlap tissue-density gradients, non-clinical 66 Gy IMRT plans were generated for 11 lung patients with the target edge defined as the maximum 3D density gradient on the 0% (end of inhale) breathing phase. Post-optimization, Beams were copied to 9 breathing phases. Monte Carlo dose computed (with 2*2*2 mm{sup 3} resolution) on all 10 breathing phases was deformably mapped to phase 0% using the Monte Carlo energy-transfer method with congruent mass-mapping (EMCM); an externally implemented tri-linear interpolation method with voxel sub-division; Pinnacle’s internal (tri-linear) method; and a post-processing energy-mass voxel-warping method (dTransform). All methods used the same base displacement-vector-field (or it’s pseudo-inverse as appropriate) for the dose mapping. Mapping was also performed at 4*4*4 mm{sup 3} by merging adjacent dose voxels. Results: Using EMCM as the reference standard, no clinically significant (>1 Gy) DMEs were found for the mean lung dose (MLD), lung V20Gy, or esophagus dose-volume indices, although MLD and V20Gy were statistically different (2*2*2 mm{sup 3}). Pinnacle-to-EMCM target D98% DMEs of 4.4 and 1.2 Gy were observed ( 2*2*2 mm{sup 3}). However dTransform, which like EMCM conserves integral dose, had DME >1 Gy for one case. The root mean square RMS of the DME for the tri-linear-to- EMCM methods was lower for the smaller voxel volume for the tumor 4D-D98%, lung V20Gy, and cord D1%. Conclusion: When tissue gradients overlap with dose gradients, organs-at-risk DME was statistically significant but not clinically significant. Target-D98%-DME was deemed clinically significant for 2/11 patients (2*2*2 mm{sup 3}). Since tri-linear RMS-DME between EMCM and tri-linear was reduced at 2*2*2 mm{sup 3}, use of this resolution is recommended for dose mapping. Interpolative dose methods are sufficiently accurate for the majority of cases. J.V. Siebers receives funding support from Varian Medical Systems.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2018

Mass, momentum, and energy transfer in supersonic aerosol deposition processes

Aerosol deposition (AD) can lead to the formation of dense coatings via deposition of particles from a gas flow; in AD the aerosol is passed through a converging-diverging nozzle, facilitating inertial particle impaction on a desired substrate at supersonic particle velocities. Unique from thermal spray methods, AD can be applied near room temperature and unique from cold spray, in AD the aerosol is typically at atmospheric pressure upstream of the nozzle. Though AD has been successfully demonstrated previously, a number of aspects related to particle motion in AD systems remain poorly understood. In this work, we simulated compressible flow field profiles and particle trajectories for typical AD working conditions for a slit type converging-diverging nozzle with a planar substrate. In examining the fluid flow profile, we show that the velocity and pressure profiles, as well as the shock structure are sensitive to the upstream and downstream operating pressures of the nozzle. These ultimately affect particle impaction speed. Importantly, in AD, the particle drag regime is dynamic; both particle Knudsen numbers and Mach numbers can vary by orders of magnitude. To aid particle trajectory simulations, we trained a neural network to predict the drag force on the particles based on existing experimental data, theoretical limits, and new direct simulation Monte Carlo (DMSC) results. The neural network based drag law, which depends upon both Mach and Knudsen numbers, shows better agreement with the DSMC simulation data than pre-existing correlations. With it, particle trajectory simulation results reveal that for a given particle density, there exists an optimal particle diameter to maximize particle impaction speed. We also find that in AD particles undergo size dependent inertial focusing, i.e. there is a particular particle diameter where the particle deposition linewidth is minimized. Particles smaller than this diameter are underfocused, and particles larger than this are overfocused, and hence have larger deposition linewidths in both cases. Using trajectory simulations, we additionally develop a framework that can be used to evaluate the position-dependent mass, momentum and kinetic energy fluxes to the deposition substrate for any aerosol size distribution function upstream the nozzle. It is shown for typical aerosol concentrations achievable in the laboratory, the kinetic energy flux can approach a magnitude normally observed in convective heat transfer with phase change, hence translational kinetic energy to thermal energy transfer in AD is likely the key contributor to the formation of dense coatings.

42 ENGINEERING↗
DOE OSTI2016

Tailoring Kinetics on a Topological Insulator Surface by Defect-Induced Strain: Pb Mobility on Bi 2 Te 3

Heteroepitaxial structures based on Bi2Te3-type topological insulators (TIs) exhibit exotic quantum phenomena. For optimal characterization of these phenomena, it is desirable to control the interface structure during film growth on such TIs. In this process, adatom mobility is a key factor. We demonstrate that Pb mobility on the Bi2Te3(111) surface can be modified by the engineering local strain, ε, which is induced around the point-like defects intrinsically forming in the Bi2Te3(111) thin film grown on a Si(111)-7 × 7 substrate. Scanning tunneling microscopy observations of Pb adatom and cluster distributions and first-principles density functional theory (DFT) analyses of the adsorption energy and diffusion barrier Ed of Pb adatom on Bi2Te3(111) surface show a significant influence of ε. Surprisingly, Ed reveals a cusp-like dependence on ε due to a bifurcation in the position of the stable adsorption site at the critical tensile strain εc ≈ 0.8%. This constitutes a very different strain-dependence of diffusivity from all previous studies focusing on conventional metal or semiconductor surfaces. Kinetic Monte Carlo simulations of Pb deposition, diffusion, and irreversible aggregation incorporating the DFT results reveal adatom and cluster distributions compatible with our experimental observations.

Chemistry↗
DOE OSTI2025

On-lattice kinetic Monte Carlo approaches for modeling molecular anisotropy in resveratrol crystallization

Stilbenes are a class of organic compounds with broad-ranging pharmaceutical and agricultural applications, which are typically isolated and purified through recrystallization. We are motivated by reducing experimental waste and optimizing yield via developing predictive simulations for processing-dependent crystal morphologies. Using resveratrol as a model stilbene system, we have developed an approach for simulating crystallization with molecular resolution using on-lattice kinetic Monte Carlo. In this work, we highlight modifications to the Stochastic Parallel PARticle Kinetic Simulator (SPPARKS) software package, which were essential to this application. Key enhancements include the incorporation of non-orthogonal cell shapes and monomer anisotropy approximations using bound hard spheres. This new SPPARKS application has been applied to resveratrol with attachment energy libraries obtained from density functional theory, resulting in excellent agreement with experimental morphology prediction.

crystallization↗
DOE OSTI2015

Crystal orbital study on the double walls made of nanotubes encapsulated inside zigzag carbon nanotubes

The structure stabilities and electronic properties are investigated by using ab initio self-consistent-field crystal orbital method based on density functional theory for the one-dimensional (1D) double-wall nanotubes made of n-gon SiO{sub 2} nanotubes encapsulated inside zigzag carbon nanotubes. It is found that formation of the combined systems is energetically favorable when the distance between the two constituents is around the Van der Waals scope. The obtained band structures show that all the combined systems are semiconductors with nonzero energy gaps. The frontier energy bands (the highest occupied band and the lowest unoccupied band) of double-wall nanotubes are mainly derived from the corresponding carbon nanotubes. The mobilities of charge carriers are calculated to be within the range of 10{sup 2}–10{sup 4} cm{sup 2} V{sup −1} s{sup −1} for the hybrid double-wall nanotubes. Young’s moduli are also calculated for the combined systems. For the comparison, geometrical and electronic properties of n-gon SiO{sub 2} nanotubes are also calculated and discussed. - Graphical abstract: Structures and band structures of the optimum 1D Double walls nanotubes. The optimized structures are 3-gon SiO2@(15,0), 5-gon SiO2@(17,0), 6-gon SiO2@(18,0) and 7-gon SiO2@(19,0). - Highlights: • The structure and electronic properties of the 1D n-gon SiO{sub 2}@(m,0)s are studied using SCF-CO method. • The encapsulation of 1D n-gon SiO{sub 2} tubes inside zigzag carbon nanotubes can be energetically favorable. • The 1D n-gon SiO{sub 2}@(m,0)s are all semiconductors. • The mobility of charge carriers and Young’s moduli are calculated.

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

NbV_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 niobium-vanadium (Nb-V). 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 Nb and V. 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 NERSC-Perlmutter and OLCF-Summit using the VASP 6.3.2. The VASP calculations for every atomic structure have been performed in 2 main steps: 1. Starting from an ideal body-centered-cubic (BCC) structure, geometry optimization with low precision has been executed to perform a preliminary optimization of the atomic structure. The output for this calculations is available in the files 0.CONTCAR, 0.OUTCAR, rlx1.out. 2. Using the atomic structure resulting from the preliminary geometry optimization, a second geometry optimization has been performed using normal precision. The output for this calculations is available in the files CONTCAR, OUTCAR, rlx2.out, vaspout.h5, and vasprun.xml. Cases 1-10 have been run without generating the file 'vaspout.h5'. Every chemical composition sampled across the composition range in the dataset has its own directory. The convention used to name the directories for binary alloys is AXBY, where A and B refer to the constituents, whereas X and Y are positive integers that represent the number of atoms for each constituent and their values still sum up to 128. Each atomic structure associated with a specific chemical composition has its own sub-directory within the directory of the corresponding chemical composition. The sub-directories for each atomic structure for each chemical composition are named 'case-*', where * is a positive integer that spans all the values from 1 through 100, extremes included. The files contained in each sub-directory 'case-*' for each atomic structure are as follows: FILES contained in each subdirectory with name case-N where N ranges between 11 and 100, extremes included: 1. INCAR: input file that contains various parameters and settings for controlling the behavior of the electronic structure calculations 2. KPOINTS: input file that specifies the Bloch vectors (k points) used to sample the Brillouin zone 3. 0.POSCAR: input file that defines the atomic structure of a system 4. 0.CONTCAR: output file that provides the atomic positions and cell parameters after the first geometry optimization has been run with the precision variable set to PREC=Low in the INCAR file 5. 0.OUTCAR: output file that contains detailed information about the progress of a calculation after the first geometry optimization has been run with the precision variable set to PREC=Low in the INCAR file 6. rlx1.out: file with diagnostic information about the execution of the first geometry optimization with precision variable set to PREC=Low in the INCAR file 7. POSCAR: input file that defines the atomic structure of a system after the first geometry optimization has been run at low precision. This represents the input for the second geometry optimization run with the precision variable set to PREC=Normal in the INCAR file 8. CONTCAR: output file that provides the atomic positions and cell parameters after the second geometry optimization has been run with the precision variable set to PREC=Normal in the INCAR file 9. OUTCAR: output file that contains detailed information about the progress of a calculation after the second geometry optimization has been run with the precision variable set to PREC=Normal in the INCAR file 10. rlx2.out: file with diagnostic information about the execution of the second geometry optimization with precision variable set to PREC=Normal in the INCAR file 11. vaspout.h5: hierarchical HDF5 file containing the inputs and outputs of a VASP calculation. To analyze the data in this file we recommend using py4vasp. This file is only produced if the VASP version used is compiled with HDF5 support 12. vasprun.xml: contains similar information to OUTCAR, but in an xml format. Subdirectories with name case-N, where N ranges between 1 and 10 (extremes included) contain all the files listed above except 'vaspout.h5'. Subdirectories with name case-N, where N ranges between 41 and 60 (extremes included), contain a duplicate copy of the files listed above except for KPOINTS. The names of the duplicate files end with -bis, and correspond to a second VASP calculation that has converged to a different optimized geometry. This research is sponsored by the Artificial Intelligence Initiative as part of the Laboratory Directed Research and Development (LDRD) Program of Oak Ridge National Laboratory, managed by UT-Battelle, LLC, for the US Department of Energy under contract DE-AC05-00OR22725. This work used resources of the Oak Ridge Leadership Computing Facility, which is supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC05-00OR22725, under Directorate Discretionary awards MAT025 (Materials Science) and LRN026 (Machine Learning), and INCITE award MAT201. This work also used resources of the National Energy Research Scientific Computing Center, which is supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC02-05CH11231, under award ERCAP0025216. REFERENCES (1) Kresse, G. and Hafner, J. Ab initio molecular dynamics for liquid metals. Phys. review B 47, 558 (1993). (2) Kresse, G. and Hafner, J. Ab initio molecular-dynamics simulation of the liquid-metal-amorphous-semiconductor transition in germanium. Phys. Rev. B 49, 14251 (1994) (3) Kresse, G. and Furthmüller, J. Efficiency of ab-initio total energy calculations for metals and semiconductors using a plane-wave basis set. Comput. materials science 6, 15-50 (1996) (4) Kresse, G. and Furthmüller, J. Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. Phys. review B 54, 11169 (1996) (5) Kresse, G. and Joubert, D. From ultrasoft pseudopotentials to the projector augmented-wave method. Phys. review b 59, 1758 (1999)

36 MATERIALS SCIENCE↗
DOE OSTI2026

Physics-informed machine learning exploration of Na storage mechanisms in disordered carbon

Sodium-ion batteries are a cost-effective, sustainable alternative to lithium-ion systems for large-scale energy storage. However, optimizing sodium storage in carbon-based anodes with microstructural complexity and atomic disorder remains a major challenge. The intrinsic inhomogeneity of these materials produces diverse local environments, making it difficult for conventional methods to predict and control ion dynamics. Hard carbon (HC) anodes, composed of ranges of ordered-to-disordered graphitic and amorphous nanodomains, offer tunable ion storage and rate capacity, yet rationale design remains a challenge due to poorly understood correlation between local atomic feature and ion transport mechanism. Here, to address this challenge, we introduce a data-driven framework that integrates validated machine-learned interatomic potentials, large-scale molecular dynamics simulations, and machine learning to elucidate sodium transport mechanisms as a function of carbon and sodium loading densities. By computing per-ion structural descriptors and applying unsupervised learning, we identify distinct diffusion modes governed by microscopic features. Supervised analysis and correlation mapping then establish quantitative links between these transport regimes and processing variables such as bulk carbon density and sodium content. This physics-informed approach establishes quantitative structure–transport relationships and offers actionable design principles for engineering high-performance HC anodes.

Data-driven framework↗
DOE OSTI2020

Size effect on the growth and pulverization behavior of Si nanodomains in SiO anode

Silicon monoxide (SiO) as a promising silicon-based anode electrode of lithium-ion batteries (LIBs) has been the subject of many recent investigations. However, the complex microstructural features of SiO hinder attempts to link the microstructure changes in SiO with its electrochemical performances upon electrochemical cycling. This study employs X-ray pair distribution function (PDF), with a high-resolution transmission electron microscope (HRTEM), to investigate the microstructure of Si nanodomains in SiO and its structural evolution over the electrochemical cycling process. The experiment results reveal the size evolution of Si nanodomains in SiO during the electrochemical cycling and find it highly dependent on the initial Si domain size in SiO. If the initial Si domain size is too large, the average size of Si nanodomains would increase at first and then decrease after a certain number of cycles, which indicates the pulverization of Si domains after the electrochemical-driven growth of Si nanodomains reaches a critical point. Furthermore, these results suggest that an optimal initial Si nanodomain size of 4-6 nm for SiO anode materials is essential to retard the growth and subsequent pulverization process of Si nanodomains and thus to improve the cyclability of SiO during long-term cycling.

25 ENERGY STORAGE↗
DOE OSTI2020

Towards defect engineering in hexagonal MoS 2 nanosheets for tuning hydrogen evolution and nitrogen reduction reactions

Combined computational and experimental approaches were used to evaluate defective 2H-MoS 2 nanosheets for their activity and selectivity for hydrogen evolution reaction (HER) and nitrogen reduction reaction (NRR). Density functional theory calculations were used to understand the relationship between HER and NRR activity on the ideal basal MoS 2 plane, seven grain-boundaries, ten single-/few-atom vacancies and anti-sites, and zigzag and armchair edge sites. The results confirm that 2H-MoS 2 should contain several defects with high activity for HER: armchair and zigzag edges, VS vacancy, MoS 2 anti-site, and S-S and Mo-Mo grain boundaries. Considering Gibbs free energy change for all the steps in the NRR mechanism and kinetic barriers for a key NRR step, we have found that activation of perspective NRR selective sites in 2H-MoS 2 , namely V MoS6 and clusters of S-vacancies, would require large overpotential, conditions at which HER dominates. The DFT conclusions are supported by the electrochemical studies of NRR activity and selectivity under aqueous conditions, which show an increase in NRR activity but a decrease in Faradaic efficiency as applied cell potential becomes more negative. The results of this work therefore highlight the challenges in activating natural 2H-MoS 2 for NRR, which would require additional material engineering or reaction condition optimization as a way to suppress HER, decrease the NRR overpotential or preferentially both.

2H-MoS2↗
DOE OSTI2021

Actinide Molten Salts: A Machine-Learning Potential Molecular Dynamics Study

We know that actinide molten salts represent a class of important materials in nuclear energy. Understanding them at a molecular level is critical to proper and optimal design of relevant technological applications. Yet, owing to the complexity of electronic structure due to the 5f orbitals, computational studies of heavy elements in condensed phases using ab initio potentials to study the structure and dynamics of these elements embedded in molten salts are difficult. This lack of efficient computational protocols makes it difficult to obtain information on properties that require extensive statistical sampling like transport. To tackle this problem, we adopted a machine-learning approach to study ThCl 4 -NaCl and UCl 3 -NaCl binary systems. The machine-learning potential, with the density functional theory accuracy, allows us to obtain long molecular dynamics trajectories (ns) for large systems (10 3 atoms) at a considerably low computing cost, thereby efficiently gaining information about their bonding structures, thermodynamics, and dynamics at a range of temperature. We observed a considerable change in the coordination environments of actinide elements and their characteristic coordination-sphere lifetime. Our study also suggests that actinides in molten salts may not follow well known entropy-scaling laws.

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

Defect Chemistry and Hydrogen Transport in La/Sr-Based Oxyhydrides

Oxyhydrides in the series La 2 – y Sr y LiH 1 + y O 3 – y are attractive materials for solid-state hydride electrolytes. However, their chemical stability is poor, and thus far only La 2 LiHO 3 (y = 0) has been demonstrated in a working fuel cell. Using a first-principles approach to study defect chemistry and stability, we find that La 2 LiHO 3 can be stabilized over a range of H-rich, O-poor chemical conditions, while Sr 2 LiH 3 O has low stability. Here, defects are critical for ionic transport and for stability, particularly in Sr 2 LiH 3 O, in which anion antisites and H vacancies readily form. In fact, we show that the ground-state structure of Sr 2 LiH 3 O is inherently disordered. La 2 LiHO 3 has low conductivity, which we connect to the high formation energy of defects such as VO and OH that are necessary for hydride conduction. When Sr 2 LiH 3 O is grown under carefully selected conditions, V H + defects can be made prevalent; these defects give rise to much higher ionic conductivity than can be obtained in La 2 LiHO 3 . In general, we show that the choice of synthesis conditions is vital when seeking to optimize the stability and ionic conductivity of these oxyhydrides, thereby tailoring them for use in solid-state hydrogen fuel cells.

08 HYDROGEN↗
DOE OSTI2014

Modeling Excited States in TiO2 Nanoparticles: On the Accuracy of a TD-DFT Based Description

We have investigated the suitability of Time-Dependent Density Functional Theory (TD-DFT) to describe vertical low-energy excitations in naked and hydrated titanium dioxide nanoparticles through a comparison with results from Equation-of-Motion Coupled Cluster (EOM-CC) quantum chemistry methods. We demonstrate that for most TiO2 nanoparticles TD-DFT calculations with commonly used exchange-correlation (XC-)potentials (e.g. B3LYP) and EOM-CC methods give qualitatively similar results. Importantly, however, we also show that for an important subset of structures, TD-DFT gives qualitatively different results depending upon the XC-potential used and that in this case only TD-CAM-B3LYP and TD-BHLYP calculations yield results that are consistent with those obtained using EOM-CC theory. Moreover, we demonstrate that the discrepancies for such structures arise from a particular combination of defects, excitations involving which are charge-transfer excitations and hence are poorly described by XC-potentials that contain no or low fractions of Hartree-Fock like exchange. Finally, we discuss that such defects are readily healed in the presence of ubiquitously present water and that as a result the description of vertical low-energy excitations for hydrated TiO2 nanoparticles is hence non-problematic.

calculations↗