Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “benchmark correlations”

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

Accelerating resonant spectroscopy simulations using multishifted biconjugate gradient

Resonant spectroscopies, which involve intermediate states with finite lifetimes, provide important insights into collective excitations in quantum materials that are otherwise inaccessible. However, theoretical understanding in this area is often limited by the numerical challenges of solving Kramers-Heisenberg-type response functions for large-scale systems. To address this, we introduce a multishifted biconjugate gradient algorithm that exploits the shared structure of Krylov subspaces across spectra with varying incident energies, effectively reducing the computational complexity to that of linear spectroscopies. Both mathematical proofs and numerical benchmarks confirm that this algorithm substantially accelerates spectral simulations, achieving constant complexity independent of the number of incident energies, while ensuring accuracy and stability. This development provides a scalable, versatile framework for simulating advanced spectroscopies in quantum materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Benchmarking Coordination Number Prediction Algorithms on Inorganic Crystal Structures

Coordination numbers and geometries form a theoretical framework for understanding and predicting materials properties. Algorithms to determine coordination numbers automatically are increasingly used for machine learning and automatic structural analysis. In this work, we introduce MaterialsCoord, a benchmark suite containing 56 experimentally-derived crystal structures (spanning elements, binaries, and ternary compounds) and their corresponding coordination environments as described in the research literature. We also describe CrystalNN, a novel algorithm for determining near neighbors. We compare CrystalNN against 7 existing near-neighbor algorithms on the MaterialsCoord benchmark, finding CrystalNN to perform similarly to several well-established algorithms. For each algorithm, we also assess computational demand and sensitivity towards small perturbations that mimic thermal motion. Finally, we investigate the similarity between bonding algorithms when applied to the Materials Project database. Finally, we expect that this work will aid the development of coordination prediction algorithms as well as improve structural descriptors for machine learning and other applications.

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 OSTI2026

A comparison of wrought and powder metallurgical FeCrAl claddings under simulated LWR accident transients

Iron-chromium-aluminum (FeCrAl) alloys are potential accident tolerant fuel (ATF) cladding candidates for light-water reactors but are difficult to fabricate as thin-walled tubes via conventional cast-and-wrought routes. Powder metallurgy (PM) offers a manufacturing alternative with improved compositional control, but its transient accident performance has not been directly benchmarked against wrought variants. This study evaluates the burst behavior of commercially developed PM-processed FeCrAl alloys, PM-C26M (Fe-12Cr-6Al-2Mo) and the high precipitate density FA-SMT (Fe-22Cr-5Al-3Mo), under simulated light-water reactor accident transient conditions. Burst testing was conducted using the Severe Accident Test Station with heating rates of 5 °C/s and 50 °C/s and internal pressures ranging from 25 MPa to 100 MPa. PM-C26M reproduced wrought C26M burst behavior within 7–37 °C across the stress range, indicating that PM processing does not compromise transient strength. FA-SMT exhibited markedly higher burst temperatures and reduced heating-rate sensitivity, consistent with its engineered precipitate strengthening. FA-SMT rupture exhibited axial "unzipping" rather than the lateral tearing characteristic of PM- and wrought C26M. Post-test EBSD and fractography indicate that this behavior is strongly correlated with strain-gated intergranular void nucleation associated with the dense precipitate architecture of FA-SMT, a response absent in the comparatively clean PM-C26M matrix and consistent with rupture morphologies reported for oxide-dispersion strengthened (ODS) FeCrAl of similar base-matrix chemistry to PM-C26M. These findings highlight the potential of powder metallurgy as a viable fabrication route for ATF claddings from an accident performance standpoint.

Bell, Sam [ORNL] (ORCID:0000000251905657)↗
DOE OSTI2026

Definitive Assessment of the Accuracy, Variationality, and Convergence of Relativistic Coupled Cluster and Density Matrix Renormalization Group in 100-Orbital Space

Accuracy, variationality, and convergence underpin the reliability of modern electronic structure methods, yet definitive benchmarks in the relativistic regime remain elusive due to the absence of numerically exact full configuration interaction (CI) references. Recent algorithmic advances in the CI framework, enabled by the small-tensor-product (STP) decomposition approach, have dramatically extended the tractable size of the configuration space, making numerically exact CI calculations feasible in large active spaces previously beyond reach. In this paper, we employ the recently developed STP-CI framework to perform large-scale numerically exact CI calculations and directly benchmark relativistic coupled cluster and density matrix renormalization group methods. Definitive benchmarking of approximate relativistic electronic structure methods is ensured through the application of the gap theorem, which provides rigorous error bounds on the CI reference and establishes a controlled standard for assessing accuracy, variationality, and convergence.

Chemical calculations↗
DOE OSTI2020

Emergence of quasiparticles in a doped Mott insulator

How a Mott insulator develops into a weakly coupled metal upon doping is a central question to understanding various emergent correlated phenomena. To analyze this evolution and its connection to the high- T c cuprates, we study the single-particle spectrum for the doped Hubbard model using cluster perturbation theory on superclusters. Starting from extremely low doping, we identify a heavily renormalized quasiparticle dispersion that immediately develops across the Fermi level, and a weakening polaronic side band at higher binding energy. The quasiparticle spectral weight roughly grows at twice the rate of doping in the low doping regime, but this rate is halved at optimal doping. In the heavily doped regime, we find both strong electron-hole asymmetry and a persistent presence of Mott spectral features. Finally, we discuss the applicability of the single-band Hubbard model to describe the evolution of nodal spectra measured by angle-resolved photoemission spectroscopy (ARPES) on the single-layer cuprate La 2– x Sr x CuO 4 (0 ≤ x ≤ 0.15). This work benchmarks the predictive power of the Hubbard model for electronic properties of high- T c cuprates.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRSDate not supplied

Implementation of Two Local Correlation-Based Transition Models in OVERFLOW 2.3e

Accurate prediction of the laminar-to-turbulent boundary layer transition is highlighted as a vital component of NASA's CFD Vision 2030 Study. With an increased emphasis on greener air transports and sustainable aviation, transition modeling is anticipated to have an added significance, particularly in the applications related to laminar flow technology. However, unmanned aerial vehicles, crewed reentry vehicles, and ground-to-flight extrapolation all benefit from transition models. It is useful to incorporate a variety of models in CFD solvers because there isn't a single transition model that is ideal for the complete spectrum of applications. While the Langtry-Menter γ–Reθt model has been widely used for CFD predictions of flows with laminar, transitional, and turbulent boundary layers, it does not meet the criteria for Galilean invariance, a desirable attribute for rotorcraft applications. To help overcome that limitation, we have implemented Menter's SST-based γ transition model within NASA’s OVERFLOW Overset CFD code (version 2.3e). In addition, an SA-based implementation of the Langtry-Menter γ–Reθt transition model has also been implemented in view of the increased robustness of the SA-turbulence model for external aerodynamic applications. An initial assessment of the implemented models has been carried out using benchmark 2D test cases involving flat plates and the NLF-0416 airfoil, focusing on bypass, Tollmien-Schlichting, and laminar separation bubble-induced transition scenarios. In comparison to the Langtry-Menter γ–Reθt model, both models were found to produce increased computation efficiency, principally as a result of the decrease from four to three auxiliary transport equations. Results from the initial analysis will be presented, along with additional details pertaining to the deployment of these models.

CFD modeling↗
DOE OSTI2024

A general Bayesian algorithm for the autonomous alignment of beamlines

Autonomous methods to align beamlines can decrease the amount of time spent on diagnostics, and also uncover better global optima leading to better beam quality. The alignment of these beamlines is a high-dimensional expensive-to-sample optimization problem involving the simultaneous treatment of many optical elements with correlated and nonlinear dynamics. Bayesian optimization is a strategy of efficient global optimization that has proved successful in similar regimes in a wide variety of beamline alignment applications, though it has typically been implemented for particular beamlines and optimization tasks. In this paper, we present a basic formulation of Bayesian inference and Gaussian process models as they relate to multi-objective Bayesian optimization, as well as the practical challenges presented by beamline alignment. We show that the same general implementation of Bayesian optimization with special consideration for beamline alignment can quickly learn the dynamics of particular beamlines in an online fashion through hyperparameter fitting with no prior information. We present the implementation of a concise software framework for beamline alignment and test it on four different optimization problems for experiments on X-ray beamlines at the National Synchrotron Light Source II and the Advanced Light Source, and an electron beam at the Accelerator Test Facility, along with benchmarking on a simulated digital twin. We discuss new applications of the framework, and the potential for a unified approach to beamline alignment at synchrotron facilities.

47 OTHER INSTRUMENTATION↗
DOE OSTI2020

A coupled cluster framework for electrons and phonons

Here, we describe a coupled cluster framework for coupled systems of electrons and harmonic phonons. Neutral and charged excitations are accessed via the equation-of-motion version of the theory. Benchmarks on the Hubbard–Holstein model allow us to assess the strengths and weaknesses of different coupled cluster approximations, which generally perform well for weak to moderate coupling. Finally, we report progress toward an implementation for ab initio calculations on solids and present some preliminary results on finite-size models of diamond with a linear electron–phonon coupling. We also report the implementation of electron–phonon coupling matrix elements from crystalline Gaussian type orbitals within the PySCF program package.

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

Improving deep learning-based protein distance prediction in CASP14

Abstract Motivation Accurate prediction of residue–residue distances is important for protein structure prediction. We developed several protein distance predictors based on a deep learning distance prediction method and blindly tested them in the 14th Critical Assessment of Protein Structure Prediction (CASP14). The prediction method uses deep residual neural networks with the channel-wise attention mechanism to classify the distance between every two residues into multiple distance intervals. The input features for the deep learning method include co-evolutionary features as well as other sequence-based features derived from multiple sequence alignments (MSAs). Three alignment methods are used with multiple protein sequence/profile databases to generate MSAs for input feature generation. Based on different configurations and training strategies of the deep learning method, five MULTICOM distance predictors were created to participate in the CASP14 experiment. Results Benchmarked on 37 hard CASP14 domains, the best performing MULTICOM predictor is ranked 5th out of 30 automated CASP14 distance prediction servers in terms of precision of top L/5 long-range contact predictions [i.e. classifying distances between two residues into two categories: in contact (<8 Angstrom) and not in contact otherwise] and performs better than the best CASP13 distance prediction method. The best performing MULTICOM predictor is also ranked 6th among automated server predictors in classifying inter-residue distances into 10 distance intervals defined by CASP14 according to the precision of distance classification. The results show that the quality and depth of MSAs depend on alignment methods and sequence databases and have a significant impact on the accuracy of distance prediction. Using larger training datasets and multiple complementary features improves prediction accuracy. However, the number of effective sequences in MSAs is only a weak indicator of the quality of MSAs and the accuracy of predicted distance maps. In contrast, there is a strong correlation between the accuracy of contact/distance predictions and the average probability of the predicted contacts, which can therefore be more effectively used to estimate the confidence of distance predictions and select predicted distance maps. Availability and implementation The software package, source code and data of DeepDist2 are freely available at https://github.com/multicom-toolbox/deepdist and https://zenodo.org/record/4712084#.YIIM13VKhQM. Supplementary information Supplementary data are available at Bioinformatics online.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRSDate not supplied

Flow field Reconstruction for Inhomogeneous Turbulence using Data and Physics Driven Models

A methodology combining Large Eddy Simulation (LES) trained data and a physics driven wave packet model to obtain a reduced order reconstruction for broadband, three-dimensional, temporally stationary but spatially inhomogeneous, incompressible turbulence. Wake turbulence generated by an axisymmetric dragging disk with a turbulent co-flow serves as the benchmark test case. We begin by studying the proper-orthogonal decomposition of the turbulent fluctuations taken from a high-resolution LES to first identify whether the fields demonstrate a low-rank character. It is argued that the presence of the turbulent co-flow results in a largely broadband character lacking any tonal properties. This is especially true for Strouhal numbers greater than 1 and only a small fraction of energy is contained in the leading order Kelvin-Helmholtz modes. As such reconstructions and reduced order modeling purely relying on data from LES does not appear to be a lucrative solution - contrary to problems with strongly tonal character. To supplement the missing energy from a low order truncated mode expansion, we utilize a physics based super-resolution (enrichment) algorithm that relies on spatio-temporally localized Gabor wave packets whose time evolution is described using a set of ordinary differential equations. The reconstructed flow has single- and two-point correlations that are consistent with the reference high resolution simulation data.

SLS↗
DOE OSTI2023

Comparison of state-interaction and spinor-representation calculations of spin-orbit coupling within exact two-component coupled-cluster theories

Here a benchmark study of state-interaction and spinor-representation calculations of spin-orbit coupling using the exact two-component Hamiltonians with atomic mean-field integrals (the X2CAMF schemes) at the equation-of-motion coupled-cluster singles and doubles level is reported. We adopt a version of the X2CAMF scheme with spin-orbit integrals correct to first order for the state-interaction calculations and a version correct to infinite order for the spinor-representation calculations. The differences between the state-interaction calculations with minimum active spaces to account for (quasi-)degeneracy and the spinor-representation calculations thus correspond to higher-order spin-orbit contributions. This state-interaction approach with scalar-relativistic effects accurately included in the reference functions and in the spin-orbit integrals is shown to exhibit robust performance for elements across the periodic table. On the other hand, the more rigorous spinor representation shows more rapid convergence with respect to the number of correlated electrons and is the preferred choice for accurate calculations for heavy elements.

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

Experimental data for damage mechanics simulation challenge

While there are many computational approaches for simulating damage in rock and other materials, few have been ground truth tested with either known experimental data or with blind data sets. Here, in this work, we present a bench-mark laboratory data set for a damage mechanics challenge to compare computational approaches on damage evolution in brittle-ductile materials. The samples were fabricated through additive manufacturing to produce repeatable specimens designed to fail in controlled ways. The failure was induced in the samples using a 3-point bending test to produce different Modes such as Mode I and mixed Modes including I-II, I-III and I-II-III Modes to generate a calibration data set and a blind challenge data set. Data collected included spatial and temporal measurements from traditional digital load–displacement sensors, 2D digital image correlation measurement to map surface deformations, 3D X-ray microscopy to ground-truth the crack-failure geometry, and laser profilometry to capture surface roughness. The data sets are available, on a data repository, to the community to advance computational models to improve our ability to predict damage in brittle-ductile materials.

3-point bending↗
DOE OSTI2026

Monitoring of Liquid Metal Reactor Heater Zones with Recurrent Neural Network Learning of Temperature Time Series

Advanced high-temperature fluid reactors (ARs), such as sodium fast reactors (SFRs) and molten salt cooled reactors (MSCRs) utilize high-temperature fluids at ambient pressure. To melt the fluid during reactor startup and prevent fluid freezing during cooldown, the thermal–hydraulic systems of such ARs include heater zones consisting of specific heaters with controllers, temperature sensors, and thermal insulation. The failure of heater zones due to insulation material degradation or improper installation, resulting in parasitic heat losses, can lead to fluid freezing. The detection of faults using a heat-transfer model is difficult because of a lack of knowledge of the experimental details. Data-driven machine learning of heater zone temperature time series offers a viable alternative. In this study, we benchmarked the performance of recurrent neural networks (RNNs) in an analysis of heat-up transient temperature time series of heater zones installed on a liquid sodium vessel. The RNN models include long short-term memory (LSTM) and gated recurrent unit (GRU) networks, as well as their bi-directional variants, BiLSTM and BiGRU. Anomalous temperature points were designated using a percentile-based threshold applied to residual fluctuations in the detrended temperature time series. Additionally, the impact of the exponentially weighted moving average (EWMA) method on detection accuracy was examined. The RNN models’ performance was assessed using precision, recall, and F 1 score metrics. Results demonstrated that RNN models effectively detect anomalies in temperature time series with the best models for each heater zone achieving F 1 scores of over 93%. To explain the variations in RNN model performance across different heater zones, we used Kullback–Leibler (KL) divergence to quantify the relative entropy between training and testing data, and the Detrended Fluctuation Analysis (DFA) to assess long-range temporal correlations. For datasets with strong long-range correlations and minimal relative entropy between training and testing data, GRU is the best-performing model. When the data exhibits weaker long-term correlations and a significant relative entropy between training and testing distributions, BiGRU shows the best performance. For the data sets with intermediate values of both KL divergence and DFA, the best performance is obtained with LSTM and BiLSTM, respectively.

gated recurrent unit↗
DOE OSTI2021

Identification and quantification of distinct active sites in Hf-Beta zeolites for transfer hydrogenation catalysis

Despite the significant progress made in characterizing different framework heteroatom sites that exist in Lewis acidic zeolites with probe molecule adsorption and spectroscopy, methods to reliably quantify site counts remain indefinite and have been primarily limited to Sn and Ti Lewis acid sites. Here, methods to quantify framework Lewis acidic Hf 4+ sites in zeolite Beta (Hf-Beta) with two Lewis base titrants (pyridine, deuterated acetonitrile) were developed using infrared (IR) spectroscopy. Lewis acid site counts for Hf-Beta zeolites were validated by measuring integrated molar extinction coefficients (IMECs; ε, cm μmol –1 ) on Sn-Beta zeolites using identical Lewis base titrants to benchmark site counts with established literature procedures to quantify Lewis acid sites in Beta zeolites from IR spectra. Highlighting the importance of benchmarking active site counts against well-established experimental protocols, IMECs of CD 3 CN bound to open (ε(Sn; 2316 cm –1 ): 1.80 ± 0.25) and closed (ε(Sn; 2308 cm –1 ): 3.76 ± 0.33) Sn sites were ~1.8x larger than those previously reported while total Lewis acid site counts agreed with those measured by pyridine (ε(Sn; 1451 cm –1 ): 1.58 ± 0.16) on six different Sn-Beta zeolites. IMECs measured for IR peaks reflecting pyridine bound to Lewis acidic Hf sites (ε(Hf; 1448 cm –1 ): 1.54 ± 0.21) and CD 3 CN bound to open (ε(Hf; 2313 cm –1 ): 2.40 ± 0.22) and closed (ε(Hf; 2307 cm –1 ): 3.55 ± 0.41) Hf sites, gave similar counts for the total number of Lewis acidic sites across six Hf-Beta zeolites (Si/Hf = 100–413). Consistent with previous reports with Sn-Beta catalysts where open Sn sites are responsible for catalytic turnover, apparent first and zero-order MPVO rate constants (0.01–1 M cyclohexanone in 2-butanol; per total Hf, 373 K) correlated with the total number of open Hf sites, per total Hf, but not with the total number of closed Hf sites or total Lewis acid site counts. Measured initial MPVO rates (0.1 M cyclohexanone in 2-butanol, per open Hf, 373 K) were ~25x higher on hydrophobic Hf-Beta-F than on hydrophilic Hf-Beta-OH zeolites. Overall, the apparent first-order MPVO rate constants (2-butanol solvent, per open Hf, 373 K) were ~6x higher on Hf-Beta-F than on Hf-Beta-OH zeolites. The characterization methods reported here enable normalization of MPVO turnover rates on Sn- and Hf-Beta zeolites by their number of open sites. Finally, this enables performing quantitative rate comparisons across Lewis acid zeolites of varying active site identity, solvation, and pore topology used in liquid-phase catalysis.

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

Rapidly convergent quantum Monte Carlo using a Chebyshev projector

The multireference coupled-cluster Monte Carlo (MR-CCMC) algorithm is a determinant-based quantum Monte Carlo (QMC) algorithm that is conceptually similar to Full Configuration Interaction QMC (FCIQMC). It has been shown to offer a balanced treatment of both static and dynamic correlation while retaining polynomial scaling, although application to large systems with significant strong correlation remained impractical. In this paper, we document recent algorithmic advances that enable rapid convergence and a more black-box approach to the multireference problem. These include a logarithmically scaling metric-tree-based excitation acceptance algorithm to search for determinants connected to the reference space at the desired excitation level and a symmetry-screening procedure for the reference space. We show that, for moderately sized reference spaces, the new search algorithm brings about an approximately 8-fold acceleration of one MR-CCMC iteration, while the symmetry screening procedure reduces the number of active reference space determinants with essentially no loss of accuracy. We also introduce a stochastic implementation of an approximate wall projector, which is the infinite imaginary time limit of the exponential projector, using a truncated expansion of the wall function in Chebyshev polynomials. Notably, this wall-Chebyshev projector can be used to accelerate any projector-based QMC algorithm. We show that it requires significantly fewer applications of the Hamiltonian to achieve the same statistical convergence. We benchmark these acceleration methods on the beryllium and carbon dimers, using initiator FCIQMC and MR-CCMC with basis sets up to cc-pVQZ quality.

Zhao, Zijun↗
DOE OSTI2023

Automatic Crack Segmentation and Feature Extraction in Electroluminescence Images of Solar Modules

The effect of cracks in solar cells on the long-term degradation of photovoltaic (PV) modules remains to be determined. To investigate this effect in future studies, it is necessary to quantitatively describe the crack features (e.g., length) and correlate them with module power loss. Electroluminescence (EL) imaging is a common technique for identifying cracks. However, it is currently challenging and time-consuming to identify cracks in a large number of EL images and quantify complex crack features by human inspection. This article introduces a fast semantic segmentation method (~0.18 s/cell) to automatically segment cracks from EL images and algorithms to extract crack features. Here we fine-tuned a UNet neural network model using pretrained VGG16 as the encoder and obtained an average F1 score of 0.875 and an intersection over union score of 0.782 on the testing set. With cracks and busbars segmented, we developed algorithms for extracting crack features, including the crack-isolated area, the brightness inside the isolated area, and the crack length. We also developed an automatic preprocessing tool for cropping individual cell images from EL images of PV modules (~0.72 s/module). Our codes are published as open-source an software, and our annotated dataset composed of various types of cells is published as a benchmark for crack segmentation in EL images.

14 SOLAR ENERGY↗
DOE OSTI2025

Polaron Delocalization and Transport in Doped Graphene Nanoribbon Thin Films

Graphene nanoribbons (GNRs) are quantum-confined π-conjugated monolayer semiconductors with attractive properties for optoelectronic applications. However, the ground- and excited-state properties of charge carriers in GNRs are still poorly understood, particularly with regards to the coupling between charges and the GNR lattice and the degree to which this coupling impacts local and macroscopic charge transport. To address this issue, we systematically correlate carrier density-dependent charge transport with spectroscopic modulations in chemically doped thin films of armchair graphene nanoribbons (9-aGNRs). This study combines Fourier transform infrared (FTIR) and ultraviolet−visible−near-infrared (UV−vis−NIR) spectroscopy with both local and macroscopic conductivity measurements to arrive at a full and self-consistent picture of transport in doped GNR thin films. Using three different molecular p-type dopants (i.e., oxidants), we demonstrate that hole polarons are the dominant quasi-particle determining charge transport in GNRs and that the degree of polaron delocalization depends sensitively on the dopant and the hole density. For all three dopants, the local conductivity probed by microwave spectroscopy substantially exceeds the long-range conductivity obtained by four-point probe measurements. Interestingly, the dopant size substantially influences charge transport at high hole densities. We ascribe this effect to different propensities for forming bipolarons with lower mobilities than polarons. Comparison of GNR transport and spectral properties to other prototypical π-conjugated semiconductors (e.g., semiconducting polymers or carbon nanotubes) benchmark the charge transport properties of GNR thin films for optoelectronic devices and applications.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗