Search NASASearch

SEARCH · Search NASA

Results for “electrode geometry”

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 19 records

Two-phase flow numerical analysis of electrode geometry for alkaline water electrolyzers

Hydrogen is a promising component of a future energy-secure and efficient economy, but its competitiveness depends on reducing production costs. One strategy is to operate alkaline water electrolyzers at higher current densities to increase output. However, this intensifies performance losses due to gas bubble accumulation, which blocks transport pathways and deactivates electrochemically active surfaces. Enhancing bubble evacuation through electrode design is therefore essential. Previous studies have explored various approaches — such as modifying surface morphology, applying sonication or pressure modulation, and introducing surfactants — but these efforts have addressed a limited range of conditions due to the complexity of two-phase flow and electrode geometries. Experiments have also largely been focused on either cell level improvements, which lack the information necessary to isolate each contributing factor, or on modified geometries that are not relevant to practical cell operation. From a modeling perspective, conventional Eulerian multiphase models do not track the complex gas–liquid interfacial dynamics and often neglect surface tension and contact angle effects, reducing their predictive accuracy. To provide insights on the effects of different electrode geometries on the performance of alklaine water electrolyzers this work employs an immersed boundary volume-of-fluid method to simulate bubble behavior in 3D porous electrodes. Multiple base electrode geometries, typically used in practice, with varying porosity are evaluated under a constant surface gas generation rate. Simulation data is analyzed to quantify electrode gas coverage, bubble size dynamics and other relevant metrics. Results show that porosity strongly influences bubble accumulation on electrode surfaces, with higher porosity reducing gas coverage, and its not strictly dependent on the electrode geometry. However, the electrode’s base geometry significantly affects gas accumulation at the separator gap, independent of porosity. A foam electrode geometry resulted in the lowest gas coverage of all electrodes with a median volumetric gas coverage of 11%, but at the cost of a 70% reduction in active area compared with the largest surface area electrode, while gyroid electrodes showed the best trade-off between gas coverage, particularly at the separator surface, and electrochemically active area. In conclusion, the results highlight the need for holistic electrode design strategies.

Alkaline

Shape Matters: Understanding the Effect of Electrode Geometry on Cell Resistance and Chemo-Mechanical Stress

Rechargeable batteries that incorporate shaped three-dimensional electrodes have been shown to have increased power and energy densities when compared to a conventional geometry, i.e. a planar cathode and anode that sandwich an electrolyte. Electrodes can be shaped to enable a higher active material loading, while keeping ion transport distances small. However, the relationship between electrical and mechanical performance of shaped electrodes remains poorly understood. Many electrode designs have been explored, where the electrodes are individually shaped or intertwined, and advances in manufacturing and shape/topology optimization have made such designs a reality. Here, we explore sinusoidal half cells and interdigitated full cells. First, we use a simple electrostatics model to understand the cell resistance as a function of shape. Here we focus on low-temperature conditions, where the electrolyte conductivity decreases relative to that of the electrode; here, LiPF 6 EC:DMC electrolyte and MnO 2 electrode are considered. Next, we use a chemo-mechanics model to examine the stress that arises due to intercalation-driven volume expansion. We show that shaped electrodes provide a significant reduction in resistance in low-temperature conditions, however, they exhibit unfavorable stress concentrations. Overall, we find that the fully interdigitated electrodes may provide the best balance with respect to this resistance-stress trade-off.

25 ENERGY STORAGE

Proximity Ferroelectricity in Compositionally Graded Structures

Proximity ferroelectricity is a novel paradigm for inducing ferroelectricity in a non-ferroelectric polar material, such as AlN or ZnO that are typically unswitchable with an external field below their dielectric breakdown field. When placed in direct contact with a thin switchable ferroelectric layer (such as Al 1-x Sc x N or Zn 1-x Mg x O), they become a practically switchable ferroelectric. Using the thermodynamic Landau-Ginzburg-Devonshire theory, in this work, we perform the finite element modeling of the polarization switching in the compositionally graded AlN-Al 1-x Sc x N, ZnO-Zn 1-x Mg x O, and MgO-Zn 1-x Mg x O structures sandwiched in both a parallel-plate capacitor geometry as well as in a sharp probe-planar electrode geometry. We reveal that the compositionally graded structure allows the simultaneous switching of spontaneous polarization in the whole system by a coercive field significantly lower than the electric breakdown field of unswitchable polar materials. The physical mechanism is the depolarization electric field determined by the gradient of chemical composition “x”. The field lowers the steepness of the switching barrier in the otherwise unswitchable parts of the compositionally graded AlN-Al 1-x Sc x N and ZnO-Zn 1-x Mg x O structures. In the MgO-like regions of the compositionally graded MgO-Zn 1-x Mg x O structure, a shallow double-well free energy potential emerges. Proximity ferroelectric switching of the compositionally graded structures placed in the probe-electrode geometry occurs due to nanodomain formation under the tip. We predict that a gradient of chemical composition “x” significantly lowers effective coercive fields of the compositionally graded AlN-Al 1-x Sc x N and ZnO-Zn 1-x Mg x O structures compared to the coercive fields of the corresponding multilayers with a uniform chemical composition in each layer. A tip-induced switching further lowers the coercive field, enabling control of ferroelectric domains in otherwise unswitchable compositionally graded structures, which can provide nanoscale domain control for memory, actuation, sensing, and optical applications.

36 MATERIALS SCIENCE

A Free-Standing Boron-Doped Diamond Grid Electrode for Fundamental Spectroelectrochemistry

Spectroelectrochemistry (SEC) is a powerful technique that enables a variety of redox properties to be studied, including formal potential (E o ), thermodynamic values (ΔG, ΔH, ΔS), diffusion coefficient (D), electron transfer stoichiometry (n), and others. SEC requires an electrode which light can pass through while maintaining sufficient electrical conductivity. This has been traditionally composed of metal or metal oxide films atop transparent substrates like glass, quartz, or metallic mesh. Robust electrode materials like boron-doped diamond (BDD) could help expand the environments in which SEC can be performed, but most designs are limited to thin films (~100–200 nm) on transparent substrates less resilient than free-standing BDD. Here, this work presents a free-standing BDD grid electrode (G-BDD) for fundamental SEC measurements, using the well-characterized Fe(CN) 6 3–/4– redox couple as proof-of-concept. With a combination of cyclic voltammetry (CV), thin-layer SEC, and chronoabsorptometry, several of the redox properties mentioned above were calculated and compared. For E o' , n, and D, similar results were obtained when comparing the CV [E o' = +0.279 (±0.002) V vs Ag/AgCl; n = 0.97; D = 4.1 × 10 –6 cm 2 ·s –1 ] and SEC [E o' = +0.278 (±0.001) V vs Ag/AgCl; n = 0.91; D = 5.2 × 10 –6 cm 2 ·s –1 ] techniques. Both values align with what has been previously reported. To calculate D from the SEC data, modification of the classical equation used in chronoabsorptometry was required to accommodate the G-BDD electrode geometry. Overall, this work expands on the applicability of SEC techniques and BDD as a versatile electrode material.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH

Acoustophoretic Additive Manufacturing for Scalable 3D Battery Electrodes

This project focused on investigating two acoustic-based processing methods: a nozzle-based printhead and a chamber that map to two different battery electrode architectures: (1) a line-pattern electrode and (2) a grid-pattern electrode. These two parallel manufacturing and electrode geometry explorations were proposed for the project to understand the process space of acoustic-based manufacturing methods to fabricate patterned battery electrodes. This two-path exploration also allowed us to de-risk the overall project and not rely on a single process for creating patterned electrodes. This project consisted of six high-level tasks aimed at transitioning the concept of acoustic focusing for battery electrodes from a technology readiness level (TRL) of 1 to 3 by project conclusion. Overall, we believe we have developed a practical and high-impact processing method that is chemistry agnostic and suitable for large-area fabrication of both 3D LIBs and other functional material systems where structuring on the scale of tens of microns has the potential to break conventional bulk material trade-offs in performance. In the case of batteries for electric vehicles, structuring 3D LIBs with our acoustic process breaks traditional energy and power trade-offs observed with conventional flat battery packs.

25 ENERGY STORAGE

The Beta-decay Paul Trap Mk IV: Design and commissioning

Here, the Beta-decay Paul Trap is an open-geometry, linear trap used to measure the decays of 8 Li and 8 B to search for a tensor contribution to the weak interaction. In the latest 8 Li measurement of Burkey et al. (2022), β scattering was the dominant experimental systematic uncertainty. The Beta-decay Paul Trap Mk IV reduces the prevalence of β scattering by a factor of 4 through a redesigned electrode geometry and the use of glassy carbon and graphite as electrode materials. The trap has been constructed and successfully commissioned with 8 Li in a new data campaign that collected 2.6 million triple coincidence events, an increase in statistics by 30% with 4 times less β scattering compared to the previous 8 Li data set.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND

Electrical Resistivity Tomography data from 2016 to 2018 at the Lower Montane site in the East River Watershed, Colorado

This dataset contains time-lapse Electrical Resistivity Tomography (ERT) data along a transect located on the northeast-facing hillslope at the lower montane site (Pumphouse site) in the upper East River Watershed. The monitoring dataset covers the period from November 2, 2016, to August 6, 2018. In addition, the archive also contains a baseline dataset from October 9, 2016. The ERT transect consisted of 128 electrodes with an electrode spacing of 1.25 m. The acquisition system was located in the middle of the transect, about 50 m on one side, and included an MPT (Multi-Phase Technologies) ERT system, a mini computer, and batteries with solar panels. Acquisition occurred daily under normal circumstances. The first 16 electrodes (from the upper end of the transect) could not be used after the cable was damaged during the 2017–2018 winter. Also, due to multiple failures in the power system, the temporal resolution of the data is much lower in 2018 compared to 2016 and 2017. The data have been processed and used in Dafflon et al., 2023, and the baseline dataset was used in Falco et al., 2019 (see reference list). This archive contains the measurements (ER.zip containing csv files) for each of the 326 acquisition times and a filtered version where only electrodes 17 to 128 are included (ERT_sm.zip containing csv files). The archive also contains the baseline dataset and two acquisitions with full reciprocals (ERT_RB.zip containing csv files), as well as all the raw MPT files (ERT_raw_MTP.zip). The geometry (electrode position and elevation) is provided in Universal Transverse Mercator (UTM) 13N Geoid2012AB in the file named ERT_Location.csv. The archive contains 1 *.csv data files, four *.zip files, and three metadata *.csv files. This work was supported by the Watershed Function Science Focus Area at Lawrence Berkeley National Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-05CH11231.

54 ENVIRONMENTAL SCIENCES

Interplay of intercalation dynamics and lithium plating in monolithic and architectured graphite anodes during fast charging

Fast charging of high-capacity anodes is challenging due to lithium plating reactions, which lead to poor cycling performance and safety concerns. Thus, accurate predictions of plating onset and an understanding of this electrochemical process are crucial for robust battery design. However, the most commonly used models, based on porous electrode theory (e.g., the pseudo-2D model), are notoriously difficult to calibrate due to their complexity, limiting their predictive power. This work studies the process of lithium plating during fast charging of (small-particle) graphite half-cells by measuring local reaction progression and plating behavior using optical operando techniques. These experiments employ a realistic 1D graphite electrode geometry with commercially-relevant mass loading charged at fast charge rates. It is demonstrated that the local reaction progression and plating onset can not only be predicted accurately with a p2D numerical model, but that these processes follow a simple scaling law. Remarkably, the entire reaction histories of different electrodes charged at different rates (e.g., 160 μm thickness at 0.5C, 111 μm at 1C or 66 μm at 4C) were observed to have self-similar intercalation profiles. It is demonstrated that plating onset is in turn governed by the reaction profile which explains why both processes exhibit the same scaling behavior. Finally, operando measurements of local reaction dynamics are conducted for the first time in electrodes with channeled architectures, quantitatively determining how channels affect reaction uniformity and plating onset. Together, these results reveal underlying simplicity in the complex electrochemical environment of fast charging and lithium plating, improving understanding of this process. These fundamental insights are broadly applicable for design processes, modeling and experimental evaluation of lithium ion batteries.

25 ENERGY STORAGE

Scaling of energy delivered through an electrostatic discharge to a small series load

We study the energy delivered through a small-resistance series “victim” load during electrostatic discharge events in air. For gap lengths over 1 mm, the fraction of the stored energy delivered is mostly gap-length independent, with a slight decrease at larger gaps due to electrode geometry. The energy to the victim scales linearly with circuit capacitance and victim load resistance but is not strongly dependent on circuit inductance. This scaling leads to a simple approach to predicting the maximum energy that will be delivered to a series resistance for the case where the victim load resistance is lower than the spark resistance.

42 ENGINEERING

Rarefied xenon flow in orificed hollow cathodes

A parametric study is conducted to quantify the effect of the keeper electrode geometry on the xenon neutral flow quantities within orificed hollow cathodes, prior to ignition. The keeper impinges directly on the flow out of the cathode orifice and its geometry influences the product between the pressure in the orifice–keeper region and the cathode-to-keeper distance. A representative cathode is simulated using the Direct Simulation Monte Carlo method. The numerical model is first validated with computational results from the literature. A parametric study is then conducted. Parameters include the cathode pressure–diameter in the range of 1–5 Torr cm and the following geometric ratios (and ranges): cathode orifice-to-inner radii (0.1–0.7), keeper orifice-to-cathode orifice radii (1–5), and keeper distance-to-cathode-orifice diameter (0.5–10). It is found that, if both keeper and cathode have identical orifice radii, the flow remains subsonic in the orifice-to-keeper region. In most cases, however, the flow becomes underexpanded and supersonic, and the static pressure within the orifice-to-keeper region is, on average, 4% that of the upstream pressure value. The orifice–keeper region pressure increases with either a decrease in the keeper orifice diameter or an increase in the distance between cathode and keeper, in agreement with literature data. Both trends are explained through conservation laws. A statistical study of numerical results reveals that the ratio of ignition-to-nominal mass flow rates has a most probable value of 50, which suggests that heaterless cathode ignition at a minimum DC voltage may be achieved by increasing the input mass flow rate by a factor of 50.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY

Photocurrent Nanoscopy of Quantum Hall Bulk

Understanding nanoscale electronic and thermal transport of two-dimensional (2D) electron systems in the quantum Hall regime, particularly in the bulk insulating state, poses considerable challenges. One of the primary difficulties arises from the presence of chiral edge channels, whose transport behavior obscures the investigation of the insulating bulk. Using near-field optical and photocurrent nanoscopy, we probe real-space variations of the optical and thermal dynamics of graphene in the quantum Hall regime without relying on complex sample or electrode geometries. Near the charge-neutrality point, we detect strong optical and photothermal signals from resonant inter-Landau-level (LL) magnetoexciton excitations between the zeroth and ± 1 st LLs, which gradually weaken with increasing doping due to Pauli blocking. Interestingly, at higher doping levels and full-integer LL fillings, photothermal signals reappear across the entire sample over an approximately 10 − μ m scale, indicating unexpectedly long cooling lengths and nonlocal photothermal heating through the insulating bulk. This observation suggests that thermal conductivity persists for the localized states even as electronic transport is suppressed—a clear violation of the Wiedemann-Franz law. Our experiments provide novel insights into nanoscale thermal and electronic transport in incompressible 2D gases, highlighting the roles of magnetoexcitons and chiral edge states in the thermo-optoelectric dynamics of the Dirac quantum Hall state. Published by the American Physical Society 2025

36 MATERIALS SCIENCE

Tip-based proximity ferroelectric switching and piezoelectric response in wurtzite multilayers

Proximity ferroelectricity is a paradigm for inducing ferroelectricity when a nonferroelectric polar material (such as Al⁢ N), which is unswitchable with an external field below the dielectric breakdown field, becomes a practically switchable ferroelectric in direct contact with a thin switchable ferroelectric layer (such as Al 1−𝑥 ⁢Sc 𝑥 ⁢N). Here, we develop a Landau-Ginzburg-Devonshire approach to study the proximity effect of local piezoelectric response and polarization reversal in wurtzite ferroelectric multilayers under a sharp electrically biased tip. Using finite-element modeling, we analyze the probe-induced nucleation of nanodomains, the features of local polarization hysteresis loops and coercive fields in the Al 1−𝑥 ⁢Sc 𝑥 ⁢N/Al⁢ N bilayers and three-layers. Similar to the wurtzite multilayers sandwiched between two parallel electrodes, the regimes of “proximity switching” (when all layers collectively switch) and the regime of “proximity suppression” (when they collectively do not switch) are the only two possible regimes in the probe-electrode geometry. However, the parameters and asymmetry of the local piezoresponse and polarization hysteresis loops depend significantly on the sequence of the layers with respect to the probe. The physical mechanism of proximity ferroelectricity in the local probe geometry is a depolarizing electric field determined by the polarization of the layers and their relative thickness. The field, whose direction is opposite to the polarization vector in the layer(s) with the larger spontaneous polarization (such as Al⁢ N), renormalizes the double-well ferroelectric potential to lower the steepness of the switching barrier in the “otherwise unswitchable” polar layers. Tip-based control of domains in otherwise nonferroelectric layers using proximity ferroelectricity can provide nanoscale control of domain reversal in memory, actuation, sensing, and optical applications. The ability of the tip-induced proximity switching to differentially switch multilayers, based on the order of the layers, provides a powerful tool for selective domain engineering.

36 MATERIALS SCIENCE

The Use of Machine Learning Models for Predicting the Dielectric Strength of Gases

Technological advancements in high voltage systems have pushed sulfur hexafluoride (SF6) to its operational limits. Furthermore, this gas has other drawbacks including a high liquefaction temperature and a high global warming potential. Therefore, there has been an urgent need to find alternative gases with high dielectric strength (DS). In this work, density functional theory (DFT) is used to calculate molecular descriptors that are fed into an artificial neural network (ANN) and a random forest (RF). These machine learning (ML) models are then used to predict the DS for hundreds of molecules. A finite element model (FEM) is also used to calculate the electric field profile of multiple simple electrode geometries as the applied voltage to the system is increased. Results indicate that the random forest model has better generalization to unseen data than the neural network. The highest DS value predicted by the RF was 2.16 relative to the experimental DS of SF6. The results also demonstrate how choosing a gas with a higher DS and a geometry with minimal edges and corners can significantly increase the operating voltage of an electrical system. Due to its superior generalization, the RF represents the most promising path toward an accurate DS predictor once sufficient experimental data are available.

Mileski, Matthew [AFIT]

Modeling Structured Electrodes and Graded Porosity for Improving Discharge Rate Capability in Ultra-Thick Graphite|LiNi 0.6 Mn 0.2 Co 0.2 O 2 Batteries

Long-range electric vehicles (EVs) require high-energy-density batteries that also meet the power demands of high current charge and discharge. Ultra-thick (>100 μm) Lithium-ion battery electrodes are critical to enable this need, but slow ion transport in conventional uniform electrodes (UEs) reduces battery capacity at increasing charge/discharge rates. We present a 3D computational analysis on the impact of structured electrode (SE) and graded electrode (GE) geometries on the discharge rate capability of ultra-thick graphite|LiNi 0.6 Mn 0.2 Co 0.2 O 2 (NMC-622) battery cells based on the footprint of a commercial EV pouch cell. SE cathodes with either a “grid” or “line” geometry and GEs with two layers of porosity were modeled. Based on the results of 230 models, we found that the electrolyte volume fraction is a key parameter that impacts capacity improvements in UEs, GEs, and SEs at 2 C–6 C discharge rates. SEs have the greatest discharge rate capability, outperforming GEs and UEs due to reduced Lithium-ion concentration gradients across the electrode thickness, which mitigates electrolyte depletion at high rates. The best SE model has a “grid” geometry with gravimetric and volumetric energy density improvements of 0.9%–4% at C/2–2 C and 18%–24% at 4 C–6 C relative to UEs.

25 ENERGY STORAGE

Design and fabrication of a four-rod radio frequency quadrupole with a variable aperture

Radiofrequency quadrupole (RFQ) linear accelerators with variable aperture cell parameters offer an effective means of accelerating high-intensity ion beams exceeding 100 mA and have been implemented in various RFQ designs. However, in a four-rod RFQ, the use of dielectric perturbation methods to measure the voltage distribution behind the rod electrodes is not applicable when implementing variable aperture geometries, because changes in the rod shape influence the electric field at the location of the perturbation element. Therefore, it was necessary to design the rod profile for each cell such that the electric field distribution behind the rods remains unaffected by changes in the aperture. In this study, we present a method to adjust the voltage distribution by modifying the geometry of the rod electrodes according to each cell's aperture parameters. This design approach enables the integration of variable aperture cells into a four-rod RFQ linac, significantly enhancing the achievable beam current. To ensure uniform capacitance and maintain a constant rod-electrode voltage despite variations in the average aperture radius, we developed a systematic method for shaping the rod electrodes. This method was refined using two-dimensional (2D) and three-dimensional (3D) electromagnetic field simulations to optimize the electrode cross-sectional profiles. The fabricated rod electrodes were then installed in a four-rod RFQ, and the electric field distribution was measured using the perturbation method and tuned by adjusting the base plates. As a result, the measured variation of ±1.1 % confirms that this design methodology can produce a sufficiently uniform voltage distribution for high-current beam acceleration. These results demonstrate the practical feasibility of a variable-aperture four-rod RFQ linac.

43 PARTICLE ACCELERATORS

An Exploration into Elemental Additions to Tungsten Electrodes

The scope of this study is to evaluate the arc start capabilities of various tungsten electrode alloys—2% Thoriated, 2% Ceriated, 2% Lanthanated, and a proprietary Tri-Mix WS2 on Stainless Steel 304 and Molybdenum. Each electrode maintained consistent geometry: 0.040" diameter, 2.000" length, and a 0.0010" flat tip with a 40° tip angle. The testing procedure aimed to determine the minimum starting amperage required to achieve 10 consecutive successful arc starts per electrode type.

36 MATERIALS SCIENCE

A tetrahedral probe constellation approach for measuring canonical momentum in self-organized laboratory plasma

To examine momentum redistribution processes and study generalized helicities during plasma relaxation in Madison Symmetric Torus, MST, reversed field pinch plasma, a new probe is being tested to measure the full 3D plasma ion flow and magnetic field vectors at four spatial locations arranged in a tetrahedral shape reminiscent of a satellite measurement constellation. These measurements permit calculation of ∇ x $\vec{u}$ and canonical momentum, where $\vec{u}$ is the plasma ion flow vector. The probe consists of four probe heads arranged in a tetrahedral pattern, with an overall probe diameter of ∼31.75 mm. The probe head diameter is ∼1.0 cm, which is of the order of the ion Larmor radius. Each head has four molybdenum electrodes, also arranged in tetrahedral geometry, which are biased relative to a common return electrode, using four power supplies (one for each head), to measure the local ion flow. Additionally, each head has three orthogonal magnetic pickup coils within it to measure equilibrium and fluctuating magnetic fields.

Physics - Plasma physics

Characterization of Electronic Stress-Induced Changes in Multilayer MoS 2

Transition metal dichalcogenides like molybdenum disulfide (MoS 2 ) are compelling for next-generation electronic devices. In this work, we investigate the impact of electronic stress on MoS 2 to illustrate that observational and phenomenological information on multiple devices can be useful to describe changes in the device, and caution against the rationalization of paltry results as representative or correlative to device behavior. Here, we stress MoS 2 by applying a sustained 20 V DC bias to study the material’s response. Post-stress electronic characterization revealed nonuniform shifts in current–voltage (I–V) behavior alongside microscale changes. Complementary mechanical, spectroscopic, and scanning microwave impedance measurements showed that stress-induced features locally modulate stiffness, surface potential, Raman intensity, and charge carrier density. We correlated I–V behavior with morphological features (wrinkles, tears, folds, height) and device-level geometry (MoS 2 overlap with electrodes, channel area, contact length) on 50 test structures across five chips to move beyond anecdotal conclusions. We found no universal correlations before DC stress. However, device-level geometry was correlated with I–V behavior after DC stress, suggesting that electrode contacts play a more dominant role than morphology in determining performance. Delamination and thinning induced by DC stress led to localized reductions in charge carrier density within the affected regions. Further, delamination and thinning appear to map to I–V device performance in a few samples, but the correlation is lost when a larger sample size is considered. This suggests significant sample-to-sample variability in surface electronic states of the test structures. We also discuss how environmental factors introduced during fabrication may contribute to the observed heterogeneous device response. Progress will require high-resolution, multimodal analysis across many samples constructed under controlled, clean conditions. By building data sets that capture variability, we can better identify the true drivers of performance.

36 MATERIALS SCIENCE