Results for “Thermodynamic States and Processes” — Page 68

DOE OSTI2024

Eliminating chemo-mechanical degradation of lithium solid-state battery cathodes during >4.5 V cycling using amorphous Nb2O5 coatings

Abstract Lithium solid-state batteries offer improved safety and energy density. However, the limited stability of solid electrolytes (SEs), as well as irreversible structural and chemical changes in the cathode active material, can result in inferior electrochemical performance, particularly during high-voltage cycling (>4.3 V vs Li/Li + ). Therefore, new materials and strategies are needed to stabilize the cathode/SE interface and preserve the cathode material structure during high-voltage cycling. Here, we introduce a thin (~5 nm) conformal coating of amorphous Nb 2 O 5 on single-crystal LiNi 0.5 Mn 0.3 Co 0.2 O 2 cathode particles using rotary-bed atomic layer deposition (ALD). Full cells with Li 4 Ti 5 O 12 anodes and Nb 2 O 5 -coated cathodes demonstrate a higher initial Coulombic efficiency of 91.6% ± 0.5% compared to 82.2% ± 0.3% for the uncoated samples, along with improved rate capability (10x higher accessible capacity at 2C rate) and remarkable capacity retention during extended cycling (99.4% after 500 cycles at 4.7 V vs Li/Li + ). These improvements are associated with reduced cell polarization and interfacial impedance for the coated samples. Post-cycling electron microscopy analysis reveals that the Nb 2 O 5 coating remains intact and prevents the formation of spinel and rock-salt phases, which eliminates intra-particle cracking of the single-crystal cathode material. These findings demonstrate a potential pathway towards stable and high-performance solid-state batteries during high-voltage operation.

Science & Technology - Other Topics↗

NASA NTRS2026

Thermochemically-Closed Sonic-Flow Inversion for Enthalpy and Temperature in Multispecies Arc-Jet Flows

A thermochemically-closed sonic-flow inversion framework (TSIF) is developed to infer bulk enthalpy and total temperature upstream of a choked nozzle in arc-jet flows. The formulation recasts a pressure-rise total enthalpy quantification technique as an inverse problem in characteristic-velocity c * space using measured mass flow rate, upstream total pressure, gas composition, and nozzle throat geometry as inputs. Unlike calorimetric energy-balance approaches or optical diagnostics, the method relies primarily on routinely measured facility quantities combined with explicit thermochemical closure. Thermochemical states are obtained using NASA’s open-source Chemical Equilibrium with Applications (CEA) code, enabling construction of a chemistry-consistent relation between characteristic velocity, total enthalpy, and total temperature under equilibrium or frozen assumptions. A discharge coefficient is self-calibrated using cold-flow (arc-off) operation data and applied to hot-flow (arc-on) measurements, enabling upstream losses to be accounted for without empirical correlations. The framework is applied to air, N 2 , and CO 2 –N 2 arc-jet flows and demonstrates expected trends for the inferred thermochemical states as function of arc power, specific energy input, mass-flow, heater configuration, and test gas. In the air limit, under equilibrium assumptions, the method recovers the classical high-enthalpy asymptotic correlation of Winovich with a mean residual of 4.4%, demonstrating compatibility with established sonic-flow scaling, while extending applicability to arbitrary multi-species mixtures and non-equilibrium chemistry. The framework provides a mixture-flexible methodology for determining bulk thermochemical states in modern arc-jet environments using routine facility pressure, mass-flow, gas-composition, and nozzle-geometry information together with a cold-flow calibration.

stagnation heat flux↗

NASA NTRS2026

Thermochemically-Closed Sonic-Flow Inversion for Enthalpy and Temperature in Multispecies Arc-Jet Flows

A thermochemically-closed sonic-flow inversion framework (TSIF) is developed to infer bulk enthalpy and total temperature upstream of a choked nozzle in arc-jet flows. The formulation recasts a pressure-rise total enthalpy quantification technique as an inverse problem in characteristic-velocity c * space using measured mass flow rate, upstream total pressure, gas composition, and nozzle throat geometry as inputs. Unlike calorimetric energy-balance approaches or optical diagnostics, the method relies primarily on routinely measured facility quantities combined with explicit thermochemical closure. Thermochemical states are obtained using NASA’s open-source Chemical Equilibrium with Applications (CEA) code, enabling construction of a chemistry-consistent relation between characteristic velocity, total enthalpy, and total temperature under equilibrium or frozen assumptions. A discharge coefficient is self-calibrated using cold-flow (arc-off) operation data and applied to hot-flow (arc-on) measurements, enabling upstream losses to be accounted for without empirical correlations. The framework is applied to air, N 2 , and CO 2 –N 2 arc-jet flows and demonstrates expected trends for the inferred thermochemical states as function of arc power, specific energy input, mass-flow, heater configuration, and test gas. In the air limit, under equilibrium assumptions, the method recovers the classical high-enthalpy asymptotic correlation of Winovich with a mean residual of 4.4%, demonstrating compatibility with established sonic-flow scaling, while extending applicability to arbitrary multi-species mixtures and non-equilibrium chemistry. The framework provides a mixture-flexible methodology for determining bulk thermochemical states in modern arc-jet environments using routine facility pressure, mass-flow, gas-composition, and nozzle-geometry information together with a cold-flow calibration.

inviscid theory↗

NASA NTRS2026

Modeling Large Dust Aerosols in the Community Earth System Model Version 2 (CESM2)

Dust aerosols have a wide size distribution from less than 0.1 to over 100 μm and dominate Earth's atmospheric aerosol mass. However, most Earth system models (ESMs) inadequately represent dust aerosols larger than 10 μm in diameter, limiting the accuracy of the simulated dust cycle and climate impacts. Here, we introduce a new modeling framework that captures the full observed size distribution of dust aerosols, incorporating recent advances into a mineral-resolved version of the Community ESM, while addressing known issues in previous versions. Comprehensive evaluation against diverse observations of bulk dust and component minerals demonstrates that the model reproduces the observed dust cycle across particle sizes. Incorporating the previously unrepresented large-dust fractions substantially alters dust budget estimates, highlighting potential changes in simulated climate impacts and underscoring the importance of comprehensive size-resolved dust modeling. Despite these advancements, uncertainties persist. Our results indicate that a size-dependent reduction in settling velocity is required to reproduce the observed dust size distribution downwind of source regions. Specifically, in the new model, the gravitational settling velocity of dust particles larger than 10 μm in diameter must be reduced by as much as 85% to achieve agreement with observations. This empirical reduction serves as a constraint on physics-based models of dust settling. Future developments should address misrepresented physical processes that hinder accurate modeling of the large dust aerosol transport. Expanding observational data sets covering the full-size distribution is also essential to better constrain the dust cycle and improve the representation of dust optical properties and climate effects.

mineral dust↗

DOE OSTI2026

Lunar accelerometer network gravitational observatory (LANGO)

With ground-based interferometers detecting hundreds of gravitational-wave (GW) events, GW astronomy has continued to blossom. U.S. and European scientists are developing plans to construct third-generation ground-based interferometers. ESA has proceeded through the mission formulation phase of space-based interferometer in a lower-frequency band, 10 –4 –0.1 Hz. Despite all these exciting developments, there is still a missing frequency band, 0.1–10 Hz. This mid-frequency band is rich with interesting astrophysical events. Coalescence and merger of intermediate-mass black holes (IMBHs) will occur in this frequency band. Coalescing stellar-mass BHs will pass through this frequency band days before they reach the frequency band of LIGO and Virgo. Detection of such signals would enable a mid-frequency detector to issue an advance notice to the high-frequency GW detectors, as well as to optical, x-ray and γ-ray telescopes. We propose Lunar Accelerometer Network Gravitational Observatory (LANGO) to detect GWs in this frequency band. In the first phase (LANGO 1), we propose to deploy four ambient-temperature (250 K) accelerometers in the tetrahedral or in a square configuration on the hemisphere facing the Earth. After successful operation at 250 K, LANGO would be upgraded to a cryogenic (4 K) version with over two orders of magnitude increased sensitivity (LANGO 2), with coherent rejection of seismic noise implemented. LANGO is a full-tensor detector, capable of determining the source direction and wave polarization. Each test mass (TM) is suspended as a pendulum with resonance frequency ∼ 0.01 Hz, thus is only weakly coupled to the lunar surface horizontally. LANGO is designed to detect the relative motion of globally separated, nearly free, TMs by using the Moon as a large quiet platform. LANGO 1 and 2 accelerometers aim at sensitivities ⩽ 10 –11 m s –2 Hz –1/2 and ⩽ 10 –13 m s –2 Hz –1/2 in the horizontal axes over the frequency band of 1 mHz–10 Hz, which yield GW sensitivities 1.1 x 10 -21 Hz -1/2 and 3.8 x 10 -24 Hz -1/2 at 1 Hz, respectively. The LANGO accelerometers will be 10 3 –10 5 times more sensitive than Apollo seismometers. With such sensitivity, LANGO will also make great contribution to the advancement of lunar geophysics.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2026

A New Shutdown Dose Rate Benchmark Problem for Representative Fusion Applications

Here, this work introduces a new benchmark problem for calculating shutdown dose rates (SDDRs) aimed at fusion reactor applications. The model is designed to represent a simplified version of a typical ITER port plug. The responses of interest include neutron flux, gamma flux, and gamma SDDR at 12 different locations scattered throughout the port. This article outlines the geometry specifications of the problem, provides material definitions for the components, specifies the required responses to be calculated, and presents the source definition information. The need for this benchmark arises from the limited availability of publicly accessible references, with only one benchmark representing the typical dimensions and materials found in fusion systems. This existing benchmark has been cited extensively, reflecting the demand within the scientific community to test both established and novel workflows for SDDR calculations. However, since its presentation at a conference in 2011, the results have become increasingly well known. Moreover, the absence of formal publication and peer review has led to the details of this benchmark being extracted from secondary sources, such as subsequent studies that reference it. As a result, analysts are left with significant flexibility in interpreting the key parameters, which can be adjusted to account for unknown systematic errors, ultimately reproducing the already well-known responses. This new benchmark serves as an updated version of that earlier work, with the aim of providing a more reliable description of the materials and their impurities, which is crucial for assessing activation and subsequent gamma emission. Additionally, it seeks to provide a geometry that more closely represents an ITER port plug. The improvements in the problem definition will lead to a more reproducible benchmark problem, while also presenting the radiation transport community with a completely new challenge. The results will be published in a future article to allow analysts adequate time to analyze this problem independently.

Benchmark↗

NASA NTRSDate not supplied

Laser Beam Welding Advancements for In-Space Servicing, Assembly, and Manufacturing

NASA is currently working to develop in-space servicing, assembly, and manufacturing (ISAM) capabilities for low Earth orbit and the lunar surface. One crucial technology for this effort is laser beam welding. Laser systems can perform joining, cleaning, cutting, and repair activities, which will enable the construction of large in-space structures that could not fit on a single launch vehicle, such as trusses for solar panels, radiators, or communications infrastructure. Multiple projects studying laser welding for space applications are currently underway at NASA Marshall Space Flight Center. One of these, the DIsk-Shaped Configurable and Modular vAcuum uNit (DISCMAN), is a compact vacuum chamber designed to support parameter development for laser welding in microgravity. It contains a rotating platen with weld samples made from aluminum, steel, and titanium, a high-power infrared laser, and integrated pumps for pulling vacuum inside the sample cartridge. The DISCMAN payload is planned to launch to the International Space Station, where welds will be performed under sustained microgravity inside the Bishop Airlock. Another effort underway at Marshall is the Lunar Assembly and Servicing by Autonomous Robotics (LASAR) initiative. This project uses a space-rated robotic arm equipped with a laser weld head, wire feeder, and multiple cameras to perform welds in a thermal vacuum chamber simulating the lunar surface environment. Some of these welds are done on snowflake joints, which are specially designed to slot together to join segments of trussN structures, allowing for the construction of tall surface infrastructure. DISCMAN, LASAR, and other projects are being carried out to advance the technological maturity of in-space laser beam welding, collect data to inform computational models, and learn reliable processes for creating weld joints in space. This work supports NASA’s greater goals to expand humanity’s presence in low Earth orbit, establish a permanent moon base, and eventually send crewed missions to Mars and beyond.

Robotics↗

NASA NTRS2026

Performance Optimization for Lunar Extravehicular Activity Readiness (POLAR) Study: Methods Paper

To better understand which aspects of physical fitness may be most related to performance during Lunar surface operations and thereby help to inform the current NASA fitness standards, much can be learned from fields encompassing the “tactical athlete.” Other physically demanding professions such as law enforcement, military, or rescue professionals often require candidates to meet occupationally-relevant fitness standards. The determination of such standards is a multistep process, including both objective and subjective measures, to determine tasks essential to occupational performance and identify the minimal fitness profile needed to meet physical demands of the job. Notably, fitness is only one component which may contribute to the demands of astronaut selection, flight assignment, and occupational performance. Utilizing a framework to systematically determine which domains of fitness most contribute to relevant job tasks can aid in the refinement of current NASA-STD-3001 fitness standards. Therefore, NASA’s Exercise Physiology & Countermeasures Laboratory conducted the Performance Optimization for Lunar Extravehicular Activity Readiness (POLAR) study to identify and examine a comprehensive list of fitness tests (including NASA-STD-3001 assessments: 1-Repetition Maximum [1-RM] bench press and deadlift) and determine preliminary relationships between identified fitness parameters and novel Artemis-relevant tasks to help inform future investigations for the continued development of aerobic and muscular fitness standards for surface EVAs. This was accomplished through 1) a review of the literature relating fitness assessments to simulated or real EVA performance to identify fitness tests that are most correlated with simulated EVA task performance; 2) a task analysis following a modified framework for physical employment standards development to down select mission critical tasks; and 3) development and pilot testing of a novel, portable Artemis-relevant EVA task circuit to relate to a battery of fitness assessments. The current report describes the methodology used to complete the task analysis, EVA task circuit development, and the pilot study.

Nicole C Strock↗

NASA NTRS2020

The Slow Demise of the Long-Lived Sn 2005ip

The Type IIn supernova (SN IIn) 2005ip is one of the most well-studied and long-lasting examples of an SN interacting with its circumstellar environment. The optical light curve plateaued at a nearly constant level for more than five years, suggesting ongoing shock interaction with an extended and clumpy circumstellar medium (CSM). Here, we present continued observations of the SN from ∼1000 to 5000 d post-explosion at all wavelengths, including X-ray, ultraviolet, near-infrared (NIR), and mid-infrared. The UV spectra probe the pre-explosion mass loss and show evidence for CNO processing. From the bolometric light curve, we find that the total radiated energy is in excess of 10 50 erg, the progenitor star’s pre-explosion mass-loss rate was >~ 1 × 10 −2 Mꙩ yr −1 , and the total mass lost shortly before explosion was >~1Mꙩ, though the mass lost could have been considerably larger depending on the efficiency for the conversion of kinetic energy to radiation. The ultraviolet through NIR spectrum is characterized by two high-density components, one with narrow high-ionization lines, and one with broader low-ionization H I, He I, [OI], Mg II, and Fe II lines. The rich Fe II spectrum is strongly affected by Lyα fluorescence, consistent with spectral modelling. Both the Balmer and He I lines indicate a decreasing CSM density during the late interaction period. We find similarities to SN 1988Z, which shows a comparable change in spectrum at around the same time during its very slow decline. These results suggest that, at long last, the shock interaction in SN 2005ip may finally be on the decline.

stars↗

DOE OSTI2024

Complex orders and chirality in the classical Kitaev-Γ model

It is well recognized that the low-energy physics of many Kitaev materials is governed by two dominant energy scales, the Ising-type Kitaev coupling 𝐾 and the symmetric off-diagonal Γ coupling. An understanding of the interplay between these two scales is therefore the natural starting point toward a quantitative description that includes subdominant perturbations that are inevitably present in real materials. This study focuses on the classical 𝐾−Γ model on the honeycomb lattice, with a specific emphasis on the region 𝐾< 0 and Γ > 0 , which is the most relevant for the available materials and which remains enigmatic in both quantum and classical limits, despite much effort. We employ large-scale Monte Carlo simulations on specially designed finite-size clusters and unravel the presence of a complex multisublattice magnetic order in a wide region of the phase diagram, whose structure is characterized in detail. We show that this order can be quantified in terms of a coarse-grained scalar-chirality order, featuring a counterrotating modulation on the two spin sublattices. Here, we also provide a comparison to previous studies and discuss the impact of quantum fluctuations on the phase diagram.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2025

Sandia National Labs - Access Delay Overview

This is a presentation discussing some of the recent activities of the Sandia Access Delay department

Hammond, Elijah David [Sandia National Laboratorie↗

DOE OSTI2025

Technology & Design Innovations to Maximize the Reduction Effect on DCFC Unit Cost Economics (Max‐REDUCE)

With the high cost and reliability issues of DCFC infrastructure the project specifically addressed the needs of innovating on the AC to DC power conversion. We proved the integration of a new power module into existing charger units, and that theoretical educational work on single stage energy conversion is commercially viable, and a reduced cost compared to existing methods. Our A sample shown 3 PH 480VAC 60hz to 350VDC 8.5kW power transfer, and our software and simulation proved efficiency greater than 97% from 20kW to > 60kW, over the DC output voltage range 200 – 1200V. The project was stopped before we could implement B sample and execute full power testing, validation and integration into a demonstration charger (budget period 3).

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2025

Complex-Valued Intermolecular Coupling Enables Directional Exciton Transport in Excitonic Circuits

Molecular systems capable of directing the flow of excitons are key to the development and optimization of optoelectronic materials. The transport of excitons across multiple molecules is governed by an intermolecular electronic coupling network. In this article, we consider the effects of complex-valued intermolecular electronic coupling on exciton transport. Here, we present a molecular motif capable of generating complex-valued coupling under excitation with circularly polarized light. We use theoretical modeling and simulation to illustrate how complex coupling can be leveraged to drive the rotational flux of excitons in cyclic molecular networks and direct the exciton population in branched molecular networks.

chromophores↗

DOE OSTI2026

Performance of Heterostructural TaC/AlGaN Schottky Diodes Based on First Principles Electronic Structure Properties

Advances in ultra-wide bandgap materials, such as high Al-content AlxGa1-xN (AlGaN), are essential for next generation power electronics, but the requirement for lattice matched substrates is currently a significant obstacle. Recently, conductive TaC has emerged as a promising virtual substrate for AlGaN heteroepitaxy, with wurtzite (0001) AlxGa1-xN lattice-matched to rocksalt (111) TaC at x ~ 0.5. Thus, understanding and controlling the electronic properties of the TaC/AlGaN interface is key for developing technological applications based on TaC/AlGaN devices. Using density functional theory and electronic structure calculations, we here investigate TaC/Al0.5Ga0.5N interfaces, where we include explicit alloy models in the slab calculations. We predict the Schottky barrier height and the electric field discontinuity resulting from interface charges. Considering all possible combinations of (Ta, C) substrate termination, (Al/Ga, N) nucleation, and (Al/Ga, N) polarity, we construct a chemical potential phase diagram to identify the stable interfaces that can be accessed through variation of the synthesis conditions. The predicted interface electronic properties are implemented in device performance simulations to demonstrate a practical design for a strain-free, high-efficiency TaC/AlGaN Schottky diode with a low barrier height and without interface charges, underscoring the potential of TaC as a substrate for ultra-wide bandgap devices.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2025

Rethinking 𝛼−RuCl 3 : Parameters, models, and phase diagram

RuCl 3 was likely the first ever deliberately synthesized ruthenium compound, following the discovery of the 44 Ru element in 1844. For a long time it was known as an oxidation catalyst, with its physical properties being discrepant and confusing, until a decade ago when its allotropic form 𝛼−RuCl 3 rose to exceptional prominence. This “rediscovery” of 𝛼−RuCl 3 has not only reshaped the hunt for a material manifestation of the Kitaev spin liquid, but it has opened the floodgates of theoretical and experimental research in the many unusual phases and excitations that the anisotropic-exchange magnets as a class of compounds have to offer. Given its importance for the field of Kitaev materials, it is astonishing that the low-energy spin model that describes this compound and its possible proximity to the much-desired spin-liquid state is still a subject of significant debate ten years later. In the present study, we argue that the existing key phenomenological observations put strong natural constraints on the effective microscopic spin model of 𝛼−RuCl 3 , and specifically on its spin-orbit-induced anisotropic-exchange parameters that are responsible for the nontrivial physical properties of this material. These constraints allow one to focus on the relevant region of the multidimensional phase diagram of the 𝛼−RuCl 3 model, suggest an intuitive description of it via a different parametrization of the exchange matrix, offer a unifying view on the earlier assessments of its parameters, and bring closer together several approaches to the derivation of anisotropic-exchange models. We explore extended phase diagrams relevant to the 𝛼−RuCl 3 parameter space using quasiclassical, Luttinger-Tisza, exact diagonalization, and density-matrix renormalization-group methods, demonstrating a remarkably close quantitative accord between them on the general structure and hierarchy of the phases, with the zigzag, ferromagnetic, and incommensurate phases that are proximate to each other. As a result, one of the highlights is the detailed agreement on the nature of the incommensurate phases that realize two distinct counterrotating helical states.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

NASA NTRSDate not supplied

Laser Beam Welding Advancements for In-Space Servicing, Assembly, and Manufacturing

NASA is currently working to develop in-space servicing, assembly, and manufacturing (ISAM) capabilities for low Earth orbit and the lunar surface. One crucial technology for this effort is laser beam welding. Laser systems can perform joining, cleaning, cutting, and repair activities, which will enable the construction of large in-space structures that could not fit on a single launch vehicle, such as trusses for solar panels, radiators, or communications infrastructure. Multiple projects studying laser welding for space applications are currently underway at NASA Marshall Space Flight Center. One of these, the DIsk-Shaped Configurable and Modular vAcuum uNit (DISCMAN), is a compact vacuum chamber designed to support parameter development for laser welding in microgravity. It contains a rotating platen with weld samples made from aluminum, steel, and titanium, a high-power infrared laser, and integrated pumps for pulling vacuum inside the sample cartridge. The DISCMAN payload is planned to launch to the International Space Station, where welds will be performed under sustained microgravity inside the Bishop Airlock. Another effort underway at Marshall is the Lunar Assembly and Servicing by Autonomous Robotics (LASAR) initiative. This project uses a space-rated robotic arm equipped with a laser weld head, wire feeder, and multiple cameras to perform welds in a thermal vacuum chamber simulating the lunar surface environment. Some of these welds are done on snowflake joints, which are specially designed to slot together to join segments of trussN structures, allowing for the construction of tall surface infrastructure. DISCMAN, LASAR, and other projects are being carried out to advance the technological maturity of in-space laser beam welding, collect data to inform computational models, and learn reliable processes for creating weld joints in space. This work supports NASA’s greater goals to expand humanity’s presence in low Earth orbit, establish a permanent moon base, and eventually send crewed missions to Mars and beyond.

Manufacturing↗

DOE OSTI2026

Optimizing Porous Transport Layer Porosity for Proton Exchange Membrane Water Electrolysis

An empirical model is presented that describes anode-side losses related to porous transport layer (PTL) morphology in proton exchange membrane water electrolysis (PEMWE). The model is based on an advanced voltage breakdown analysis that links various overpotentials to PTL morphology. Custom Ti PTLs, spanning uncommonly low porosities (22 - 31%), were fabricated and analyzed with X-ray CT to obtain pore and particle size distributions. Particle size distributions were consistent across samples with an average particle diameter of 12.0?..mu..m, whereas average pore diameters ranged from 6.0 to 7.0?..mu..m. The PTLs were tested in standard PEMWE cell assemblies with anode catalyst loadings of 0.1 mgIr cm-2 to obtain polarization curves, electrochemical impedance spectra, and augmented Tafel analysis. The PTL-dependent anode side losses were deconvoluted and assigned to excess utilization, concentration, ion transport resistance, and electrical contact resistance overpotentials. The data and model reveal an optimal 20 - 28% PTL porosity region where utilization and contact resistance overpotentials are minimized without triggering concentration and ion transport losses related to water deprivation. The optimal PTL porosity depends on the operating current density and is demonstrated at realistic PEMWE water flow rates to establish PTL design guidance for operation at scale.

08 HYDROGEN↗

DOE OSTI2026

Ferroelectric Fractals: Switching Mechanism of Wurtzite AlN

The advent of wurtzite ferroelectrics is enabling new ferroelectric devices for computer memory that have the potential to bypass the von Neumann bottleneck due to their robust polarization and silicon compatibility. However, the atomistic switching mechanism of wurtzites is still undetermined due to the limitations of density functional theory simulation size and experimental temporal and spatial resolution. Thus, physics-informed materials engineering to reduce coercive field and breakdown in these devices has been limited. In this work, the atomistic mechanism of domain wall migration and domain growth in aluminum nitride-based wurtzites is uncovered using molecular dynamics and Monte Carlo simulations. We reveal the anomalous switching mechanism of fast 1D single columns of atoms propagating from a slow-moving 2D fractallike domain wall. We find that the critical nucleus is a single aluminum ion that breaks its bond with one nitrogen and bonds to another nitrogen; this creates a cascade that flips atoms directly only in the same column, due to the extreme locality (sharpness) of the domain walls in wurtzites. We further show how the fractallike shape of the domain wall in the 2D plane breaks assumptions in the Kolmogorov, Avrami, and Ishibashi (KAI) model and leads to the anomalously fast switching in wurtzite structured ferroelectrics.

36 MATERIALS SCIENCE↗