Search NASA⌕ Search

SEARCH · Search NASA

Results for “strongly correlated systems”

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.

1,408 records · Page 9
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Boundary Layer Analysis of Shock Tube Flows

Shock tubes offer a controlled environment to reproduce kinetic and radiative phenomena characteristic of atmospheric entry flows under ground-test conditions. The boundary layer developing behind the incident shock wave determines the available test time, influences particle residence times important for similarity scaling, and can affect radiative energy transport. In this study, we couple a quasi-1D space marcher with the compressible boundary-layer equations to numerically compute the post-shock flow in a shock frame of reference for various test gas mixtures representative of different planetary atmospheres. The solvers are individually verified against CFD simulations and analytical correlations available in the literature. Coupled solutions are computed for a finite-rate chemistry in the boundary layer and a non-catalytic isothermal wall. Results yield refined estimates of the maximum separation distance, as well as insights into concentration profiles of relevant species within the boundary layer.

Andrea Fagnani↗
NASA NTRS2026

Micromechanical Design of Carbon Nanotube Ribbon Reinforced Polymer Composite Materials

Lightweight materials are an important component of the design of aerospace structures. Carbon nanotube materials have been considered for this purpose due to the strength and stiffness of individual nanotubes, and the commercial availability of bulk formats such as fibers. These fibers can have a ribbon cross section which results in a different design space for their composites relative to traditional reinforcements which have a round cross section. This work applies brick-and-mortar micromechanical models and classical lamination theory with an inverse approach to investigate the design space of these composites. Using this approach, the influences of fiber geometry and axial and transverse mechanical properties are mapped. Finally, a sensitivity study is performed and the relative impacts of ±10% variations in the constituent material and geometric properties are ranked. Lamina axial moduli were found to range from a maximum of 3x to 1x minimum relative to a target quasi-isotropic laminate modulus depending on the anisotropy and shear modulus of the lamina. The fiber targets depended strongly on the fiber volume fraction in the lamina and the fiber axial modulus target was found to range from 4.8x to 3.2x the quasi-isotropic laminate target. The sensitivity analysis found that the largest driver of performance was the volume fraction, followed by the fiber axial modulus. While bio-based brick-and-mortar composites, such as nacre, can benefit from reinforcement aspect ratios above 10, for carbon nanotube ribbon (or carbon fiber)/polymer composites the sensitivity study indicated that the optimal cross-sectional aspect ratio was relatively smaller, potentially less than three.

Micromechanics↗
DOE OSTI2023

2,5‐Dimercapto‐1,3,4‐Thiadiazole (DMCT)‐Based Polymers for Rechargeable Metal–Sulfur Batteries

Organosulfur materials are a sustainable alternative to the present-day layered oxide cathodes in lithium-based batteries. One such organosulfur material that was intensely explored from the 1990s to early 2010s is 2,5-dimercapto-1,3,4-thiadiazole (DMCT). However, research interest declined as the electrode reactions with DMCT were assumed to be too sluggish to be practical. Armed with the advances in metal–sulfur batteries, we revisit DMCT-based materials in the form of poly[tetrathio-2,5-(1,3,4-thiadiazole)], referred to as pDMCT-S. With an appropriate choice of electrode design and electrolyte, pDMCT-S cathode paired with a Li-metal anode shows a capacity of 715 mA h g −1 and a Coulombic efficiency of 97.7% at a C/10 rate, thus quelling the concerns of sluggish reactions. Surprisingly, pDMCT-S shows significantly improved long-term cyclability compared to a sulfur cathode. Investigations into the origin of the stability reveals that the discharge product Li-DMCT in its mesomeric form can strongly bind to polysulfides, preventing their dissolution into the electrolyte and shuttling. This unique mechanism solves a critical problem faced by sulfur cathodes. Encouragingly, this mechanism results in a stable performance of pDMCT-S with Na-metal cells as well. In conclusion, this study opens the potential for exploring other organic materials that have inherent polysulfide sequestering capabilities, enabling long-life metal–sulfur batteries.

2,5-dimercapto-1,3,4-thiadiazole↗
DOE OSTI2026

Healable Coatings as a Mechanism to Repair Leading Edge Erosion in Wind Energy

Wind turbine blades are highly engineered structures designed to face temperature extremes and high winds. However, erosion of the blade's leading edge and subsequent repair remains a significant and costly challenge for the wind energy industry. Repair of these leading edges can lead to large amounts of downtime for the turbine and significant operational inefficiencies. In this work, the strength of adhesion and healing ability of a commercially available vitrimer (Mallinda's VITRIMAX) was compared to that of a thermoplastic resin, which has previously been demonstrated in wind energy applications (Arkema's Elium) to evaluate their efficacy as surface coatings for wind turbine blades, particularly their leading edges. Vitrimers are a class of inherently reprocessable thermosets, and it was theorized that vitrimer-based leading edge coatings could enable more robust and efficient wind turbine blades with decreased operational downtime and safer maintenance practices. It was found that the VITRIMAX adhered better to the wind blades' surfaces than both the manufacturer's paint and Elium, with increases in pull-off strength of adhesion ranging from 24% to 83% above that of the original paint. Furthermore, the VITRIMAX adhered strongly to the underlying composite of each blade with strength of adhesion values increasing in ranges from 42% to 97% above that of the original paint. Finally, the vitrimer coating showed an 88% decrease in surface roughness compared to end-of-life blade materials, and initial healing demonstrations in which coatings were manually scratched and subsequently healed exhibited an ~84.5% decrease in scratch depths.

Hubbard, Amber [ORNL]↗
NASA NTRSDate not supplied

Performance Assessment of LunaNet’s Augmented Forward Signal

LunaNet provides a common set of interoperable specifications for communication and position, navigation and time (PNT) services and interfaces soon to be implemented in lunar vicinity. The LunaNet Interoperability Specification (LNIS) provides the design for the GNSS-like Augmented Forward Signal (AFS), which enables orbiting and surface users in lunar space, such as Artemis, to estimate their position, velocity, and time. The specification of AFS defines two orthogonal signal components on a single carrier: the in-phase component (AFS-I), a lower-chip-rate data channel tailored for applications where low SWaP (Size, Weight, and Power) is critical (e.g., IoT devices or search and rescue), and the quadrature component (AFS-Q), a high-chip-rate data-less pilot signal for high-precision, robust lunar navigation and positioning applications. An initial description of AFS was provided in LNIS 2023, with initial analysis results shown in Dafesh 2024 and Dafesh 2025, and the current signal in space description provided in LNIS 2025. As part of NASA's Lunar Communication Relay and Navigation Systems (LCRNS) project, this work expands upon the initial analysis results and proposes a new expanded set of AFS-Q spreading codes that exceed the cross-correlation and autocorrelation sidelobe performance of L1C and other GNSS signals, while providing additional expansion capabilities for future service satellites. A set of 420 codes was selected from a Weil-based code derived from the prime number 10247, which is larger than the 10243 prime number used to derive BeiDou’s B1C Weil sequences. Both the initial set of 210 codes and the expanded set of 420 codes are shown to provide the best cross-correlation of any 10230-chip satellite navigation codes. The performance is demonstrated for hierarchical sets of spreading codes optimized and organized in sets of 30 codes. The work also compares LunaNet’s AFS to terrestrial GNSS signals in terms of acquisition, tracking, and data demodulation performance. Performance is evaluated for receivers that only track the 1.023 MCPS data channel spreading code for low SWaP IoT use cases, as well as for receivers that track both the 1.023 MCPS data channel and the 5.115 MCPS pilot channel spreading code for high-performance use cases. Performance is assessed in the presence of interference and thermal noise. The analysis is performed in terms of expected operating conditions on the lunar surface. Several unique flexibility aspects of the augmented forward signal are described, including the use of the Q channel’s secondary and tertiary codes to enable variable coherent integrations during acquisition. This is compared to GNSS signals such as L5/E5 and MBOC in terms of achievable processing gain for interference mitigation versus acquisition complexity. The work details acquisition and tracking techniques used to optimally acquire and track the primary, secondary, and tertiary codes on the Q channel, as well as acquisition of the I channel spreading code. Acquisition of the 8 ms, Q channel spreading code is also compared to joint acquisition of the I and Q channel primary codes in noise and interference environments

LCRNS↗
DOE OSTI2026

Process control-enabled mitigation of microstructural and plastic heterogeneities in additively manufactured Grade 91 steel

Synergizing wire arc-directed energy deposition (WA-DED) additive manufacturing (AM) with particle-strengthened creep strength-enhanced ferritic (CSEF) steels enables fabrication and repair of critical power-plant components. Investigations focused on fusion-welded particle-strengthened CSEF steels, such as Grade 91 steel, have linked microstructurally heterogeneous regions—forming due to heat affected zones (HAZ)—with premature failure during elevated temperature service. Fusion-based AM, including WA-DED, likewise generates microstructurally and plastically heterogeneous regions due to spatiotemporally varying thermokinetics during deposition. However, works investigating such microstructural heterogeneities, their implications for mechanical behavior, and strategies to mitigate their formation remain scarce. This work identifies microstructurally and plastically heterogeneous regions within the WA-DED-processed Grade 91 steel. Spatial microhardness variations in the as-fabricated specimen correlate with the variation in the attributes of grain, martensitic microstructure, and precipitates across the fusion zone and HAZ. Digital image correlation-enabled tensile tests performed at 500 °C revealed pronounced deformation localization and a wave-like strain distribution, with wavelength close to the melt pool depth, indicating susceptibility of the as-fabricated components to premature creep failure. Such heterogeneity in microstructural and mechanical behavior was attributed to recurring solid-state phase transformations. Subsequently, an interlayer temperature control strategy was implemented, wherein maintaining interlayer temperature above the martensitic start temperature mitigated the heterogeneous microstructural and plastic response in the as-fabricated condition. Findings open pathways to achieving deformation-localization- and creep-resistant microstructures in WA-DED fabricated particle-strengthened CSEF steel components, reducing reliance on post-welding heat treatments—conventionally required to enhance creep resistance—and enabling on-demand, short lead-time fabrication of next-generation power-plant components.

Heat affected zones↗
DOE OSTI2026

SAVY-4000 Finite-Element Drop Test Analysis

PFE Auxiliary Systems conducted drop testing on SAVY-4000 containers to evaluate structural response under 12-foot drop conditions. In support of that effort, a finite-element modeling capability was developed to simulate drop response across multiple container sizes and impact orientations. The purpose of this work was to provide a consistent analysis framework that could support interpretation of testing, compare response trends across multiple configurations, and generate quantities of interest for later comparison with experimental data. More broadly, the analysis and testing were intended to assess whether the containers continued to perform their primary function after a 12-foot drop, namely maintaining structural integrity and containment of the contents. The modeling approach combined an implicit preload analysis with an explicit drop simulation so that each drop event began from a mechanically realistic assembled condition, including compression of the silicone O-ring. Separate models were developed for 2-quart, 5-quart, 12-quart, and 10-gallon containers. The results were evaluated in terms of strain-gauge response, collar-lid gap behavior, and accumulated plastic strain. In addition, parametric studies were performed on the 2-quart container to assess sensitivity to O-ring stiffness, friction, canister thickness, geometry tolerance, and mesh density. The simulations showed that predicted drop responses depended strongly on both container size and drop orientation. Gap metrics identified cases in which the predicted collar-lid opening exceeded the nominal O-ring cross-section threshold, while plastic strain metrics identified localized regions of elevated permanent deformation. Parametric studies showed that the predicted response was especially sensitive to the assumed O-ring stiffness and contact friction, while the geometry tolerance study produced smaller changes in the cases examined. The main value of this work was that it established a repeatable modeling and simulation workflow to support drop-test implementation, evaluate effects of future configuration changes, and understand modeling assumptions that most influenced predicted response. At the current stage, the results were viewed as preliminary model predictions rather than validated predictions. The next step would be to compare drop-test data to the model so that predictive values of the workflow could be refined and used with greater confidence to assess whether the containers maintained structural integrity and containment of the contents after a 12-foot drop.

42 ENGINEERING↗
DOE OSTI2026

Block-Type Antiferromagnetism in Single Chain Quasi-One-Dimensional K 3 ⁢Fe 2 ⁢Se 4

One-dimensional (1D) structures provide a unique platform to study the correlated quantum interactions and phase transitions such as unconventional magnetism and superconducting states. Here, we report that iron chalcogenide K 3 ⁢Fe 2 ⁢Se 4 exhibits an unusual block-type canted antiferromagnetic (AFM) order with a clear single chain quasi-1D structure, which is structurally different from the two-leg ladder BaFe 2 ⁢Se 3 , through both experimental measurements and density matrix renormalization group (DMRG) calculations. The narrow bandgap semiconductor K 3 ⁢Fe 2 ⁢Se 4 has a quasi-1D edge-shared FeSe4 tetrahedra chain structure and orders antiferromagnetically below 110 K. The magnetic moments couple antiferromagnetically along the quasi-1D chain direction of the 𝑏 axis and form an up-down-down-up (↑−↓−↓−↑)–like spin structure with a commensurate propagation vector 𝒌=⁢(0,0,0), where block-type spin ↑−↑ or ↓−↓ coupling are between the longer Fe-Fe bonds of the quasi-1D chain. DMRG results show that block antiferromagnetic state is stable in K 3 ⁢Fe 2 ⁢Se 4 and reveal that the block-ordered arrangement of Fe 2.5+ ions spins arise from the competition between ferromagnetic and AFM interaction in the presence of strong electronic correlation. Our research results not only report the discovery of a clear block-type canted antiferromagnetic structure in a real quasi-1D chain material but also provide a theoretical approach to understand the block-type antiferromagnetism in quasi-1D iron chalcogenides.

antiferromagnetism↗
DOE OSTI2026

Universal Relationship between Limiting Current and Electrochemical Transport Properties in Malonate-Based Polymer Electrolytes

There is considerable interest in developing high-performance electrolytes for rechargeable lithium batteries. For practical applications, the electrolyte must support large dc currents. However, the parameters most often reported in the literature, conductivity, κ, and current fraction, ρ + , reflect ion transport in the limit of infinitesimal currents. In this limit, the efficacy of an electrolyte is given by the product κρ + . The limiting current density, i lim , is the maximum current density that can be applied across an electrolyte; the cell voltage diverges if the applied current density exceeds ilim. This parameter reflects ion transport in the limit of large dc currents and is therefore of practical interest. It would therefore be convenient if i lim could be predicted from measurements of κρ + . In order to explore this possibility, we studied six malonate-based polymers and PEO at a fixed salt concentration (r = 0.08) and temperature (90°C) using symmetric cells with planar electrodes. Unfortunately, there is no correlation between ilim and κρ + . When the applied current density, i, is less than ilim, the cell voltage approaches a stable plateau, ϕ plateau . Here, we found a linear dependence between i and thickness-normalized plateau potential, ϕ plateau L –1 , irrespective of the magnitude of the applied current. In all seven polymer electrolytes, we found a linear correlation between ilim and the slopes of these lines, σ. In other words, measurements of σ can be used to predict the limiting current.

Jana, Rounak [Lawrence Berkeley National Laborator↗
DOE OSTI2025

Cooperative and bifunctional Ga-Ca-Cr 2 O 3 @CaO structured monoliths as versatile platform for reactive capture of CO 2 and its subsequent conversion to ethylene

Cooperative and bifunctional materials (BFMs) that integrate adsorbents and catalysts offer a promising strategy for the reactive capture of CO 2 to produce valuable fuels and chemicals. In this study, we developed structured BFMs via 3D printing that combine CaO as an adsorbent with Ga–Ca–Cr 2 O 3 metal oxides as the catalyst for the reactive capture of CO 2 and its subsequent conversion to C 2 H 4 via the oxidative dehydrogenation of C 2 H 6 (CO 2 -ODHE). Three different Ga–Ca compositions were used to modify the catalyst surface characteristics and enhance C 2 H 4 selectivity. In these formulations, Ga ions stabilize the oxygen lattice of the BFM, while Ca ions interact strongly with Cr to form CaCrO 4 , thereby altering the oxygen species and enhancing the material’s basic properties. Under adsorption–reaction conditions at 600–650 °C, the optimal BFM achieved an excellent C 2 H 4 selectivity of 96.4 %, attributed to a balanced redox process and improved basicity that facilitate efficient C 2 H 6 conversion and rapid desorption of C 2 H 4 without excessive oxidation. Overall, this work provides new insights into the formulation of BFMs monoliths and highlights the critical role of catalytic surface modification in enhancing C 2 H 4 selectivity in the CO 2 -ODHE reactive capture process.

C2H4 production↗
DOE OSTI2026

Multi-purpose quantum laboratories from superconducting circuits

Superconducting circuits (SCs) are the cornerstone of modern quantum technology, enabling scalable computing through coherent control of macroscopic quantum states. Through a legacy that predates modern quantum computing, SCs have emerged as high-precision instruments for discovery. In this review, we highlight the role of SCs as general-purpose quantum laboratories, outlining the emerging landscape of correlated matter-circuit science. We review and unify the capabilities of superconducting quantum hardware across condensed matter, high energy and quantum information sciences. We trace the technical evolution of these architectures, illustrating how their foundational development has culminated in a toolkit for resolving the complexities of macroscopic quantum states.

Arora, Arpit [UCLA, Los Angeles (main); UCLA; Haim↗
DOE OSTI2026

The impact of metastability on the high-pressure behavior of cerium

The structures adopted by solids under pressure are often assumed to reflect thermodynamic equilibrium, yet in many materials phase selection is strongly influenced by kinetic pathways and microstructural inheritance. Elemental cerium (Ce) exemplifies this challenge, with decades of conflicting reports describing different high-pressure crystal structures emerging under nominally identical conditions. Here we use neutron diffraction from large ( ~ 60 mm 3 ) sample volumes to follow the structural evolution in ultra-high-purity Ce during controlled pressure-temperature cycling between 85 and 295 K and up to 8 GPa. We find that the crystal structure formed at high pressure depends on the compression pathway: slow compression ( ~ 0.25 GPa hr −1 ) at room temperature favors an orthorhombic phase (α'), whereas slow ( ~ 0.25 GPa hr −1 ) and also moderately faster ( ~ 0.5 GPa hr −1 ) compression at low temperature stabilizes a pure monoclinic phase (α"). The low-temperature phase persists metastably over a wide temperature range but transforms irreversibly upon heating above ~ 280 K, or modest pressure cycling. Remarkably, the lower-pressure γ phase remains trapped far beyond the expected stability field, persisting to the highest pressures studied. These observations show that phase selection in Ce is governed by kinetics and microstructural memory rather than equilibrium thermodynamics or sample morphology alone, establishing path dependence as a defining feature of its high-pressure behavior.

Ridley, Christopher J. [Oak Ridge National Laborat↗
NASA NTRS1995

Analysis of Developing Laminar Flows in Circular Pipes Using a Higher-Order Finite-Difference Technique

A higher-order finite-difference technique is developed to calculate the developing-flow field of steady incompressible laminar flows in the entrance regions of circular pipes. Navier-Stokes equations governing the motion of such a flow field are solved by using this new finite-difference scheme. This new technique can increase the accuracy of the finite-difference approximation, while also providing the option of using unevenly spaced clustered nodes for computation such that relatively fine grids can be adopted for regions with large velocity gradients. The velocity profile at the entrance of the pipe is assumed to be uniform for the computation. The velocity distribution and the surface pressure drop of the developing flow then are calculated and compared to existing experimental measurements reported in the literature. Computational results obtained are found to be in good agreement with existing experimental correlations and therefore, the reliability of the new technique has been successfully tested.

Herbert J Gladden↗
DOE OSTI2026

Microstructure-dependent particulate populations in ODS steels revealed by single-particle ICP-TOF-MS

In this work, we evaluated the ability of single-particle inductively coupled plasma time-of-flight mass spectrometry (SP-ICP-TOF-MS) to characterize particulates generated from 300 series stainless steels and oxide dispersion–strengthened (ODS) steels. Particulates were produced using two different rotary tools (cutting disk and grinding bit) to assess the influence of particle generation on effluent particulates. Across all materials and preparation techniques, SP-ICP-TOF-MS detected a substantial population of particles containing the expected constituent elements, with average compositions closely matching the nominal alloy values. A secondary population of Y-bearing particles was evident in ODS steels; their formation is strongly favored in particulates generated using the grinding bit. Specifically, in PM2000, Y-bearing particles were predominantly colocated with Al, whereas in 14YWT, Y-bearing particles were less frequently Fe-free but showed stronger colocation with Ti. These trends are consistent with the known microstructural differences between Y–Al–O oxide particles in PM2000 and finer Y–Ti–O nanoclusters in 14YWT. Together, these results demonstrate that SP-ICP-TOF-MS can resolve distinct particulate populations in complex, multicomponent steels at the femtogram mass scale, enabling particle-resolved insight into microstructure-driven phase separation and preparation-dependent particle generation.

Brubaker, Zach E. [Oak Ridge National Laboratory (↗
DOE OSTI2025

A Cross‐Linked Flexible Metaferroelectrolyte Regulated by 2D/2D Perovskite Heterostructures for High‐Performance Compact Solid‐State Sodium Batteries

Abstract To address the issues of limited ionic conductivity and poor interface stability at room and low temperatures in solid‐state electrolytes, a robust intrinsic ferroelectrolyte or nanoferroelectrolyte strategy for engineering solid‐state flexible ferroelectric composite electrolytes utilizing strongly coupled intrinsic ion conducting 2D/2D sodium‐rich anti‐perovskite (NaRAP)/ferroelectric perovskite heterostructures is introduced. Herein, highly scalable PVDF‐based metaferroelectrolytes with Na 2.99 Ba 0.005 OCl/Ca 2 Na 2 Nb 5 O 16 − (CNNO − ) nanosheets into a ferroelectric poly(vinylidene fluoride‐co‐hexafluoropropylene) (PVDF‐HFP) matrix, through an in situ cross‐linking and spontaneous bridging method, for compact solid‐state sodium batteries (SSBs), are reported. Benefiting from unique well‐dispersed 3D ferroelectric coupled network and the Na 2.99 Ba 0.005 OCl/CNNO − ‐induced PVDF‐HFP ferroelectric β phase, the Na + flux is regulated, thereby inhibiting Na dendrite growth at the interface. Notably, the optimized PH‐5% NC metaferroelectrolyte exhibits rapid ion transport (1.11 × 10 −4 S cm −1 at 25 °C), a wide electrochemical window (> 4.8V), superior conformal mechanical compatibility, improved flexibility, good elasticity and flame retardancy. The solid‐state Na 3 V 2 (PO 4 ) 3 /PH‐5% NC/Na batteries present a stable cycling performance (remaining 56.4 mAh g −1 after 500 cycles at 1 C) even at 0 °C, potential for cost‐effective, safe, stable and compact SSB energy storage over 600 Wh L −1 , vastly surpassing 365 Wh L −1 of the current commercial sodium‐ion liquid‐electrolyte batteries.

Chemistry↗
NASA NTRS1964

Thermal Design

An essential phase of the Explorer XIII development was a thermal study and design to provide an acceptable temperature environment for the electronics components and external surfaces of the satellite during ascent and in orbit. It was necessary to establish by preflight analysis and tests that tolerable temperatures could be maintained in three regimes of flight. In the first regime (from launch until release of the heat shield) the satellite was subjected to radiative and conductive heat from the shield. In the second regime, after release of the shield at 350,000 feet, the satellite was heated by free molecular flow. In the third regime - the orbit phase - it was necessary that temperatures within limits for a 1-year lifetime be established. In addition to these three regimes, a study was made to investigate the effect of elevated rocket-motor temperatures during launch and after burnout. This chapter will deal with some pertinent preflight estimates and correlation of these estimates with flight data.

Explorer XIII Satellite↗
NASA NTRS1965

Manufacturing Research in Support of Saturn V

Research achievements applicable to the production of Saturn V components a r e described in this report. Leakage at the mechanical connections in high-pressure lines has been minimized or eliminated through the use of precisely flared tubes. Flaring to close tolerance is done with a n orbital flaring adapter, and the forming operation is precisely controlled with an electronic control console. A technique of “aging” or heat-treatment hardening of aluminum materials has been used to form cylindrical-tank segments and double-contoured gore segments precisely, and to remove weld distortions). The control of weld porosity has been improved as a result of research which demonstrated the close correlation between electrical energy input per inch of weld and weld porosity and metallurgical quality. Based on this information, the gas metallic arc method has been used successfully for two-pass welding, and thinner gage aluminum has been welded (LOX tunnel to fuel bulkhead) with an MSFC designed gas metallic arc unit.

COMPOSITE MATERIAL↗
DOE OSTI2026

Relativistic Magnetohydrodynamic Simulations of Giant Magnetar Bursts

Gradual crustal deformation can generate strongly twisted magnetic fields around magnetars, potentially triggering giant flares with total energies exceeding 10 44 erg. In this letter, we present the first relativistic magnetohydrodynamic simulation of a surface shear-driven magnetar eruption, capturing reconnection-driven plasma heating, the ejection of relativistically hot plasma, and the formation of a hot fireball confined within the inner magnetosphere. We find that magnetic reconnection in the equatorial current sheet launches a hot trailing outflow capable of powering the initial spike observed in giant flares, while simultaneously leaving behind a thermally stratified fireball with sufficient thermal energy to produce the pulsating, decaying tail. Together, these features provide a self-consistent physical framework for understanding the observed energetics of magnetar giant flares. The eruption also expels a magnetically dominated giant plasmoid carrying up to ∼9% of the magnetosphere’s total magnetic energy. Furthermore, our simulation demonstrates how the plasmoid drives the formation of a blast wave—an important ingredient in models linking magnetar eruptions to fast radio bursts.

79 ASTRONOMY AND ASTROPHYSICS↗