Adsorption Site Screening on a PGM-Free Electrocatalyst: Insights from Grand Canonical Density Functional Theory
Not Available
SEARCH · Engineering Papers
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
A chromium triiodide (CrI 3 ) monolayer has an interesting ferromagnetic ground-state. In this work, we calculate band structures and magnetic moments of tensile-strained and bent zigzag CrI 3 nanoribbons with density functional theory. The edge iodine atoms form flat low-lying conduction bands and couple with chromium atoms ferromagnetically, while the non-edge iodine atoms weakly couple antiferromagnetically. Narrow CrI 3 nanoribbons have two locally stable magnetic moment orientations, namely out-of-plane and in-plane (along the nanoribbon periodic direction) configurations. This enables four magnetization states in CrI 3 nanoribbons, including two out-of-plane ones (up and down) and two in-plane ones (forward and backward along the nanoribbon periodical direction), increasing the operating controllability. Based on the one-dimensional Ising spin chain model, the spin correlation length of the narrow CrI 3 nanoribbon is estimated as about 10 Å at its estimated Curie temperature of 27 Kelvin, which is lower than the measured 45 Kelvin of the monolayer CrI 3 . The optical absorption and magneto-optical properties of CrI 3 nanoribbons are investigated with many-body perturbation GW-BSE (Bethe-Salpeter equation), including magnetic dichroism, Faraday and magneto-optical Kerr effects. Here, the low-energy dark excitons are mainly from transitions between electrons and holes with unlike spins and are non-Frenkel-like, while the bright excitons have mixed spin configurations. The intrinsic lifetime of excitons can be over one nanosecond, suitable for quantum information processes. Tensile strains and bending manifestly modulate the absorption spectra and magneto-optical properties of CrI 3 nanoribbons within a technologically important photon energy-range of ~1.0-2.0 eV. The CrI 3 nanoribbons can be used in 1D or 2D magnetic storage nanodevices, tunable magnetic optoelectronics, and spin-based quantum information controls.
Charge generation and separation are two key processes for semiconductor photocatalysis. Here, we use Pt as a single atom catalyst to systematically examine the facet-dependent electronic band structure and light absorption of the layered perovskite-type wide-gap semiconductor, lanthanum titanium oxide (La 2 Ti 2 O 7 , LTO) by means of density functional theory simulations. In this study, it is found that single Pt atom doping of different LTO surfaces (here, (100), (101) and (001)) can not only create states in the bandgap that would promote the formation of recombination centers, but also shift the optical absorption edge to the visible region. Interestingly, the Pt doping forms a heterojunction, with the valence band maximum consisting of Pt states and the conduction band minimum consisting of LTO states, respectively. To study the facet dependent surface activity, acetic acid (CH 3 COOH) was used as a model molecule to investigate the adsorption and charge transfer on the (101), (100) and (001) Pt-LTO surface facets. The results show that the (101) facet could enable stronger adsorption of CH 3 COOH by promoting more electron transfer during the interfacial interaction. Our theoretical findings aim to promote the design and optimization of the single atom catalysts (SACs) for photocatalytic applications and other broad catalysis systems.
Accurately modeling photochemical reactions is difficult due to the presence of conical intersections and locally avoided crossings, as well as the inherently multiconfigurational character of excited states. As such, one needs a multistate method that incorporates state interaction in order to accurately model the potential energy surface at all nuclear coordinates. The recently developed linearized pair-density functional theory (L-PDFT) is a multistate extension of multiconfiguration PDFT, and it has been shown to be a cost-effective post-MCSCF method (as compared to more traditional and expensive multireference many-body perturbation methods or multireference configuration interaction methods) that can accurately model potential energy surfaces in regions of strong nuclear–electronic coupling in addition to accurately predicting Franck–Condon vertical excitations. Here, in this paper, we report the derivation of analytic gradients for L-PDFT and their implementation in the PySCF-forge software, and we illustrate the utility of these gradients for predicting ground- and excited-state equilibrium geometries and adiabatic excitation energies for formaldehyde, s-trans-butadiene, phenol, and cytosine.
The number and properties of the neutrons and photons emitted in nuclear fission are directly related to the excitation energy of the fission fragments when they are formed at scission. Though not observable experimentally because of the extremely short timescales, the excitation energy of fission fragments can be predicted by microscopic theory based on time-dependent density functional theory (TDDFT). Initial results on the value of the total kinetic energy of fission reactions were very promising, but could not probe all possible fragmentations. Here, in this work, we perform large-scale TDDFT calculations in 240 Pu enabled by the development of a new TDDFT solver. We obtain TDDFT trajectories covering nearly all possible fragmentations. We find that the total kinetic energy is close to experimental values only for the most likely fission while it is severely underestimated at both small and large asymmetries. This conclusion seems rather independent of the parametrization of the energy functional, both in its particle-hole and particle-particle channels.
Digital fringe projection is a surface-profiling technique used for highly accurate non-contact measurements. As with any measurement technique, a variety of sources degrade to the measurement accuracy of the method. Here, this paper presents an analytically-derived probability density function that explicitly models the surface height measurement error due to inevitable phase measurement error, and it includes the specific case of pixel noise inducing the phase measurement error that ultimately leads to the height estimation error. The accuracy of the model was validated through Monte-Carlo simulations of resultant height distributions subject to arbitrarily correlated pixel intensity noise and experimental digital fringe projection measurements where the pixel-by-pixel height uncertainty estimations were compared to the predictions of the derived model.
A computationally efficient and transferable approach for modeling reactions at metal/water interfaces could significantly accelerate our understanding and ultimately the development of new catalytic transformations, particularly in the context of the emerging field of biomass conversion. Here, we present a parametrization of Pt–X (X = H, O, C) density-functional tight-binding (DFTB) for addressing this need. We first constructed Pt–H, Pt–O, and Pt–C repulsive potential splines. These pairwise parameters were then augmented to include many-body interactions using the Chebyshev Interaction Model for Efficient Simulation (ChIMES). We compare the geometrical and energetic performances of both DFTB and DFTB/ChIMES methods with DFT reference data across a variety of organic molecules at the platinum surface from nanoparticles to single-crystal surfaces. DFTB shows limited transferability between extended crystal surfaces and small nanoparticles. This transferability is significantly improved through the introduction of three-body interactions with Pt in DFTB/ChIMES, which provides consistent results across various systems, with reductions in the RMSD from around 30 kcal/mol in DFTB to around 10 kcal/mol. We demonstrate the stability and reliability of the obtained parameters by performing metadynamic simulations for the adsorption of phenol on Pt(111). We observe that DFTB itself is undersolvating the surface, leading to only one or two chemisorbed water molecules in a c(4 × 6) unit cell. In contrast, DFTB/ChIMES leads to a coverage of about 0.5 ML and successfully captures the chemisorbed mode of phenol at both the solid/liquid and the solid/gas interfaces. Furthermore, in agreement with experimental measurements, the adsorption at the solid/liquid interface is significantly weaker than that at the solid/gas interface. As a result, we highlight that even with DFTB, where we can accumulate dynamics for more than 1 ns for a given system, the simulations are not fully converged.
In this work, we extend the damped Zaremba-Kohn model (dZK) for long-range dispersion interaction between a molecule and a planar surface to molecules adsorbed on a curved cylindrical surface, and employ this extended model as an additive correction to the semilocal density functionals PBE (Perdew-Burke-Ernzerhof) and SCAN (strongly constrained and appropriately normed). The resulting PBE+vdW (van der Waals)-dZK and SCAN+vdW-dZK are applied to two systems, NH3 and NO2 molecules adsorbed on a single-wall carbon nanotube (CNT), for calculations of binding energies and equilibrium distances. For comparison, the results from vdW nonlocal functionals, such as SCAN+rVV10 and PBE+rVV10, are also presented. The binding energies from PBE+rVV10 (Vydrov and Van Voorhis), SCAN+rVV10, PBE+vdW-dZK, and SCAN+vdW-dZK are about 70–115 meV for the system of CNT + NH 3 and 300–500 meV for the system of CNT + NO 2 . The results from PBE+vdW-dZK and SCAN+vdW-dZK are closer to each other than those from PBE+rVV10 and SCAN+rVV10 are. The relatively closer results from PBE+vdW-dZK and SCAN+vdW-dZK indicate the consistency of our developed vdW–dZK model for cylindrical surfaces. All methods, including PBE, SCAN, PBE+rVV10, SCAN+rVV10, PBE+vdW-dZK, and SCAN+vdW-dZK, give approximately the same binding energy differences between two adsorption configurations (types I and II) for the two systems. This implies that the two adsorption sites have approximately the same adsorption stability. The exponent of the vdW interaction power law from our vdW-dZK model for the two systems is about 0 at short distance, largely due to the damping factor, and tends slowly to –4 to –4.5 at distances D about 20–50 Å. At even larger distances, the vdW power-law exponent approaches –5. This feature is very similar to the one calculated with random-phase approximation and renormalization group approaches, supporting the applicability of our methods. Our developed vdW-dZK method provides a highly efficient and reliable method for large systems with cylindrical surfaces, such as vdW interactions with nanotubes.
Developing efficient and low-cost urea oxidation reaction (UOR) catalysts is a promising but still challenging task for environment and energy conversion technologies such as wastewater remediation and urea electrolysis. In this work, NiO nanoparticles that incorporated graphene as the NiO@Graphene composite were constructed to study the UOR process in terms of density functional theory. The single-atom model, which differed from the previous heterojunction model, was employed for the adsorption/desorption of urea and CO 2 in the alkaline media. As demonstrated from the calculated results, NiO@Graphene prefers to adsorb the hydroxyl group than urea in the initial stage due to the stronger adsorption energy of the hydroxyl group. After NiOOH@Graphene was formed in the alkaline electrolyte, it presents excellent desorption energy of CO 2 in the rate-determining step. Electronic density difference and the d band center diagram further confirmed that the Ni(III) species is the most favorable site for urea oxidation while facilitating charge transfer between urea and NiO@Graphene. Moreover, graphene provides a large surface for the incorporation of NiO nanoparticles, enhancing the electron transfer between NiOOH and graphene and promoting the mass transport in the alkaline electrolyte. Notably, this work provides theoretical guidance for the electrochemical urea oxidation work.
Engineering the electrolyte microenvironment represents an attractive route to tuning the selectivity of electrocatalytic reactions beyond catalyst composition and morphology. However, harnessing the full potential of this approach requires understanding the interplay between voltage, electrolyte composition, and adsorbate binding within the electrical double layer, which is absent from the usual theoretical approaches. In this study, we apply a recently developed density functional theory (DFT)–continuum approach based on the effective screening medium method and reference interaction site model (ESM–RISM) to explore electrolyte effects with an enhanced description of the electrochemical interface. Furthermore, applying this method to the binding of CO adsorbates in potassium-containing electrolytes on copper, a problem of direct relevance to CO 2 electroreduction to value-added products, we show that the interdependence of voltage and pH leads to an unexpected change in adsorption site preference on Cu(001) terraces. Our findings highlight the often-overlooked importance of the electrical double-layer structure for predicting catalyst operation.
Fe 2 O 3 hematite is a technologically important material with applications in energy storage as well as being the key phase formed in the rust of iron-based materials. Despite this central importance, there is much that is still unknown regarding the properties of defects in Fe 2 O 3 and consequently oxide growth. In this work, using screened hybrid density functional theory (HSE06), we consider the thermodynamics of vacancies in Fe 2 O 3 , considering the effects of ionic and electronic chemical potentials on both iron and oxygen vacancy formation. We find that, in the oxygen-rich limit, iron vacancies are easier to form, though the difference in formation energy between the two vacancies is only about 1 eV. In contrast, in the Fe-rich limit, oxygen vacancies have an extremely low formation energy, only 0.07 eV at mid-gap, and would spontaneously form as the Fermi level is reduced, while Fe vacancies require over 5 eV to form. Consistent with experiment, this indicates that Fe 2 O 3 is relatively easily reduced but not oxidized. However, the theoretical picture is very different when considering other exchange–correlation functionals (GGA + U or SCAN), emphasizing the critical role of the exchange–correlation functional in describing this system accurately.
The radiation-induced damages in bio-molecules are ubiquitous processes in radiotherapy and radio-biology, and critical to space projects. In this study, we present a precise quantification of the fragmentation mechanisms of deoxyribonucleic acid (DNA) and the molecules surrounding DNA such as oxygen and water under non-equilibrium conditions using the first-principle calculations based on density functional theory (DFT). Our results reveal the structural stability of DNA bases and backbone that withstand up to a combined threshold of charge and hydrogen abstraction owing to simultaneously direct and indirect ionization processes. We show the hydrogen contents of the molecules significantly control the stability in the presence of radiation. This study provides comprehensive information on the impact of the direct and indirect induced bond dissociations and DNA damage and introduces a systematic methodology for fine-tuning the input parameters necessary for the large-scale Monte Carlo simulations of radio-biological responses and mitigation of detrimental effects of ionizing radiation.
Abstract Chlorinated organics have been frequently detected in groundwaters, threatening the quality of drinking water supplies worldwide. A promising method for groundwater remediation involves reductive dechlorination (RD), in which chlorine atoms are sequentially removed and substituted by hydrogen, producing less harmful byproducts. In this paper, for the first time, RD research is reviewed in light of the growing incorporation of density functional theory (DFT) as a research tool. DFT has been used to uncover a variety of reaction properties for a range of relevant groundwater pollutants, including 1,2,3‐trichloropropane, hexachlorobenzene, and various dioxins. DFT models have revealed the role of surface interactions in driving the kinetics of catalytically driven RD. Mechanisms involved with biologically mediated RD have also been elucidated with insights gleaned from DFT. Issues and challenges for future research are also discussed.
Accurate hydrodynamic modeling for laser-direct-drive (LDD) inertial-confinement-fusion (ICF) relies on precise calculations of the electron thermal conduction in all target materials. The nonlocal stopping range of electrons in ICF plasmas directly influences thermal conduction; yet, no first principles model exists for the electron mean free path in the conduction-zone regime. This work utilized time-dependent stochastic density-functional theory (TD-sDFT) to calculate the electron stopping power in deuterium-tritium (DT) plasmas at (ρ, T) conditions relevant to the conduction zone and the compressed shell in ICF. Using a combination of our TD-sDFT data and already established analytical models, we developed and trained an artificial neural network to create a global model for the nonlocal electron deposition range, λ E . We compared our machine-learning (ML) based model for λ E to the currently-used modified-Lee-More model in LDD radiation-hydrodynamic codes, such as lilac, and saw an overall decrease in the deposition range. To understand the effects of λ E on LDD ICF implosion dynamics, we implemented the ML-based model into lilac; specifically, we looked at designs consistent with a current experiment on the OMEGA laser and for a newly designed LDD-ICF target for the future OMEGA-Next facility. In both cases, we saw an overall drop in predicted ablation pressure, peak areal density, and neutron yield due to the reduced thermal conduction (smaller λ E ) in DT plasmas. Comparisons with the experiment on OMEGA are also made.
High-nickel-content layered oxides are among the most promising electric vehicle battery cathode materials. However, their interfacial reactivity with electrolytes and tendency toward oxygen release (possibly yielding reactive 1 O 2 ) remain degradation concerns. Elucidating the most relevant (i.e., fastest) interfacial degradation mechanism will facilitate future mitigation strategies. We apply screened hybrid density functional (HSE06) calculations to compare the reaction kinetics of Li x NiO 2 surfaces with ethylene carbonate (EC) with those of O 2 release. Further, on both the (001) and (104) facets, EC oxidative decomposition exhibits lower activation energies than O 2 release. Our calculations, coupled with previously computed liquid-phase reaction rates of 1 O 2 with EC, strongly question the role of “reactive 1 O 2 ” species in electrolyte oxidative degradation. The possible role of other oxygen species is discussed. To deal with the challenges of modeling Li x NiO 2 surface reactivity, we emphasize a “local structure” approach instead of pursuing the global energy minimum.
The mineral schreibersite, e.g., Fe 3 P, is commonly found in iron-rich meteorites and could have served as an abiotic phosphorus source for prebiotic chemistry. However, atomistic calculations of its degradation chemistry generally require quantum simulation approaches, which can be too computationally cumbersome to study sufficient time and length scales for this process. In this regard, we have created a computationally efficient semiempirical quantum density functional tight binding (DFTB) model for iron and phosphorus-containing materials by adopting an existing semiautomated workflow that represents many-body interactions by linear combinations of Chebyshev polynomials. Here, we have utilized a relatively small training set to optimize a DFTB model that is accurate for schreibersite physical and chemical properties, including its bulk properties, surface energies, and water absorption. We then show that our model shows strong transferability to several iron phosphide solids as well as multiple allotropes of iron metal. Our resulting DFTB parametrization will allow us to interrogate schreibersite aqueous decomposition at longer time and length scales than standard quantum approaches, providing for more detailed investigations of its role in prebiotic chemistry on early Earth.
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 present a full experimental and theoretical study of an interesting series of polycyclic aromatic hydrocarbons, the oligorylenes. The absorption spectra of perylene, terrylene and quaterrylene in neutral, cationic and anionic charge states are obtained by matrix-isolation spectroscopy in Ne. The experimental spectra are dominated by a bright state that red shifts with growing molecular size. Excitation energies and state symmetry assignments are supported by calculations using time dependent density functional theory methods. These calculations also provide new insight into the observed trends in oscillator strength and excitation energy for the bright states: the oscillator strength per unit mass of carbon increases along the series.