Search NASA⌕ Search

SEARCH · Search NASA

Results for “Molecular Hydrogen Generation”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Water adsorption on vanadium oxide thin films in ambient relative humidity

Here in this work, ambient pressure x-ray photoelectron spectroscopy (APXPS) is used to study the initial stages of water adsorption on vanadium oxide surfaces. V 2p, O 1s, C 1s, and valence band XPS spectra were collected as a function of relative humidity in a series of isotherm and isobar experiments. Experiments were carried out on two VO 2 thin films on TiO 2 (100) substrates, prepared with different surface cleaning procedures. Hydroxyl and molecular water surface species were identified, with up to 0.5 ML hydroxide present at the minimum relative humidity, and a consistent molecular water adsorption onset occurring around 0.01% relative humidity. The work function was found to increase with increasing relative humidity, suggesting that surface water and hydroxyl species are oriented with the hydrogen atoms directed away from the surface. Changes in the valence band were also observed as a function of relative humidity. The results were similar to those observed in APXPS experiments on other transition metal oxide surfaces, suggesting that H 2 O-OH and H 2 O-H 2 O surface complex formation plays an important role in the oxide wetting process and water dissociation. Compared to polycrystalline vanadium metal, these vanadium oxide films generate less hydroxide and appear to be more favorable for molecular water adsorption.

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

Hydrogen-induced degradation dynamics in silicon heterojunction solar cells via machine learning

Abstract Among silicon-based solar cells, heterojunction cells hold the world efficiency record. However, their market acceptance is hindered by an initial 0.5% per year degradation of their open circuit voltage which doubles the overall cell degradation rate. Here, we study the performance degradation of crystalline-Si/amorphous-Si:H heterojunction stacks. First, we experimentally measure the interface defect density over a year, the primary driver of the degradation. Second, we develop SolDeg, a multiscale, hierarchical simulator to analyze this degradation by combining Machine Learning, Molecular Dynamics, Density Functional Theory, and Nudged Elastic Band methods with analytical modeling. We discover that the chemical potential for mobile hydrogen develops a gradient, forcing the hydrogen to drift from the interface, leaving behind recombination-active defects. We find quantitative correspondence between the calculated and experimentally determined defect generation dynamics. Finally, we propose a reversed Si-density gradient architecture for the amorphous-Si:H layer that promises to reduce the initial open circuit voltage degradation from 0.5% per year to 0.1% per year.

14 SOLAR ENERGY↗
DOE OSTI2025

Superionic-like diffusion in yttrium dihydride

For the next-generation high temperature microreactors, yttrium dihydride (YH 2 ) is an attractive solid state neutron moderator. Despite a number of recent investigations, the mechanism of hydrogen transport remains poorly understood. Experimental evaluations of diffusivity are inconclusive with large variations in diffusivities and activation energies. In this work, we perform ab initio molecular dynamics (AIMD) simulations on YH 2 for temperatures spanning 300 K to 1200 K. Our main finding is that YH 2 shows a superionic-like behavior with hydrogen atoms hopping from one native site to another above a characteristic temperature of 800 K. This correlated motion results in quasi-one-dimensional string-like displacements that enable the hydrogen atoms to diffuse rapidly. We confirm that the octahedral sites are mostly unoccupied, although channeling through them is the most favored pathway between lattice hops above 800 K. At the highest temperature of 1200 K, the string relaxation time is merely of the order of a few picoseconds, which indicates a liquid-like diffusive behavior. Based on the formation of spontaneous thermal vacancies, an order-disorder crossover temperature T α ~ 800 K is established for YH 2 with an activation energy of 0.83 eV for hydrogen diffusion in the superionic-like state.

Superionic-like Diffusion↗
DOE OSTI2024

Considerations for Defining G-Values for Aluminum-Clad Spent Nuclear Fuel

Sealed-canister dry storage of aluminum-clad spent nuclear fuel (ASNF) generated by research reactors is an alternative to current storage and disposition pathways as directed by the U.S. Department of Energy. The major challenge faced for this storage approach is radiolytic H 2 generation, including from the aluminum (oxy)hydroxide layers on the surface of ASNF. Experimental and modeling activities have been carried out to characterize the radiolytic yield as part of a DOE-sponsored research program to develop the technical basis for ASNF dry storage. The G-value is a commonly way to report results of radiolysis testing and is defined as the radiolytic yield of a species (e.g. molecular hydrogen) per unit radiation energy deposited into the material system. An independent technical review of the ASNF dry storage technical basis performed by Pacific Northwest National Laboratory raised questions about differences in G-value definitions used for experiments on ASNF surrogates consisting of aluminum samples with adherent (oxy)hydroxides compared to G-values reported in prior literature and how the magnitudes compared between different studies. Material systems resembling ASNF pose complications for measuring/defining G-values to predict the evolution of H 2 in a sealed canister, including i) accounting for radiolytic yields potentially arising from multiple sources, i.e., residual free (vapor), physisorbed, and chemisorbed/chemically bound waters; ii) deciding what portions of the multi-material system to include in the absorbed energy (radiation dose) calculation, considering possible energy exchange between materials as well as measurement limitations, and iii) capturing variations in G-value associated with non-linear yield vs. dose curves and/or dependence on the cover gas. This report summarizes previous literature information on radiolytic H 2 generation and associated G-values from mixed-material systems (generally oxides in contact with water or organic compounds) and from (oxy)hydroxides/hydrates to compare with the definitions and values for ASNF surrogate samples containing adherent aluminum (oxy)hydroxides.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

Syngas Production at Si Hybrid Photoelectrodes Modified with Re(I) and Mn(I) Tricarbonyl Phenanthroline Complexes Containing Reactive Aryl Azide Groups

We installed molecular CO 2 reduction (CO 2 R) catalysts directly onto Si (photo)electrodes. The highly reactive M(5-azido-1,10-phenanthroline)(CO) 3 X (where M = Mn or Re, X = Br or Cl) complexes readily bubbled when dissolved in polar organic solvents, in both the presence and absence of an ultraviolet light source. When placed on hydrogen-terminated Si (H-Si) and native silicon oxide (SiOx), similar amounts of the complex were attached to the surface under illumination (367 nm, 50–200 mW/cm 2 ) or in the dark. Surprisingly, these films revealed submonolayer coverages instead of the multilayered structures we expected. DFT analyses support monolayer formation, showing that the triplet-state nitrene of the complex is more energetically favorable than the singlet state. Using controlled-potential electrolysis experiments, we showed that Re- and Mn-containing films on pSi photoelectrodes generated small amounts of CO when exposed to 1 atm of CO 2 and 1 sun illumination. These amounts of CO were an order of magnitude greater than control surfaces, producing 5.59 × 10 –7 mol CO/h for Re(az-phen) and 7.83 × 10 –7 mol CO/h for Mn(az-phen) films. Much of the charge passed at the pSi electrodes was consumed by the competing hydrogen evolution reaction, which we attribute to the low molecular coverage and the presence of native oxide on the electrode surface after attachment. Finally, this work demonstrates the feasibility of reacting azide-containing ligands with Si surfaces. Still, it highlights the need for alternative ligand structures and reaction conditions to form multilayer films.

azides↗
DOE OSTI2022

Catalytic Chemical Recycling of Post-Consumer Polyethylene

Among commercial plastics, polyolefins are the most widely produced worldwide but have limited recyclability. Here, we report a chemical recycling route for the conversion of post-consumer high-density polyethylene (HDPE) into telechelic macromonomers suitable for circular reprocessing. Unsaturation was introduced into HDPE by catalytic dehydrogenation using an Ir-POCOP catalyst without an alkene acceptor. Cross-metathesis with 2-hydroxyethyl acrylate followed by hydrogenation transformed the partially unsaturated HDPE into telechelic macromonomers. The direct repolymerization of the macromonomers gave a brittle material due to the low overall weight-average molecular weight. Aminolysis of telechelic macromonomers with a small amount of diethanolamine increased the overall functionality. The resulting macromonomers were repolymerized through transesterification to generate a polymer with comparable mechanical properties to the starting post-consumer HDPE waste. Finally, depolymerization of the repolymerized material catalyzed by an organic base regenerated the telechelic macromonomers, thereby allowing waste polyethylene materials to enter a chemical recycling pathway.

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

Solvent molecules form surface redox mediators in situ and cocatalyze O 2 reduction on Pd

Solvent molecules influence the reactions of molecular hydrogen and oxygen on palladium nanoparticles. Organic solvents activate to form reactive surface intermediates that mediate oxygen reduction through pathways distinct from reactions in pure water. Kinetic measurements and ab initio quantum chemical calculations indicate that methanol and water cocatalyze oxygen reduction by facilitating proton-electron transfer reactions. Methanol generates hydroxymethyl intermediates on palladium surfaces that efficiently transfer protons and electrons to oxygen to form hydrogen peroxide and formaldehyde. Formaldehyde subsequently oxidizes hydrogen to regenerate hydroxymethyl. Water, on the other hand, heterolytically oxidizes hydrogen to produce hydronium ions and electrons that reduce oxygen. These findings suggest that reactions of solvent molecules at solid-liquid interfaces can generate redox mediators in situ and provide opportunities to substantially increase rates and selectivities for catalytic reactions.

Adams, Jason S.↗
DOE OSTI2025

Epoxy-based vitrimeric semi-interpenetrating network/MXene nanocomposites for hydrogen gas barrier applications

Herein, we report MXene-filled epoxy-based vitrimeric nanocomposites featuring a semi-interpenetrating network (S-IPN) to develop a hydrogen gas (H 2 ) barrier coating with self-healing characteristics for compressed H 2 storage applications. The reversible epoxy network was formed by synthesizing linear epoxy chains with pendent bis-hydroxyl groups using amino diol, which were then crosslinked with 1,4-benzenediboronic acid to generate dynamic boronic ester linkages. To achieve the S-IPN-type molecular arrangement, the epoxy chains were in situ crosslinked in the presence of poly(ethylene-co-vinyl alcohol) (EVOH), giving rise to a self-healing network (EEP) with a healing efficiency of 87%. Into the S-IPN vitrimer (EEP), a 2D platelet-type nanofiller MXene was incorporated to introduce a tortuous path for H 2 gas diffusion along with improved mechanical properties. The nanocomposite coating was applied to nylon 6 liner material, which is conventionally used in all-composite H 2 storage vessels. Further, the application of a 2 wt% MXene/EEP nanocomposite coating showed a permeability coefficient of 0.062 cm 3 mm m -2 d -1 atm -1 exhibiting ~96% reduction in gas permeability compared to uncoated nylon 6. The same nanocomposite exhibited a healing efficiency of 79%. Increasing the MXene loading to 10 wt% further reduced the permeability coefficient to 0.002 cm 3 mm m -2 d -1 atm -1 ; however, the healing efficiency decreased due to restricted chain mobility. In essence, the current work highlights the potential of vitrimeric S-IPN nanocomposite coatings for H 2 gas-barrier applications, enhancing safety and performance.

42 ENGINEERING↗
DOE OSTI2023

ORNL_AISD-Ex: Quantum chemical prediction of UV/Vis absorption spectra for over 10 million organic molecules

We performed calculations of electronic excitation energies and associated oscillator strengths based on the time-dependent density-functional tight-binding (TD-DFTB) method [1]. The SMILES (Simplified molecular-input line-entry system) strings of the molecules from the AISD HOMO-LUMO database [2] were converted to a 3D atomistic structure and stored in a PDB file after preliminary geometry optimization using the Merck Molecular Force Field (MMFF94) in RDKit [3,4]. The primary information stored in the PDB file archive consists of Cartesian coordinates for each atom of the molecule in their 3D location in space, along with summary information about the structure, sequence, and experiment. We then performed molecular geometry optimization using the density-functional tight-binding (DFTB) method [5] in the electronic ground state, followed by single-point excited states calculations, as described below. We note that, since RDKit employs a random choice for the generation of molecular conformers, the molecular geometries obtained in this dataset could be different from the ones that were generated when the AISD HOMO-LUMO dataset was generated. The computed excitation energies and associated oscillator strengths can be converted to predict UV/Vis absorption spectra, where excitation energies correspond to absorption peak positions, and oscillator strengths are a good measure of the probability of absorption of visible or UV light in transitions between electronic ground and excited states. The conversion of SMILES strings to 3D Cartesian coordinates of fully DFTB-optimized molecules was successful for 10,502,904 out of 10,502,917 molecules. For these molecules, both geometry optimizations and excited states calculations were successful. The DFTB calculations did not complete for 13 molecules of the original AISD HOMO-LUMO dataset. We still provide information about the geometry of these molecules. The molecules are diverse for chemical compositions (which span 5 non-hydrogen elements: oxygen, carbon, nitrogen, fluorine, sulfur) and molecular size (the smallest molecule contains 5 non-hydrogen atoms, and the largest molecule contains 71 non-hydrogen atoms). The DFTB method [5] is an approximation to density functional theory (DFT), utilizing a minimal basis set in conjunction with a two-center approximation to the electronic Hamiltonian and overlap matrix elements. The DFTB total energy is the sum of an electronic and a repulsive energy contribution, and their calculation requires optimized electronic parameters and diatomic repulsive potential energy functions. All DFTB calculations were performed using the DFTB+ code [6] (version 21.2) and the wrapper for DFTB+ in the Atomic Simulation Environment (ASE) (version 3.22.1) [7], which performed an internal conversion of Cartesian coordinates from PDB to the .gen file format. For the geometry optimizations on the electronic ground state potential energy surface of the molecules, we have chosen the third-order DFTB (DFTB3) method [5c] and employed the matching 3ob set of electronic parameters and repulsive potentials [8]. The empirical γ-damping for hydrogen bond correction, and Grimme's D3 empirical dispersion correction with Becke-Johnson damping (D3(BJ)) [9] dispersion correction was included to improve the description of non-covalent interactions. For excited states single-point energy calculations, we employed the TD-DFTB method in conjunction with the DFTB2 method [5b] and the matching mio [5b,10] and halorg [11] parameter sets. We opted to request the simultaneous calculation of 50 excited states for singlet transition to investigate sufficient number of excited states, based on linear response theory using the Casida equation [Ref: T. A. Niehaus, S. Suhai, F. Della Sala, P Lugli, M. Elstner, G. Seifert, and Th. Frauenheim. Tight-binding approach to time-dependent density-functional response theory. Phys. Rev. B, 63:085108, 2001] and the ARPACK diagonalizer [R. B. Lehoucq, D. C. Sorensen, and C. Yang. Arpack users guide: Solution of large-scale eigenvalue problems by implicitly restarted arnoldi methods, 1997. 46, 51]. The dataset contains 1001 tar.gz files. Tar files are named as “ornl_aisd_ex_1.tar.gz†through “ornl_aisd_ex_1000.tar.gzâ€. Additionally, the 13 failed molecules are in “ornl_aisd_ex_unprocessed.tar.gzâ€. Except for the tar files listed below, each tar file contains 10,500 molecules. Tar files numbered 34, 121, 128, 352, 360, 429, 495, 509, 518, 627, 676, 668, and 862 contain 10,499 molecules each. The last tar file numbered 1000 contains 13,417 molecules. The total size of the uncompressed dataset is over 283 Gigabytes. The code for calculating the electronic excitation energies and statistical analysis of the dataset is provided at the following GitLab repository: https://github.com/ORNL/Analysis-of-Large-Scale-Molecular-Datasets-with-Python Calculating the UV spectrum of a molecule requires performing 3 main operations: 1. Converting the smiles string representation of a molecule into a geometric structure where each atom is assigned XYZ coordinates. The geometric structure is written to the file smiles.pdb. 2. Using smiles.pdb to compute the relaxed geometry of the molecule, which corresponds with the position of the atoms at the position of equilibrium at the ground state. This generates the files band.out, detailed.out, and geo_end.gen. 3. Using geo_end.gen to calculate the UV spectrum of the molecule which is written into the file EXC.DAT. Every molecule in the dataset has its own directory. The files contained in each molecule directory are as follows: 1. geo_end.gen 2. detailed.out 3. band.out 4. EXC.DAT 5. smiles.pdb REFERENCES [1] Niehaus, T. A.; Suhai, S.; Della Salla, F.; Lugli, P.; Elstner, M.; Seifert, G.; Frauenheim, Th. Tight-binding approach to time-dependent density-functional response theory. Phys. Rev. B, 2001, 63, 085108/1-9. [2] Blanchard, A.; Gounley, J.; Metha, K.; Yoo, P.; Irle, S. AISD HOMO-LUMO. DOI: 10.13139/ORNLNCCS/1869409 [3] RDKit: Cheminformatics and Machine Learning Software. 2013, [http://www.rdkit.org] [4] Tosco, P.; Stiefl, N. and Landrum, G. Bringing the MMFF force field to the RDKit: implementation and validation. J Cheminform. 2014, 6, 1–4. [5] a) Porezag, D.; Frauenheim, T.; Kohler, T.; Seifert, G.; Kaschner, Construction of tight-binding-like potentials on the basis of density-functional theory: Application to carbon, R. Phys. Rev. B 1995, 51, 12947-12957; b) Elstner, M.; Porezag, D.; Jungnickel, G.; Elsner, J.; Haugk, M.; Frauenheim, Th.; Suhai, S.; Seifert, G.; Phys. Rev. B 1998, 58, 7260-7268; c) Gaus, M.; Cui, Q.; Elstner, M. DFTB3: Extension of the Self-Consistent-Charge Density-Functional Tight-Binding Method (SCC-DFTB), J. Chem. Theory Comput. 2011, 7, 931-948; d) Cui, Q.; Elstner, M. Density functional tight binding: values of semi-empirical methods in an ab initio era, Phys. Chem. Chem. Phys. 2014, 16, 14368-14377. [6] Hourahine, B. et al. DFTB+, a software package for efficient approximate density functional theory based atomistic simulations, J. Chem. Phys. 2020, 152, 124101/1-19. [7] Larsen, A. H. et al. The atomic simulation environment—a Python library for working with atoms. J. Phys.: Cond. Matter 2017, 29, 273002. [8] Kubillus, M.; Kubar, T.; Gaus, M.; Rezac, J.; Elstner, M. Parameterization of the DFTB3 Method for Br, Ca, Cl, F, I, K, and Na in Organic and Biological Systems, J. Chem. Theory Comput. 2015, 11, 332-342. [9] Brandenburg, J. G.; Grimme, S. Accurate Modeling of Organic Molecular Crystals by Dispersion-Corrected Density Functional Tight Binding (DFTB), J. Phys. Chem. Lett. 2014, 5, 1785−1789. [10] a) Niehaus, T. A.; Elstner, M.; Frauenheim, Th.; Suhai, S. Application of an approximate density-functional method to sulfur containing compounds. J. Mol. Struct.: THEOCHEM 2001, 541, 185-94; b) Elstner, M.; Hobza, P.; Frauenheim, Th.; Suhai, S.; Kaxiras, E. Hydrogen bonding and stacking interactions of nucleic acid base pairs: A density-functional-theory based treatment. J. Chem. Phys. 2001, 114, 5149-55. [11] Kubar, T.; Bodrog, Z.; Gaus, M.; Köhler, C.; Aradi, B.; Frauenheim, Th.; Elstner, M. Parametrization of the SCC-DFTB Method for Halogens. J. Chem. Theory Comput. 2013, 9, 2939-49.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Low Temperature Oxidation of Methylpropyl Ether

Ethers are potential high-performance fuels that can be produced from renewable carbon sources such as biomass, but their combustion properties are not clearly understood. The cetane numbers (CN) for ethers, a measure of autoignition tendency, can range from less than 10 to greater than 100, and predictive tools based upon the molecular structures are desirable in order to focus fuel development research. In this study, we used molecular modeling, chemical kinetic modeling, and experimental measurements to understand the important reactions that lead to the autoignition of methylpropyl ether (MPE), a model for ethers containing alkyl chains long enough for intramolecular hydrogen transfer. These reactions are known to be important in the autoignition of alkanes and ethers and we propose to investigate the relationship between alkyl and ether reaction pathways. A reaction mechanism for MPE was created using MIT's Reaction Mechanism Generator (RMG) and G4 quantum calculation results and models were tested using the ChemKin Pro suite of software. Results from atmospheric pressure, flow tube reactor experiments and rapid compression machine (RCM) experiments were used to test the kinetic model and suggest refinements.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2021

Thermal gradient effect on helium and self-interstitial transport in tungsten

First-wall materials in a fusion reactor are expected to withstand harsh conditions, with high heat and particle fluxes that modify the materials microstructure. These fluxes will create strong gradients of temperature and concentration of diverse species. Besides the He ash and the hydrogenic species, neutron particles generated in the fusion reaction will collide with the material creating intrinsic defects, such as vacancies, self-interstitials atoms (SIAs), and clusters of such point defects. These defects and the He atoms will then migrate in the presence of the aforementioned gradients. In this study, we use nonequilibrium molecular dynamics to analyze the transport of He and SIAs in the presence of a thermal gradient in tungsten. We observe that, in all cases, the defects and impurity atoms tend to migrate toward the hot regions of the tungsten sample. The resulting species concentration profiles are exponential distributions, rising toward the hot regions of the sample, in agreement with irreversible thermodynamics analysis. For both He atoms and SIAs, we find that the resulting species flux is directed opposite to the heat flux, indicating that species transport is governed by a Soret effect (thermal-gradient-driven diffusion) characterized by a negative heat of transport that drives species diffusion uphill (from the cooler to the hot regions of the sample). Here, we demonstrate that the steady-state species profiles obtained accounting for the Soret effect vary significantly from those where temperature-gradient-driven transport is not considered and discuss the implications of such a Soret effect on the response to plasma exposure of plasma-facing tungsten.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Synthesis of TiO{sub 2} nanoparticles decorated with gold nanoclusters using pulsed laser ablation in liquid

Amorphous spherical titanium dioxide (TiO{sub 2}) nanoparticles, decorated gold nanoparticles, are shown to be formed under irradiation of commercial TiO{sub 2} nanopowders dispersed in aqueous solution of hydrogen tetrachloroaurate (HAuCl{sub 4}) by second-harmonic nanosecond pulses of an Nd : YAG laser. It is found that these hybrid nanostructures are formed during laser-induced remelting of the initial nanoparticles, stimulated by gold nanoclusters (providing more efficient absorption of visible laser radiation by titanium dioxide) recovered on their surface. The morphology and chemical composition of the newly formed hybrid nanomaterials are investigated in detail by electron microscopy, Raman spectroscopy, and energy-dispersive analysis. It is shown that the average size and number of gold nanoparticles recovered on the TiO{sub 2} surface can be controlled by varying the concentration ratio of the initial nanomaterial and HAuCl{sub 4}. A spectroscopic analysis of light scattering by single hybrid nanoparticles and the results of numerical calculation of the structure of electromagnetic fields near their surface indicate good prospects of these hybrid nanomaterials for various applications in modern optics, optoelectronics, and nanophotonics; for example, for designing chemo- and biosensor platforms and new-generation solar cells. (nanoparticles)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Mesoscale Frank–Kasper Crystal Structures from Dendron Assembly by Controlling Core Apex Interactions

Single-component polymeric materials open up a great potential for self-assembly into mesoscale complex crystal structures that are known as Frank-Kasper (FK) phases. Predicting the packing structures of the soft-matter spheres, however, remains a challenge even when the molecular design is precisely known. Here, we investigate the role of the molecules' enthalpic interaction in determining the low-symmetry crystal structures. To this end, we synthesize architecturally asymmetric dendrons by varying their apex functionalities and examine the packing structures of the second-generation (G2) dendritic wedges. Our work shows that weakening the hydrogen bonding of the dendron apex makes the particles softer and smaller, and leads to the formation of various FK structures at lower temperatures, including the new observation of a FK C14 phase in the cone-shaped dendron systems. As a consequence of the free energy balance between the particle's interfacial tension and the chain's stretching, various packing structures are mainly tuned by designing the hydrogen bonding interaction.

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

Prediction of chlorine and fluorine crystal structures at high pressure using symmetry driven structure search with geometric constraints

The high-pressure properties of fluorine and chlorine are not yet well understood because both are highly reactive and volatile elements, which have made conducting diamond anvil cell and x-ray diffraction experiments a challenge. Here, we use ab initio methods to search for stable crystal structures of both elements at megabar pressures. We demonstrate how symmetry and geometric constraints can be combined to efficiently generate crystal structures that are composed of diatomic molecules. Furthermore, our algorithm extends the symmetry driven structure search method by adding constraints for the bond length and the number of atoms in a molecule while still maintaining generality. As a method of validation, we have tested our approach for dense hydrogen and reproduced the known molecular structures of Cmca-12 and Cmca-4. We apply our algorithm to study chlorine and fluorine in the pressure range of 10 GPa–4000 GPa while considering crystal structures with up to 40 atoms per unit cell. We predict chlorine to follow the same series of phase transformations as elemental iodine from Cmca to Immm to Fm3¯m, but at substantially higher pressures. We predict fluorine to transition from a C2/c to Cmca structure at 70 GPa, to a novel orthorhombic and metallic structure with P42/mmc symmetry at 2500 GPa, and finally to its cubic analog form with Pm3¯n symmetry at 3000 GPa.

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

Buried Interfaces in Organic Photocathodes for H 2 Evolution: Fermi-Level Pinning and Recombination

Herein, we demonstrate how Fermi-level pinning at buried contacts impacts solar fuel generation in all-polymer photocathodes by systematically comparing the effects of work function, hydroxyl coverage, and hydrogen evolution using chemically modified indium tin oxide (ITO) supports. Photovoltages and net photocathode performance are improved when the ITO is passivated using phosphonic acids, independent of work function, suggesting that the passivation reduces Fermi-level pinning at the buried interface arising from blended heterojunction interactions with surface metal hydroxyls. Transient photovoltage decay reveals differences in recombination mechanisms, supported by light intensity-dependent measurements. Briefly, nonpassivated, hydrophilic contacts exhibit trap-assisted recombination, while passivated, hydrophobic contacts follow bimolecular recombination. We then investigate changes in electroactivity of hole-transfer processes as a function of scan rate and repetitive cycling using a diffusion-controlled molecular redox probe, analogous to a hole-only device achieved via the electrolyte. The nonpassivated buried contacts exhibit higher overpotentials for oxidation, indicative of hole injection/extraction barriers. We observe irreversible electron transfer via the hole-transport level of the blended heterojunction and a strong cycle dependence, consistent with changes in the hole trap state density. Passivation results in more reversible redox behaviors, consistent with more Ohmic-like contacts. Collectively, these results provide context toward the realization of durable organic photoelectrodes with optimized photovoltages and net solar-to-hydrogen conversion efficiencies via fundamental understanding of the rates of carrier generation, recombination, and transport in high-dielectric aqueous environments and opportunities to characterize buried interfaces under device-relevant electric fields.

buried interfaces↗
DOE OSTI2022

Improving Rare-Earth Mineral Separation with Insights from Molecular Recognition: Functionalized Hydroxamic Acid Adsorption onto Bastnäsite and Calcite

Enhancing the separation of rare-earth elements (REEs) from gangue materials in mined ores requires an understanding of the fundamental interactions driving the adsorption of collector ligands onto mineral interfaces. In this work, we examine five functionalized hydroxamic acid ligands as potential collectors for the REE-containing bastnäsite mineral in froth flotation using density functional theory calculations and a suite of surface-sensitive analytical spectroscopies. These include vibrational sum frequency generation, attenuated total reflectance Fourier transform infrared, Raman, and X-ray photoelectron spectroscopies. Differences in the chemical makeup of these ligands on well-defined bastnäsite and calcite surfaces allow for a systematic relationship connecting the structure to adsorption activity to be framed in the context of interfacial molecular recognition. Here we show how the intramolecular hydrogen bonding of adsorbed ligands requires the inclusion of explicit water solvent molecules to correctly map energetic and structural trends measured by experiments. We anticipate that the results and insights from this work will motivate and inform the design of improved flotation collectors for REE ores.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗