Search NASA⌕ Search

SEARCH · Search NASA

Results for “Optimized energy density functionals”

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 OSTI2025

Bayesian optimization of entropy-stabilized metal fluoride conversion cathodes and their synthesis

Fluoride-based conversion cathodes are promising for next-generation Li ion batteries because of their high voltage and energy densities. High entropy earth abundant fluoride cathodes are attractive for sustainable battery because of the potentially higher cycling stability. While Ni and Co containing equimolar alloys stabilized by high entropy were considered previously, Ni and Co free non-equimolar metal alloy compositions might offer new opportunities to increase materials stability and average voltage. This opportunity was pursued using Bayesian optimization of average voltage for the computational prediction of (Fe x1 Cu x2 Zn x3 Mn x4 Mg x5 )F 2 cathode chemistries. Theoretical voltage, room temperature free energy, and electronic structure of potentially high voltage compositions were computed using density functional theory calculations. Here, the experimental realization of non-equimolar quinary Co and Ni free fluorides with stable rutile single phase and improved conversion reaction voltage by up to 28 % compared to the equimolar alloy demonstrated the power of Bayesian optimization methods for the computationally guided experimental synthesis of fluoride cathodes. Finally, non-equimolar cathode chemistries were found with cycling stability superior to known multi-element fluorides.

25 ENERGY STORAGE↗
DOE OSTI2023

High-throughput ab initio reaction mechanism exploration in the cloud with automated multi-reference validation

Quantum chemical calculations on atomistic systems have evolved into a standard approach to studying molecular matter. These calculations often involve a significant amount of manual input and expertise, although most of this effort could be automated, which would alleviate the need for expertise in software and hardware accessibility. Here, we present the AutoRXN workflow, an automated workflow for exploratory high-throughput electronic structure calculations of molecular systems, in which (i) density functional theory methods are exploited to deliver minimum and transition-state structures and corresponding energies and properties, (ii) coupled cluster calculations are then launched for optimized structures to provide more accurate energy and property estimates, and (iii) multi-reference diagnostics are evaluated to back check the coupled cluster results and subject them to automated multi-configurational calculations for potential multi-configurational cases. All calculations are carried out in a cloud environment and support massive computational campaigns. Key features of all components of the AutoRXN workflow are autonomy, stability, and minimum operator interference. We highlight the AutoRXN workflow with the example of an autonomous reaction mechanism exploration of the mode of action of a homogeneous catalyst for the asymmetric reduction of ketones.

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

Shadow energy functionals and potentials in Born–Oppenheimer molecular dynamics

In Born–Oppenheimer molecular dynamics (BOMD) simulations based on the density functional theory (DFT), the potential energy and the interatomic forces are calculated from an electronic ground state density that is determined by an iterative self-consistent field optimization procedure, which, in practice, never is fully converged. The calculated energies and forces are, therefore, only approximate, which may lead to an unphysical energy drift and instabilities. Here, we discuss an alternative shadow BOMD approach that is based on backward error analysis. Instead of calculating approximate solutions for an underlying exact regular Born–Oppenheimer potential, we do the opposite. Instead, we calculate the exact electron density, energies, and forces, but for an underlying approximate shadow Born–Oppenheimer potential energy surface. In this way, the calculated forces are conservative with respect to the approximate shadow potential and generate accurate molecular trajectories with long-term energy stabilities. We show how such shadow Born–Oppenheimer potentials can be constructed at different levels of accuracy as a function of the integration time step, δt, from the constrained minimization of a sequence of systematically improvable, but approximate, shadow energy density functionals. For each energy functional, there is a corresponding ground state Born–Oppenheimer potential. These pairs of shadow energy functionals and potentials are higher-level generalizations of the original “zeroth-level” shadow energy functionals and potentials used in extended Lagrangian BOMD. The proposed shadow energy functionals and potentials are useful only within this extended dynamical framework, where also the electronic degrees of freedom are propagated as dynamical field variables together with the atomic positions and velocities. The theory is quite general and can be applied to MD simulations using approximate DFT, Hartree–Fock, or semi-empirical methods, as well as to coarse-grained flexible charge models.

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

Active Learning Accelerated Discovery of Stable Iridium Oxide Polymorphs for the Oxygen Evolution Reaction

The discovery of high-performing and stable materials for sustainable energy applications is a pressing goal in catalysis and materials science. Understanding the relationship between a material’s structure and functionality is an important step in the process, such that viable polymorphs for a given chemical composition need to be identified. Machine-learning-based surrogate models have the potential to accelerate the search for polymorphs that target specific applications. Herein, we report a readily generalizable active-learning (AL) accelerated algorithm for identification of electrochemically stable iridium oxide polymorphs of IrO 2 and IrO 3 . The search is coupled to a subsequent analysis of the electrochemical stability of the discovered structures for the acidic oxygen evolution reaction (OER). Structural candidates are generated by identifying all 956 structurally unique AB2 and AB3 prototypes in existing materials databases (more than 38000). Next, using an active learning approach, we find 196 IrO 2 polymorphs within the thermodynamic amorphous synthesizability limit and reaffirm the global stability of the rutile structure. We find 75 synthesizable IrO 3 polymorphs and report a previously unknown FeF 3 -type structure as the most stable, termed α-IrO 3 . To test the algorithms performance, we compare to a random search of the candidate space and report at least a 2-fold increase in the rate of discovery. Additionally, the AL approach can acquire the most stable polymorphs of IrO 2 and IrO 3 with fewer than 30 density functional theory optimizations. Analysis of the structural properties of the discovered polymorphs reveals that octahedral local coordination environments are preferred for nearly all low-energy structures. Subsequent Pourbaix Ir–H 2 O analysis shows that α-IrO3 is the globally stable solid phase under acidic OER conditions and supersedes the stability of rutile IrO 2 . Calculation of theoretical OER surface activities reveal ideal weaker binding of the OER intermediates on α-IrO 3 than on any other considered iridium oxide. We emphasize that the proposed AL algorithm can be easily generalized to search for any binary metal oxide structure with a defined stoichiometry.

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

High-shade dryland agrivoltaic conditions enhanced carbon uptake and water-use efficiency in zucchini ( Cucurbita pepo )

Introduction: The increasing global demand for food and energy is intensifying land-use competition. Agrivoltaic systems are a multifunctional land-use approach that vertically integrates the production of agricultural crops and solar power on the same land area. Most food crops are adapted to full-sun conditions, and the physiological responses of these crops to the novel microclimate under solar panels remain poorly understood. We hypothesized that the microclimate beneath the high-density photovoltaic system would influence carbon uptake, water use, and yield outcomes of zucchini summer squash.Methods: We conducted a field experiment in a hot, semi-arid climate on zucchini (Cucurbita pepo). Plants were grown under an agrivoltaic system with a 75% ground cover ratio (GCR) and in a full-sun control plot, each with two irrigation regimes (100 and 50%). We measured leaf-level photosynthesis, microclimate variables, and fruit yield at plant maturity and throughout the growing season.Results: The agrivoltaic array reduced photosynthetically active radiation (PAR) by ~79%, resulting in a cooler (−1.1 °C), more humid environment with higher soil moisture. These microclimatic conditions enhanced midday photosynthesis and daily cumulative carbon uptake. However, fruit yield was consistently lower under the panels, indicating a shift in carbon allocation toward vegetative growth. Photosynthesis was primarily driven by PAR across treatments, while soil moisture significantly influenced photosynthesis only in the control plots, suggesting water limitation was alleviated under the panels.Discussion: These findings highlight a trade-off between improved physiological performance and reduced yield under high-density agrivoltaics. While the system buffered heat and drought stress and improved overall plant function, excessive shade reduced reproductive output. Optimizing panel density or selecting crops cultivated for non-fruit yields will be essential for balancing food production and energy generation in dryland agrivoltaic settings.

14 SOLAR ENERGY↗
DOE OSTI2020

A fundamental look at electrocatalytic sulfur reduction reaction

The fundamental kinetics of the electrocatalytic sulfur reduction reaction (SRR), a complex 16-electron conversion process in lithium–sulfur batteries, is so far insufficiently explored. In this study, by directly profiling the activation energies in the multistep SRR, we reveal that the initial reduction of sulfur to the soluble polysulfides is relatively easy owing to the low activation energy, whereas the subsequent conversion of the polysulfides into the insoluble Li 2 S 2 /Li 2 S has a much higher activation energy, contributing to the accumulation of polysulfides and exacerbating the polysulfide shuttling effect. We use heteroatom-doped graphene as a model system to explore electrocatalytic SRR. We show that nitrogen and sulfur dual-doped graphene considerably reduces the activation energy to improve SRR kinetics. Density functional calculations confirm that the doping tunes the p-band centre of the active carbons for an optimal adsorption strength of intermediates and electroactivity. This study establishes electrocatalysis as a promising pathway to tackle the fundamental challenges facing lithium–sulfur batteries.

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

Combined Experimental and Computational Study of Molybdenum and Niobium for Nuclear Sensor Application

Due to their novel electromagnetic and thermal properties, molybdenum (Mo) and niobium (Nb) become optimal temperature sensor materials for nuclear energy applications. We leveraged voltage recorded during a heat ramp to tune a computational method to predict the Seebeck electromotive force (EMF) of Mo and Nb. Using a combined Density Functional Theory (DFT) and Boltzmann Transport Equations (BTE) method the voltage was predicted but did not include the effects of temperature on atomic structure. Combining Ab Initio Molecular Dynamics (AIMD) and BTE included temperature effects on structure optimization and yielded voltages in a good agreement with experiment. Lanthanum (La) and Phosphorus (P) additives in Mo and Nb, respectively, could increase the EMF compared to those of the pure metals. The presence of oxygen (O) in Mo increases the EMF while O in Nb slightly reduces the EMF. Furthermore, our studies suggested that heat treatment-induced structural changes that lead to a reduction in voltage occur not only at the mesoscale as previously understood but also at the atomic scale.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2026

Structured for success: conjugated polymer binders with tailored composition and architecture for lithium-ion batteries

Conjugated polymer binders are replacing conventional binders in lithium-ion batteries. Herein, we examine how molecular engineering and hierarchical nanostructuring govern binder functionality and electrochemical performance. Lithium-ion batteries (LIBs) are the leading energy storage technology, yet enhancing their energy density and cycle life remains critical. Significant progress has been made in high-capacity anodes and high-voltage cathodes, but their performance is hindered by electrode degradation, where it is related to the behaviors of binders at the surface and interface. Conventional non-conductive binders like poly(vinylidene difluoride) (PVDF), combined with conductive additives, often fail to maintain electrical pathways under repeated volume changes. Alternatively, conjugated polymer binders have emerged as a superior alternative, simultaneously offering intrinsic conductivity, mechanical flexibility, and strong adhesion through π-conjugated backbones and functional groups. Their tunable molecular structure enables efficient electron/ion transport while mitigating electrode cracking. Additionally, the development of hierarchically ordered nanostructures in conjugated polymer binder can further enhance their electrochemical performance. This review examines the design principles of conjugated polymer binders, focusing on molecular engineering and nanostructural control to optimize their performance in high-loading electrodes, such as silicon-based anodes. By addressing key challenges in binder functionality, these advanced materials pave the way for next-generation high-energy-density LIBs.

Jin, Xiuyu↗
DOE OSTI2023

Density functional modeling of the binding energies between aluminosilicate oligomers and different metal cations

Interactions between negatively charged aluminosilicate species and positively charged metal cations are critical to many important engineering processes and applications, including sustainable cements and aluminosilicate glasses. In an effort to probe these interactions, here we have calculated the pair-wise interaction energies (i.e., binding energies) between aluminosilicate dimer/trimer and 17 different metal cations M n+ (M n+ = Li + , Na + , K + , Cu + , Cu 2+ , Co 2+ , Zn 2+ , Ni 2+ , Mg 2+ , Ca 2+ , Ti 2+ , Fe 2+ , Fe 3+ , Co 3+ , Cr 3+ , Ti 4+ and Cr 6+ ) using a density functional theory (DFT) approach. Analysis of the DFT-optimized structural representations for the clusters (dimer/trimer + M n+ ) shows that their structural attributes (e.g., interatomic distances) are generally consistent with literature observations on aluminosilicate glasses. The DFT-derived binding energies are seen to vary considerably depending on the type of cations (i.e., charge and ionic radii) and aluminosilicate species (i.e., dimer or trimer). A survey of the literature reveals that the difference in the calculated binding energies between different M n+ can be used to explain many literature observations associated with the impact of metal cations on materials properties (e.g., glass corrosion, mineral dissolution, and ionic transport). Analysis of all the DFT-derived binding energies reveals that the correlation between these energy values and the ionic potential and field strength of the metal cations are well captured by 2nd order polynomial functions ( R 2 values of 0.99–1.00 are achieved for regressions). Given that the ionic potential and field strength of a given metal cation can be readily estimated using well-tabulated ionic radii available in the literature, these simple polynomial functions would enable rapid estimation of the binding energies of a much wider range of cations with the aluminosilicate dimer/trimer, providing guidance on the design and optimization of sustainable cements and aluminosilicate glasses and their associated applications. Finally, the limitations associated with using these simple model systems to model complex interactions are also discussed.

Gong, Kai↗
DOE OSTI2021

Nitromethane Decomposition via Automated Reaction Discovery and an Ab Initio Corrected Kinetic Model

In the explore the systematic construction of kinetic models from in silico reaction data for the decomposition of nitromethane. Our models are constructed in a computationally affordable manner by using reactions discovered through accelerated molecular dynamics simulations using the ReaxFF reactive force field. The reaction paths are then optimized to determine reaction rate parameters. We introduce a reaction barrier correction scheme that combines accurate thermochemical data from density functional theory with ReaxFF minimal energy paths. We validate our models across different thermodynamic regimes, showing predictions of gas phase CO and NO concentrations and high-pressure induction times that are similar to experimental data. The kinetic models are analyzed to find fundamental decomposition reactions in different thermodynamic regimes.

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

Binding of Sulfates and Water to Monovalent Cations

The binding of the sulfate ligand group to monovalent cations in the presence of water is important for many systems. To understand the structure and energetics of sulfate complexes, we use density functional theory to study ethyl sulfate binding to the monovalent cations Li + , Na + , and K + , and to water. The free energies of binding and optimal structures are calculated for a range of the number of ethyl sulfates and waters. Without water, the most optimal structure for all the cations is bidentate binding by two ethyl sulfates, yielding a 4-fold coordination. With water, the lowest free energy structures also have two ethyl sulfates, but the coordination varies with cations. For complexes with water, the four oxygen atoms in the sulfate group enable multiple binding geometries for the cations and for hydrogen bonding with water. Many of these geometries differ in free energy by only a small amount (1–2 kcal/mol), meaning there will be multiple binding configurations in bulk solution. In comparison to the optimal structures for binding to the carboxylate group, there is more variation for binding to the sulfate group as a function of cation type and the number of waters. Further, the polarization of the atoms is significant and varies among the sulfate oxygen atoms. The water oxygen charge is often larger than that of sulfate oxygen, which plays a role in the preference for monodentate ligand binding to cations in the presence of water.

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

Simulating excited-state complex ensembles: Fluorescence and solvatochromism in amine-arene exciplexes

Exciplexes are excited-state complexes formed as a result of partial charge transfer from the donor to the acceptor species when one moiety of the donor–acceptor pair is electronically excited. The arene–amine exciplex formed between oligo-(p-phenylene) (OPP) and triethylamine (TEA) is of interest in the catalytic photoreduction of CO2 because it can compete with complete electron transfer to the OPP catalyst. Therefore, formation of the exciplex can hinder the generation of a radical anion OPP·− necessary for subsequent CO2 reduction. We report an implementation of a workflow automating quantum-chemistry calculations that generate and characterize an ensemble of structures to represent this exciplex state. We use FireWorks, Pymatgen, and Custodian Python packages for high-throughput ensemble generation. The workflow includes time-dependent density functional theory optimization, verification of excited-state minima, and exciplex characterization with natural transition orbitals, exciton analysis, excited-state Mulliken charges, and energy decomposition analysis. Fluorescence spectra computed for these ensembles using Boltzmann-weighted contributions of each structure agree better with experiment than our previous calculations based on a single representative exciplex structure [Kron et al., J. Phys. Chem. A 126, 2319–2329 (2022)]. The ensemble description of the exciplex state also reproduces an experimentally observed red shift of the emission spectrum of [OPP-4–TEA]* relative to [OPP-3–TEA]*. The workflow developed here streamlines otherwise labor-intensive calculations that would require significant user involvement and intervention.

Chemistry↗
DOE OSTI2025

Global stellarator coil optimization with quadratic constraints and objectives

Most present stellarator designs are produced by costly two-stage optimization: the first for an optimized equilibrium, and the second for a coil design reproducing its magnetic configuration. Few proxies for coil complexity and forces exist at the equilibrium stage. Rapid initial state finding for both stages is a topic of active research. Most present convex coil optimization codes use the least square winding surface method by Merkel (NESCOIL), with recent improvements in conditioning, regularization, sparsity, and physics objectives. While elegant, the method is limited to modeling the norms of linear functions in coil current. We present QUADCOIL, a global coil optimization method that targets combinations of linear and quadratic functions of the current. It can directly constrain and/or minimize a wide range of physics objectives unavailable in NESCOIL and REGCOIL, including the Lorentz force, magnetic energy, curvature, field-current alignment, and the maximum density of a dipole array. QUADCOIL requires no initial guess and runs nearly $10$ 2 x faster than filament optimization. Integrating it in the equilibrium optimization stage can potentially exclude equilibria with difficult-to-design coils, without significantly increasing the computation time per iteration. QUADCOIL finds the exact, global minimum in a large parameter space when possible, and otherwise finds a well-performing approximate global minimum. It supports most regularization techniques developed for NESCOIL and REGCOIL. We demonstrate QUADCOIL’s effectiveness in coil topology control, minimizing non-convex penalties, and predicting filament coil complexity with three numerical examples.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Analytical gradients for molecular-orbital-based machine learning

We report molecular-orbital-based machine learning (MOB-ML) enables the prediction of accurate correlation energies at the cost of obtaining molecular orbitals. Here, we present the derivation, implementation, and numerical demonstration of MOB-ML analytical nuclear gradients, which are formulated in a general Lagrangian framework to enforce orthogonality, localization, and Brillouin constraints on the molecular orbitals. The MOB-ML gradient framework is general with respect to the regression technique (e.g., Gaussian process regression or neural networks) and the MOB feature design. We show that MOB-ML gradients are highly accurate compared to other ML methods on the ISO17 dataset while only being trained on energies for hundreds of molecules compared to energies and gradients for hundreds of thousands of molecules for the other ML methods. The MOB-ML gradients are also shown to yield accurate optimized structures at a computational cost for the gradient evaluation that is comparable to a density-corrected density functional theory calculation.

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

Development of a Reactive Force Field for Simulating Photoinitiated Acrylate Polymerization

Light-driven and photo-curable polymer based additive manufacturing (AM) has enormous potential due to its excellent resolution and precision. Acrylated radical chain-growth polymerized resins are widely used in photopolymer AM due to their fast kinetics, and often serve as a departure point for developing other resin materials for photopolymer-based AM technologies. For successful control of the photopolymer resins, the molecular basis of the acrylate free-radical polymerization has to be understood in detail. We present an optimized reactive force field (ReaxFF) for molecular dynamics (MD) simulations of acrylate polymer resins that captures radical polymerization thermodynamics and kinetics. The force field is trained against an extensive training set including density functional theory (DFT) calculations of reaction pathways along the radical polymerization from methyl acrylate to methyl butyrate, bond dissociation energies, and structures and partial charges of several molecules and radicals. We also found that it was critical to train the force field against an incorrect, nonphysical reaction pathway observed in simulations that used parameters not optimized for acrylate polymerization. As a result, the parameterization process utilizes a parallelized search algorithm, and the resulting model can describe polymer resin formation, crosslinking density, conversion rate, and residual monomers of the complex acrylate mixtures.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Hydrogen Dissociation Controls 1-Hexyne Selective Hydrogenation on Dilute Pd-in-Au Catalysts

Increasing the selectivity of the catalytic hydrogenation of alkynes to alkenes is of major importance for the processing of petrochemicals and the production of fine chemicals. Achieving high selectivity for alkene formation at high conversions, however, remains a long-standing challenge in heterogeneous catalysis. Here, the mechanism and origin of the high selectivity of dilute Pd-in-Au catalysts has been studied by a combination of first-principles calculations, microkinetic simulations, and isotopic exchange hydrogenation experiments using Pd 0.04 Au 0.96 nanoparticles embedded in raspberry colloid-templated silica. The Pd is predominantly in the form of isolated atoms, only surrounded by Au atoms, based on prior studies. The simulations indicate that the rate-limiting process for 1-hexyne hydrogenation on Pd monomers in Au(111) is H 2 dissociation, which has a large free energy barrier of 0.86 eV at 363 K and 0.2 bar of H 2 . The C–H bond formation steps, on the other hand, proceed with lower barriers, which contrasts with previous studies of extended Pd catalysts. The microkinetic simulations identify the sizable H 2 dissociation barrier and the small barrier for the hydrogenation of 1-hexyne as key factors that lead to a high selectivity for the production of 1-hexene from 1-hexyne, even at high conversion. The unconventional H 2 dissociation limiting process in combination with the low coverage of weakly bound hydrocarbon intermediates explains the near-zero order of 1-hexyne found experimentally. Furthermore, the partial hydrogenation of 1-hexyne to form 1-hexene is shown to be an irreversible process from our isotopic exchange hydrogenation experiments and is explained by the strongly exothermic nature of the reaction. Diluting active species, like Pd, in a less active host metal, like Au, hence appears promising as a means of tuning the binding energy of reactants and altering reaction profiles, leading to distinct kinetic behavior for an optimal catalytic activity and selectivity. Furthermore, the combination of microkinetic modeling, density functional theory calculations, and isotopic exchange experiments is thus demonstrated to be an effective approach to modeling important catalytic phenomena.

1-hexyne↗
DOE OSTI2021

3D Printed, Low Tortuosity Garnet Framework For Beyond 500 Wh/kg Batteries

In this project, we developed LLZO garnet ink recipes and processes for 3D-printing highly ordered ionically conductive garnet porous structures on dense garnet separators. Using this technique, we are able to fabricate controlled architecture LLZO garnet solid-state electrolyte (SSE) trilayers for application in solid-state lithium batteries. The trilayer comprises a thin dense center layer sandwiched between a 3D-printed patterned porous layer and a random porous layer. The dense layer functions as the ionic separator between the anode and cathode. The random porous layer hosts the lithium-metal anode and provides the structural support. The 3D-printed SSE patterned porous layer hosts the cathode, providing continuous, low tortuosity pathways for fast 3D Li+ transport through the cell while increasing the electrode/electrolyte interface area to decrease the interfacial resistance. Compared to the random porous structure, this ordered patterned structure possesses more vacant space for higher cathode loading without sacrificing ionically conducting capability, thus potentially greatly increasing the cell energy density. For demonstration purposes, we developed two patterns for the 3D-printed porous layer: grids and columns, for hosting sulfur and NMC cathode, respectively. The corresponding two types of cells were fabricated and tested, and have demonstrated achievement of theoretical discharge capacity without cathode calendaring. In addition, we developed a fundamental solid-state ionic and electronic transport model to optimize the 3D-printed structures for maximum energy and power density. The model was validated by experiment and provides the critical design criteria for achieving the >500 Wh/kg energy goal as function of C-rate.

25 ENERGY STORAGE↗
DOE OSTI2025

Interligand Coupling Drives Fast Triplet Energy Transfer Routes in PbS/Tetracene Quantum Dot Hybrids

The binding of photoactive organic ligands to inorganic quantum dots (QDs) creates a versatile hybrid architecture that allows access to photophysical processes such as efficient triplet exciton generation with near-infrared radiation. Here we report the subnanosecond generation of a hybrid triplet state with mixed ligand-QD character by replacing native oleate ligands on small PbS QDs with 5,12-tetracenepropiolic acid, a bifunctional ligand with two carboxylic acids that tends to lie face-on with the QD surface at low loadings. The face-on geometry engenders a regime of strong electronic coupling that is evident in steady-state absorption and hastens triplet energy flow by several orders of magnitude compared with more typical tetracene-based ligands exhibiting weak coupling. We further determined via Fourier transform infrared (FTIR) and supported by density functional theory (DFT)-based geometry optimizations that high ligand loading causes a shift in QD-ligand mutual disposition toward an edge-on geometry that instigates the formation of intermolecular excited states characterized by triplet excimer-like features in photoluminescence and transient absorption. Our results demonstrate the ability to control strongly coupled ligand-QD systems toward ultrafast generation of photophysically relevant species such as triplets that are valuable for photon upconversion and catalysis.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗