Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Chemical initiation”

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 559 records · Page 31
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Why Perovskite Thermal Stress is Unaffected by Thin Contact Layers

Metal halide perovskite photovoltaics have emerged as a high efficiency, low-cost alternative that can potentially rival or enhance conventional silicon technology. Despite exceptional initial power conversion efficiencies, achieving compliance with international standards and widespread adoption requires further enhancements to their operational stability. Notably, addressing mechanical strain and stress in brittle perovskites has emerged as a pivotal approach to mitigate chemical degradation and improve reliability during thermal cycling. Here, in this study, a popularized strain engineering strategy is investigated in which a high coefficient of thermal expansion (CTE) hole transport layer (i.e., PDCBT) is cast onto inorganic perovskite (CsPbI 2 Br) at 100 °C. Contrary to previously published results, the X-ray diffraction (XRD):Sin 2 ψ and substrate curvature measurement techniques show that the hole transport layer has no discernible impact on perovskite strain. The accuracy of the XRD:Sin 2 ψ method for measuring strain is highlighted in contrast to an analysis based on shifts of single XRD peaks which can be influenced by multiple artifacts. The findings in this study are in accordance with mechanics theory: thin layers are unable to induce significant strain changes in perovskite thin films as the force they apply is negligible compared to that applied by a thick and stiff substrate.

14 SOLAR ENERGY↗
DOE OSTI2021

Non-fullerene acceptor organic photovoltaics with intrinsic operational lifetimes over 30 years

Abstract Organic photovoltaic cells (OPVs) have the potential of becoming a productive renewable energy technology if the requirements of low cost, high efficiency and prolonged lifetime are simultaneously fulfilled. So far, the remaining unfulfilled promise of this technology is its inadequate operational lifetime. Here, we demonstrate that the instability of NFA solar cells arises primarily from chemical changes at organic/inorganic interfaces bounding the bulk heterojunction active region. Encapsulated devices stabilized by additional protective buffer layers as well as the integration of a simple solution processed ultraviolet filtering layer, maintain 94% of their initial efficiency under simulated, 1 sun intensity, AM1.5 G irradiation for 1900 hours at 55 °C. Accelerated aging is also induced by exposure of light illumination intensities up to 27 suns, and operation temperatures as high as 65 °C. An extrapolated intrinsic lifetime of > 5.6 × 10 4 h is obtained, which is equivalent to 30 years outdoor exposure.

14 SOLAR ENERGY↗
DOE OSTI2023

Microwave-assisted synthesis of iron sulfide motifs for electrochemical applications

Abstract The syntheses of FeS 2 and Fe 3 S 4 nanomaterials were optimized using a novel facile, surfactant-free, and microwave-assisted, one-pot synthesis method, run under ambient and reasonably mild reaction conditions. Synthetic parameters, such as metal precursor salt identity, reaction time, reaction temperature, metal:sulfur molar ratios, and solvent combinations, were all systematically investigated and optimized. A series of FeS 2 (pyrite) samples was initially fabricated using thioacetamide (TAA) as the sulfur precursor to generate a distinctive, uniform octahedra-based morphology. Switching the sulfur precursor from TAA to L-cysteine resulted in a corresponding transformation in not only chemical composition from FeS 2 to an iron thiospinel structure, Fe 3 S 4 (otherwise known as greigite), but also an associated morphological evolution from octahedra to nanosheet aggregates. The study of these materials has enabled crucial insights into the formation mechanisms of these materials under a relatively non-conventional microwave-assisted setting. Furthermore, in separate experiments, multi-walled carbon nanotubes (MWNTs) and graphene were added in with underlying metal sulfide species to create conductive Fe–S/MWNT composites and Fe–S/graphene composites, respectively. The method of addition of either MWNTs or graphene was also explored, wherein an ‘ ex-situ ’ synthetic procedure was found to be the least disruptive means of attachment and immobilization onto iron sulfide co-reagents as a means of preserving the latter’s inherent composition and morphology. The redox acidity for the parent material and associated composites demonstrates the utility of our as-developed synthetic methods for creating motifs relevant for electrochemical applications, such as energy storage.

25 ENERGY STORAGE↗
DOE OSTI2022

Organic Evaporation and Oxidation Testing in Support of Hanford Sample-and-Send

The Hanford site has approximately 56 million gallons of radioactive mixed waste stored in 177 underground storage tanks. The Hanford Waste Treatment and Immobilization Plant (WTP) is being built to treat and immobilize the tank waste. The baseline method for immobilization of Low Activity Waste (LAW) through the WTP is vitrification, but additional immobilization capacity is needed to supplement the initial LAW melters. An alternative cementitious waste form is being investigated for that future immobilization method to supplement vitrification. However, one impediment to a cementitious waste form is the presence of Land Disposal Restricted (LDR) organic chemicals in tank waste. This work evaluates potential avenues to eliminate that impediment to permit possible use of a cementitious waste form and work towards a decision whether additional LDR organic pretreatment would be required. Savannah River National Laboratory (SRNL) performed testing using simulants to examine evaporation as a method to remove some prevalent organics from LAW. Spiking the caustic LAW simulant with selected regulated organic chemicals found one that clearly decomposes because of caustic instability. Oxidation testing of other organic chemicals found some LDR organics degrade as desired and others are stable in the presence of peroxide and permanganate. In addition to studies with simulants, a literature review was performed to evaluate radiological stability of LDR organics. Descriptions of the experimental details, equipment, and results are included in this report. Evaporation testing consisted of preparing the LAW simulant, spiking that simulant with organic chemicals, and evaporating the mixture via differential distillation. The apparatus was a laboratory-scale vacuum evaporator operated at 60 ±5 torr absolute (vacuum evaporation) and also at atmospheric pressure. The LAW simulant represented the liquid expected to be retrieved from the Hanford tank farms at approximately 4.0 M [Na + ] total sodium ion concentration. The concentration of the organic chemicals added was significantly higher than typically found in the tank waste samples since the higher levels were necessary to assist in analytical measurement and tracking of the spiked species. Organic chemicals were chosen for the work with a consideration of how their volatility compares with that of methanol. This was done by comparing the ratio of the pure water Henry’s law coefficient (K h ) of methanol to that of the compound in question (hereafter termed the K h ratio), where ratios above unity indicated less volatility than methanol. Methanol was chosen because it is a common regulated chemical with relatively low volatility but which has been removed by evaporation in previous laboratory work. While organic separation results depend on evaporator design, laboratory experiments verified that organic partitioning to the overhead condensate stream by evaporation is a practical process. The work reported here found difficulties in quantitative analysis of the organic chemicals in aqueous samples. Most of the time there was insufficient analysis to close a mass balance for evaporator runs, but qualitative evidence of carryover was obtained. The methods were also able to show whether organic chemicals were susceptible or resistant to solution oxidation in permanganate or hydrogen peroxide tests.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2020

A Micelle Electrolyte Enabled by Fluorinated Ether Additives for Polysulfide Suppression and Li Metal Stabilization in Li-S Battery

The Li-S battery is a promising next-generation technology due to its high theoretical energy density (2600 Wh kg -1 ) and low active material cost. However, poor cycling stability and coulombic efficiency caused by polysulfide dissolution have proven to be major obstacles for a practical Li-S battery implementation. In this work, we develop a novel strategy to suppress polysulfide dissolution using hydrofluoroethers (HFEs) with bi-functional, amphiphlic surfactant-like design: a polar lithiophilic "head" attached to a fluorinated lithiophobic "tail." A unique solvation mechanism is proposed for these solvents whereby dissociated lithium ions are readily coordinated with lithiophilic "head" to induce self-assembly into micelle-like complex structures. Complex formation is verified experimentally by changing the additive structure and concentration using small angle X-ray scattering (SAXS). These HFE-based electrolytes are found to prevent polysulfide dissolution and to have excellent chemical compatibility with lithium metal: Li||Cu stripping/plating tests reveal high coulombic efficiency (>99.5%), modest polarization, and smooth surface morphology of the uniformly deposited lithium. Li-S cells are demonstrated with 1395 mAh g-1 initial capacity and 71.9% retention over 100 cycles at >99.5% efficiency-evidence that the micelle structure of the amphiphilic additives in HFEs can prohibit polysulfide dissolution while enabling facile Li + transport and anode passivation.

25 ENERGY STORAGE↗
DOE OSTI2026

A fresh perspective on the role of band bending, and related contributors, in light-driven production of electricity and chemicals

It is widely known that semiconductor-based solar energy conversion could power our planet. This is in part because high-quality semiconductor structures are unrivalled in their ability to separate photogenerated electrons and holes. One effective approach to achieving this photoinduced charge separation relies on a phenomenon known as “band bending”. But details to justify why band bending results in photoinduced charge separation are more complex than often appreciated. This underappreciation is an impediment to the rational, hypothesis-driven design of next-generation approaches to solar energy conversion. Herein we show, by means of derivations rooted in physical chemistry, that several phenomena – not just band bending – can facilitate photoinduced charge separation, and that each is influenced by nonequilibrium species concentration and a parameter, such as diffusion coefficient or rate coefficient, that introduces dynamics. To help visualize the impact of each phenomenon, we introduce plots that depict their contributions as free energy, force, flux, force constant, and rate. We reveal that spatial dopant distributions that define band bending are predictors of initial photogenerated species transport rates. But charge separation alone does not guarantee high-efficiency operation. A photogenerated change in energy that is freely available to do useful work is also essential, and is strongly dependent on semiconductor optical properties and reaction kinetics. Notably, this information reveals that specificity of interfacial chemical reactions – even when they are not preceded by charge separation elsewhere – can result in efficient solar energy conversion. We expect that this tutorial will guide researchers in their pursuit to uncover new mechanisms for light to perform useful work.

14 SOLAR ENERGY↗
DOE OSTI2024

NbZr_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-zirconium (Nb-Zr). 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 Zr. 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. 13. CHGCAR: contains the charge density data of a system. This data is crucial for analyzing electronic structures, calculating electrostatic potential, and studying the distribution of charge in a crystal or molecular system FILES contained in each subdirectory with name "case-N" where N ranges between 1 and 10, 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. {ID}.POSCAR: input file that defines the atomic structure of a system at the beginning of ID execution of geometry optimization with PREC=LOW 4. {ID}.CONTCAR: output file that provides the atomic positions and cell parameters at the end of ID execution of geometry optimization with PREC=LOW in the INCAR file 5. {ID}.OUTCAR: output file that contains detailed information about the progress of a calculation after the ID execution of geometry optimization that has been run with the precision variable set to PREC=Low in the INCAR file 6. rlx1.{ID}.out: file with diagnostic information about the execution of the ID execution of the geometry optimization with precision variable set to PREC=Low in the INCAR file 7. N{ID}.POSCAR: input file that defines the atomic structure of a system after the geometry optimization run at low precision. This represents the input for the ID execution of the second geometry optimization run with the precision variable set to PREC=Normal in the INCAR file 8. N{ID}.CONTCAR: output file that provides the atomic positions and cell parameters after the ID execution of the geometry optimization run with the precision variable set to PREC=Normal in the INCAR file 9. N{ID}.OUTCAR: output file that contains detailed information about the progress of a calculation after the ID execution of the geometry optimization run with the precision variable set to PREC=Normal in the INCAR file 10. rlx2.{ID}.out: file with diagnostic information about the ID execution of 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. 13. CHGCAR: contains the charge density data of a system. This data is crucial for analyzing electronic structures, calculating electrostatic potential, and studying the distribution of charge in a crystal or molecular system 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. & Hafner, J. Ab initio molecular dynamics for liquid metals. Phys. review B 47, 558 (1993). (2) Kresse, G. & 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. & 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. & 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. & Joubert, D. From ultrasoft pseudopotentials to the projector augmented-wave method. Phys. review b 59, 1758 (1999)

36 MATERIALS SCIENCE↗
DOE OSTI2024

TaZr_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 tantalum-zirconium (Ta-Zr). 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 Ta and Zr. 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 80, 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. 13. CHGCAR: contains the charge density data of a system. This data is crucial for analyzing electronic structures, calculating electrostatic potential, and studying the distribution of charge in a crystal or molecular system FILES contained in each subdirectory with name "case-N" where N ranges between 1 and 10 and between 81 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. {ID}.POSCAR: input file that defines the atomic structure of a system at the beginning of ID execution of geometry optimization with PREC=LOW 4. {ID}.CONTCAR: output file that provides the atomic positions and cell parameters at the end of ID execution of geometry optimization with PREC=LOW in the INCAR file 5. {ID}.OUTCAR: output file that contains detailed information about the progress of a calculation after the ID execution of geometry optimization that has been run with the precision variable set to PREC=Low in the INCAR file 6. rlx1.{ID}.out: file with diagnostic information about the execution of the ID execution of the geometry optimization with precision variable set to PREC=Low in the INCAR file 7. N{ID}.POSCAR: input file that defines the atomic structure of a system after the geometry optimization run at low precision. This represents the input for the ID execution of the second geometry optimization run with the precision variable set to PREC=Normal in the INCAR file 8. N{ID}.CONTCAR: output file that provides the atomic positions and cell parameters after the ID execution of the geometry optimization run with the precision variable set to PREC=Normal in the INCAR file 9. N{ID}.OUTCAR: output file that contains detailed information about the progress of a calculation after the ID execution of the geometry optimization run with the precision variable set to PREC=Normal in the INCAR file 10. rlx2.{ID}.out: file with diagnostic information about the ID execution of 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. 13. CHGCAR: contains the charge density data of a system. This data is crucial for analyzing electronic structures, calculating electrostatic potential, and studying the distribution of charge in a crystal or molecular system 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. & Hafner, J. Ab initio molecular dynamics for liquid metals. Phys. review B 47, 558 (1993). (2) Kresse, G. & 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. & 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. & 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. & Joubert, D. From ultrasoft pseudopotentials to the projector augmented-wave method. Phys. review b 59, 1758 (1999)

36 MATERIALS SCIENCE↗
DOE OSTI2024

Mechanochemically accelerated deconstruction of chemically recyclable plastics

Plastics redesign for circularity has primarily focused on monomer chemistries enabling faster deconstruction rates concomitant with high monomer yields. Yet, during deconstruction, polymer chains interact with their reaction medium, which remains underexplored in polymer reactivity. Here, we show that, when plastics are deconstructed in reaction media that promote swelling, initial rates are accelerated by over sixfold beyond those in small-molecule analogs. This unexpected acceleration is primarily tied to mechanochemical activation of strained polymer chains; however, changes in the activity of water under polymer confinement and bond activation in solvent-separated ion pairs are also important. Together, deconstruction times can be shortened by seven times by codesigning plastics and their deconstruction processes.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Air Carbon for Algae Production (AirCAP) – Expanding algae resource potential via direct (in-pond) air-CO2 capture.

Microalgae biofuel production potential in the United States is constrained by a lack of suitable cultivation locations within a close proximity to CO2 point-sources. Initial estimates suggest that removing the requirement of a concentrated CO2 source expands resource potential ten-fold, allowing the projection of production levels that approach levels of current U.S. diesel consumption. Decoupling cultivation locations from concentrated CO2 point-sources requires direct capture of air-CO2 into the pond bulk at a rate sufficient to meet maximum expected carbon assimilation rates (10 g C/m2-day and higher). Enhancement of the ingassing rate via chemical reaction, the direct reaction of dissolved CO2 with hydroxide ions, has been proposed as a mechanism to increase the ingassing rate of air-CO2. The central objectives of this project are to validate whether chemically enhanced air-CO2 ingassing rates in large raceway ponds can meet or exceed the level required to support high levels of algae productivity, then subsequently identify alkaliphilic strains able to withstand such conditions while maintaining a biomass productivity that enables process economics. Full scale ponds (~1 acre, unlined) will be used to validate ingassing rates as a function of pH expected at scale in a series of abiotic experiments. To date, mass-transfer dynamics have been characterized in ponds ranging from 1 to 1,000 m2. Initial results will be presented describing the mass transfer coefficient in 1-acre (~4,000 m2), unlined ponds. The maximum expected carbon ingassing rate will be measured during ingassing trials, in which pond pH is displaced from equilibrium via the addition of a strong base, with the subsequent decrease in pH (and increase in dissolved inorganic carbon) corresponding to the air-CO2 absorption rate. A mass-transfer model, validated with experimental results, will be used to project ingassing rates expected at pond scales (10 acres) required for production of commodity chemicals.

09 BIOMASS FUELS↗
DOE OSTI2022

Observation of the surface layer of lithium metal using in situ spectroscopy

In this work, we have investigated the surface of lithium metal using x-ray photoemission spectroscopy and optical spectroscopic ellipsometry. Even if we prepare the surface of lithium metal rigorously by chemical cleaning and mechanical polishing inside a glovebox, both spectroscopic investigations show the existence of a few tens of nanometer-thick surface layers, consisting of lithium oxides and lithium carbonates. When lithium metal is exposed to room air (~50% moisture), in situ real-time monitoring of optical spectra indicates that the surface layer grows at a rate of approximately 24 nm/min, presumably driven by an interface-controlled process. Our results hint that surface-layer-free lithium metals are formidable to achieve by a simple cleaning/polishing method, suggesting that the initial interface between lithium metal electrodes and solid-state electrolytes in fabricated lithium metal batteries can differ from an ideal lithium/electrolyte contact.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Preface: special issues on photosystem II

Here, the appearance of oxygenic photosynthesis 3–4 billion years ago is arguably one of the most important evolutionary developments which radically transformed both the biotic and abiotic landscapes of Earth. The presence of molecular oxygen enabled the evolution of complex lifeforms requiring large amounts of chemical energy which could only be produced by aerobic respiration. Additionally, molecular oxygen remodeled the geochemistry of the planet during its transition from a reducing to an oxidizing atmosphere during the Great Oxidation Event. Central to these processes is Photosystem II (PS II), the light energy-driven, water-plastoquinone oxidoreductase which initially evolved in cyanobacteria and subsequently was incorporated, along with Photosystem I (PS I), the cytochrome b 6 f complex and the Ndh complex, into all algae and higher plants by the process of endosymbiosis. This membrane protein complex releases molecular oxygen and protons as the byproducts of water splitting and produces plastoquinol, which provides reducing equivalents that are used sequentially by the cytochrome b 6 /f complex, PS I, and ultimately, as NADPH + H + , by the Calvin–Benson–Bassham cycle. Additionally, the protons released by PS II (and those pumped by the cytochrome b 6 /f complex and the NDH complex (during cyclic electron transfer)) form the proton motive force that drives ATP synthesis via the ATP Synthase, thus providing the chemical energy source for carbon fixation.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Evolution of PTCDA-derived seeds prior to graphene nanoribbon growth on Ge(001)

The synthesis of graphene nanoribbons (GNRs) can be realized via CH4 chemical vapor deposition (CVD) on substrates such as Ge(001) that promote highly anisotropic growth kinetics. Small polycyclic aromatic hydrocarbons (PAHs) are first sublimed onto Ge near room temperature to form graphene-like seeds that subsequently initiate GNR growth with CH4 exposure at 1173 K. The behaviors of PAHs in the time in between their sublimation onto Ge and GNR growth are unclear. Here, we study an archetypical PAH, perylene-3,4,9,10-tetracarboxyl acid dianhydride (PTCDA), on Ge(001) using both scanning tunneling microscopy (STM) and density functional theory (DFT) to characterize PAH conformation, surface diffusivity, and clustering, versus temperature. PTCDA becomes mobile above 673 K (consistent with a DFT calculated diffusion barrier of 2.01 eV). The mobile PTCDA molecules meet, cluster, and fuse at Ge step edges. These clusters have an apparent height of 0.4 nm, similar to small graphene islands, and grow laterally with increasing temperature – reaching 1.2 – 2.1 nm in width at 1173 K, consistent with the seed size extrapolated from CVD experiments. These results provide a plausible picture for how PTCDA forms seeds for anisotropic GNR CVD and show that PAHs with reduced surface diffusivity and inter-molecular reactivity are needed to enable more monodisperse PAH-seeded GNR synthesis.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Swift heavy ion irradiation effects in zirconium and hafnium carbides

The behavior of microcrystalline zirconium carbide (ZrC) and hafnium carbide (HfC) was studied under highly ionizing irradiation conditions at room temperature. The induced structural modifications were characterized via synchrotron-based X-ray diffraction experiments. Unit-cell expansion and buildup of microstrain were determined across a wide fluence range and linked to chemical compositions of the target material. The observed swelling resulting from irradiation with 198 MeV Xe ions in both carbide materials is characterized by two distinct mechanisms that operate within different fluence regimes. Unit-cell expansion initially proceeds by a direct-impact behavior that reaches saturation, followed at higher fluences by a second, linear swelling regime. The overall behavior, particularly the direct-impact regime, is similar for ZrC and HfC, with a more pronounced second defect accumulation process in HfC. Swelling in ZrC shows the same two distinct mechanisms upon irradiation with 198 MeV Xe ions and 946 MeV Au ions, but swelling induced by the lower-energy ions is greater across the entire fluence series. Accounting for the difference in energy deposition density between the two irradiation conditions reveals that the first swelling mechanism (direct-impact behavior) is likely related to the formation of more simple defects. In contrast, the second damaging mechanism at higher fluences (linear increase) cannot be fully explained by the induced energy density, and swelling remains somewhat higher for the low-velocity Xe irradiation. This may suggest that more complex defects and defect clusters are responsible for this swelling regime, with either their size and/or morphology modified at different energy densities.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Sub-picosecond Production of Solute Radical Cations in Tetrahydrofuran after Radiolysis

Ultrafast hole transfer from solvent radical cations produced by radiolysis with ~10 ps, 9 MeV electron pulses to solutes was investigated in tetrahydrofuran (THF). Because of rapid fragmentation of initially produced THF +• , solute radical cations are not expected and have not previously been reported. Using 9,9-Dihexyl-2,7-dibromofluorene (Br 2 F) at 5 to 1000 mM, Br 2 F +• was observed with radiation chemical yields up to G = 2.23 / 100 eV absorbed. While more than half of this was the result of direct solute ionization, the results highlight the importance of capture of holes from THF +• prior to solvation and fragmentation. The observed data show a time-resolution limited (15 ps) rise in transient absorption of Br 2 F +• , identical in form to reports of pre-solvated or dry electron capture in water and a few organic liquids, including THF. The results were thus interpreted with a similar formalism, finding C 37 = 1.7 M, the concentration at which 37% of holes escape capture. The yield of solvent hole capture can be accounted for by the formation of solvent holes adjacent to solute molecules reacting faster than they can fragment, however mechanisms such as delocalized holes or rapid hopping may play a role. Finally, low temperature results find over two times more capture, supporting the speculation that if THF +• was longer lived, the yield of capture in under 15 ps would have been at least 2x larger at 1 M Br 2 F, possibly capturing nearly all available holes from the solvent.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2025

Low-energy electron driven reactions in 2-bromo-5-nitrothiazole

Thiazole derivatives are biologically relevant molecules, used also in pharmaceutical applications. Herein, we report results for electron attachment to 2-bromo-5-nitrothiazole (BNT) in the gas phase. Employing two crossed electron–molecule beam experiments, we determined the efficiency curves of various fragment anions as a function of the initial electron energy between about 0 and 10 eV as well as the emission angle and kinetic energy distributions of Br− and NO2− ions formed from a resonance near 4 eV. The experiments were supported by quantum chemical calculations, exploring possible dissociation pathways along with their reaction energies. We also compared the electron attachment characteristics of BNT with those of the native thiazole molecule by performing electron attachment experiments and calculations for this molecule as well. Compared to thiazole, which is primarily degraded only by electrons with kinetic energies between about 5 and 10 eV, BNT is susceptible to low-energy electrons near 0 eV with enhanced cross section for (dissociative) electron attachment. However, although BNT offers two localization sites with high electron affinity (Br and NO2 moieties), we do not find the corresponding anions as the dominant negatively charged species formed upon electron attachment. Instead, the reaction channels with an abstraction of Br and NO2 as neutral radicals prevail, accompanied by the opening of the thiazole ring due to the relatively weak C–S bond.

Chen, Jiakuan (ORCID:0009000729608641)↗
DOE OSTI2021

Hard x-ray photoelectron spectroscopy: a snapshot of the state-of-the-art in 2020

Hard X-ray photoelectron spectroscopy (HAXPES) is establishing itself as an essential technique for the characterisation of materials. The number of specialised photoelectron spectroscopy techniques making use of hard X-rays is steadily increasing and ever more complex experimental designs enable truly transformative insights into the chemical, electronic, magnetic, and structural nature of materials. This paper begins with a short historic perspective of HAXPES and spans from developments in the early days of photoelectron spectroscopy to provide an understanding of the origin and initial development of the technique to state-of-the-art instrumentation and experimental capabilities. The main motivation for and focus of this paper is to provide a picture of the technique in 2020, including a detailed overview of available experimental systems worldwide and insights into a range of specific measurement modi and approaches. Herein, we also aim to provide a glimpse into the future of the technique including possible developments and opportunities.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Valence and magnetism in EuPd 3 S 4 and (Y, La) x Eu 1-x Pd 3 S 4

151 Eu Mössbauer spectroscopy shows that yttrium substitution in mixed-valent EuPd 3 S 4 drives the initial 50:50 mix of Eu 3+ and E 2+ towards pure Eu 2+ , whereas lanthanum substitution has the opposite effect, but only for substitution levels above 50%. Here, we find that total valence electron count and chemical pressure effects cannot account for the observed behavior; however, conserving the cell volume provides a consistent description of the changes in the Eu 2+ :Eu 3+ ratio. Remarkably, lanthanum substitution also leads to a clear transition from static mixed-valent behavior at lower temperatures to dynamic mixed-valent behavior at higher temperatures, with the onset temperature monotonically increasing with Eu content and extrapolating to a value of ~ 340 K for the pure EuPd 3 S 4 compound. Magnetic order persists at least as far as x=0.875 in both series, despite the drastic reduction in the amount of moment-carrying Eu 2+ ions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗