Generalizable density functional theory based photoemission model for the accelerated development of photocathodes and other photoemissive devices
Not Available
SEARCH · Search NASA
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.
Not Available
Not Available
From the lightest Hydrogen isotopes up to the recently synthesized Oganesson (Z = 118), it is estimated that as many as about 8,000 atomic nuclei could exist in nature. Most of these nuclei are too short-lived to be occurring on Earth, but they play an essential role in astrophysical events such as supernova explosions or neutron star mergers that are presumed to be at the origin of most heavy elements in the Universe. Understanding the structure, reactions, and decays of nuclei across the entire chart of nuclides is an enormous challenge because of the experimental difficulties in measuring properties of interest in such fleeting objects and the theoretical and computational issues of simulating strongly-interacting quantum many-body systems. Nuclear density functional theory (DFT) is a fully microscopic theoretical framework which has the potential of providing such a quantitatively accurate description of nuclear properties for every nucleus in the chart of nuclides. Thanks to high-performance computing facilities, it has already been successfully applied to predict nuclear masses, global patterns of radioactive decay like β or γ decay, and several aspects of the nuclear fission process such as, e.g., spontaneous fission half-lives. Yet, predictive simulations of nuclear spectroscopy—the low-lying excited states and transitions between them—or of nuclear fission, or the quantification of theoretical uncertainties and their propagation to basic or applied nuclear science applications, would require several orders of magnitude more calculations than currently possible. However, most of this computational effort would be spent into generating a suitable basis of DFT wavefunctions. Such a task could potentially be considerably accelerated by borrowing tools from the field of machine learning and artificial intelligence. In this paper, we review different approaches to applying supervised and unsupervised learning techniques to nuclear DFT.
Ni coarsening in solid oxide cell (SOC) electrodes is known to be significantly faster under a humid atmosphere. The underlying mechanisms, though, have not been fully understood. In this work, we examine the surface diffusion of Ni(OH) x by a combination of density-functional theory and phase-field modeling. Here, density-functional theory is used to evaluate the adsorption and diffusion of the Ni(OH) x species on Ni (111) surface, and the results are used to obtain the effective surface diffusivity of Ni versus steam and hydrogen gas partial pressures. This diffusivity is used as an input for a phase-field model to investigate the Ni coarsening in SOC electrodes. It is found that Ni(OH) x formed through chemical reaction cannot accelerate coarsening unless an extremely high steam to hydrogen ratio is reached. However, if Ni(OH) x is formed through electrochemical reactions near triple phase boundaries (TPBs), surface diffusion of Ni(OH) x may cause faster Ni coarsening in fuel cell mode under a large overpotential. Specifically, surface diffusion of Ni(OH) may be comparable to or faster than that of Ni under an overpotential that is large but still possible under fuel cell operating conditions.
The successful commercialization of electrochemical energy-conversion systems hinges on a deeper understanding of electrocatalytic reaction kinetics. Despite extensive research, a key descriptor that characterizes electrolyte effects on reaction kinetics remains elusive. Here, surface charge in electrical double layers (EDLs) is introduced as a descriptor for electrolyte-dependent kinetics. The surface charge is calculated with a continuum EDL model parameterized by density-functional theory. The model is validated by reproducing the anomalously low slope of Pt(111) in Parsons-Zobel plots. Strong correlations are observed between calculated surface charge and experimental kinetic currents for hydrogen evolution, oxygen reduction, and CO 2 -reduction reactions across various pH levels and cationic species. These correlations can be either promotional or inhibitory, depending on solute-intermediate interactions. In acidic media, incorporating adsorbate charge captures specific adsorption effects in oxygen reduction reaction. In conclusion, these findings establish surface charge density as a key descriptor for electrolyte-dependent kinetics, which will guide the design of the electrode/electrolyte interface.
Ordering of different chalcogens, S, Se, and Te, on Au(111) exhibit broad similarities but also some distinct features, which must reflect subtle differences in relative values of the long-range pair and many-body lateral interactions between adatoms. We develop lattice-gas (LG) models within a cluster expansion framework, which includes about 50 interaction parameters. These LG models are developed based on density functional theory (DFT) analysis of the energetics of key adlayer configurations in combination with the Monte Carlo (MC) simulation of the LG models to identify statistically relevant adlayer motifs, i.e., model development is based entirely on theoretical considerations. The MC simulation guides additional DFT analysis and iterative model refinement. Given their complexity, development of optimal models is also aided by strategies from supervised machine learning. The model for S successfully captures ordering motifs over a broader range of coverage than achieved by previous models, and models for Se and Te capture the features of ordering, which are distinct from those for S. More specifically, the modeling for all three chalcogens successfully explains the linear adatom rows (also subtle differences between them) observed at low coverages of ~0.1 monolayer. The model for S also leads to a new possible explanation for the experimentally observed phase with a (5 × 5)-type low energy electron diffraction (LEED) pattern at 0.28 ML and to predictions for LEED patterns that would be observed with Se and Te at this coverage.
We report extended Lagrangian Born–Oppenheimer molecular dynamics (XL-BOMD) is formulated for orbital-free Hohenberg–Kohn density-functional theory and for charge equilibration and polarizable force-field models that can be derived from the same orbital-free framework. The purpose is to introduce the most recent features of orbital-based XL-BOMD to molecular dynamics simulations based on charge equilibration and polarizable force-field models. These features include a metric tensor generalization of the extended harmonic potential, preconditioners, and the ability to use only a single Coulomb summation to determine the fully equilibrated charges and the interatomic forces in each time step for the shadow Born–Oppenheimer potential energy surface. The orbital-free formulation has a charge-dependent, short-range energy term that is separate from long-range Coulomb interactions. This enables local parameterizations of the short-range energy term, while the long-range electrostatic interactions can be treated separately. The theory is illustrated for molecular dynamics simulations of an atomistic system described by a charge equilibration model with periodic boundary conditions. The system of linear equations that determines the equilibrated charges and the forces is diagonal, and only a single Ewald summation is needed in each time step. The simulations exhibit the same features in accuracy, convergence, and stability as are expected from orbital-based XL-BOMD.
Understanding ion adsorption in porous carbons is crucial for a range of technologies, including water desalination and energy storage. In this work, we combined density functional theory with a continuum solvation model to investigate thermodynamics and kinetics of the adsorption process of alkali metal ions from aqueous solutions into carbon nanopores with different sizes and geometries. We found that cations with a larger ionic radius are more favorable to enter the nanopores because of a lower energy penalty of dehydration. In addition, the pore size and geometry were found to have a significant impact on the ion–pore interaction under confinement and cation selectivity. Our study highlights a complex interplay among nanopore geometry, ion size, and hydration on the cation adsorption selectivity, suggesting that tuning the porosity could represent a general strategy for improving ion separations.
Strain-engineering of bimetallic nanomaterials is an important design strategy for developing new catalysts. In this paper, we introduce an approach for including strain effects into a recently introduced, density functional theory (DFT)-based alloy stability model. The model predicts adsorption site stabilities in nanoparticles and connects these site stabilities with catalytic reactivity and selectivity. Strain-based dependencies will increase the model’s accuracy for nanoparticles affected by finite-size effects. In addition to the stability of small nanoparticles, strain also influences the heat of adsorption of epitaxially grown metal-on-metal adlayers. In this respect, we successfully benchmark the strain-including alloy stability model with previous experimentally determined trends in the heats of adsorption of Au and Cu adlayers on Pt (111). For these systems, our model predicts stronger bimetallic interactions in the first monolayer than monometallic interactions in the second monolayer. We explicitly quantify the interplay between destabilizing strain effects and the energy gained by forming new metal–metal bonds. While tensile strain in the first Cu monolayer significantly destabilizes the adsorption strength, compressive strain in the first Au monolayer has a minimal impact on the heat of adsorption. Hence, this study introduces and, by comparison with previous experiments, validates an efficient DFT-based approach for strain-engineering the stability, and, in turn, the catalytic performance, of active sites in bimetallic alloys with atomic level resolution.
Here we present a dynamic density functional theory for modeling the effects of applied electric fields on the local structure of polymers with added salt (polymer electrolytes). Time-dependent equations for the local electrostatic potential and volume fractions of polymer, cation, and anion of added salt are developed using the principles of linear irreversible thermodynamics. For such a development, a field theoretic description of the free energy of polymer melts doped with salts is used, which captures the effects of local variations in the dielectric function. Connections of the dynamic density functional theory with experiments are established by relating the three phenomenological Onsager’s transport coefficients of the theory to the mutual diffusion of electrolyte, ionic conductivity, and transference number of one of the ions. The theory is connected with a statistical mechanical model developed by Bearman and Kirkwood [J. Chem. Phys. 28, 136 (1958)] after relating the three transport coefficients to friction coefficients. The steady-state limit of the dynamic density functional theory is used to understand the effects of dielectric inhomogeneity on the phase separation in polymer electrolytes. The theory developed here provides not only a way to connect with experiments but also to develop multi-scale models for studying connections between local structure and ion transport in polymer electrolytes.
Recent activity in the area of chiroptical phenomena has been focused on the connection between structural asymmetry, electron spin configuration and light/matter interactions in chiral semiconductors. In these systems, spin‐splitting phenomena emerge due to inversion symmetry breaking and the presence of extended electronic states, yet the connection to chiroptical phenomena is lacking. Here, in this study, we develop an analytical effective mass model of chiral excitons, parameterized by density functional theory. The model accounts for parity mixing of the band edge Bloch functions resulting from polar distortions, resulting in allowed magnetic dipole transitions. Through the study of a prototypical chiral 2D hybrid perovskite semiconductor, we show that circular dichroism of the chiral exciton and its interband continuum emerges from spin‐splitting via cross‐coupling of Rashba‐like and chiral/helical spin‐texture components. To demonstrate the generality of our approach, and as a counterpoint, we apply our model to describe chiroptical properties of three‐dimensionally confined excitons in perovskite nanocrystals that occur without chiral lattice distortions.
In 1954, Salpeter calculated nuclear reaction rates would increase via weak and strong electron screening in stars. Astrophysical observations and laboratory astrophysics experiments have confirmed that electron screening reduces the Coulomb Barrier between charged particles. Fermi degenerate strong screening occurs when the electron density exceeds 10e23 electrons/cm3 in the cores of gas giant planets, White Dwarf Stars, the conduction bands of metals, deuterated LENR materials and Inertial Confinement Fusion targets at maximum compression. We have modelled these quantum mechanical effects with Density Functional Theory codes.
As state-of-the-art (SOA) lithium-ion (Li-ion) batteries approach their specific energy limit (~250 Wh kg -1 ), layer-structured, nickel-rich (Ni-rich) lithium transition metal oxide-based cathode materials, e. g., LiNi 0.8 Mn 0.1 Co 0.1 O 2 (NMC811), have attracted great interest owing to their practical high specific capacities (>200 mAhg -1 ). Coupled with their high average discharge voltages (~4 V vs Li/Li + ), Ni-rich cathode-based lithium batteries possess a great potential to achieve much higher specific energies (>350 Wh kg -1 at the cell level) than the SOA Li-ion counterparts. In addition, Ni-rich oxides are low-cost battery cathode materials due to their low cobalt contents. However, Ni-rich cathode-based lithium batteries suffer quick capacity degradations upon cycling, particularly at high upper cutoff voltages (e.g., ≥= 4.5 V vs Li/Li + ), due to crystal structure changes of the active cathode materials and parasitic side reactions at the electrolyte/electrode interfaces. In this study, a fluorinated-solvent-based, high-voltage stable electrolyte (HVE), i.e., 1 M Li bis(trifluoromethanesulfonyl)imide (LiTFSI) in fluoroethylene carbonate (FEC), bis(2,2,2-trifluoroethyl) carbonate (FDEC), and methyl (2,2,2-trifluoroethyl) carbonate (FEMC) with Li difluoro(oxalate)borate (LiDFOB) additive, was formulated and evaluated in Li/NMC811 battery cells. Further, to the best of our knowledge, this class of electrolyte has not been investigated for Ni-rich cathode-based lithium batteries. Li/NMC811 cells with HVE exhibited a superior long-term cycle performance stability, maintaining ~80% capacity after ~500 cycles at a high cutoff voltage of 4.5 V (vs Li/Li + ) than a baseline carbonate-solvent-based electrolyte (BE). The superior cycle stability of the Li/NMC811 cells is attributed to the inherently high-voltage stability of HVE, supported by the physical and electrochemical analyses. This conclusion is supported by our density functional theory (DFT) modeling where HVE shows a less tendency of deprotonation/oxidation than BE, leading to the observed cycle stability. The findings in this study are important to help tackle the technical challenges facing Ni-rich cathode-based lithium batteries to realize their high energy density potentials with a long cycle life.
Chloride poisoning is a serious problem for the electrocatalytic reduction of aqueous nitrate (NO 3 – ) and improved electrocatalysts are needed. Here we study the electrocatalytic activity of rhodium sulfide supported on carbon (Rh x S y /C) for the reduction of nitrate and compare it against Pt/C and Rh/C in the presence of chloride. Between 0.05–0.15 V vs. RHE, Rh x S y /C has a steady-state nitrate reduction current density in 1 M H 2 SO 4 + 1 M NaNO 3 that is 1.6–5.6 times greater than Rh/C (the most active metal electrocatalyst) and 10–24 times greater than Pt/C. Current densities are decreased by 37% for Rh x S y /C, 62% for Rh/C, and 40% for Pt/C at 0.1 V vs. RHE in the presence of 1 mM chloride. The decrease in nitrate reduction activity for Pt, Rh, and Rh x S y is due to the competitive adsorption of chloride and nitrate on the surface. Density functional theory (DFT) modeling predicts that chloride poisoning will persistently inhibit nitrate reduction on metals due to linear adsorbate scaling relations between nitrate and chloride. DFT calculations and microkinetic modeling of our experimental measurements predict that nitrate converts to nitrite via an H-assisted dissociation mechanism on Pt and direct nitrate dissociation on Rh and Rh x S y . Pristine Rh x S y (i.e., Rh 3 S 4 , Rh 2 S 3 , and Rh 17 S 15 ) terraces are predicted to be inactive toward nitrate reduction. In contrast, sulfur vacancies in Rh 3 S 4 terraces are predicted to be active for nitrate reduction, but also bind chloride strongly. Furthermore, sulfur-defected Rh 3 S 4 rationalize the experimentally observed high activity but moderate chloride poison-resistance of Rh x S y /C for nitrate reduction.
Nonconventional NMR spin-coupling constants were investigated to determine their potential as conformational constraints in MA’AT modeling of the O-glycosidic linkages of oligosaccharides. Four ( 1 J C1',H1' , 1 J C1',C2' , 2 J C1',H2' , and 2 J C2',H1' ) and eight ( 1 J C4,H4 , 1 J C3,C4 , 1 J C4,C5 , 2 J C3,H4 , 2 J C4,H3 , 2 J C5,H4 , 2 J C4,H5 , and 2 J C3,C5 ) spin-couplings in methyl β-d-galactopyranosyl-(1→4)-β-d-glucopyranoside (methyl β-lactoside) were calculated using density functional theory (DFT) to determine their dependencies on O-glycosidic linkage C–O torsion angles, $\phi$ and ψ, respectively. Long-range 4 J H1',H4 was also examined as a potential conformational constraint of either $\phi$ or ψ. Secondary effects of exocyclic (hydroxyl) C–O bond rotation within or proximal to these coupling pathways were investigated. Based on the findings of methyl β-lactoside, analogous J-couplings were studied in five additional two-bond O-glycosidic linkages [βGlcNAc-(1→4)-βMan, 2-deoxy-βGlc-(1→4)-βGlc, αMan-(1→3)-βMan, αMan-(1→2)-αMan, and βGlcNAc(1→2)-αMan] to determine whether the coupling behaviors observed in methyl β-lactoside were more broadly observed. Of the 13 nonconventional J-couplings studied, 7 exhibit properties that may be useful in future MA’AT modeling of O-glycosidic linkages, none of which involve coupling pathways that include the linkage C–O bonds. The findings also provide new insights into the general effects of exocyclic C–O bond conformation on the magnitude of experimental spin-couplings in saccharides and other hydroxyl-containing molecules.
Modern density functional theory (DFT) is a powerful tool for accurately predicting self-consistent material properties such as equations of state, transport coefficients and opacities in high energy density plasmas, but it is generally restricted to conditions of local thermodynamic equilibrium (LTE) and produces only averaged electronic states instead of detailed configurations. We propose a simple modification to the bound-state occupation factor of a DFT-based average-atom model that captures essential non-LTE effects in plasmas—including autoionization and dielectronic recombination—thus extending DFT-based models to new regimes. Finally, we then expand the self-consistent electronic orbitals of the non-LTE DFT-AA model to generate multi-configuration electronic structure and detailed opacity spectra. This article is part of the theme issue ‘Dynamic and transient processes in warm dense matter’.