Search NASA⌕ Search

SEARCH · Search NASA

Results for “Modeling and Simulation”

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.

3,221 records · Page 107
All sourcesNASA NTRSDOE OSTI
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

Molecular origin of anisotropic shear elastoplasticity in chitin

Chitin nano- and mesoscale structures present in the exoskeleton of crustaceans exhibit exceptional longitudinal stiffness and toughness, rivaling or even exceeding that of many synthetic polymer architectures. Here, we reveal the origin of the asymmetric shear response in chitin multiscale architectures, marked by pronounced anisotropy in deformation. Under shear aligned with the molecular axis, chitin accommodates strain through coherent atomic rearrangements that enable elastic recovery. In contrast, shear applied perpendicular to this axis induces liquid-like plasticity via localized sliding. These results demonstrate the intrinsic mechanical anisotropy of chitin, underpinning a dual function in resisting repetitive loading while dissipating internal stress during high-strain events. Our findings establish the molecular basis of shear elastoplasticity in multiscale chitin structures, wherein axial elasticity supports energy storage in load-bearing regions, whereas transverse plasticity enables controlled energy dissipation. These atomic-scale insights lay a foundation for the predictive design of chitin-based materials with tunable strength-toughness profiles.

Wan, Zhangmin [University of British Columbia, Van↗
DOE OSTI2026

Aliovalent gallium dopants remove Ti 3+ defects and improve photocatalytic and photoelectrochemical water oxidation properties of LaTiO 2 N

LaTiO 2 N is a promising semiconductor for the water splitting reaction due to its 2.1 eV band gap and stability against corrosion. However, its solar energy conversion is limited by Ti 3+ recombination defects introduced during ammonolysis. Here we show for the first time that Ti 3+ defects in LaTiO 2 N can be suppressed with incorporation of 2, 5, and 10% aliovalent gallium (Ga 3+ ) dopants during synthesis via the layered La 2 Ti 2 O 7 intermediate. Electron paramagnetic resonance (EPR) spectroscopy on the solid powders confirms a reduction in the Ti 3+ donor density from 2.97 × 10 17 cm −3 for the non-doped material to ∼6.24 × 10 16 cm −3 for 5% Ga-doped LaTiO 2 N. The remaining Ti 3+ defects are concentrated near the LaTiO 2 N surface, according to X-ray photoelectron spectroscopy. The defect reduction shifts the optical absorption edge from 2.02 to 2.09 eV and eliminates a broad absorption band at 1050 nm from the optical absorption spectra. It also removes a 1.0–1.7 eV sub-band gap photovoltage signal from surface photovoltage spectra. This suggests that empty Ti 3+ d-orbitals are located 1.0–1.7 eV above the LaTiO 2 N valence band edge. Removing these recombination states with increasing Ga 3+ content enhances the photoconversion efficiency of LaTiO 2 N during water oxidation. The optimized 2 wt% CoO x -loaded 5% Ga-doped LaTi O2 N material has a 16% AQE (400 nm) for O 2 production from aqueous silver nitrate solution and a ∼2.1 mA cm −2 water oxidation photocurrent at 1.23 V under 100 mW cm −2 Xe arc lamp illumination. The water oxidation photocurrent is stable during a 50 min test, and the Faraday efficiency for O 2 generation is 97%, confirming short-term corrosion stability of LaTiO 2 N. Altogether, these results provide a better understanding of the distribution, concentration, and impact of Ti 3+ defects on the optical, photovoltage, and photoelectrochemical properties of LaTiO 2 N. In combination with other defect control strategies, aliovalent Ga 3+ doping can help bring the solar energy conversion efficiency of LaTiO 2 N closer to the theoretical limit.

Wang, Li [University of California, Davis, CA (Uni↗
DOE OSTI2026

Evolution of magnetic bubble domains in the uniaxial ferromaget CeRu 2 Ga 2 B inferred from the Hall effect and ac magnetic susceptibility

We study the Hall effect, AC magnetic susceptibility (χ ac ), and magnetic force microscopy of the uniaxial ferromagnet CeRu 2 Ga 2 B with a centrosymmetric crystal structure. We observe a finite topological Hall effect (THE) within the ordered phase, before the magnetization is polarized by applied field. By comparing the field dependences of the area fraction of the magnetic bubbles, the derivative of χ ac , and the THE signal, we deduce that the magnetic bubbles in CeRu 2 Ga 2 B evolve from the trivial to topological spin texture with field. Our findings enable the expansion of the search for magnetic materials hosting topological spin textures to include uniaxial ferromagnets and open a new possibility to tailor the topological spin texture.

AC magnetic susceptibility↗
DOE OSTI2026

Hierarchical Epoxy Structures via Tunable Polymerization-Induced Phase Separation Combined with Additive Manufacturing

Polymerization-induced phase separation (PIPS) allows for the control of thermoset morphologies and properties, enabling the tuning of domain sizes and thermomechanical response. However, its use in generating substructural features in additively manufactured materials has been limited. In this work, we combine epoxy PIPS with UV curable acrylate and rheological modifiers to print nano- to macro-phase separating materials via a two-step, dual-cure approach. This method enables direct ink write printing of hierarchical structures with both controlled morphologies through phase separation and macroscale architecture through print design. We find that formulations for phase-separating materials require judicious incorporation of additives to enable printability and to provide sufficient green strength. Atomic force microscopy-nano infrared mapping reveals tunable, reticulated nano- to micron-scale domains of the resultant multiphase materials and their morphology changes due to additives, resulting in alterations to thermomechanical and tensile properties. Shape memory behavior is also demonstrated through multimaterial additive manufacturing of epoxies with functionally graded internal morphology using active mixing techniques, highlighting this method’s ability to fabricate complex architectures with controlled morphologies and thermomechanical response.

Van Meter, Kylie E [Organic Materials Science, San↗
NASA NTRS2026

Experimental Study of Aeropropulsive Interactions for Electric Ducted Fans on a Wing

The Adaptable Distributed Electric Propulsion Testbed (ADEPT)Wing was tested in the NASA Langley 12-Foot Low-Speed Tunnel to investigate electric propulsion technologies for commercial transport aircraft. The purpose of the test was to study the aeropropulsive effects across angle of attack and fan speed variations for several propulsor placements and array sizes, and to investigate how distributed electric propulsion (DEP) can be leveraged for augmented lift and vehicle control. Six electric ducted fan (EDF) array placements were tested by positioning the fans along the leading and trailing edges on the top and bottom surfaces, as well as centered on the chord line in front of and behind the wing. Each fan placement configuration was tested with an array of one, three, and five EDFs to characterize an aeropropulsive extrapolation effect, and the bare wing was tested separately to provide baseline aerodynamic measurements. The results compare the longitudinal forces and moments for each configuration by analyzing the combined aeropropulsive interaction effects, and offer valuable insight into considerations for fan placement along a wing to enable the benefits of DEP integration.

Rose Weinstein↗
NASA NTRS2026

The Relevance of the VIPER Mission to NASA’s Artemis Human Exploration of the Moon

NASA has ambitious plans to send astronaut crews to the south polar region of the Moon and explore this lunar terrain with humans for the first time through a series of Artemis missions. Artemis will provide economic benefits, drive technological advancement, and inspire the next generation of explorers. Simultaneously, Artemis presents a unique opportunity for humans to conduct high-priority planetary science in situ in the south polar region of the Moon. Prior to human exploration, landed robotic precursor investigations can provide valuable information to reduce risk and maximize discovery and productivity of subsequent crewed missions. NASA’s Volatiles Investigating Polar Exploration Rover (VIPER) is a robotic mission designed to explore and characterize the lunar south polar region prior to crewed Artemis landed missions. VIPER is designed to explore multiple thermal regimes to characterize the lunar polar environment and regolith properties, and to explore for surface and subsurface ices. Each of these objectives are key for both scientific investigation and to detect and assess volatile deposits to support in situ resource utilization (ISRU) and a sustained human presence on the Moon. VIPER would also operate with a real-time mission operations architecture with relevance and feed-forward to Artemis crewed operations on the lunar surface. VIPER would provide key information regarding the lunar environment, scientific exploration, technologies, and real-time operations to optimize the valuable surface exploration time of Artemis crew members while simultaneously reducing risk to crew and increasing health and human safety on the lunar surface.

Jennifer L Heldmann↗
DOE OSTI2026

Advanced LiFSI-LiPF6 Electrolyte for Wide-Temperature and Thermally Stable Lithium-Ion Batteries

A new optimized LiFSI–LiPF6 dual-salt controlled-solvation electrolyte (E-DS) is demonstrated to enable practical graphite||LiNi0.8Mn0.1Co0.1O2 cells (˜4.0 mAh cm?²) to achieve exceptional performance and safety under extreme conditions. By optimizing anion coordination with the smaller, more dissociating FSI? anion, the E-DS forms ultrathin, dense, and inorganic-rich electrode/electrolyte interphases that dramatically suppress solvent decomposition, transition-metal dissolution, and surface reconstruction compared to the conventional LiPF6/carbonate electrolyte. Consequently, E-DS cells deliver >78% capacity retention after 300 cycles at 60 °C, retain fast discharging capacity at 30 °C, and operate effectively at -20 °C. Most strikingly, fully charged full cells with E-DS, even under overcharging to 4.8 V, show a lower heat evolution in stable formulations — transforming a traditionally unstable high-voltage/high-temperature configuration into an intrinsically safe state. This work establishes a new benchmark for carbonate-containing electrolytes, simultaneously achieving high energy density, fast-discharging capability, wide-temperature operation (-20 to 60 °C), and outstanding thermal safety in nickel-rich lithium-ion batteries.

electrode/electrolyte interphase↗
NASA NTRS2026

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

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

stagnation heat flux↗
NASA NTRS2026

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

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

inviscid theory↗
NASA NTRS1992

Long Duration Exposure Facility experiment M0003 deintegration observation database

The four trays of the M0003 materials experiment on LDEF contained 1274 samples from 19 subexperiments. The complete sample complement represented a broad range of materials, including thin film optical coatings, paints, polymer sheets and tapes, adhesives, and composites, for use in various spacecraft applications including thermal control, structures, optics, and solar power. Since some subexperiments included duplicate samples exposed on the leading and trailing edge trays, for 9 week, 19 week, 40 week, and 68 month durations on both the leading and trailing edges, comparisons provided a valuable time history of degradation. During the deintegration of the subexperiments from the M0003 trays, each sample was examined using bright field, dark field, and Normarski light microscopy techniques. Over 3,000 color macrographs and photomicrographs were made of the condition of M0003 trays and the individual samples. Records of the condition of the samples, photographic records, and all available references of publications on postflight analyses on each sample were collected in a database. The photographs provide a visual comparison of the response of materials common to various subexperiments in different LDEF environments. Aerospace Corporation is distributing the database records and photographs and written records on materials of interest from the database is given in this presentation.

S. R. Gyetvay↗
NASA NTRS1993

Evidence for the Interaction of the IRS 16 Wind With the Ionized and Molecular Gas at the Galactic Center

We present a number of high-resolution radio images showing evidence for the dynamical interaction of the outflow arising from the IRS 16 complex with the ionized gas associated with the Northern Arm of Sgr A West, and with the northwestern segment of the circumnuclear molecular disk which engulfs the inner few parsecs of the Galactic center. These interactions offer an opportunity to explain some important aspects of the complex morphology of ionized and molecular gas at the Galactic center. In particular we suggest that the wind disturbs the dynamics of the Northern Arm within 0.1 pc of the center, is responsible for the waviness of the arm at larger distances and is collimated by Sgr A West and the circumnuclear disk. The waviness is discussed in terms of the Rayleigh-Taylor instability induced by the ram pressure of the wind incident on the surface of the Northern Arm. Another consequence of this interaction is the strong mid IR polarization of the Northern Arm in the vicinity of the IRS 16 complex which is explained as a result of the ram pressure of the wind compressing the gas and the magnetic field. The IRS 16 complex must therefore lie close to the Northern Arm. On a large scale we argue that the wind has swept up the diffuse material in the inner parsec of the Galactic center and that the interaction with the inner part of the northwestern part of the circumnuclear ring is responsible for the H 2 v = 1-0 emission and the ionized streamers. Lastly, we discuss the possibility that the outflow from IRS 16 could be responsible for eating away the inner edge of the circumnuclear disk and provide friction needed for the infall of the material in the arms of Sgr A West.

IRS 16↗
DOE OSTI2026

Flame stabilization in DME spray flames under engine-relevant conditions characterized by OH* chemiluminescence and formaldehyde laser-induced fluorescence

The transient and quasi-steady flame structures of Dimethyl Ether (DME) fuel sprays, produced by a single-hole injector (Spray D), were investigated using Planar Laser-Induced Fluorescence (PLIF) and chemiluminescence imaging in a constant-volume chamber under Engine Combustion Network (ECN) Spray A conditions (900 K ambient temperature, 60 bar ambient pressure, 1500 bar injection pressure, and 22.8 kg/m 3 ambient density). Low-temperature chemical reaction zones were visualized using formaldehyde (CH 2 O) PLIF with 355 nm excitation, while high-temperature flame regions were captured via chemiluminescence imaging of excited-state hydroxyl radicals (OH*). Both transient and quasi-steady flame structures clearly show the transition from CH 2 O to OH*, highlighting the progression from low- to high-temperature combustion, while the position of the flame is displaced for DME compared to reference hydrocarbon n-dodecane. Homogeneous reactor calculations with detailed chemistry and using adiabatic mixing for initial temperature show that CH 2 O peaks are significantly higher for DME at the same equivalence ratio, with a higher heat-release during the cool-flame regime with respect to the fuel heating value. Thus, the cool-flame dynamic as a precursor to high-temperature combustion and flame stabilization exhibit distinct behavior for DME relative to conventional hydrocarbons, and these phenomena are effectively resolved through the soot-free nature of DME and the high-speed, time-resolved diagnostics.

CH2O laser-induced fluorescence↗
DOE OSTI2024

A thermal evaporation–trapping strategy to synthesize flexible and robust oxygen electrocatalysts for rechargeable zinc–air batteries

Great efforts have been devoted to the development of bifunctional electrocatalysts to accelerate the sluggish kinetics of cathodic oxygen reduction/evolution reactions (ORR/OER) in zinc–air batteries (ZABs). Here we report a thermal evaporation–trapping synergistic strategy to fabricate a bifunctional electrocatalyst of flexible N-doped carbon fiber cloth loaded with both CoFe-oxide nanoparticles and single-atom Co/Fe–Nx sites, in which the thermal evaporation process functions in both downsizing CoFe-oxide nanoparticles and trapping the evaporated Co/Fe species to generate Co/Fe–Nx sites. The obtained flexible electrocatalyst, directly serving as an oxygen electrode, displays a small potential gap of 0.542 V for the OER/ORR, large peak power densities (liquid-state ZAB: 237.4 mW cm2 ; solid-state ZAB: 141.1 mW cm2 ), and excellent charge–discharge cycling stability without decay after working more than 770 hours. Furthermore, in situ Raman spectroscopy characterization and theoretical calculations reveal that CoFe2O4 species is responsible for the OER while atomic Fe/Co sites play a key role in the ORR

Zhang, Hong-Bo↗
NASA NTRS2026

NASA Hybrid Thermally Efficient Core (HyTEC) Phase 2 Q3’GFY25 to Q2’GFY26 Unlimited Rights Annual Report for CLIN 1 and 2

The HyTEC Phase 2 Project focuses on the development and demonstration of high-power density, small engine core gas turbine engine technologies. The focus of HyTEC on these small engine core technologies will provide direct benefits to the next single-aisle class aircraft in terms of thermal efficiency, as well as integrate with other technologies, through increased hybridization, as they mature to provide substantial continuous fuel burn reductions during the aircraft lifecycle. These advances will strengthen the US position in the commercial aviation engine market and enable cost benefits in commercial aviation. In Government fiscal year 26 the HyTEC Project was closed and contract moved to the Subsonic Vehicle Technologies and Tools Project. GE Aerospace is executing two major work plans as part of the HyTEC Phase 2 contract. The first is Contract Line-Item Number (CLIN) 001 and its objective is to mature High Pressure Turbine (HPT) aerodynamics technology to TRL 5 through multiple rigs so that it may be incorporated into a complete TRL 6 evaluation and unlock learnings for the engine core demonstration as applicable. CLIN 002, the second work plan, is to design, procure and test a cost shared portion of technologies on the TRL6 Compact Core test vehicle. The compact core demonstrated through HyTEC Phase 2 along with the HPT rig TRL advancement directly matures the technologies being developed as part of the CFM RISE demonstrator program which is targeting a more than 20% fuel burn reduction at the engine level versus today’s state of the art. The RISE program’s objective is to advance both the novel Open Fan architecture, compact core technologies and advanced systems such as hybrid electric needed to achieve the desired fuel burn reduction to TRL6 ahead of a new product introduction in the single-aisle class aircraft in the 2030s timeframe.

gas turbine engine↗
DOE OSTI2026

An Automated Probabilistic Asteroid Prediscovery Pipeline

We present an automated and probabilistic method to make prediscovery detections of near-Earth asteroids (NEAs) in archival survey images, with the goal of reducing orbital uncertainty immediately after discovery. We refit the Minor Planet Center's astrometry and propagate the full six-parameter covariance to survey epochs to define search regions. We build low-threshold source catalogs for viable images and evaluate every detected source in a search region as a candidate prediscovery. We eliminate false positives by refitting a new orbit to each candidate and probabilistically linking detections across images using a likelihood ratio. Applied to the Zwicky Transient Facility's (ZTF) imaging, we identify approximately 3000 recently discovered NEAs with prediscovery potential, including a doubling of the observational arc for about 500. We use archival ZTF imaging to make prediscovery detections of the potentially hazardous asteroid 2021 DG1, extending its arc by 2.5 yr and reducing future apparition sky plane uncertainty from many degrees to arcseconds. We also recover 2025 FU24 nearly 7 yr before its first known observation, when its sky plane uncertainty covers hundreds of square degrees across thousands of ZTF images. The method is survey agnostic and scalable, enabling rapid orbit refinement for new discoveries from Rubin, NEO Surveyor, and NEOMIR.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Robust Heat-Flux Sensors for Coal-Fired Boiler Extreme Environments

In this project, robust heat-flux measurement systems were developed. The heat-flux sensors utilize thermoelectric effects to directly transduce the heat-flux inputs to analog electrical voltage signals. They were constructed from dedicated materials that can withstand temperatures of at least 1000°C and maintain adequate performance at these conditions for prolonged periods of time. The proposed approaches took into account numerous considerations, including system cost, sensor head resilience, sensor footprint, data accuracy, response time, and maintenance requirements. Through modern thermoelectric materials design, methodical materials selection and rigorous testing in materials characterization labs and medium-scale fire research facilities, we have demonstrated functioning laboratory prototypes, upon which one could base industrial heat-flux sensing platforms capable of operating in the challenging high-temperature, corrosive environments of the boilers of coal-fired power plants. A distributed sensor array for heat-flux measurements throughout the furnace water-wall, the superheater area and the economizer coils can provide critical data for the power plant control systems to increase efficiency, improve safety and reduce down times. For example, the combined heat-flux sensor/control systems can contribute to the optimization of burner and boiler operations under flexible loads, the optimization of heat-exchange conditions and overall reduction of heat rate and emissions, the prediction of imminent overheating conditions, and the optimization of the soot-blowing protocols.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2024

Unravelling ultralow thermal conductivity in perovskite Cs2AgBiBr6: dominant wave-like phonon tunnelling and strong anharmonicity

Abstract Understanding the lattice dynamics and heat transport physics in the lead-free halide double perovskites remains an outstanding challenge due to their lattice dynamical instability and strong anharmonicity. In this work, we investigate the microscopic mechanisms of anharmonic lattice dynamics and thermal transport in lead-free halide double perovskite Cs 2 AgBiBr 6 from first principles. We combine self-consistent phonon calculations with bubble diagram correction and a unified theory of lattice thermal transport that considers both the particle-like phonon propagation and wave-like tunnelling of phonons. An ultra-low thermal conductivity at room temperature (~0.21 Wm −1 K −1 ) is predicted with weak temperature dependence( ~ T −0.34 ), in sharp contrast to the conventional ~T −1 dependence. Particularly, the vibrational properties of Cs 2 AgBiBr 6 are featured by strong anharmonicity and wave-like tunnelling of phonons. Anharmonic phonon renormalization from both the cubic and quartic anharmonicities are found essential in precisely predicting the phase transition temperature in Cs 2 AgBiBr 6 while the negative phonon energy shifts induced by cubic anharmonicity has a significant influence on particle-like phonon propagation. Further, the contribution of the wave-like tunnelling to the total thermal conductivity surpasses that of the particle-like propagation above around 310 K, indicating the breakdown of the phonon gas picture conventionally used in the Peierls-Boltzmann Transport Equation. Importantly, further including four-phonon scatterings is required in achieving the dominance of wave-like tunnelling, as compared to the dominant particle-like propagation channel when considering only three-phonon scatterings. Our work highlights the importance of lattice anharmonicity and wave-like tunnelling of phonons in the thermal transport in lead-free halide double perovskites.

Chemistry↗