SEARCH · Search NASA

Results for “interface kinetics”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 91 records · Page 5

All sources NASA NTRS DOE OSTI

DOE OSTI2025

An Atomistic Study of Reactivity in Solid-State Electrolyte Interphase Formation for Li/Li7P3S11

Lithium metal batteries offer superior volumetric and gravimetric specific capacities compared to those based on traditional graphite anodes. Although advancements in solid-state electrolytes address safety concerns, challenges remain, particularly regarding interphase formation in lithium metal anodes. This work presents a computational framework based on high-throughput first-principles density functional theory and machine-learning interatomic potentials (MLIPs) including automated iterative, active learning to enable robust computational exploration of interphase formation between lithium metal anodes and an inorganic solid-state electrolyte. As a demonstration, we apply the framework to a Li/Li7P3S11 interface and find that it accurately identifies the experimentally observed, thermodynamically stable interphase products as well as their overall spatial arrangement within a heterogeneous, amorphous layered structure, with Li2S domains of nanocrystallinity. Our simulations show two stages, a fast and slow diffusion reaction regime, that corroborate the relative phase formation rate of Li x P, Li2S, and Li3P. Using the Onsager transport theory, we capture time-dependent ionic diffusion within the reacting interface, including cross-correlation effects. We found that cross-correlation effects between Li-P and P-S ionic motion significantly influence P-ion diffusion, making it highly sensitive to the local environment and potentially leading to "kinetic trapping" of Li-P phases. The passivation of the interface is shown as the ionic fluxes all approach zero, effectively halting interphase growth.

Diffusion↗

DOE OSTI2024

TChem-atm v1.0

SAND2024-11300O TChem-atm is a software library that was developed to solve complex kinetic models for atmospheric chemistry applications. TChem-atm interface employs a hierarchical parallelism design to exploit the massive parallelism available from modern computing platforms. It also supports gas atmospheric chemistry applications, e.g., the energy exascale earth system model. TChem can be used as a box model or coupled with a climate model to compute the time evolution of gas tracer species. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

Safta, Cosmin↗

DOE OSTI2024

A multiscale model to understand the interface chemistry, contacts, and dynamics during lithium stripping

A reversible Li-metal electrode, paired with a solid electrolyte, is critical for attaining higher energy density and safer batteries beyond the current lithium-ion cells. A stable stripping process may be even harder to attain as the stripping process will remove Li-atoms from the surface, and naturally reduce surface contact area, if not self-corrected by other mechanisms, such as diffusion and plastic deformation under an applied external stack pressure. Here, we capture these mechanisms occurring at multiple length- and time- scales, i.e., interface interactions, vacancy hopping, and plastic deformation, by integrating density functional theory (DFT) simulations, kinetic Monte Carlo (KMC), and continuum finite element method (FEM). By assuming the self-affine nature of multiscale contacts, we predict the steady-state contact area as a function of stripping current density, interface wettability, and stack pressure. We further estimate the exponential increase of overpotential due to contact area loss to maintain the same stripping current density. We demonstrate that a lithiophilic interface requires less stack pressure to reach the same steady-state contact area fraction than a lithiophobic interface. A “tolerable steady-state” contact area loss for maintaining stable stripping is estimated at 20 %, corresponding to a 10 % increase in overpotential. To constrain contact loss within the tolerance, the required stack pressure is 0.1, 0.5, and 2 times the yield strength of lithium metal for three distinct interfaces, lithiophilic Li/lithium oxide(Li2O), Li/lithium lanthanum zirconium oxide(LLZO), and lithiophoblic Li/lithium fluoride(LiF), respectively. The modeling results agree with experiments on the impact of the stack pressure quantitatively, while the discrepancy in stripping rate sensitivity is attributed to the simplifying interface interaction in our simulations. Overall, this multiscale simulation framework demonstrates the importance of electrochemical-mechanical coupling in understanding the dynamics of the Li/SE interface during stripping.

Feng, Min↗

DOE OSTI2025

How Solvent Mediates the Catalytic Fate of 2-Methoxyphenol in Hydrodeoxygenation Catalysis?

Kinetic and isotopic assessments combined with theoretical calculations shed light on the mechanistic differences in terms of the number of hydrogen addition steps and the electronic charges of reactive hydrogen in the catalytic sequence leading to the Ar−OCH 3 scission of 2-methoxyphenol occurring at various interfaces of vapor−Ru, cyclohexane−Ru, and water−Ru. At vapor−Ru and cyclohexane−Ru interfaces, the initial hydrogen adatom (H*) additions on the phenyl ring disrupt its aromaticity, as a step required to occur before kinetically relevant Ar−OCH 3 scission on an ensemble of Ru sites. The number of H* addition events prerequisite to the formation of the kinetically relevant Ar−OCH 3 scission transition state, however, is larger at the cyclohexane−Ru interface than those at the vapor−Ru and water−Ru interfaces. The preferential solvation of 2-methoxyphenol derived reaction precursors and H* weakens their adsorption on Ru sites, which leads to a decrease in the activation barrier for the third H* addition step and in turn a larger number of H* additions. At the water−Ru interface, water as a polar protic solvent promotes the Ar−OCH 3 scission by enabling a mechanistically distinct catalytic route an adjacent water molecule from the solvent layer assists with the intramolecular proton shuttling from the hydroxyl group (−OH) to the vicinal methoxy group (−OCH 3 ) of 2-methoxyphenol, evolving a highly charged transition state that is further stabilized by the water solvent layer, thereby leading to Ar−OCH 3 scission rate constants that are much higher than those at the vapor−Ru interface and, when comparing at the same H 2 chemical potential, higher Ar−OCH 3 scission turnovers. Furthermore, these distinct mechanistic details of the Ar−OCH 3 scission at various reaction interfaces, inferred from experiment and theory, illustrate how a solvent either solvates reaction intermediates and transition states or directly participates in the formation of kinetically relevant transition states, leading to changes in the free energy landscape, thus rewriting the catalytic fates in hydrodeoxygenation reactions.

09 BIOMASS FUELS↗

DOE OSTI2025

Off-Equilibrium Reactivity of Boron-Enriched Metal Diboride Surfaces in Electroreduction Conditions

Boron-based materials, featuring B-dependent reactivity and diverse phases, are emerging as promising catalyst systems. However, the catalytic mechanism on many borides remains poorly understood due to complex surface reconstructions under reaction conditions. Here, we investigate the MoB 2 surface in conditions of hydrogen evolution reaction in acidic media, using grand canonical global optimization, grand canonical density functional theory, ab initio molecular dynamics, free energy surface sampling, and an analytical model for electrochemical barrier evaluation. We propose a boron-enrichment strategy to tune the surface reactivity of the hexagonal face of MoB 2 . We reveal the dynamic nature of the B-enriched surface under H coverage and kinetic trapping of the system in the metastable regime with an extensive examination of the deactivation pathways. The metastable center B site on B-enriched surfaces, featuring buckled-up configuration and a usual relaxation effect, is found to be highly active toward HER via the Volmer–Heyrovsky mechanism. In conclusion, this work demonstrates how off-equilibrium behaviors can arise from the interplay between adsorbate coverage and surface reconstruction on a seemingly simple surface, and we present a theoretical framework and computational workflows to address these behaviors, along with other realistic complexities, in kinetics simulations.

Adsorption↗

DOE OSTI2025

Measuring the flatband potential in 2D semiconductors: Pitfalls and a possible SECCM solution

The flatband potential (V fb ) is a critical parameter in semiconductor electrochemistry, defining the potential at which no excess charge exists at the semiconductor/electrolyte interface. It serves as a key reference for interpreting charge transfer kinetics and current–voltage behavior. However, conventional methods like Mott–Schottky analysis fail for atomically thin 2D materials due to the breakdown of the depletion approximation. This perspective examines the limitations of traditional V fb measurements for 2D semiconductors and the experimental challenges that arise. To address these issues, we propose using scanning electrochemical cell microscopy (SECCM) to spatially resolve the potential of zero charge (V pzc ), equivalent to V fb . This approach mitigates sample heterogeneity issues, such as pinholes or multilayer defects, and offers a pathway to more accurate electrochemical characterization. Ultimately, this method will enhance understanding of current–potential behavior in 2D materials, supporting the design of advanced systems for photoelectrocatalysis, energy conversion, and sensing.

2D semiconductors↗

DOE OSTI2026

Temperature dependent interfacial microstructure in HIP-bonded Al6061/Zr joints

Pressing (HIP) at 450 °C and 560 °C under 103 MPa for 90 minutes to evaluate interfacial diffusion and microstructural evolution. Sound metallurgical bonding was achieved at both temperatures, with Zr serving as an effective diffusion barrier. Microstructural analysis revealed the the formation of a (Al,Si)3Zr intermetallic layer at the Al/Zr interface, whose thickness increased from ~100 nm at 450 °C to ~2100 nm at 560 °C, consistent with temperature-driven diffusion kinetics. The thin and continuous reaction layer observed at 450 °C indicates controlled interdiffusion and minimal reaction zone growth, favorable for maintaining interface stability and ductility. At 560 °C, enhanced atomic mobility led to pronounced intermetallic development, confirming diffusion-controlled growth behavior. These results provide valuable insights for optimizing Al/Zr interface design and diffusion barrier performance in advanced structural and energy systems.

Al/Zr interface↗

DOE OSTI2025

Interfacial Cation Arrangement Controls Electrocatalytic Kinetics in CO 2 Reduction

The identity of electrolyte cations is known to strongly influence electrocatalytic activity, but the relationship between their interfacial arrangement and observed performance remains poorly understood. Organic cations, with their molecular tunability, provide a powerful platform for systematically probing these effects. Here, we leverage phosphonium-based geminal dications to control interfacial cation arrangement and identify the variables that most strongly influence catalytic rates. As a case study, we examine CO 2 reduction to CO over polycrystalline silver electrodes in dry aprotic acetonitrile. Through a combination of rotating disk electrode measurements, electrochemical impedance spectroscopy, and molecular dynamics simulations, we decouple the effects of cation–electrode distance and interfacial cation density on catalytic rates. We find that smaller, more densely packed cations induce stronger interfacial electric fields, which lower the activation barrier for CO 2 adsorption and increase reaction rates. Using geminal phosphonium dications [C n (P mmm ) 2 ][ClO 4 ] 2 , we demonstrate that both the vertical and lateral positioning of organic cations within the electrical double layer independently affect reactivity. These results demonstrate that electrolyte cation identity primarily influences catalytic kinetics by determining how efficiently charge can be arranged at electrochemical interfaces. Altogether, our findings support an electrostatic view of cation effects in catalysis and provide design principles for next-generation electrolytes.

Cations↗

DOE OSTI2025

Dual-function lignin monomers enable high-performance graphene electrodes via interface confinement and proton transfer enhancement

Graphene oxide (GO)-based energy storage faces dual bottlenecks: unsustainable reduction methods and sluggish proton transfer kinetics. Here, we introduce a groundbreaking green strategy using lignin-derived vanillyl alcohol (VA) as a dual-function monomer to simultaneously address these challenges. By thermally annealing GO/VA films at mild temperatures (<100 °C), VA triggers an interface-confined reduction of GO while self-polymerizing into redox-active oligomers (P-VA) that intercalate between graphene layers. This dual role-reducing agent and proton highway enables a 3D conductive network with minimized graphene restacking, abundant redox sites, and rapid H + transport pathways. Density Functional Theory (DFT) reveals how P-VA optimizes proton dynamics, while the resulting rGO-P-VA4-T90 electrode achieves a record volumetric capacitance of 311.1F/cm 3 (777.8F/cm 2 ) and retains 87.8 % capacity after 10,000 cycles. Flexible solid-state supercapacitors deliver 94.2 μWh/cm 2 energy density at 63.8 μW/cm 2 , rivaling state-of-the-art devices. Furthermore, this work redefines sustainable graphene engineering, merging biomass valorization with high-performance energy storage in a scalable, eco-friendly paradigm.

Graphene oxide reduction↗

DOE OSTI2025

Probing the Impact of Vacancy Diffusion on Void Dynamics at the Lithium Metal–Solid Electrolyte Interface

Lithium (Li) metal-based solid-state batteries (SSBs) are considered promising candidates for next-generation energy storage due to their superior energy density and enhanced safety compared to conventional Li-ion systems. However, their practical application is limited by challenges such as void formation at the Li-solid electrolyte (SE) interface, which disrupts ion transport and accelerates interfacial degradation. This work investigates how the coupled effects of electro-dissolution kinetics and surface diffusion at the Li metal surface govern the evolution of interfacial morphology during stripping. This work examines the influence of three distinct surface diffusion modes, which are terrace diffusion, step diffusion, and interlayer diffusion, on maintaining interfacial stability. In addition, how the dominant surface diffusion mechanism can overcome the contact loss due to high reaction kinetics is explored. Furthermore, the roughness of the Li metal anode surface is quantified, and the influence of different diffusion mechanisms on the evolution of the dynamic solid–solid interface is examined. The critical role of temperature in enhancing Li surface diffusivity and expanding the regime of stable contact is highlighted. By identifying distinct regimes of interface stability, this study analyzes how non-uniform electrochemical dynamics dictate void morphology evolution and interfacial contact. These insights offer guiding principles for engineering robust Li–SE interfaces in SSBs.

lithium metal anode↗

DOE OSTI2024

Real-Time KMC Simulation of Vacancy-Mediated Intermixing in Au@Ag Octahedral Core–Cubic Shell Nanocrystals with Ab Initio-Guided Kinetics

Utilization of core–shell rather than monometallic nanocrystals (NCs) facilitates fine-tuning of NC properties for applications. However, compositional evolution via intermixing can degrade these properties prompting recent experimental studies. We develop an atomistic-level stochastic model for vacancy-mediated intermixing exploiting a formalism which allows incorporation at an ab initio density functional theory level of not just the thermodynamics of vacancy formation, but also relevant diffusion barriers for a vast number of possible local environments (in the core and in the shell, at the interface, and in the intermixed phase). This facilitates a predictive treatment and comprehensive understanding of intermixing on the relevant time scale (e.g., 10 1 –10 3 s). In contrast, previous modeling at the atomistic level utilized only unrealistic generic prescriptions of barriers or employed simplified continuum treatments. For Au@Ag octahedral core–cubic shell NCs, our modeling not only captures the experimentally observed rate or time scale for intermixing of ~100 s at 450 °C for 60 nm NCs, but also elucidates the underlying rate controlling processes and the effective intermixing barrier.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

Engineering Assembly Kinetics and Line Roughness in Solvent Vapor-Annealed Block Copolymer/Homopolymer Blends

Block copolymer (BCP) directed self-assembly (DSA) is a promising route to enhance lithography resolution by multiplying nanopattern density and reducing feature roughness. Eliminating kinetically trapped self-assembly defects requires fast self-assembly. However, acceleration strategies like solvent vapor annealing or homopolymer blending broaden domain interfaces, implying a trade-off in increased feature roughness. In this work, we experimentally investigate this apparent dilemma between self-assembly kinetics and line roughness for solvent vapor-annealed thin films of a lamellar poly(styrene-block-2-vinylpyridine) (PS-b-P2VP) BCP blended with PS and P2VP homopolymers. Binary blends with PS or P2VP homopolymers and ternary blends incorporating both in equal weight fractions were solvent vapor annealed using acetone, a near-neutral solvent for PS and P2VP, followed by P2VP-selective vapor-phase infiltration with alumina (AlOx) and polymer etching. Binary blends with P2VP exhibit a modest kinetic enhancement but also higher line-edge and -width roughness due to the increased frequency of P2VP protrusions and bridge defects in the alumina line patterns. In contrast, binary blends with PS self-assemble noticeably faster, while domain asymmetry from the added PS homopolymer reduces roughness by curbing the number of alumina protrusions and bridge defects. Ternary blends maintain DSA line patterns across a wider composition window and, at higher homopolymer loadings, reduce roughness at length scales near the lamellar period, consistent with a reduced impact of intradomain compositional fluctuations. These findings provide important insights for codesigning blend compositions and process flows to achieve high-resolution, defect free patterns with minimal roughness through BCP DSA.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Metal–Oxide Interface Sites Created Using Atomic Layer Deposition and Tested for CO Oxidation

The performance of catalysts made out of Pt supported on TiO 2 thin films grown on SBA-15 (a silica mesoporous material) by atomic layer deposition (ALD) was characterized systematically by combining in situ infrared absorption spectroscopy (IR) with other techniques including electron microscopy and adsorption−desorption isothermal measurements. The titania films in the resulting high-surface-area catalysts were evenly distributed throughout the inner surface of the SBA-15 mesopores, and their thickness could be controlled at a submonolayer level, with 3 to 4 TiO 2 ALD cycles needed for the complete coverage of the silica sites. The titania films could be deposited either before or after adding the metal (Pt), which was dispersed in the form of small nanoparticles (NPs) approximately 4−6 nm in diameter, in order to exert some control on the density and nature of the Pt/TiO 2 interface sites. One important lesson deriving from this work is that such an order of deposition leads to significantly different catalysts in spite of the fact that most of their structural properties are similar. If the Pt is deposited on the titania films, the resulting metal NPs are slightly smaller than those grown on silica and display CO adsorption sites with lower surface Pt coordination numbers. On the other hand, when TiO 2 is deposited on the Pt/SBA-15 starting material, some titania grows on the metal and partially blocks its surface while also creating new interface sites where CO binds more weakly and displays lower C−O stretching frequencies. In terms of catalytic performance, the results from in situ IR CO site titration and kinetic measurements combined suggest a mechanism where CO first adsorbs on Pt atop sites and then migrates to Pt/TiO 2 interface sites, where oxidation takes place. Both types of sites appear to be similar in all the catalysts tested, but catalytic performance could be optimized by tuning their surface densities. Maximum catalytic activity was obtained when the TiO2 films were deposited first and with TiO 2 coverages of at least half a monolayer, that is, after at least 2 ALD cycles.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Enhancing the Reaction Kinetics and Stability of Co‐Free Li‐Rich Cathode Materials via a Multifunctional Strategy

Co-free Li-rich layered oxides (CFLLOs) with anionic redox activity are among the most promising cathode materials for high-energy-density and low-cost lithium-ion batteries (LIBs). However, irreversible oxygen release often causes severe structural deterioration, electrolyte decomposition, and the formation of unstable cathode-electrolyte interface (CEI) film with high impedance. Additionally, the elimination of cobalt elements further deteriorates the reaction kinetics, leading to reduced capacity and poor rate performance. Here, in this study, a multifunctional strategy is proposed, incorporating Li 2 MnO 3 phase content regulation, micro-nano structure design, and heteroatom substitution. The increased content of Li 2 MnO 3 phase enhances the capacity through oxygen redox. The smaller nanoscale primary particles induce greater tensile strain and introduce more grain boundaries, thereby improving the reaction kinetics and reactivity, while the larger micron-sized secondary particles help to reduce interfacial side reactions. Furthermore, Na⁺ doping modulates the local coordination environment of oxygen, stabilizing both the anion framework and the crystal structure. As a result, the designed cathode exhibits enhanced rate performance, delivering a capacity of 158 mAh g⁻¹ at 5.0 C and improved cyclic stability, with a high capacity retention of 99% after 400 cycles at 1.0 C. This multifunctional strategy holds great promise for advancing the practical application of CFLLOs in next-generation LIBs.

Co-free Li-rich layered oxide↗

DOE OSTI2025

Investigation of Corrosion Behavior of Carbide Phase Strengthened Nickel‐Based Alloy in Molten Fluoride Salt

The corrosion behavior of a high creep strength carbide-phase strengthened Ni-based alloy in molten FLiNaK (LiF-NaF-KF: 46.5-11.5-42 mol%) salt in the temperature range of 700°C–750°C has been investigated as a part of the development of structural alloys for fluoride salt-based molten salt reactors (MSRs). The alloy composition was designed based on Hastelloy-N, but with the goal of improving creep strength. Cr depletion depth, a measure of corrosion, was observed to be single micrometers after several hundred hours of corrosion testing. Sequential corrosion testing involving testing of pre-tested samples in fresh salt coupled with SEM-EDS, scanning transmission electron microscopy (STEM) examinations and thermodynamic and kinetic modeling, suggest that the corrosion rate at the alloy-salt interface is governed by diffusion of elements from the alloy bulk to the surface. The carbide phases in the corrosion-tested sample microstructure were identified to be largely M 6 C-type Mo-rich carbides and MC-type mixed carbides. Atom probe tomography (APT) showed some partitioning of Cr and Ti to the carbide phase and showed the carbide phases to be stable at the salt-facing surface.

Gordon, Ryan [Univ. of Wisconsin, Madison, WI (Uni↗

DOE OSTI2024

Generative Physics-Informed Neural Network Solving Multi-Scale and Multi-Phase Plasma Chemical Flow Field

Low-temperature plasmas (LTPs) are non-equilibrium systems with near-room-temperature gas and highly energetic electrons. This makes them ideal for delicate applications in biomedicine and semiconductor manufacturing, enabling processes like wound healing, sterilization, etching, and plasma-enhanced chemical vapor deposition without thermal damage. However, LTPs involve complex chemistries, with hundreds of species and thousands of reactions, complicating their diagnosis, prediction, and control. Conventional diagnostics, such as Fourier-transform infrared spectroscopy (FTIR), laser-induced fluorescence (LIF), and optical emission spectroscopy (OES), offer limited species detection, while mass spectrometry (MS) struggles with low-sensitivity species. Additionally, LTP simulations face multi-scale challenges, as macroscopic fluid dynamics and microscopic particle collisions operate on vastly different timescales. To address these issues, we developed an artificial intelligence (AI) based diagnostic system: a generative physics-informed neural network (PINN-Gen) that can predict spatially resolved species concentrations and temperatures in LTPs by integrating experimental data from planar LIF with microscopic plasma chemical kinetics and macroscopic fluid mechanics, including plasma-liquid interactions at the interface between two phases. PINN-Gen solves no equations but checks the errors of physical laws by substituting the output from neural network, and the comparison with the experimental results. Thus, it naturally avoids the multi-scale difficulty of numerical simulations and predicts the results of conventionally unsolvable multi-scale and multi-phase problems. The real-time prediction will be robust due to the physical information used in the training of such a neural network, and only very limited input of condition required due to its generative feature.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2026

Roadmap on advanced and / real-time characterisation of solid state materials and devices for energy applications

A strong societal and political drive is motivating the development and optimization of novel energy conversion and storage systems for decarbonization. The successful implementation of solid state devices such as fuel cells and secondary batteries depends, however, on achieving ambitious targets in terms of performance, reliability and cost competitiveness. Research and technology are addressing these needs through a holistic approach including exploration of new materials and nanoarchitectures, as well as system engineering. These significant efforts require the support of appropriate characterization tools capable of assessing nanometer-scale phenomena such as concentration profiles of ionic and electronic charges, local chemical compositions and their evolution over time across interfaces. This roadmap provides an overview of selected advanced characterization techniques for energy materials and devices. Specific focus is put on in situ/operando methods for probing electrochemical phenomena in real-time under realistic working conditions. Experts in the field provide an extensive review of the current state of the art in 2025 and the current and future challenges for the characterization of local chemistry and kinetics in the bulk of the material, in nanoarchitectures (e.g. thin films) and at the interfaces (e.g. grain boundaries, phase contacts, solid/liquid and solid/gas interfaces) . The aim is to provide a detailed guide to the techniques, describing opportunities and bottlenecks for their practical deployment and examples of successful applications.

25 ENERGY STORAGE↗

DOE OSTI2025

Polymorphs of the n–Type Polymer P(NDI2OD–T2): A Comprehensive Description of the Impact of Processing on Crystalline Morphology and Charge Transport

A systematic study of the polymorphs emerging in P(NDI2OD-T2) (also commercially known as N2200), a prototypical organic semiconducting n-type polymer, is presented. Using a tightly integrated experimental and computational approach, detailed atomistic-level descriptions are provided investigating the three known P(NDI2OD-T2) polymorphs observed at room temperature as a function of thin-film processing. Importantly, over the course of the work, a missing link is uncovered, a fourth polymorph referred to here as Form I-β; this new form is a morphological intermediary observed upon thermal annealing, which evolves from Form I but tends to disappear upon full polymer chain melting. The computationally derived polymorph structures show excellent agreement with experimental X-ray scattering characterization. The relative stabilities of each polymorph are calculated in terms of both the bulk material and the polymorph-air interface. An energy landscape is then constructed to qualitatively compare the thermodynamic versus kinetic origins of each polymorph, and the factors driving (supra)assembly and associated transformations among polymorphs using an approach generalizable to other organic semiconducting polymers. Lastly, the relationships among preferential polymorphic crystallinity, relative chain orientations, and directional charge transport properties in P(NDI2OD-T2) are explored. Altogether, this work provides unprecedented insights into complex structure-processing-transport relationships in a representative semiconducting organic polymer.

36 MATERIALS SCIENCE↗