Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “stability regime”

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 271 records · Page 15

A discontinuous piecewise polynomial generalized moving least squares scheme for robust finite element analysis on arbitrary grids

A variational approach is developed with a meshless discretization to enable accurate and robust numerical simulation of partial differential equations for meshes that are of poor quality. Traditional finite element methods use the mesh to both discretize the geometric domain and to define the finite element shape functions. The latter creates a dependence between the quality of the mesh and the properties of the finite element basis that may adversely affect the accuracy of the discretized problem. Here, we propose a new approach for defining finite element shape functions that breaks this dependence and separates mesh quality from the discretization quality, which we call discontinuous piecewise polynomial generalized moving least squares (DPP-GMLS). At the core of the approach is a meshless definition of the shape functions, which limits the purpose of the mesh to representing the geometric domain and integrating the basis functions without having any role in their approximation quality. The resulting non-conforming space can be utilized within a standard discontinuous Galerkin framework, providing a rigorous foundation for solving partial differential equations on low-quality meshes. We present a collection of numerical experiments demonstrating our approach in a wide range of settings: strongly coercive elliptic problems, linear elasticity in the compressible regime, and the stationary Stokes problem. We demonstrate convergence for all problems and stability for element pairs for problems which usually require inf-sup compatibility for conforming methods, also referring to a minor modification possible through the symmetric interior penalty Galerkin framework for stabilizing element pairs that would otherwise be traditionally unstable. Mesh robustness is particularly critical for elasticity, and we provide an example that our approach provides a greater than 5 x improvement in accuracy and allows for taking an 8 x larger stable timestep for a highly deformed mesh, compared to the continuous Galerkin finite element method.

97 MATHEMATICS AND COMPUTING↗

The effect of shear flow on the resistive tearing instability

We develop a new scaling theory for the resistive tearing mode instability of a current sheet with a strong shear flow across the layer. The growth rate decreases with increasing flow shear and is completely stabilized as the shear flow becomes Alfvénic: both in the constant-Ψ regime, as in previous results, but we also show that the growth rate is in fact suppressed more strongly in the nonconstant-Ψ regime. As a consequence, for sufficiently large flow shear, the maximum of the growth rate is always affected by the shear suppression, and the wavenumber at which this maximum growth rate is attained is an increasing function of the strength of the flow shear. These results may be important for the onset of reconnection in imbalanced MHD turbulence.

astrophysical plasmas↗

High-Performance Near-Infrared Quantum Emission from Color Centers in hBN

Color centers hosted in hexagonal boron nitride (hBN) have emerged as a highly promising platform for single-photon emission and spin-photon technologies relevant to quantum communication and quantum networking. As a wide bandgap van der Waals material, hBN can host optically active quantum defects across a broad spectral range. Here, we demonstrate a simple and scalable oxygen-plasma process that reproducibly creates single quantum emitters in hBN with blinking-free zero-phonon lines (ZPLs) spanning near-infrared (NIR) from 700 up to 971 nm. These emitters combine MHz-level brightness, single-photon purity up to 99.9%, and ultranarrow cryogenic line widths down to 2.7 GHz under quasi-resonant excitation, placing them in a particularly attractive regime for quantum photonics. Photostability measurements further reveal resistance to photobleaching, subnanometer spectral stability over long time scales, and near-shot-noise-limited intensity fluctuations. Analysis of the phonon sidebands shows weak vibronic coupling and ZPL-dominated emission, with Debye–Waller factors approaching 50%. Control experiments together with elemental mapping support oxygen incorporation as a necessary ingredient in activating the NIR emitter population, while first-principles calculations identify O N V N and O N V N H as the leading defect candidates. These results establish a high-performance NIR quantum-emitter platform in hBN for free-space quantum networking and future integrated quantum-photonic architectures.

2D materials↗

Prospects and challenges of colloidal quantum dot laser diodes

Colloidal quantum dots (QDs) combine the superior light-emission characteristics of quantum-confined semiconductors with the chemical flexibility of molecular systems. We find that these properties could, in principle, enable solution-processable laser diodes with an ultrawide range of accessible colours. However, the realization of such devices has been hampered by fast optical gain decay due to non-radiative Auger recombination and poor stability of QD solids at the high current densities required for the lasing regime. Recently, these problems have been resolved, which resulted in the development of electrically pumped optical gain devices operating at ultrahigh current densities of around 1,000 A cm -2 . The next step is the realization of a QD laser diode (QLD). Here we assess the status of the QD lasing field, examine the remaining challenges on the path to a QLD and discuss practical strategies for attaining electrically pumped QD lasing.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Regression Based Approach for Robust Finite Element Analysis on Arbitrary Grids. LDRD Final Report

This report summarizes the work performed under a one-year LDRD project aiming to enable accurate and robust numerical simulation of partial differential equations for meshes that are of poor quality. Traditional finite element methods use the mesh to both discretize the geometric domain and to define the finite element shape functions. The latter creates a dependence between the quality of the mesh and the properties of the finite element basis that may adversely affect the accuracy of the discretized problem. In this project, we propose a new approach for defining finite element shape functions that breaks this dependence and separates mesh quality from the discretization quality. At the core of the approach is a meshless definition of the shape functions, which limits the purpose of the mesh to representing the geometric domain and integrating the basis functions without having any role in their approximation quality. The resulting non-conforming space can be utilized within a standard discontinuous Galerkin framework providing a rigorous foundation for solving partial differential equations on low-quality meshes. We present a collection of numerical experiments demonstrating our approach in a wide range of settings: strongly coercive elliptic problems, linear elasticity in the compressible regime, and the stationary Stokes problem. We demonstrate convergence for all problems and stability for element pairs for problems which usually require inf-sup compatibility for conforming methods, also referring to a minor modification possible through the symmetric interior penalty Galerkin framework for stabilizing element pairs that would otherwise be traditionally unstable. Mesh robustness is particularly critical for elasticity, and we provide an example that our approach provides a greater than 5x improvement in accuracy and allows for taking an 8x larger stable timestep for a highly deformed mesh, compared to the continuous Galerkin finite element method. The report concludes with a brief summary of ongoing projects and collaborations that utilize or extend the products of this work.

97 MATHEMATICS AND COMPUTING↗

Spin orbit torque in ferromagnet/topological-quantum-material heterostructures

The main objective of this proposal was to investigate the spin-orbital driven phenomena in topological insulator (TI)/ferromagnet (FM) or transition metal (TM) heterostructures. In particularly, we tried to utilize the spin-momentum locking properties on the surface of TIs to achieve ideal spin-orbit torques (SOTs). In practice, we found that the surface properties have been greatly modified, resulting in much more complicated magneto-transport behaviors than what is expected in SOT-driven magnetotransport. Through systematic and extended investigations, we start to unravel that the nonlinear magnetotransport properties, which was, in the past, mainly attributed to SOT-driven phenomena, can be used to infer the modification of the electronic band structures of TIs and spin textures in FMs which are closely related with charge-spin interconversion phenomena. Regarding the band structures In TIs, we find under electrostatic perturbations imposed by deposition of nonmagnetic TMs, the quantum well states formed on the surface of TIs develop substantial Rashba spin-splitting, which is detected as a bilinear magnetoresistance (BMR). Under magnetic perturbations, the nonlinear planar Hall effect (NPHE) in the topological surface states is greatly enhanced, possibly reflecting the modified surface band structure due to time-reversal symmetry breaking. In a double HM/FM heterostructure where the magnetic anisotropy can be robustly tuned, we observed second-harmonic Hall voltage responses corresponding to domain wall structures with fixed chirality, stabilized by interfacial Dzyaloshinskii–Moriya interaction (DMI). In addition, in the linear transport regime, we utilized inverse-SHE-induced anomalous Hall effect in Pt to probe the exchange bias on the surface of a doped magnetic insulator Tm 3 Fe 5 O 12 .

36 MATERIALS SCIENCE↗

An experimental evaluation of the Sternberg task as a workload metric for helicopter Flight Handling Qualities (FHQ) research

The objective was to determine whether the Sternberg item-recognition task, employed as a secondary task measure of spare mental capacity for flight handling qualities (FHQ) simulation research, could help to differentiate between different flight-control conditions. FHQ evaluations were conducted on the Vertical Motion Simulator at Ames Research Center to investigate different primary flight-control configurations, and selected stability and control augmentation levels for helicopters engaged in low-level flight regimes. The Sternberg task was superimposed upon the primary flight-control task in a balanced experimental design. The results of parametric statistical analysis of Sternberg secondary task data failed to support the continued use of this task as a measure of pilot workload. In addition to the secondary task, subjects provided Cooper-Harper pilot ratings (CHPR) and responded to workload questionnaire. The CHPR data also failed to provide reliable statistical discrimination between FHQ treatment conditions; some insight into the behavior of the secondary task was gained from the workload questionnaire data.

Hemingway, J. C.↗

The Sternberg Task as a Workload Metric in Flight Handling Qualities Research

The objective of this research was to determine whether the Sternberg item-recognition task, employed as a secondary task measure of spare mental capacity for flight handling qualities (FHQ) simulation research, could help to differentiate between different flight-control conditions. FHQ evaluations were conducted on the Vertical Motion Simulator at Ames Research Center to investigate different primary flight-control configurations, and selected stability and control augmentation levels for helicopers engaged in low-level flight regimes. The Sternberg task was superimposed upon the primary flight-control task in a balanced experimental design. The results of parametric statistical analysis of Sternberg secondary task data failed to support the continued use of this task as a measure of pilot workload. In addition to the secondary task, subjects provided Cooper-Harper pilot ratings (CHPR) and responded to a workload questionnaire. The CHPR data also failed to provide reliable statistical discrimination between FHQ treatment conditions; some insight into the behavior of the secondary task was gained from the workload questionnaire data.

Hemingway, J. C.↗

Experimental evaluation of active-member control of precision structures

The results of closed loop experiments that use piezoelectric active-members to control the flexible motion of a precision truss structure are described. These experiments are directed toward the development of high-performance structural systems as part of the Control/Structure Interaction (CSI) program at JPL. The focus of CSI activity at JPL is to develop the technology necessary to accurately control both the shape and vibration levels in the precision structures from which proposed large space-based observatories will be built. Structural error budgets for these types of structures will likely be in the sub-micron regime; optical tolerances will be even tighter. In order to achieve system level stability and local positioning at this level, it is generally expected that some form of active control will be required.

Fanson, James↗

NSTX-U National Research Program: White Paper in Response to Call from FESAC Sub-Committee

Both scientific and technical innovation is needed for the realization of an attractive engineering solution for a timely and cost-effective Pilot Plant, the design and construction of which is the overarching recommendation of the FESAC Long Range Plan, and the 2021 NASEM Pilot Plant reports, which underpin the Bold Decadal Vision. The two most significant plasma physics gaps to close for a Compact Pilot Plant (CPP) are core confinement improvement and heat flux mitigation, neither of which have been closed in an integrated fashion for any planned fusion power production device. High core confinement and stability are essential for producing majority self-driven plasmas in CPPs with reduced size and auxiliary heating power requirements, with an improvement in confinement being the major driver for cost reduction of a CPP. The National Spherical Tokamak Experiment - Upgrade (NSTX-U) is a unique low aspect ratio research facility that will address the fundamental challenge of developing the science and technology basis for a CPP design that integrates high core and edge confinement with the ability to mitigate very high incident heat fluxes. NSTX-U capabilities will enable the high performance, already achieved on NSTX, to extend into physics regimes much closer to those anticipated in Spherical Tokamak (ST)-based CPPs. These confinement and stability properties will be assessed by a full complement of diagnostics and analysis tools, which will also aid in the development of the underlying theory and predictive models needed for further optimization. Both conventional and transformative heat flux mitigation methods, such as liquid lithium plasma-facing components, will be developed and tested in-situ in NSTX-U at incident heat fluxes of ~100 MW/m 2 , and will inform plans and reduce risk for a subsequent major upgrade to the device to fully heated, high-Z wall and full liquid lithium divertor capability, a technology that potentially could then be implemented on any magnetic confinement device at any aspect ratio. NSTX-U research is fully complementary to programs performed on other STs, nationally and internationally. Furthermore, NSTX-U research has a direct connection to the private sector by informing design choices for future power production facilities being developed by these companies. The NSTX-U program will operate as a national User Facility, with collaborating researchers, engineers, and graduate students from 19 outside institutions, and open to participation and experiments led by researchers from both public and private entities. The research program will advance workforce development through training of young scientists, engineers, and technicians, and it will also serve for further diagnostic innovation, especially for high heat flux and high-Z wall environments, and implementation of advanced artificial intelligence (AI) for plasma and heat flux control.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Capability in Theory, Modeling, and Validation for a Range of Innovative Fusion Concepts using High-Fidelity Moment-Kinetic Models

A computational modeling capability is created and available to the fusion community to understand and design lower-cost and innovative fusion concepts. The approach uses high- fidelity kinetic, moment-kinetic, and moment models and includes sophisticated plasma- boundary interactions. A majority of fusion-relevant simulations are performed with magnetohydrodynamic models and hybrid particle-in-cell codes, with limited-fidelity electron and kinetic physics. However, in fusion configurations like Z-pinches, field-reversed- configurations, plasma jet magneto-inertial fusion, spinning mirrors, and others, kinetic effects (both electron and ions) are critical to understand the physics and design scaling into the highly kinetic regime of a burning fusion plasma. Furthermore, as present fusion machines move towards a burning plasma regime, liquid-metal blankets are needed to handle first-wall heat- flux, reduce erosion, and eventually for energy conversion and fuel breeding. The work performed under this ARPA-E BETHE Capability Team advances the state-of-the-art in modeling and understanding plasma dynamics in fusion devices and its coupling with liquid-metal dynamics. These are critical areas of research for fusion energy to become realizable. To address these complex problems, we have leveraged and extended computational capabilities through the code, Gkeyll (developed jointly with Princeton Plasma Physics Laboratory and academic partners), for kinetic and moment modeling of fusion plasmas. The Concept Teams supported by this Capability Team include the Wisconsin High-field Axisymmetric Mirror (WHAM), Centrifugal Mirror Experiment (CFME), Plasma-Jet Magneto- Inertial Fusion (PJMIF), and solid and liquid wall plasma-material interaction studies relevant to a number of fusion concepts including Zap Energy’s Z-pinch. This software is open-source and available to the fusion community as a high-fidelity tool for the design of lower-cost fusion experiments. 3D gyrokinetic simulations of WHAM are now possible for long enough time scales to understand the evolution of interchange instabilities. 3D multi-fluid simulations of CMFE at higher Mach numbers are now possible for detailed design iterations with the goal of stability. The state-of-the-art in understanding shock formation and shock mitigation regimes in merging liners for PJMIF have been furthered by our kinetic simulations. Our novel models and frameworks studying plasma-material interaction by incorporating wall emission for various solid wall materials of relevance to pulsed and steady fusion concepts have advanced the state-of-the-art in our understanding of particle fluxes, heat fluxes, and other quantities at cathodes and anodes. The results from this work may explain discrepancies between experimental and theoretical predictions of achieved current densities in pulsed concepts such as Z-pinches. Another significant contribution of this Capability Team is the development and deployment of a novel experimental platform, LEX (Liquid Electrode eXperiment), at Virginia Tech to understand liquid metal free-surface response to electromagnetic pulses. The novel experiments along with model validation quantified the effect of different materials and sizes of liquid metal droplets on the radiative power balance of fusion plasmas for pulsed concepts. Furthermore, these experiments provided mitigation strategies for violent liquid metal response for high current pulses as would be expected in fusion regimes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

MAS: A versatile Landau-fluid eigenvalue code for plasma stability analysis in general geometry

We have developed a new global eigenvalue code, Multiscale Analysis for plasma Stabilities (MAS), for studying plasma problems with wave toroidal mode number (n) and frequency (ω) in a broad range of interest in general tokamak geometry, based on a five-field Landau-fluid description of thermal plasmas. Beyond keeping the necessary plasma fluid response, we further retain the important kinetic effects including diamagnetic drift, ion finite Larmor radius, finite parallel electric field (E||), ion and electron Landau resonances in a self-consistent and non-perturbative manner without sacrificing the attractive efficiency in computation. The physical capabilities of the code are evaluated and examined in the aspects of both theory and simulation. In theory, the comprehensive Landau-fluid model implemented in MAS can be reduced to the well-known ideal MHD model, electrostatic ion-fluid model, and drift-kinetic model in various limits, which clearly delineates the physics validity regime. In simulation, MAS has been well benchmarked with theory and other gyrokinetic and kinetic-MHD hybrid codes in a manner of adopting the unified physical and numerical framework, which covers the kinetic Alfv\'en wave (KAW), ion sound wave (ISW), low-n kink, high-n ion temperature gradient mode (ITG) and kinetic ballooning mode (KBM). Moreover, MAS is successfully applied to model the Alfv\'en eigenmode (AE) activities in DIII-D discharge #159243, which faithfully captures the frequency sweeping of reversed shear Alfv\'en eigenmode (RSAE), the tunneling damping of toroidal Alfv\'en eigenmode (TAE), as well as the polarization characteristics of kinetic beta-induced Alfv\'en eigenmode (KBAE) and beta-induced Alfv\'en-acoustic eigenmode (BAAE) being consistent with former gyrokinetic theory and simulation. With respect to the key progress contributed to the community, MAS has the advantage of combining rich physics ingredients, realistic global geometry and high computation efficiency together for plasma stability analysis in linear regime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Liquified SO 2 induced solid/cathode electrolyte interphase for lithium ion batteries

Formation of robust solid/cathode electrolyte interphases (S/CEI) is vital for long-term stability and high-performance operation of lithium-ion batteries (LIBs), particularly under high voltage regimes. However, engineering electrochemically stable S/CEIs that effectively suppress interfacial side reactions remains a key challenge. Herein, we introduce a liquefied sulfur dioxide (SO 2 )– ionic liquid complex as a fluorine-free multifunctional electrolyte additive for the first time that significantly improves the formation of sulfate/sulfite-rich S/CEI layers at both graphite and NMC811 interfaces. The unique SO 2 -N coordination with a 1,2,4-triazolide-based ionic liquid enables homogeneous SO 2 dissolution, resulting in controlled SO 2 decomposition during the initial electrochemical cycle. This decomposition yields sulfur-rich interphase species that stabilize the electrolyte-electrode interface, reduce impedance growth, and lessen electrolyte decomposition. Electrochemical tests show significantly improved cycle life, reduced polarization, and increased Coulombic efficiency for both anodes and cathodes. XPS confirms the presence of SO 2 -derived surface species that contribute to interfacial stability. In conclusion, this approach highlights a new direction for interphase engineering using liquefied gas additives and opens pathways for sulfur-based S/CEI chemistry in advanced battery systems.

Graphite↗

Asymptotic-preserving dynamical low-rank method for the stiff nonlinear Boltzmann equation

In kinetic theory, numerically solving the full Boltzmann equation is extremely expensive. This is because the Boltzmann collision operator involves a high-dimensional, nonlinear integral that must be evaluated at each spatial grid point and every time step. The challenge becomes even more pronounced in the fluid (strong collisionality) regime, where the collision operator exhibits strong stiffness, causing explicit time integrators to impose severe stability restrictions. In this paper, we propose addressing this problem through a dynamical low-rank (DLR) approximation. The resulting algorithm requires evaluating the Boltzmann collision operator only r 2 times, where r, the rank of the approximation, is much smaller than the number of spatial grid points. We propose a novel DLR integrator, called the XL integrator, which reduces the number of steps compared to the available alternatives (such as the projector splitting or basis update & Galerkin (BUG) integrator). For a class of problems including the Boltzmann collision operator which enjoys a separation property between physical and velocity space, we further propose a specialized version of the XL integrator, called the sXL integrator. This version requires solving only one differential equation to update the low-rank factors. Furthermore, the proposed low-rank schemes are asymptotic-preserving, meaning they can capture the asymptotic fluid limit in the case of strong collisionality. Our numerical experiments demonstrate the efficiency and accuracy of the proposed methods across a wide range of regimes, from non-stiff (kinetic) to stiff (fluid).

97 MATHEMATICS AND COMPUTING↗

On the stability and stationarity of the Super H-mode combined with an ion transport barrier in the core

The Super H-mode, a regime with high pedestal pressure and stored energy, is explored on DIII-D and combined with an ion transport barrier in the plasma core to increase performance. A significant improvement of ion temperatures and confinement is facilitated by favorable conditions such as high rotational shear and high ion pedestal temperatures. As a result of a rise in density and simultaneous decrease in rotation, the ion transport barrier disappears during the discharge evolution, leading to a transition from a very high confinement state at early times, to a reduced but still high confinement phase. Additionally, in many discharges, a global MHD event consistent with the coupling of a destabilized internal mode to an edge localized mode (ELM) causes a large energy loss and leads to a reorganization of the plasma into a lower temperature, higher density state. Depending on the magnitude of the global MHD event, the plasma edge collisionality can increase significantly and shift the operational boundary from the peeling to the ballooning side, which can be understood as a drop out of the Super H-mode channel into standard H-mode. Hence, in Super H-mode discharges with ion transport barriers, both the improved pedestal height and rotational shear contribute to the high stored energy. Here, at very low levels of rotation, the confinement factor for SH modes is still expected to exceed standard H-mode by 20-30 %. With their overall stationarity and high-performance levels, Super H-mode discharges provide an attractive regime for ITER and may enable a more compact design of future fusion power plants.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Refined spin-wave model and multimagnon bound states in Li 2 CuO 2

Here we report a study of the spin dynamics in the ferromagnetic chain compound Li 2 CuO 2 . Inelastic neutron scattering measurements allow for the spin Hamiltonian to be determined using a J 1 –J 2 XXZ-Heisenberg spin chain model with weak interchain interactions. The primary exchange parameters determined from our data are qualitatively consistent with those of Lorenz et al. [Europhys. Lett. 88, 37002 (2009)], and our data allow for the resolution of additional interchain exchange interactions. We also observe the formation of two- and, potentially, three-magnon bound states. The two-magnon bound state exists only in the magnetically ordered phase of this material, consistent with stabilization by the weak, Ising-like exchange anisotropy of the nearest-neighbor intrachain interaction. In contrast, the potential three-magnon state persists in a finite temperature regime above T N , indicating an unconventional character. Finally, our results establish Li 2 CuO 2 as an experimental platform for the study of exchange anisotropy-stabilized bound states in a ferromagnetic chain.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Climate Drives Modeled Forest Carbon Cycling Resistance and Resilience in the Upper Great Lakes Region, USA

Forests dominate the global terrestrial carbon budget, but their ability to continue doing so in the face of a changing climate is uncertain. A key uncertainty is how forests will respond (resistance) and recover from (resilience) to rising levels of disturbance of varying intensities. This knowledge gap can optimally be addressed by integrating manipulative field experiments with ecophysiological modeling. We used the ED-2.2 model to project carbon fluxes for a northern temperate deciduous forest subjected to a real-world disturbance severity manipulation experiment. ED-2.2 was run for 150 years, starting from near bare ground in 1900, and subjected to three disturbance treatments under an ensemble of climate conditions. Both disturbance severity and climate strongly affected carbon fluxes, and interacted with one another. We then calculated metrics of functional ecosystem stability. Modeled gross primary production exhibited a mean resistance of -0.13, -0.23, and -0.42 (unitless) in response to increasing disturbance severity of 45%, 65%, and 85% mortality, respectively; and resilience of 0.002, 0.004, 0.008 (yr-1) respectively. This pattern held true for other carbon fluxes, meaning that in general simulated GPP declined more initially but recovered more quickly from higher mortality. Notably, however, heterotrophic respiration responded more slowly to disturbance, and its highly variable response was affected by different drivers. This work provides insight into how future conditions might affect the functional stability of mature forests in this region under ongoing climate change and changing disturbance regimes.

54 ENVIRONMENTAL SCIENCES↗

Topological prethermal strong zero modes on superconducting processors

Abstract Symmetry-protected topological phases 1–4 cannot be described by any local order parameter and are beyond the conventional symmetry-breaking model 5 . They are characterized by topological boundary modes that remain stable under symmetry respecting perturbations 1–4,6–8 . In clean, gapped systems without disorder, the stability of these edge modes is restricted to the zero-temperature manifold; at finite temperatures, interactions with mobile thermal excitations lead to their decay 9–11 . Here we report the observation of a distinct type of topological edge mode 12–14 , which is protected by emergent symmetries and persists across the entire spectrum, in an array of 100 programmable superconducting qubits. Through digital quantum simulation of a one-dimensional disorder-free stabilizer Hamiltonian, we observe robust long-lived topological edge modes over up to 30 cycles for a wide range of initial states. We show that the interaction between these edge modes and bulk excitations can be suppressed by dimerizing the stabilizer strength, leading to an emergent U(1) × U(1) symmetry in the prethermal regime of the system. Furthermore, we exploit these topological edge modes as logical qubits and prepare a logical Bell state, which exhibits persistent coherence, despite the system being disorder-free and at finite temperature. Our results establish a viable digital simulation approach 15–18 to experimentally study topological matter at finite temperature and demonstrate a potential route to construct long-lived, robust boundary qubits in disorder-free systems.

Science & Technology - Other Topics↗