Search NASA⌕ Search

SEARCH · Search NASA

Results for “thermal perturbations”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Monte Carlo Perturbation Analysis of Fuel Temperature Variations in the MCNP Model of the Annular Core Research Reactor

The Annular Core Research Reactor (ACRR) Monte Carlo N-Particle (MCNP) model is used by ACRR reactor operators and experiment designers at Sandia National Laboratories for a variety of computational calculations ranging from reactor kinetics parameter estimates and safety analyses to experimental planning. To understand the dominant source of uncertainty within the MCNP model, perturbations in temperature were applied to individual ACRR MCNP fuel rods. Fuel rod temperatures were randomly sampled from a uniform distribution from operational temperatures to quantify temperature-related uncertainty effects. Stochastic mixing was used to blend the cross sections of the desired temperatures using the MCNP continuous and Thermal Neutron Scattering Treatment [S(α,β)] libraries in ENDF/B-VII.1. Furthermore, this uncertainty analysis produced a 640 row × 640 column correlation and covariance matrix of the neutron energy spectra. Positive covariance was produced around the 1-MeV region and the 0.2-eV region. Correlation was found in the thermal and fast energy regions, but no correlation was observed in the slowing-down energy region because interactions in this region are not dominated by fuel.

ACRR↗
DOE OSTI2022

Do ambient shear and thermal stratification impact wind turbine tip-vortex breakdown?

Modern wind turbines experience uneven inflow conditions across the rotor, due to the ambient flow’s shear and thermal stratification. Such conditions alter the shape and length of turbine wakes and thus impact the loads and power generation of downstream turbines. To this end, understanding the spatial evolution of the individual wakes under different atmospheric conditions is key to controlling and optimising turbine arrays. With this numerical study we aim to obtain a better understanding of the fundamental physics governing the near-wake dynamics of wind turbines under shear and thermal stability, by examining their tip-vortex breakup mechanisms. Our approach considers scale-resolving simulations of a single turbine wake under a linear shear profile as well as the application of harmonic tip perturbations to trigger flow instabilities. For the subsequent analysis we use the proper orthogonal decomposition (POD) method to extract coherent structures from the flow, and we also calculate mean kinetic energy fluxes to quantify each coherent structure’s contribution to wake recovery. The wake’s helical spiral is found to hinder wake recovery for all studied ambient flow conditions, whereas the mutual inductance instability has positive MKE flux leading to an enhanced wake recovery. Finally, the ambient shear has the largest impact on the local MKE flux with respect to downstream location by changing the shape of the curve and location of extrema, whereas thermal stratification has only a minimal impact on the magnitude of the near-wake local MKE flux distribution.

17 WIND ENERGY↗
DOE OSTI2020

Measured Thermal and Fast Neutrons Fluence Rates for ATF-1 Holders During ATR Cycle 166B

This report contains the thermal (2200 m/s) and fast (E>lMeV) neutron fluence rate data for the ATF-1 holders located in core for ATR Cycle 166B which were measured by the Radiation Measurements Laboratory (RML) as requested by the Power Reactor Programs (ATR Experiments) Radiation Measurements Work Order. This report contains measurements of the fluence rates corresponding to the particular elevations relative to the 80-ft. core elevation. The data in this report consist of (1) a table of the ATR power history and distribution, (2) a hard copy listing of all thermal and fast neutron fluence rates, and (3) plots of both the thermal and fast neutron fluence rates. All "BR" holder monitor wires for this cycle are 56.375 inches long. This length allows measurements from 31.17 inches above core midplane to 19.81 inches below core midplane. The distance from the end of the wires to the first count position was 4.25 inches for all wires counted from this cycle. The monitors were located in an outward small I-holes,1-22 and 1-23. Thermal neutron fluence was determined based on the activation of the CoAl wires and resonance corrections determined in the irradiation of Cycle 152A when the RML measured cadmium-covered as well as bare neutron monitors in the Southwest, Southeast, and H positions at the Center lobe to determine the division in neutron energy between epithermal and thermal neutrons. For more information please refer to letter “Radiation Measurements Laboratory measurements of in core ATR physics testing during cycle 152A. The measured values documented in RML procedure, “ACMM-3600, Flux Monitoring,” will be used indefinitely until further measurements supersede the current values. Fast neutron fluence was determined using the activation of nickel wires. Typically, a fission spectrum averaged cross section (E>1 MeV) is used as determined in detailed measurements described in “ATR Neutron Spectral Characterization,” which gives an effective cross section of 97 mb for the small I-holes. However, those measurements were performed with unfueled experiments. The fast fluence values reported here use a perturbed 114 mb effective>1 MeV spectrum averaged cross section for the small I-holes based on MCNP simulations for a previous fueled experiment in the1-24 position. Refer to “LWR- ATR Irradiation Measured Thermal and Fast Neutron Fluence Rates ATR CYCLE 133A 04/19/04 thru 05/17/04” for more information.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Uncertainty quantification and sensitivity analysis of a nuclear thermal propulsion reactor startup sequence

The research presented in this article describes progress in applying stochastic methods, uncertainty quantification, parametric studies, and variance-based sensitivity analysis (also known as Sobol sensitivity analysis) to a full-core model of a nuclear thermal propulsion (NTP) system simulated via the radiation transport code Griffin to simulate neutronics. Our goal is to develop a reduced-order (surrogate) model that can be rapidly sampled with perturbations to multiple input parameters. In this NTP system, reactivity and power feedback affect the rotation of control drums (CDs), which is itself controlled by a hybrid proportional-integral-derivative (PID) controller actuated by the power demand and reactivity feedback from the numerical model. This model uses reactor kinetic feedback (mean generation time [Λ] and effective delayed neutron fraction [ β eff ] from a transient Griffin simulation executed via Griffin’s improved quasi-static solver to provide the kinetic parameters) as inputs to functions that control the CD rotation angle. By investigating numerous stochastic approaches, we developed a dual-purpose surrogate model of the NTP system, using polynomial regression in the Multiphysics Object-Oriented Simulation Environment (MOOSE) Stochastic Tools Module (STM). The trained model can be rapidly sampled while simultaneously perturbing various input parameters, such as coefficients on the PID control or temperature (directly affecting the neutron cross section). The surrogate model delivers accurate (within 5%) results at speeds orders of magnitude faster (minutes, not days of computational time) than the base model. Once the surrogate model has been trained, distributions of the uncertain parameters can be changed at will to investigate the effects of perturbing multiple inputs as well as the effects of these inputs on the model output. For example, coefficients used in the PID control system may vary due to some type of physical interference, or uncertainty may exist in the temperature of the neutron cross sections in various regions of the reactor. A distribution can be placed on these parameters, and operational boundaries can be determined. The goal of this work is to support development of an advanced control system for operating CDs in a functioning NTP system. This work is a scoping study of the MOOSE STM.

21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN↗
DOE OSTI2025

Fast ion stabilization of tilt in large radius FRCs

The field reversed configuration (FRC) has been a curious case in plasma physics research in that early MHD analysis suggested FRCs should be grossly unstable, while experimental results contradicted that prediction. Later, this theory was able to resolve this contradiction by understanding that finite Larmor radius effects largely negated the MHD predictions. Similarly, previous theoretical studies of beam driven FRCs predicted that such system would be unstable to beam driven modes while, again, experimental results indicated the contradiction. In this paper, we reconcile the theoretical understanding of beam driven modes with experimental observations of stability in these systems. By self-consistently capturing fast ion generation from neutral beam injection and its impact on the plasma equilibrium, we show that low amplitude perturbations in the magnetic field, driven by betatron particles, modify the precession frequencies of the betatron particles such that the drive for compressional Alfvén waves in the thermal plasma is reduced. Finally, we are able to demonstrate, for the first time, stable beam driven FRC evolution at high S*/E in 3D kinetic simulations.

Betatrons↗
DOE OSTI2023

Finite-temperature many-body perturbation theory for anharmonic vibrations: Recursions, algebraic reduction, second-quantized reduction, diagrammatic rules, linked-diagram theorem, finite-temperature self-consistent field, and general-order algorithm

A unified theory is presented for finite-temperature many-body perturbation expansions of the anharmonic vibrational contributions to thermodynamic functions, i.e., the free energy, internal energy, and entropy. The theory is diagrammatically size-consistent at any order, as ensured by the linked-diagram theorem proved in this study, and, thus, applicable to molecular gases and solids on an equal footing. It is also a basis-set-free formalism, just like its underlying Bose–Einstein theory, capable of summing anharmonic effects over an infinite number of states analytically. It is formulated by the Rayleigh–Schrödinger-style recursions, generating sum-over-states formulas for the perturbation series, which unambiguously converges at the finite-temperature vibrational full-configuration-interaction limits. Two strategies are introduced to reduce these sum-over-states formulas into compact sum-over-modes analytical formulas. One is a purely algebraic method that factorizes each many-mode thermal average into a product of one-mode thermal averages, which are then evaluated by the thermal Born–Huang rules. Canonical forms of these rules are proposed, dramatically expediting the reduction process. The other is finite-temperature normal-ordered second quantization, which is fully developed in this study, including a proof of thermal Wick’s theorem and the derivation of a normal-ordered vibrational Hamiltonian at finite temperature. The latter naturally defines a finite-temperature extension of size-extensive vibrational self-consistent field theory. These reduced formulas can be represented graphically as Feynman diagrams with resolvent lines, which include anomalous and renormalization diagrams. Two order-by-order and one general-order algorithms of computing these perturbation corrections are implemented and applied up to the eighth order. The results show no signs of Kohn–Luttinger-type nonconvergence.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Scaling functions for eigenstate entanglement crossovers in harmonic lattices

For quantum matter, eigenstate entanglement entropies obey an area law or log-area law at low energies and small subsystem sizes and cross over to volume laws for high energies and large subsystems. This transition is captured by crossover functions, which assume a universal scaling form in quantum critical regimes. We demonstrate this for the harmonic lattice model, which describes quantized lattice vibrations and is a regularization for free scalar field theories, modeling, e.g., spin-0 bosonic particles. In one dimension, the ground-state entanglement obeys a log-area law. For dimensions d ≥ 2, it displays area laws, even at criticality. The distribution of excited-state entanglement entropies is found to be sharply peaked around subsystem entropies of corresponding thermodynamic ensembles in accordance with the eigenstate thermalization hypothesis. Numerically, we determine crossover scaling functions for the quantum critical regime of the model and do a large-deviation analysis. We show how infrared singularities of the system can be handled and how to access the thermodynamic limit using a perturbative trick for the covariance matrix. Eigenstates for quasi-free bosonic systems are not Gaussian. Here, we resolve this problem by considering appropriate squeezed states instead. For these, entanglement entropies can be evaluated efficiently.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

exoALMA. V. Gaseous Emission Surfaces and Temperature Structures

An analysis of the gaseous component in protoplanetary disks can inform us about their thermal and physical structure, chemical composition, and kinematic properties, all of which are crucial for understanding various processes within the disks. By exploiting the asymmetry of the line emission, or via line profile analysis, we can locate the emitting surfaces. Here, we present the emission surfaces of the exoALMA sources in 12 CO J = 3–2, 13 CO J = 3–2, and CS J = 7–6. We find that 12 CO traces the upper disk atmosphere, with mean $\langle$ z/r $\rangle$ values of ≈0.28, while 13 CO and CS trace lower regions of the disk with mean $\langle$ z/r $\rangle$ values of ≈0.16 and ≈0.18, respectively. We find that 12 CO $\langle$ z/r $\rangle$ and the disk mass are positively correlated with each other; this relationship offers a straightforward way to infer the disk mass. We derive 2D r – z temperature distributions of the disks. Additionally, we search for substructure in the surfaces and radial intensity profiles; we find evidence of localized substructure in the emission surfaces and peak intensity profiles of nearly every disk, with this substructure often being coincident between molecular tracers, intensity profiles, and kinematic perturbations. Four disks display evidence of potential photodesorption, implying that this effect may be common even in low far-ultraviolet star-forming regions. For most disks, we find that the physical and thermal structure is more complex than analytical models can account for, highlighting a need for more theoretical work and a better understanding of the role of projection effects on our observations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Kinetic Landau-fluid closures of non-Maxwellian distributions

New kinetic Landau-fluid closures, based on the cutoff Maxwellian distribution, are derived. A special static case is considered (the frequency ω=0). In the strongly collisional regime, our model reduces to Braginskii's heat flux model, and the transport is local. In the weak collisional regime, our model indicates that the heat flux is non-local and recovers the Hammett–Perkins model while the value of the cutoff velocity approaches to infinity. We compare the thermal transport coefficient χ of Maxwellian, cutoff Maxwellian and super-Gaussian distribution. The results show that the reduction of the high-speed tail particles leads to the corresponding reduction of the thermal transport coefficient χ across the entire range of collisionality, more reduction of the free streaming transport toward the weak collisional regime. In the collisionless limit, χ approaches to zero for the cutoff Maxwellian and the super-Gaussian distribution but remains finite for Maxwellian distribution. χ is complex if the cutoff Maxwellian distribution is asymmetric. The Im(χ) approaches to different convergent values in both collisionless and strongly collisional limit, respectively. It yields an additional streaming heat flux in comparison with the symmetric cutoff Maxwellian distribution. Furthermore, due to the asymmetric distribution, there is a background heat flux q0 though there is no perturbation. Furthermore, the derived Landau-fluid closures are general for fluid moment models, and applicable for the cutoff Maxwellian distribution in an open magnetic field line region, such as the scape-off-layer of Tokamak plasmas, in the thermal quench plasmas during a tokamak disruption, and the super-Gaussian electron distribution function due to inverse bremsstrahlung heating in laser-plasma studies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Chiral kinematic theory and converse vortical effects

Response theories in condensed matter typically describe the response of an electron fluid to external electromagnetic fields, while perturbations on neutral particles are often designed to mimic such fields. Here, we study the response of fermions to a space-time-dependent velocity field, thereby sidestepping the issue of gauge charge. First, we use a semiclassical chiral kinematic theory to obtain the local density of current and extract the orbital magnetization. The theory immediately predicts a "converse vortical effect," defined as an orbital magnetization driven by linear velocity. It receives contributions from magnetic moments on the Fermi surface and the Berry curvature of the occupied bands. Then, transcending semiclassics via a complementary Kubo formalism reveals that the uniform limit of a clean system receives only the Berry curvature contribution while other limits sense the Fermi surface magnetic moments too. We propose CoSi as a candidate material and suggest magnetometry of a sample under a thermal gradient to detect the effect. Overall, our study sheds light on the effects of a space-time-dependent velocity field on electron fluids and paves the way for exploring quantum materials using new probes and perturbations.

Chen, Kai↗
DOE OSTI2024

Generation of shear flows induced by AE / EPM in LHD plasma

The generation of shear flows (SFs) by Alfven Eigenmodes (AEs) and energetic particle modes (EPMs) have important effects on the operation of future nuclear fusion reactors, because SFs regulate the saturation of the AEs/EPMs, the transport of EPs and thermal plasma, as well as the formation of transport barriers among other consequences. The aim of this study is the analysis of SFs generation during the saturation phase of AEs and EPMs in LHD plasma. Experiments performed in the 23rd and 24th LHD experimental campaigns are dedicated to explore the destabilization of AEs/EPMs in discharges with different heating patterns, thermal plasma and magnetic field configurations. In particular, the shots 176490 and 179697 show the destabilization of MHD bursts and energetic-ion-driven resistive interchange modes (EIC), respectively. Charge exchange spectroscopy measurements in both discharges indicate that the generation of SFs by AE/EPM is uncorrelated with the perturbation induced by the neutral beam injector (NBI). Nonlinear simulations performed using the gyro-fluid code FAR3d show the generation of zonal structures, especially SFs, induced during the saturation phase of Toroidal Alfven Eigenmodes (TAEs) triggered in the MHD burst as well as by the 1/1 EIC in the bursting phase. The simulations indicate that SFs are caused by the radial electric fields powered by energy transfers from the unstable AE/EPM towards the thermal plasma. The strongest SFs are measured during the EIC bursting phase once the 1/1 EPM overlaps with nearby resonances at the plasma periphery. Likewise, the largest SFs during the MHD burst are observed once TAEs radially overlap in the inner-middle plasma region.

AE↗
DOE OSTI2021

Magnetic Field Perturbations to a Soft X-ray-Activated Fe (II) Molecular Spin State Transition

The X-ray-induced spin crossover transition of an Fe (II) molecular thin film in the presence and absence of a magnetic field has been investigated. The thermal activation energy barrier in the soft X-ray activation of the spin crossover transition for [Fe{H2B(pz)2}2(bipy)] molecular thin films is reduced in the presence of an applied magnetic field, as measured through X-ray absorption spectroscopy at various temperatures. The influence of a 1.8 T magnetic field is sufficient to cause deviations from the expected exponential spin state transition behavior which is measured in the field free case. We find that orbital moment diminishes with increasing temperature, relative to the spin moment in the vicinity of room temperature.

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

Direct Recycling of Spent NCM Cathodes through Ionothermal Lithiation

Ionic liquids (ILs) are a family of nonconventional molten salts that offer many advantages, such as negligible vapor pressures, negligible flammability, wide liquidus ranges, good thermal stability, and much synthesis flexibility. The unique solvation environment of these ILs provides new reaction or flux media for controlling formation of solid-state materials with a minimum perturbation of morphologies. A successful lithiation via ionothermal synthesis using a cost-effective Li halide as Li source and recyclable ILs as solvents is reported in this study for the direct recycling of LiNi 1/3 Co 1/3 Mn 1/3 O 2 (NCM 111) cathodes. In addition, the ionic liquids can be readily recycled and reused after ionothermal lithiation. The lithiation of spent cathodes can enable the direct recycling of spent cathode materials in lithium-ion batteries.

25 ENERGY STORAGE↗
DOE OSTI2024

North‐East Peri‐Tethyan Water Column Deoxygenation and Euxinia at the Paleocene Eocene Thermal Maximum

Abstract The Paleocene–Eocene Thermal Maximum (PETM) is associated with climatic change and biological turnover. It shares features with the Oceanic Anoxic Events (OAEs) of the Mesozoic, such as transient global warming and biogeochemical perturbations. However, the PETM experienced a more muted expansion of marine anoxia compared to the Mesozoic OAEs (especially OAE2), with geographically limited evidence for photic zone euxinia (PZE). We explore the extent and drivers of marine deoxygenation during the PETM using biomarkers for water column euxinia and anoxia as well as an intermediate complexity Earth system model (cGEnIE). These reveal that the water column in the North‐East Peri‐Tethys became anoxic, with euxinic conditions reaching the photic zone (PZE) during the PETM. Our model shows that euxinia developed due to a global increase in the ocean nutrient inventory with concomitant oxygen consumption, similar to findings for OAE2. The particularly strong regional response in the NE Peri‐Tethys appears to arise from a combination of global CO 2 ‐weathering forcing, regionally restricted circulation and upwelling of sulphidic thermocline waters. Unlike OAE2, anoxia and PZE do not become widespread in our PETM simulations, consistent with new and existing geochemical and biological proxy data. This globally muted response could result from reduced biogeochemical feedbacks to climate forcing relative to the mid‐Cretaceous climate. Our observations suggest that similar mechanisms operated in response to disparate Cenozoic (PETM) and Mesozoic (OAEs) transient global warming events, while also highlighting that background conditions are crucial in modulating the sensitivity of Earth's system to them.

Behrooz, L.↗
DOE OSTI2020

Wave function methods for canonical ensemble thermal averages in correlated many-fermion systems

We present a wave function representation for the canonical ensemble thermal density matrix by projecting the thermofield double state against the desired number of particles. Furthermore, the resulting canonical thermal state obeys an imaginary-time evolution equation. Starting with the mean-field approximation, where the canonical thermal state becomes an antisymmetrized geminal power (AGP) wave function, we explore two different schemes to add correlation: by number-projecting a correlated grand-canonical thermal state and by adding correlation to the number-projected mean-field state. As benchmark examples, we use number-projected configuration interaction and an AGP-based perturbation theory to study the hydrogen molecule in a minimal basis and the six-site Hubbard model.

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

Economic Dispatch Optimization of Multi-Unit SMR Site

With the increase of integration of variable renewable into the grid, grid resilience can be compromised without flexible electricity suppliers. Flexible operation of NPP can expand the scope of integration of new nuclear power plants in regions with high share of renewables. It also provides opportunity to increase profits through participation in ancillary services. In this research, NuScale modules are used as the case study. NuScale is an integrated pressurized water reactor with thermal power of 250 MWth. Having multiple units can provide higher flexibility in changing power levels and supplying the grid with electricity all year long. Additionally, load-following might perturbate the ordered sequence of refueling (every 24 months for NuScale). To maintain the 2 year refueling cycles, the loading patterns will need to be adjusted in a year-to-year basis, resulting in additional unnecessary costs. Another solution is to find the optimal outage schedule. This will also shift the focus of the planned outage to a more optimized predictive outages and maintenance. Results from short-term optimization (i.e., daily participation in energy and ancillary services markets) shows that average commitment to regulation down ancillary service markets during late night and early morning hours can increase revenues. Participating in upward reserves during pre-work commuting hours and early evening hours can also generate additional revenues due to high ramp up of electricity demand. Having multiple units at a single site can also help operate flexibly. Results from long-term optimization show that daily fluctuation in electric power is usually handled by few units while other units operate on base-load profile. It is possible to combine base-load operation with flexible operation to extend component remaining useful life in some units and find optimal outage schedule for the units.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2021

Statistical emulation of a perturbed basal melt ensemble of an ice sheet model to better quantify Antarctic sea level rise uncertainties

Abstract. Antarctic ice shelves are vulnerable to warming ocean temperatures, and some have already begun thinning in response to increased basal melt rates. Sea level is therefore expected to rise due to Antarctic contributions, but uncertainties in its amount and timing remain largely unquantified. In particular, there is substantial uncertainty in future basal melt rates arising from multi-model differences in thermal forcing and how melt rates depend on that thermal forcing. To facilitate uncertainty quantification in sea level rise projections, we build, validate, and demonstrate projections from a computationally efficient statistical emulator of a high-resolution (4 km) Antarctic ice sheet model, the Community Ice Sheet Model version 2.1. The emulator is trained to a large (500-member) ensemble of 200-year-long 4 km resolution transient ice sheet simulations, whereby regional basal melt rates are perturbed by idealized (yet physically informed) trajectories. The main advantage of our emulation approach is that by sampling a wide range of possible basal melt trajectories, the emulator can be used to (1) produce probabilistic sea level rise projections over much larger Monte Carlo ensembles than are possible by direct numerical simulation alone, thereby providing better statistical characterization of uncertainties, and (2) predict the simulated ice sheet response under differing assumptions about basal melt characteristics as new oceanographic studies are published, without having to run additional numerical ice sheet simulations. As a proof of concept, we propagate uncertainties about future basal melt rate trajectories, derived from regional ocean models, to generate probabilistic sea level rise estimates for 100 and 200 years into the future.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

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↗