Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Density functional calculations”

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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Selective Hydrogenation of Furfural Acetone over a Cu Catalyst: Combined Theoretical and Experimental Study

The selective hydrogenation of biomass-derived hydroxymethyl furfural (HMF)-acetone-HMF (HAH) presents an alternative route to producing highly functionalized polyesters and polyurethanes. While HAH undergoes furan ring hydrogenation over Pd, Ru, and Ni catalysts, furan ring hydrogenation is not observed over Cu. Herein, we combine reaction kinetics experiments and density functional theory calculations to elucidate the selective hydrogenation behavior of HAH over Cu catalysts. We identified furfural acetone (FAc) as a suitable surrogate for modeling HAH hydrogenation over Cu surfaces and performed reaction kinetics experiments between temperatures of 313–393 K and a H 2 partial pressure of 55 bar. Similar to the behavior of HAH, hydrogenation of FAc follows a consecutive two-step hydrogenation pathway over Cu and does not undergo furan ring hydrogenation. The apparent activation energy barriers for hydrogenation of the aliphatic double bond (0.58 eV) and carbonyl (0.43 eV) of FAc measured in a continuous flow reactor setup are consistent with those reported in prior studies for batch HAH hydrogenation. Reaction orders with respect to each reactant, including H 2 and FAc, were determined to be nearly one. Density functional theory (DFT; GGA-PBE-D3) calculations on Cu(111) showed that the hydrogenation of the aliphatic double bond of FAc is more facile than the hydrogenation of the furan ring, which displays weak interactions with the Cu surface. We determined an apparent activation energy barrier for FAc hydrogenation (0.58 eV) that agreed with our DFT predictions (highest barrier for FAc hydrogenation of 0.57 eV). Finally, our DFT calculations further show that weak interactions between the furan ring and Cu surface are responsible for the selective hydrogenation behavior.

Biomass↗
DOE OSTI2024

Machine learning-enhanced design of lead-free halide perovskite materials using density functional theory

The investigation of emerging non-toxic perovskite materials has been undertaken to advance the fabrication of environmentally sustainable lead-free perovskite solar cells. This study introduces a machine learning methodology aimed at predicting innovative halide perovskite materials that hold promise for use in photovoltaic applications. The seven newly predicted materials are as follows: CsMnCl 4 , Rb 3 Mn 2 Cl 9 , Rb 4 MnCl 6 , Rb 3 MnCl 5 , RbMn 2 Cl 7 , RbMn 4 Cl 9 , and CsIn 2 Cl 7 . The predicted compounds are first screened using a machine learning approach, and their validity is subsequently verified through density functional theory calculations. CsMnCl 4 is notable among them, displaying a bandgap of 1.37 eV, falling within the Shockley-Queisser limit, making it suitable for photovoltaic applications. Here, through the integration of machine learning and density functional theory, this study presents a methodology that is more effective and thorough for the discovery and design of materials.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Investigation of Uranyl Perchlorate Anion Complexes in the Gas Phase via Infrared Multiphoton Dissociation and Collision‐Induced Dissociation

The infrared multiphoton dissociation spectrum of the gas-phase uranyl perchlorate anion ([UO 2 (ClO 4 ) 3 ] − ) was measured. Here, the asymmetric uranyl stretch lies somewhere between 930 and 1030 cm −1 , though it could not be deconvoluted from a closely located perchlorate stretching mode. The experimentally measured spectrum was in good agreement with spectra calculated using density functional theory with the B3LYP and TPSSh functionals and Grimme's D3 dispersion corrections with Becke–Johnson damping. The calculations suggest that the asymmetric uranyl stretch lies in the range of 970–980 cm −1 , about 20–30 cm −1 higher than for the analogous uranyl trinitrato complex and consistent with perchlorate as a weaker ligand than nitrate. Like the uranyl trinitrato anion, fragmentation of [UO 2 (ClO 4 ) 3 ] − by collision-induced dissociation resulted in loss of a ClO 3 • radical to form [UO 3 (ClO 4 ) 2 ] − . Infrared multiphoton dissociation measurements of [UO 3 (ClO 4 ) 2 ] − indicate that the structure of the product ion has an O − ligand bound to the uranyl in a T-arrangement and two perchlorate ligands, matching the findings from uranyl nitrato complexes. The uranyl stretch in this complex was slightly separated from the perchlorate modes and was measured at 906 cm −1 , slightly higher than for the analogous nitrate complex, again consistent with perchlorate as a weaker ligand than nitrate in the gas phase.

37 - INORGANIC, ORGANIC, PHYSICAL AND ANALYTICAL C↗
DOE OSTI2023

Precisely doping the surface of tin-based electrocatalysts for improved CO 2 conversion to liquid chemicals

Doping tin catalysts with sulfur can improve the electrochemical CO 2 conversion into formate/formic acid, but the lack of composition-dependent activity trends hinders further catalyst development. Here, we precisely controlled the composition of sulfur-doped Sn catalysts to show that sulfur doping only improves CO 2 conversion over a very narrow composition range, achieving maximum activity at 1.4 at% S. In situ Raman spectroscopy indicted working catalysts were in a primarily metallic state (e.g. S-Sn), and we achieved some of the highest reported partial current densities in both H-cell and full-cell electrolyzer configurations. Density Functional Theory calculations predicted S atoms preferentially occupied the catalyst surface and improved CO 2 reduction by localizing charge density at the catalyst/intermediate interface, which stabilized the *OCOH intermediate and lowered the CO 2 conversion thermodynamic barrier. Our work quantifies the composition-dependent influence of S dopants on Sn-based CO 2 reduction catalysts and provides a pathway for maximizing their CO 2 conversion activity.

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

Is Cs 2 TiBr 6 a promising Pb-free perovskite for solar energy applications?

In a quest for Pb-free perovskites suitable for solar energy applications, Cs 2 TiBr 6 has recently been reported as a promising compound, with appropriate optical and electrical properties as well as high stability under environmental stresses. In this study, we pursue investigation on this compound, demonstrating phase pure Cs 2 TiBr 6 powder formation using solution synthesis and providing complementary experimental characterization and theoretical calculations. An experimental absorption onset of around 2.0 eV is extracted and a weak broad photoluminescence is measured. Density functional theory calculations predict an indirect bandgap, parity-forbidden for both the direct and indirect transitions, which explains the weak and Stokes shifted luminescence. Additionally, we highlight the strong instability of Cs 2 TiBr 6 powder in ambient atmosphere. Therefore, our experimental results supported by theoretical calculations differ from previous results and raise doubts on the suitability of Cs 2 TiBr 6 in its pristine form for solar energy applications.

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

Cu[Ni(2,3-pyrazinedithiolate) 2 ] Metal–Organic Framework for Electrocatalytic Hydrogen Evolution

The application of metal–organic frameworks (MOFs) as electrocatalysts for small molecule activation has been an emerging topic of research. Previous studies have suggested that two-dimensional (2D) dithiolene-based MOFs are among the most active for the hydrogen evolution reaction (HER). Here, a three-dimensional (3D) dithiolene-based MOF, Cu[Ni(2,3-pyrazinedithiolate) 2 ] (1), is evaluated as an electrocatalyst for the HER. In pH 1.3 aqueous electrolyte solution, 1 exhibits a catalytic onset at –0.43 V vs the reversible hydrogen electrode (RHE), an overpotential (η 10 mA/cm 2 ) of 0.53 V to reach a current density of 10 mA/cm 2 , and a Tafel slope of 69.0 mV/dec. Interestingly, under controlled potential electrolysis, 1 undergoes an activation process that results in a more active catalyst with a 200 mV reduction in the catalytic onset and η 10 mA/cm 2 . It is proposed that the activation process is a result of the cleavage of Cu–N bonds in the presence of protons and electrons. Furthermore, this hypothesis is supported by various experimental studies and density functional theory calculations.

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

GPU Acceleration of Large-Scale Full-Frequency GW Calculations

Many-body perturbation theory is a powerful method to simulate electronic excitations in molecules and materials starting from the output of density functional theory calculations. By implementing the theory efficiently so as to run at scale on the latest leadership high-performance computing systems it is possible to extend the scope of GW calculations. Here, we present a GPU acceleration study of the full-frequency GW method as implemented in the WEST code. Excellent performance is achieved through the use of (i) optimized GPU libraries, e.g., cuFFT and cuBLAS, (ii) a hierarchical parallelization strategy that minimizes CPU-CPU, CPU-GPU, and GPU-GPU data transfer operations, (iii) nonblocking MPI communications that overlap with GPU computations, and (iv) mixed precision in selected portions of the code. A series of performance benchmarks has been carried out on leadership high-performance computing systems, showing a substantial speedup of the GPU-accelerated version of WEST with respect to its CPU version. Good strong and weak scaling is demonstrated using up to 25 920 GPUs. Finally, we showcase the capability of the GPU version of WEST for large-scale, full-frequency GW calculations of realistic systems, e.g., a nanostructure, an interface, and a defect, comprising up to 10 368 valence electrons.

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

Benchmarking boron carbide equation of state using computation and experiment

Boron carbide ( B 4 C ) is of both fundamental scientific and practical interest due to its structural complexity and how it changes upon compression, as well as its many industrial uses and potential for use in inertial confinement fusion (ICF) and high-energy density physics experiments. Here, we report the results of a comprehensive computational study of the equation of state (EOS) of B 4 C in the liquid, warm dense matter, and plasma phases. Our calculations are cross-validated by comparisons with Hugoniot measurements up to 61 megabar from planar shock experiments performed at the National Ignition Facility (NIF). Our computational methods include path integral Monte Carlo, activity expansion, as well as all-electron Green's function Korringa-Kohn-Rostoker and molecular dynamics that are both based on density functional theory. We calculate the pressure-internal energy EOS of B 4 C over a broad range of temperatures ( ~ 6 × 10 3 – 5 × 10 8 K) and densities (0.025–50 g / cm 3 ). We assess that the largest discrepancies between theoretical predictions are ≲ 5 % near the compression maximum at 1– 2 × 10 6 K. This is the warm-dense state in which the K shell significantly ionizes and has posed grand challenges to theory and experiment. By comparing with different EOS models, we find a Purgatorio model (LEOS 2122) that agrees with our calculations. The maximum discrepancies in pressure between our first-principles predictions and LEOS 2122 are ~ 18 % and occur at temperatures between 6 × 10 3 – 2 × 10 5 K, which we believe originate from differences in the ion thermal term and the cold curve that are modeled in LEOS 2122 in comparison with our first-principles calculations. To account for potential differences in the ion thermal term, we have developed three new equation-of-state models that are consistent with theoretical calculations and experiment. We apply these new models to 1D hydrodynamic simulations of a polar direct-drive NIF implosion, demonstrating that these new models are now available for future ICF design studies.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Transuranium Sulfide via the Boron Chalcogen Mixture Method and Reversible Water Uptake in the NaCu T S 3 Family

The behavior of 5f electrons in soft ligand environments makes actinides, and especially transuranium chalcogenides, an intriguing class of materials for fundamental studies. Due to the affinity of actinides for oxygen, however, it is a challenge to synthesize actinide chalcogenides using non-metallic reagents. Using the Boron Chalcogen Mixture (BCM) method, we achieved the synthesis of the transuranium sulfide NaCuNpS 3 starting from the oxide reagent, NpO 2 . Via the same synthetic route, the isostructural composition of NaCuUS 3 was synthesized and the material contrasted with NaCuNpS 3 . Single crystals of the U-analog, NaCuUS 3 , were found to undergo an unexpected reversible hydration process to form NaCuUS 3 ·xH 2 O (x ≈ 1.5). Here, a large combination of techniques was used to fully characterize the structure, hydration process, and electronic structures, specifically a combination of single crystal, powder, high temperature powder X-ray diffraction, extended X-ray ab-sorption fine structure, infrared, and inductively coupled plasma spectroscopies, thermogravimetric analysis, and density functional theory calculations. The outcome of these analyses enabled us to determine the composition of NaCuUS 3 ·xH 2 O and obtain a structural model that demonstrated the retention of the local structure within the [CuUS 3 ] – layers throughout the hydration-dehydration process. Band structure, density of states, and Bader charge calculations for NaCuUS 3 , NaCu-US 3 ·xH 2 O, and NaCuNpS 3 along with X-ray absorption near edge structure, UV-vis-NIR and work function measurements on ACuUS 3 (A = Na, K, and Rb) and NaCuUS 3 ·xH 2 O samples were carried out to demonstrate that electronic properties arise from the [CuTS 3 ] – layers and show surprisingly little dependence on the interlayer distance.

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

Two-Dimensional Molecular Charge Density Waves in Single-Layer-Thick Islands of a Dirac Fermion System

Charge density waves have been intensely studied in inorganic materials such as transition metal dichalcogenides however its counterpart in organic materials has yet to be explored in detail. Here we report the finding of a robust two dimensional charge density waves in molecular layers formed by α-(BEDT-TTF) 2 -I 3 on a Ag(111) surface. Low temperature scanning tunnelling microscopy images of a multi-layer thick α-(BEDT-TTF) 2 -I 3 on Ag(111) substrate reveal coexistence of 5a0 x 5a0 and √31 a0 x √31 a0 R9° charge density wave patterns commensurate with the underlying molecular lattice at 80 K. Both charge density wave patterns remain in nano-size molecular islands with just a single constituent molecular-layer thickness at 80 K and 5 K. Local tunneling spectroscopy measurements reveal the variation of the gap from 244 meV to 288 meV between the maximum and minimum charge density wave locations. Density functional theory calculations further confirm a vertical positioning of BEDT-TTF molecules in the molecular layer. While the observed charge density wave patterns are stable for the defect sites, they can be reversely switched for one molecular lattice site by means of inelastic tunnelling electron energy transfer with the electron energies exceeding 400 meV using a scanning tunneling microscope manipulation scheme.

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

Model Parameter Development for Complex Materials: Species-Specific Diffusion Barriers in 316 Stainless Steel from Systematic DFT Calculations

Vacancy-mediated diffusion barriers in 316 stainless steel have been systematically calculated using density functional theory to provide essential parameters for mesoscale microstructure evolution models. A statistical sampling approach employing 210 nudged elastic band calculations across multiple special quasi-random structures captures the effects of local chemical environments in this concentrated alloy. The computational methodology addresses challenges specific to chemically disordered systems, including proper magnetic treatment throughout multi-step calculations and validation against experimental structural properties. The calculated activation barriers reveal clear species-dependent diffusion behavior with the hierarchy Ni >> Fe ˜ Cr >> Mo. Nickel exhibits the highest barriers (0.74–1.31 eV, mean 1.045 eV), confirming its role as the slowest-diffusing major component. Iron and chromium show similar moderate barriers averaging 0.587 eV and 0.522 eV, respectively. Remarkably, molybdenum demonstrates exceptionally low barriers (0.12–0.28 eV, mean 0.194 eV), suggesting much higher mobility than previously recognized and potentially significant implications for precipitation kinetics and microstructure evolution. The barrier ranges remain consistent across different 316 SS compositions, supporting parameter transferability for modeling applications. The overall mean barrier of 0.64 eV provides a practical approximation for phase field simulations, while species-specific values enable detailed treatments of diffusion-controlled processes. This systematic approach establishes a validated framework for generating diffusion parameters in other concentrated alloys where experimental data are limited, while providing the first systematic set of species-specific barriers for predictive modeling of 316 stainless steel microstructure evolution.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Assessment of exchange-correlation functionals on oxygen vacancy formation energies of metal oxides

Oxygen vacancies are a common defect in oxides and play vital roles in many technological applications such as oxygen separation, catalytic reactors, solid oxide fuel cells, and solar thermochemical water splitting. The oxygen vacancy formation energy is directly related to the oxide reduction enthalpy and the vacancy concentration and, therefore, is a key quantity underlying these applications. The complexity of measuring oxygen vacancy formation energies experimentally suggests the utility of calculations based on density functional theory (DFT). However, the calculated results are strongly dependent on the exchange-correlation functionals and other parameters such as the Hubbard U . In this work, we compare the performance of the meta-GGA strongly constrained and appropriately normed (SCAN) functional and the commonly used semilocal generalized gradient approximation (GGA) functionals with experimental values of structural parameters, band gaps, magnetic structures, and oxygen vacancy formation energies for six representative oxides, i.e., CaMnO 3 , SrMnO 3 , LaMnO 3 , YMnO 3 , LaFeO 3 , and CeO 2 . Further, our results show that SCAN usually has better agreement with the experimental lattice constants and band gaps and larger magnetic moments than the commonly used GGA functionals. Although SCAN overestimates the oxygen vacancy formation energies of transition metal oxides and therefore requires unusually large Hubbard U values to reproduce the experimental reduction enthalpies, it does predict the correct oxygen vacancy formation energy of CeO 2 , which challenges the commonly used GGA and hybrid functionals. Our results underscore the challenges that exist in describing these complex oxides by DFT and may shed light on developing more accurate exchange-correlation functionals for oxygen vacancy formation energy calculations.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Interplay between copper redox and transfer and support acidity and topology in low temperature NH 3 -SCR

Low-temperature standard NH 3 -SCR over copper-exchanged zeolite catalysts occurs on NH 3 -solvated Cu-ion active sites in a quasi-homogeneous manner. As key kinetically relevant reaction steps, the reaction intermediate Cu II (NH 3 ) 4 ion hydrolyzes to Cu II (OH)(NH 3 ) 3 ion to gain redox activity. The Cu II (OH)(NH 3 ) 3 ion also transfers between neighboring zeolite cages to form highly reactive reaction intermediates. Via operando electron paramagnetic resonance spectroscopy and SCR kinetic measurements and density functional theory calculations, we demonstrate here that such kinetically relevant steps become energetically more difficult with lower support Brønsted acid strength and density. Consequently, Cu/LTA displays lower Cu atomic efficiency than Cu/CHA and Cu/AEI, which can also be rationalized by considering differences in their support topology. By carrying out hydrothermal aging to eliminate support Brønsted acid sites, both Cu II (NH 3 ) 4 ion hydrolysis and Cu II (OH)(NH 3 ) 3 ion migration are hindered, leading to a marked decrease in Cu atomic efficiency for all catalysts.

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

Computational exploration of biomedical HfNbTaTiZr and Hf 0.5 Nb 0.5 Ta 0.5 Ti 1.5 Zr refractory high-entropy alloys

Refractory high entropy alloys (RHEAs) have been proven to be a potential candidate in the biomedical field due to their balanced mechanical properties and biocompatible composition. Recent experimental findings show that RHEAs like HfNbTaTiZr and Hf 0.5 Nb 0.5 Ta 0.5 Ti 1.5 Zr have good mechanical properties such as high polarization and wear resistance than others which establish them as potential materials for biomedical application. In this work, we performed first-principles density functional theory calculations on the mechanical and thermal properties of HfNbTaTiZr and Hf 0.5 Nb 0.5 Ta 0.5 Ti 1.5 Zr. The predicted lattice constant, density, Young's modulus, and Vickers hardness are consistent with the available experimental report, which verifies the accuracy of the applied model. The thermal coefficient of linear expansion of both RHEAs has been investigated by utilizing the Debye theory. The present methods could be applied to study other future RHEAs on exploration of their physical properties.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Shear strength measurements and hydrostatic compression of rhenium diboride under high pressures

Shear strength measurements have been carried out on rhenium diboride, ReB 2 , to a pressure of 74 GPa using a Radial X-ray Diffraction (R-XRD) technique in a diamond anvil cell using platinum as an internal x-ray pressure standard. The R-XRD result has provided a unique insight into the deformation of hexagonal crystal lattice under non-hydrostatic compression and variation of shear strength with increasing pressure. From R-XRD data, we have estimated hydrostatic component of compression to determine an equation of state of rhenium diboride yielding a bulk modulus of K 0 = 366 ± 25 GPa with a pressure derivative $K^{'}_{0}$ = 4.3 ± 0.5 in good agreement with hydrostatic density functional theory calculations. Here, the average lower bound of shear strength (τ) from various diffraction planes was then calculated using the measured interplanar d-spacing (d m ) and hydrostatic component of d-spacing (d p ) to be shown to approach 6.7 ± 0.4 GPa at 70 GPa. Our results show that the anisotropic compression effects observed in ReB 2 under hydrostatic compression are correlated to electronic structure changes under compression as predicted by theoretical calculations.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Validation of the Cossee–Arlman mechanism for propylene oligomerization on Ni/UiO-66

Steady state rate expressions can be derived to distinguish the Cossee–Arlman and metallacycle mechanisms postulated for propylene oligomerization on nickel-based catalysts based on product selectivities, where product selectivities for the former are a function of olefin pressure because sequential coordination and insertion steps lead to independent mechanistic pathways for different hexene isomers. In contrast, the metallacycle mechanism presents pressure-independent product selectivities due to successive coordination prior to the kinetically relevant steps in each mechanism. In this work, steady state propylene oligomerization rates and selectivities were measured in the absence of an activator on nickel functionalized UiO-66 metal organic framework (MOF), Ni/UiO-66, to validate the Cossee–Arlman mechanism for light olefin oligomerization. In situ NO titrations reveal that ~5% of nickel sites were active during the reaction, and thus, not all nickel sites are relevant for catalysis. Propylene dimerization was first order in propylene pressure from 5 to 500 kPa with an apparent activation energy of ~20 kJ mol -1 from 453 to 493 K. Calculated apparent activation energies with density functional theory (DFT) calculations on cluster models of Ni/UiO-66 are in agreement with experiment to corroborate the Cossee–Arlman mechanism. Selectivities of hexene products and the ratio of hexene product selectivities on Ni/UiO-66 are in accordance with selectivity expressions derived from the Cossee–Arlman mechanism. In conclusion, analysis of product selectivities can be used more extensively to demarcate the Cossee–Arlman and metallacycle mechanisms for olefin oligomerization on metal-based catalysts.

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

Significantly increased Raman enhancement on defect-rich O-incorporated 1T-MoS 2 nanosheets

Improving charge transfer is the key to the performance of non-noble metal semiconductor-based surface enhanced Raman scattering (SERS) substrates. In this work, the O-incorporated 1T-MoS 2 nanosheets with rich sulfur defects (ID-MoS 2 ) are obtained by simple calcination of 1T-MoS 2 nanosheets in air atmosphere. Using rhodamine 6G (R6G) as typical probe molecules, ID-MoS 2 nanosheets show ultrahigh Raman enhancement effects with an enhancement factor of 1.24 × 10 7 due to sulfur defects and O incorporation in the 1T-MoS 2 lattice. First-principle density functional theory calculations suggest that the existence of sulfur defects and O incorporation significantly increase the Fermi energy level (E f ) and electronic density of states of ID-MoS 2 . Furthermore, O incorporation can enhance the interactions between the substrate and the adsorbed molecules through electrostatic and hydrogen bonding. All these improve the charge transfer resonance and result in the remarkable SERS activity of ID-MoS 2 nanosheets. This is the first study on the increasing SERS performance of semiconductor substrates by simultaneously employing defect and dopant incorporation. This study provides an approach to optimize the performance of semiconductor-based SERS substrates.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Molecular-Scale Characterization of Photoinduced Charge Separation in Mixed-Dimensional InSe–Organic van der Waals Heterostructures

Layered indium selenide (InSe) is an emerging two-dimensional semiconductor that has shown significant promise for high-performance transistors and photodetectors. The range of optoelectronic applications for InSe can potentially be broadened by forming mixed-dimensional van der Waals heterostructures with zero-dimensional molecular systems that are widely employed in organic electronics and photovoltaics. Here, we report the spatially resolved investigation of photoinduced charge separation between InSe and two molecules (C 70 and C 8 -BTBT) using scanning tunneling microscopy combined with laser illumination. We experimentally and computationally show that InSe forms type-II and type-I heterojunctions with C 70 and C 8 -BTBT, respectively, due to an interplay of charge transfer and dielectric screening at the interface. Laser-excited scanning tunneling spectroscopy reveals a ~0.25 eV decrease in the energy of the lowest unoccupied molecular orbital of C 70 with optical illumination. Furthermore, photoluminescence spectroscopy and Kelvin probe force microscopy indicate that electron transfer from InSe to C 70 in the type-II heterojunction induces a photovoltage that quantitatively matches the observed downshift in the tunneling spectra. In contrast, no significant changes are observed upon optical illumination in the type-I heterojunction formed between InSe and C 8 -BTBT. Density functional theory calculations further show that, despite the weak coupling between the molecular species and InSe, the band alignment of these mixed-dimensional heterostructures strongly differs from the one suggested by the ionization potential and electronic affinities of the isolated components. Self-energy-corrected density functional theory indicates that these effects are the result of the combination of charge redistribution at the interface and heterogeneous dielectric screening of the electron–electron interactions in the heterostructure. In addition to providing specific insight for mixed-dimensional InSe–organic van der Waals heterostructures, this work establishes a general experimental methodology for studying localized charge transfer at the molecular scale that is applicable to other photoactive nanoscale systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗