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.
A formalism for target fragment transport is presented with application to energy loss spectra in thin silicon devices. A nuclear data base is recommended that agrees well with the measurements of McNulty et al. using surface barrier detectors. High-energy events observed by McNulty et al., which are not predicted by intranuclear cascade models, are well represented by the present work.
Observations of solar Lyman alpha have been interpreted as indicating that the polar mass flux density is lower than the equatorial average. This led Lallement et al (1986) to make a parametric study of solar wind acceleration, along the lines of earlier the Munro-Jackson study (1977), in which they concluded that uncertainties in the polar mass flux were large enough to be consistent with two extreme opposites: (1) a substantial energy supply beyond classical thermal conduction is required; or (2) classical thermal conduction is adequate to drive the flow. This ambiguity has been clarified by Ulysses observations of the polar outflow (Phillips et al, 1994). The polar mass flux density lies in the middle of the range studies by Lallement et al (1986), which suggests that extended heating is going on out to at least approximately 5 AU. Independent, purely energetic arguments can be made to estimate the required coronal source (electron) temperature that would be required to account for the observed energy flux density. An electron temperature of at least 2 x 10(exp 6) K would be required for the classical conduction flux density to be comparable to the total energy flux density; such a high temperature is thought to be unlikely in a coronal hole. These arguments strongly suggest that some extended heating or momentum transfer mechanism is required to drive the solar wind from the polar coronal hole. A number of mechanisms are discussed.
Spin-crossover (SCO) complexes are materials that exhibit changes in the spin state in response to external stimuli, with potential applications in molecular electronics. It is challenging to know a priori how to design ligands to achieve the delicate balance of entropic and enthalpic contributions needed to tailor a transition temperature close to room temperature. Here, we leverage the SCO complexes from the previously curated SCO-95 data set [Vennelakanti et al. J. Chem. Phys. 159, 024120 (2023)] to train three machine learning (ML) models for transition temperature (T 1/2 ) prediction using graph-based revised autocorrelations as features. We perform feature selection using random forest-ranked recursive feature addition (RF-RFA) to identify the features essential to model transferability. Of the ML models considered, the full feature set RF and recursive feature addition RF models perform best, achieving moderate correlation to experimental T 1/2 values. We then compare ML T 1/2 predictions to those from three previously identified best-performing density functional approximations (DFAs) which accurately predict SCO behavior across SCO-95, finding that the ML models predict T 1/2 more accurately than the best-performing DFAs. In addition, we study ML model predictions for a set of 18 SCO complexes for which only estimated T 1/2 values are available. Upon excluding outliers from this set, the RF-RFA RF model shows a strong correlation to estimated T 1/2 values with a Pearson’s r of 0.82. In contrast, DFA-predicted T 1/2 values have large errors and show no correlation to estimated T 1/2 values over the same set of complexes. Overall, our study demonstrates slightly superior performance of ML models in comparison with some of the best-performing DFAs, and we expect ML models to improve further as larger data sets of SCO complexes are curated and become available for model training.
A 3D auditory display can potentially enhance information transfer by combining directional and iconic information in a quite naturalistic representation of dynamic objects in the interface. Another aspect of auditory spatial clues is that, in conjunction with other modalities, it can act as a potentiator of information in the display. For example, visual and auditory cues together can reinforce the information content of the display and provide a greater sense of presence or realism in a manner not readily achievable by either modality alone. This phenomenon will be particularly useful in telepresence applications, such as advanced teleconferencing environments, shared electronic workspaces, and monitoring telerobotic activities in remote or hazardous situations. Thus, the combination of direct spatial cues with good principles of iconic design could provide an extremely powerful and information-rich display which is also quite easy to use. An alternative approach, recently developed at ARC, generates externalized, 3D sound cues over headphones in realtime using digital signal processing. Here, the synthesis technique involves the digital generation of stimuli using Head-Related Transfer Functions (HRTF's) measured in the two ear-canals of individual subjects. Other similar approaches include an analog system developed by Loomis, et. al., (1990) and digital systems which make use of transforms derived from normative mannikins and simulations of room acoustics. Such an interface also requires the careful psychophysical evaluation of listener's ability to accurately localize the virtual or synthetic sound sources. From an applied standpoint, measurement of each potential listener's HRTF's may not be possible in practice. For experienced listeners, localization performance was only slightly degraded compared to a subject's inherent ability. Alternatively, even inexperienced listeners may be able to adapt to a particular set of HRTF's as long as they provide adequate cues for localization. In general, these data suggest that most listeners can obtain useful directional information from an auditory display without requiring the use of individually-tailored HRTF's.
In the context of recent observations, microphysical models, and laboratory data, a photochemical model of Titan's atmosphere, including updated chemistry focusing on rate coefficients and cross sections measured under appropriate conditions, has been developed to increase understanding of these processes and improve upon previous Titan photochemical models. The model employs a two-stream discrete ordinates method to characterize the transfer of solar radiation, and the effects of electron-impact, cosmic-ray deposition, and aerosol opacities from fractal and Mie particles are analyzed. Sensitivity studies demonstrate that an eddy diffusion profile with a homopause level of 850 km and a methane stratospheric mole fraction of 2.2% provides the best fit of stratospheric and upper atmosphere observations and an improved fit over previous Titan photochemical models. Lack of fits for C3H8, HC3N, and possibly C2H3CN can be resolved with adjustments in aerosol opacity. The model presents a benzene profile consistent with its detection in Titan's stratosphere [Coustenis et al., 2003], which may play an important role in the formation of Titan hazes. An electron peak concentration of 4200 cm(exp -3) is calculated, which exceeds observations by 20%, considerably lower than previous ionosphere models. With adjustments in aerosol opacities and surface fluxes the model illustrates that reasonable fits to existing observations are possible with a single eddy diffusion profile, contrary to the conclusions of previous Titan models. These results will aid in the receipt and interpretation of data from Cassini-Huygens, which will arrive at Titan in 2004 and deploy a probe into Titan's atmosphere in January 2005.
This project focuses first and foremost on metal sulfide geochemistry and mineralogy as controlled by a complex matrix of environmental factors. The principal investigator’s group aim to illuminate the metal-sulfide reaction mechanisms, rates, and pathways through systematic experimentation and data collection and analyzing the relationships between the characteristics of the produced metal sulfide solid-phase/aqueous complexes and the environmental factors. This understanding is essential for obtaining a full picture of the complex cycling patterns of single or multi metal species in sulfidic environments ranging from deep-see basins, hydrothermal vents, inland seas, terrestrial water bodies, to engineered remediation systems. The specific goal of this past project was to illuminate the reaction mechanisms and kinetics of metal anions and sulfide in mixed metal cation-metal anion-sulfide systems under various aqueous conditions (which resembled a range of naturally occurring euxinic settings). For the period of this contract, we investigated the molybdenum-iron-sulfide system, with an emphasis on the conditions that caused solid phase formation. We focused on quantifying the mobility/sequestration of molybdenum under each experimental condition and identified the changes of valence states for each involved element (i.e., Mo, Fe, and S) in the precipitate. We also proposed pathways for the electron transfer that occurred in aqueous chemistry. The major analytical tools used for this study include UV-visible light spectroscopy, transmission electron microscopy, X-ray photoelectron spectroscopy, and synchrotron-based X-ray absorption spectroscopy (access gained through facility proposals to the Canadian Light Source). The biggest finding of this project was that besides pH, the iron-sulfur chemistry has a dominant control of the thiolation kinetics and subsequent reduction of Mo(VI), which are likely prerequisites of molybdenum sequestration in anoxic conditions. The results have been written up as manuscript by the end of this project (see Phillips et al.). The experimental results of this project may be critical for advancing our understanding of (1) basic chemistry involving transition metals and reduced sulfur species, (2) the validity of certain geochemical proxies, and (3) the stability and evolution of euxinic geochemical environments. It is noted that the basic results obtained through this project also have implications for Mo-S cluster-based catalyst development in inorganic chemistry and materials sciences.
The charge-transfer rate coefficient for the reaction Si(3+)(3s(sup 2)S) + He yields products is measured by means of a combined technique of laser ablation and ion storage. A cylindrical radio-frequency ion trap was used to store Si(3+) ions produced by laser ablation of solid silicon targets. The rate coefficient of the reaction was derived from the decay rate of the ion signal. The measured rate coefficient is 6.27(exp +0.68)(sub -0.52) x 10(exp -10)cu cm/s at T(sub equiv) = 3.9 x 10(exp 3)K. This value is about 30% higher than the Landau-Zener calculation of Butler and Dalgarno and is larger by about a factor of 3 than the recent full quantal calculation of Honvault et al.
Development of Forged Ingot Niobium Technology: Worlds science frontier programs and SRF accelerator applications demand high performance and cost-effective SRF accelerator technology [1-8]. Fine-grain (FG) and Large-grain (LG) ingot niobium technologies have been very well developed and implemented in all the present-day accelerator projects. However, forged ingot niobium technology which is the focus of this development proposal will be much more cost-effective and expected to have several technical advantages. FG niobium sheet production is very complex involving more than ten processing steps making them prone to contamination. As a result, they are very expensive to produce and require stringent QA procedure to be ready for SRF cavity production. The accelerating cavity process steps are also numerous and require strict procedures in order to achieve high accelerating gradients and quality factors needed for science frontier programs. LG niobium disc production, directly sliced from the ingot, is relatively simple and straight forward to keep surface cleanness. The disc production cost is significantly low compared to FG niobium sheet production. However, there are (some) draw backs due to non-homogeneity of the grain boundary distribution, resulting in non-uniform mechanical properties and complex cavity fabrication, although the LG cavities achieve the expected high-gradient performance goals with lower cost. Medium-grain (MG) niobium disk production may be realized with a new approach/process, the disc directly sliced from the forged ingot, involves a simpler process steps contributing major production cost reduction [9]. These discs are expected to be superior as they tend to be homogenous with uniform sub millimeter grains and mechanical properties. We are eagerly looking forward to developing the forged ingot niobium SRF accelerator technology for the benefit of the world-wide science frontier programs, green energy subcritical nuclear energy systems and a wide variety of industrial applications including the production of radio isotopes and nuclear transmutation applications. Measurement of thermal characteristics of the forged ingot niobium: Measurement of the thermal diffusivity, D, of superconducting MG niobium is important to be understood in comparison with FG sheet and LG disc, as well as that of advanced composite material of Nb3Sn film sputtered on forged ingot Nb, which will be determined using transient pump-probe thermo-modulation [10,11]. Transient thermo-modulation is based on using an ultrafast laser pulse to heat the superconducting materials by a few K, then a synchronized laser pulse probes the reflectance of the heated material. For thin films on a substrate (e.g, 100-nm Nb3Sn on Nb), it takes <100 ps for the heat to reach the substrate by diffusion. Therefore, an ultrafast method is needed to probe D of the studied material. D will be measured from room temperature to liquid helium temperatures. Fundamental Research: The increase of rf loss in SRF cavities is related to trapping of residual magnetic field during the cavity cool-down. The study suggested that the micro structure plays the role in flux trapping sensitivity. The research includes the flux trapping and expulsion study as forged ingot niobium goes through several mechanical deformations, crystallization leading to the optimal SRF cavity performance. Furthermore, the cavity made from FG, MG, and LG sheet/disc will be compared to understand the optimal re-crystallization temperature that cavity needed to be heat treated which minimize the flux trapping and increase SRF cavity performance. References: 1. S. Belomestnykh, Overview of recent SRF developments for ERLs, presented at the SRF 2015, Whistler, BC, Canada (2015). 2. A. Yamamoto, M. Yamanaka and G. Myneni, Ingot Nb based SRF Technology for the International Linear Collider, in Science and Technology of Ingot Niobium for Superconducting Radio Frequency Applications, AIP Conf. Proc. 1687, 030005-1 ? 03005-6, 2015 3. Report of the Workshop on Energy and Environmental Applications of Accelerators, DOE Workshop Report. https://science.osti.gov/-/media/hep/pdf/Reports/2020/CASM_WorkshopReport.pdf? la=en&hash=AEB0B318ED0436B1C5FF4EE0FDD6DEB84C2F15B2 4. G. Ciovati, et al., Design of a cw, low-energy, high-power superconducting linac for environmental applications, Phys. Rev. Accel. Beams 21, 091601 (2018). 5. P. Dhakal, et al., Effect of high temperature heat treatments on the quality factor of a large-grain superconducting radio-frequency niobium cavity, Phys. Rev. ST Accel. Beams 16, 042001 (2013). 6. P. Kneisel et al., Review of ingot niobium as a material for superconducting radio frequency accelerating cavities, Nuclear Instruments and Methods in Physics Research A 774, 133 (2015). 7. G. Ciovati, P. Dhakal, and G. R, Myneni, Superconducting radio-frequency cavities made from medium and low-purity niobium ingots, Supercond. Sci. Technol. 29, 064002 (2016). 8. M. Drury et al., commissioning of the prototype C75 cavities in a CEBAF cryomodule, in Proc. of IPAC 2018, Vancouver, BC, Canada (2018) 9. Feasibility of forged-ingot niobium disc and SRF cavity fabrication technology recently demonstrated in cooperation of ATI, BSCE, and KEK (2020). 10. Hani E. Elsayed-Ali, ?Measurements of heat transport in thin films by ultrafast laser-based techniques,? 3rd International Conference on Thermal Issues in Emerging Technologies Theory and Applications, Cairo, Egypt, pp. 347 ? 350 (2010). DOI: 10.1109/ThETA17616.2010 11. W. M. G. Ibrahim, H. E. Elsayed-Ali, M. Schinn, and C. A. Bonner, Jr., ?Ultrafast investigation of electron dynamics in multilayer metals,? Int. J. Heat and Mass Transfer, 47(10?11), 2261?2268 (2004).
Magnetohydrodynamics (MHD) can be used to model capillary discharge waveguides in laser-wakefield accelerators. However, the predictive capability of MHD can suffer due to poor microscopic closure models. Here, we study the impact of electron heating and thermal conduction on the capillary waveguide performance as part of an effort to understand and quantify uncertainties in modeling and designing next-generation plasma accelerators. To do so, we perform two-dimensional high-resolution MHD simulations using an argon-filled capillary discharge waveguide with three different electron transport coefficients models. The models tested include (i) Davies et al., (ii) Spitzer, and (iii) Epperlein–Haines (EH). We found that the EH model overestimates the electron temperature inside the channel by over 20% while predicting a lower azimuthal magnetic field. Furthermore, the Spitzer model, often used in MHD simulations for plasma-based accelerators, predicts a significantly higher electron temperature than the other models suggest.
The paper presents a theoretical model for the response of a tissue-equivalent proportional counter (TEPC) irradiated with charged particles. Heavy ions and iron ions in particular constitute a significant part of radiation in space. TEPCs are used for all space shuttle and International Space Station (ISS) missions to estimate the dose and radiation quality (in terms of lineal energy) inside spacecraft. The response of the tissue-equivalent proportional counters shows distortions at the wall/cavity interface. In this paper, we present microdosimetric investigation using Monte Carlo track structure calculations to simulate the response of a TEPC to charged particles of various LET (1 MeV protons, 2.4 MeV alpha particles, 46 MeV/nucleon 20Ne, 55 MeV/nucleon 20Ne, 45 MeV/nucleon 40Ar, and 1.05 GeV/nucleon 56Fe). Data are presented for energy lost and energy absorbed in the counter cavity and wall. The model calculations are in good agreement with the results of Rademacher et al. (Radiat. Res. 149, 387-389, 1998), including the study of the interface between the wall and the sensitive region of the counter. It is shown that the anomalous response observed at large event sizes in the experiment is due to an enhanced entry of secondary electrons from the wall into the gas cavity.
Power-over-Fiber (PoF) technology has been used extensively in settings where high voltages require isolation from ground and electromagnetic isolation is critical. In cryogenic environments, PoF offers a reliable power transmission technology, leveraging optical fibers to transfer power with minimal system degradation. PoF technology excels in maintaining low noise levels and isolation when delivering power to sensitive electronic systems operating in extreme temperature ranges and high voltage environments. Here, in a novel application of PoF for a HEP detector, power is provided to photon detector modules located on a surface at ~300 kV with respect to ground in the planned DUNE experiment. This summary paper of the PoF talk at the 16th PISA Meeting on Advanced Detectors highlights the R&D effort of PoF in extreme conditions and underscores its capacity to revolutionize power delivery and management in critical applications offering a dependable solution with low noise, optimal efficiency, and superior isolation. The DUNE (Abi et al., 2020) experiment will soon deploy large liquid argon (LAr) time projection chambers (TPC) to detect neutrino interactions and other particle physics phenomena. In addition to the particle tracking provided by the TPC, photon detectors, powered by a first ever PoF system, in the cryostat will leverage the high scintillation light yield of LAr to provide crucial timing and additional calorimetric information.
Previously, using an incompressible von Kármán–Howarth formalism, the behavior of cross-scale energy transfer in magnetic reconnection and turbulence was found to be essentially identical to each other, independent of an external magnetic (guide) field, in the inertial and energy-containing ranges [Adhikari et al., Phys. Plasmas 30, 082904 (2023)]. However, this description did not account for the energy transfer in the dissipation range for kinetic plasmas. In this Letter, we adopt a scale-filtering approach to investigate this previously unaccounted-for energy transfer channel in reconnection. Using kinetic particle-in-cell simulations of antiparallel and component reconnection, we show that the pressure–strain interaction becomes important at scales smaller than the ion inertial length, where the nonlinear energy transfer term drops off. Also, the presence of a guide field makes a significant difference in the morphology of the scale-filtered energy transfer. These results are consistent with kinetic turbulence simulations, suggesting that the pressure strain interaction is the dominant energy transfer channel between electron scales and ion scales.
A recombination mechanism that significantly reduces charge collection from very dense ion tracks in silicon devices was postulated by Zoutendyk et al. The theoretical analysis presented here concludes that Auger recombination is such a mechanism and is of marginal importance for higher density tracks produced by 270-MeV krypton, but of major importance for higher density tracks. The analysis shows that recombination loss is profoundly affected by track diffusion. As the track diffuses, the density and recombination rate decrease so fast that the linear density (number of electron-hole pairs per unit length) approaches a non-zero limiting value as t yields infinity. Furthermore, the linear density is very nearly equal to this limiting value in a few picoseconds or less. When Auger recombination accompanies charge transport processes that have much longer time scales, it can be simulated by assigning a reduced linear energy transfer to the ion.
Exciplexes are excited-state complexes formed as a result of partial charge transfer from the donor to the acceptor species when one moiety of the donor–acceptor pair is electronically excited. The arene–amine exciplex formed between oligo-(p-phenylene) (OPP) and triethylamine (TEA) is of interest in the catalytic photoreduction of CO2 because it can compete with complete electron transfer to the OPP catalyst. Therefore, formation of the exciplex can hinder the generation of a radical anion OPP·− necessary for subsequent CO2 reduction. We report an implementation of a workflow automating quantum-chemistry calculations that generate and characterize an ensemble of structures to represent this exciplex state. We use FireWorks, Pymatgen, and Custodian Python packages for high-throughput ensemble generation. The workflow includes time-dependent density functional theory optimization, verification of excited-state minima, and exciplex characterization with natural transition orbitals, exciton analysis, excited-state Mulliken charges, and energy decomposition analysis. Fluorescence spectra computed for these ensembles using Boltzmann-weighted contributions of each structure agree better with experiment than our previous calculations based on a single representative exciplex structure [Kron et al., J. Phys. Chem. A 126, 2319–2329 (2022)]. The ensemble description of the exciplex state also reproduces an experimentally observed red shift of the emission spectrum of [OPP-4–TEA]* relative to [OPP-3–TEA]*. The workflow developed here streamlines otherwise labor-intensive calculations that would require significant user involvement and intervention.
Abstract Energetic electron precipitation (EEP) associated with pulsating aurora can transfer greater than 30 keV electrons from the outer radiation belt region into the upper atmosphere and can deplete atmospheric ozone via collisions that produce NOx and HOx molecules. Our knowledge of exactly how EEP occurs is incomplete. Previous studies have shown that pitch angle scattering between electrons and lower‐band chorus waves can cause pulsating aurora associated with EEP and that substorms play an important role. In this work, we quantify the timescale of chorus wave decay following substorms and compare that to previously determined timescales. We find that the chorus decay e‐folding time varies based on magnetic local time (MLT), magnetic latitude, and wave frequency. The shortest timescales occur for lower‐band chorus in the 21 to 9 MLT region and compares, within uncertainty, to the energetic pulsating aurora timescale of Troyer et al. (2022, https://doi.org/10.3389/fspas.2022.1032552 ) for energetic pulsating aurora. We are able to further support this connection by modeling our findings in a quasi‐linear diffusion simulation. These results provide observations of how chorus waves behave after substorms and add additional statistical evidence linking energetic pulsating aurora to substorm driven lower‐band chorus waves.
The dynamics of molecular excitonic systems are complicated by a competition between electronic coupling (which drives delocalization) and vibrational-electronic (vibronic) interactions (which tend to encourage electronic localization). A particular challenge of molecular systems is that they typically possess a large number of independent vibrations, with frequencies often spanning the entire spectrum of relevant electronic energy gaps. Recent spectroscopic observations and numerical simulations on a water-soluble chlorophyll-binding protein (WSCP) reveal a transition between two regimes of vibronic behavior, a Redfield-like regime in which low-frequency vibrations respond to a delocalized excitonic state, and a Förster-like regime where high-frequency vibrations act as incoherent excitations on individual pigments. Although numerical simulations can reproduce these effects, there is a need for a simple, systematic theory that accurately describes the smooth transition between these two regimes in experimental spectra. Here we address this challenge by generalizing the variational polaron transform approach of [Bloemsma et al., Chem. Phys. 481, 250 (2016)] to include arbitrary bath densities for systems with or without symmetry. We benchmark this theory against both numerical matrix-diagonalization methods and experimental 77 K fluorescence spectra for two WSCP variants, obtaining quite satisfactory agreement in both cases. Here, we apply this theory to offer an explanation for the large loss in apparent electronic coupling in the WSCP Q57K mutant and to examine the likely impact of the interplay between excitonic delocalization and vibrational localization on vibrational sideband shapes and apparent coupling strengths in high-resolution optical spectra for chlorophyll-protein complexes such as WSCP.
The overarching scientific goals of our multidisciplinary grant was to further expand understanding about the key microorganisms (players), metabolic strategies (processes), and interspecies relationships (interactions) involved in the formation and oxidation of methane in the environment. This research applied novel environmental metagenomics, transcriptomics, and proteomic techniques, state-of-the-art analytical imaging, stable isotope geochemistry, and reaction-transport modeling to address these goals and develop an ‘ecosystems level’ understanding of the factors which regulate microbial methane cycling in anoxic sedimentary ecosystems. For decades, it was believed that the obligate step in methanogenesis catalyzed by methyl coenzyme M reductase (Mcr) was limited to a specific branch of the archaeal Domain, formerly known as the Euryarchaeota. Less than a decade ago co-I Tyson’s team published a surprising metagenomic-based discovery of divergent Mcr genes in a novel uncultured phylum (Bathyarchaeota), catalyzing a major shift in thinking about the diversity of microorganisms that encode the potential for methane (or higher alkane) metabolism in anoxic environments (Evans et al., 2015). In our work here, we further expand on the groups of archaea harboring the genomic potential for methane or hydrocarbon metabolism using environmental metagenomics and new gene targeted bioinformatics techniques. We additionally advanced understanding about the terminal electron acceptors and metabolic potential supporting the anaerobic oxidation of methane (AOM) in terrestrial ecosystems, specifically focused on new lineages of ANME archaea capable of respiring manganese oxides with methane presumably using large extracellular multi-heme cytochrome complexes. New details about specific syntrophic mechanisms underlying the exchange of electrons during sulfate-coupled methane oxidation between ANME-2 archaea and their sulfate-reducing bacterial partners were also elucidated as part of this funded project. Through a series of experimental ‘omics and single cell stable isotope probing studies with incubated environmental sediment samples and a cultured model electrogenic microorganism combined with model-based predictions. Combined, this work provides strong support for the hypothesis of direct interspecies electron transfer (DIET) is the dominant syntrophic mechanism controlling the anaerobic oxidation of methane with sulfate over other proposed mechanisms and additionally illustrates important spatial constraints and the underlying physico-chemical factors influencing AOM syntrophic consortia structure for DIET and extracellular metal respiration. This collaborative multi-institutional project successfully advanced several of our milestone goals including the identification of new microbial players containing methyl coenzyme M reductases hypothesized to be central to methane or hydrocarbon cycling in anoxic environments and enhancing fundamental knowledge about the role extracellular electron transfer plays in the ecophysiology of methanotrophic archaea respiring metal oxides and in the physical and metabolic structuring of syntrophic interactions in methane-rich sedimentary ecosystems.
We study the effect of self-interaction errors on the barrier heights of chemical reactions. For this purpose we use the well-known Perdew-Zunger [J. P. Perdew and A. Zunger, Phys. Rev. B,23, 5048 (1981)] self-interaction-correction(PZSIC), as well two variations of the recently developed, locally scaled self-interaction correction (LSIC) [Zope et al.,J. Chem. Phys.151, 214108 (2019)] to study the barrier heights of the BH76 benchmark dataset. Our results show that both PZSIC and especially the LSIC methods improve the barrier heights relative to the local density approximation(LDA). The version of LSIC that uses the iso-orbital indicator z as a scaling factor gives a more consistent improvement than an alternative version that uses an orbital-dependent factor based on the ratio of orbital densities to the total electron density. We show that LDA energies evaluated using the self-consistent and self-interaction-free PZSIC densities can be used to assess density-driven errors. Furthermore, the LDA reaction barrier errors for the BH76 set are found to contain significant density-driven errors for all types of reactions contained in the set, but the corrections due to adding SIC to the functional are much larger than those stemming from the density for the hydrogen transfer reactions and of roughy equal size for the non-hydrogen transfer reactions.