Engineering PapersSearch

SEARCH · Engineering Papers

Results for “Atmospheric dynamics”

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.

337 records · Page 8

Thermal Disorder‐Induced Strain and Carrier Localization Activate Reverse Halide Segregation

The reversal of halide ions is studied under various conditions. However, the underlying mechanism of heat-induced reversal remains unclear. This work finds that dynamic disorder-induced localization of self-trapped polarons and thermal disorder-induced strain (TDIS) can be co-acting drivers of reverse segregation. Localization of polarons results in an order of magnitude decrease in excess carrier density (polaron population), causing a reduced impact of the light-induced strain (LIS – responsible for segregation) on the perovskite framework. Meanwhile, exposing the lattice to TDIS exceeding the LIS can eliminate the photoexcitation-induced strain gradient, as thermal fluctuations of the lattice can mask the LIS strain. Under continuous 0.1 W cm -2 illumination (upon segregation), the strain disorder is estimated to be 0.14%, while at 80 °C under dark conditions, the strain is 0.23%. However, in situ heating of the segregated film to 80 °C under continuous illumination (upon reversal) increases the total strain disorder to 0.25%, where TDIS is likely to have a dominant contribution. Therefore, the contribution of entropy to the system's free energy is likely to dominate, respectively. Various temperature-dependent in situ measurements and simulations further support the results. These findings highlight the importance of strain homogenization for designing stable perovskites under real-world operating conditions.

36 MATERIALS SCIENCE

Transforming Windows from Energy Liabilities to Zero-Energy Assets: Next-Generation Solutions for Buildings

Windows have traditionally contributed to a building's HVAC load, but they can also become a source of net energy gain or even operate as zero-energy components. For heating applications, highly insulating windows can harness more solar heat than the energy lost through them, transforming windows from energy liabilities to assets. Dynamic glazings provide further benefits by regulating solar heat gain, reducing cooling loads in summer and heating demands in winter. This simulation study focuses on developing the next generation of zero-energy windows (ZEW) for residential new construction. Through annual energy simulations across climate zones 1-8, ZEW performance benchmarks were established based on current code-level buildings, and we've identified the regions where meeting ZEW standards are most achievable. This work evaluates both static and dynamic window technologies, assessing their effects on annual energy use and cost. Key findings demonstrate that ZEW performance is achievable across diverse climate zones, with specific regional requirements. Most climate zones from 3-8 can achieve ZEW with specific configurations, while some warm climates (1-2) appear challenging for ZEW implementation. Climate zones 4-6 consistently allow for zero energy window implementation, offering multiple pathways through either static or dynamic window technologies. Colder climate zones (7-8) ZEW products allow for higher SHGC values while requiring low U-values.

Yu, Lili

Numerical simulation of frost formation and heat transfer on fin-and-tube heat exchangers in turbulent cross-flow

Frost formation in fin-and-tube heat exchangers in turbulent cross-flow presents significant challenges in industrial refrigeration applications, affecting heat transfer efficiency and operational reliability. The purpose of this work is to investigate frost deposition and growth on a staggered bank of a fin-and-tube freezer coil under turbulent forced convection conditions. The focus here is on investigating conditions that closely replicate real-world scenarios in large walk-in industrial freezers. Using a direct numerical simulation approach, we examine the flow dynamics and thermal behaviour in the presence of frost, considering turbulent regimes characterized by a Reynolds number in the range 1050 ≤ R e D , avg ≤ 4800 , with the characteristic length being the outer diameter of the tube and the velocity being the bulk fluid velocity between the plates (fins). Computational fluid dynamics simulations are employed to resolve the interactions between turbulent airflow and the frost layer. Our approach incorporates a modified immersed boundary method and a slow-time acceleration technique to address the complex dynamic interface between the continuously evolving frost layer and the flowing air stream. Our findings indicate that frost forms more on the sides of the finned surfaces (plates) and less on the tubes themselves. This article is part of the theme issue ‘Heat and mass transfer in frost and ice’.

Science & Technology - Other Topics

Purpose, Design, and Analysis of the NSTX-U Centerstack Casing Lateral Support Shims

The National Spherical Torus Experiment Upgrade (NSTX-U) centerstack casing surrounds the central magnet core of NSTXU. It provides the inner vacuum boundary of the plasma chamber and supports the radially inboard plasma facing components lining the plasma chamber. The casing is supported like a cantilever, fixed at the lower flange by a connection to the TF Bundle through a G-10 ring. It is flexibly supported at the top by bellows, and insulated shims intended to react the largest part of the inventory of the disruption loads at the upper support, thereby protecting the upper and lower connections from excessive loads. A large part of the qualification of the casing and tiles requires addressing and quantifying disruption loads. These have been estimated based on extensions of the behavior of NSTX during its run period roughly 15 years ago. Structural qualification of the Upgrade used these estimates of the loading as input to calculations to demonstrate adequacy for a “most probable” disruption. The estimates are included in a requirements document to guide conservative design of all components loaded by disruptions. However, the changes made for the upgrade made the extrapolation of disruption behavior uncertain. Disruptions include axisymmetric loading from axisymmetric currents induced in the vessel shells. In addition, more complex currents called Halo currents are a result of plasma poloidal currents impinging on the casing. The halo currents add another layer of uncertain behavior and net lateral loading to the casing. As a result, NSTX-U is equipped with a number of instruments intended to assess the nature and severity of disruptions. The casing shims are required for lateral support of the casing, but these have also been converted to load cells to provide another measure of the magnitude of symmetric and asymmetric loading. Some of the logic in the development of the requirements document are presented. Electromagnetic analysis of axisymmetric and halo loading is described. Global simulations of the casing transient dynamic response, and local analysis and testing of the shim/load cell are included. The shim load cells and other instrumentation will be an important guide in planning NSTX-U experiments beyond its next run period expected in late 2026 or early 2027.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY

Nuclear Safety [Vol. 32, No. 4, October-December 1991]

Nuclear Safety is a review journal that covers significant developments in the field of nuclear safety. Its scope includes the analysis and control of hazards associated with nuclear energy, operations involving fissionable materials, and the products of nuclear fission and their effects on the environment. Primary emphasis is on safety in reactor design, construction, and operation; however, the safety aspects of the entire fuel cycle, including fuel fabrication, spent-fuel processing, nuclear waste disposal, handling of radioisotopes, and environmental effects of these operations, are also treated. Table of Contents for this issue follows. GENERAL SAFETY CONSIDERATIONS: 477 Report on the International Symposium on the Use of Probabilistic Safety Assessment for Operational Safety—PSA '91, S. Chakraborty and M. Khatib-Rahbar; 488 Good Relationships Are Pivotal in Nuclear Data Bases, A. S. Heger and B. V. Koen; 494 Technical Note: The Interagency Nuclear Safety Review Panel's Evaluation of the Ulysses Space Mission, J. A. Sholtis, Jr., D. A. Huff, L. B. Gray, N. P. Klug, and R. O. Winchester; ACCIDENT ANALYSIS: 502 The Severe Accident Analysis Program for the Savannah River Nuclear Production Reactors, M. L. Hyder; CONTROL AND INSTRUMENTATION: 511 A Framework for Selecting Suitable Control Technologies for Nuclear Power Plant Systems, R. A. Kisner; DESIGN FEATURES: 521 Containments for Gas-Cooled Power Reactors History and Status, P.W. Williams; ENVIRONMENTAL EFFECTS: 537 Indoor Radon: A Natural Risk, N. H. Harley and J.H. Harley; On the Importance of the Atmospheric Parameters in the Fission Products Distribution of a Severe Reactor Accident, M. C. Barla and A. R. Bayulken; Book Review: Health Effects of Exposure to Low Levels of Ionizing Radiation BEIR V, C. R. Richmond; WASTE AND SPENT FUEL MANAGEMENT: 555 Activities Related to Waste and Spent Fuel Management, M. D. Muhlheim and E. G. Silver OPERATING EXPERIENCES: 567 Effects of Component Aging on the Westhinghouse Control Rod Drive System, K. Sullivan and W. Gunther; 577 Reactor Shutdown Experience, Compiled by J. W. Cletcher, 580 Selected Safety-Related Events, Compiled by G. A. Murphy; 582 Operating U 8 Power Reactors, Compiled by M. D. Muhlheim and E. G. Silver, RECENT DEVELOPMENTS: 596 General Administrative Activities, Compiled by M. D. Muhlheim and E. G. Silver; 610 Reports, Standards, and Safety Guides, D. S. Queener; 614 Proposed Rule Changes as of June 30, 1991; ANNOUNCEMENTS: 520 Workshop on PC PRAISE: A Probabilistic Fracture Mechanics Personal Computer Code for Nuclear Power Plant Piping Reliability Assessment; 536 Fifth Workshop on Nuclear Power Plant Containment Integrity; 624 New OECD "International Information System on Occupational Exposure (ISOE)"; 625 Harvard School of Public Health Announces Short Courses; 625 Eighth Power Plant Dynamics, Control and Testing Symposium; 626 Second Training Course on Off-Site Emergency Planning and Response for Nuclear Accidents; 618 The Authors; 622 Reviewers of Nuclear Safety, Vol 32.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS

Observation of Extraordinary Vibration Scatterings Induced by Strong Anharmonicity in Lead‐Free Halide Double Perovskites

Abstract Lead‐free halide double perovskites provide a promising solution for the long‐standing issues of lead‐containing halide perovskites, i.e., the toxicity of Pb and the low stability under ambient conditions and high‐intensity illumination. Their light‐to‐electricity or thermal‐to‐electricity conversion is strongly determined by the dynamics of the corresponding lattice vibrations. Here, the measurement of lattice dynamics is presented in a prototypical lead‐free halide double perovskite(Cs 2 NaInCl 6 ). The quantitative measurements and first‐principles calculations show that the scatterings among lattice vibrations at room temperature are at the timescale of ≈1 ps, which stems from the extraordinarily strong anharmonicity in Cs 2 NaInCl 6 . Further the degree of anharmonicity of each type of atom is quantitatively characterized in the Cs 2 NaInCl 6 single crystal, which stems from the interatomic forces, and demonstrate that this strong anharmonicity is synergistically contributed by the bond hierarchy, the tilting of the NaCl 6 and InCl 6 octahedral units, and the rattling of Cs + ions. Consequently, the crystalline Cs 2 NaInCl 6 possesses an ultralow thermal conductivity of ≈0.43 W mK −1 at room temperature, and a weak temperature dependence ofT −0.41 . These findings uncovered the underlying mechanisms behind the dynamics of lattice vibrations in double perovskites, which can largely benefit the design of optoelectronics and thermoelectrics based on halide double perovskites.

Chemistry

SODAs: sparse optimization for the discovery of differential and algebraic equations

Differential-algebraic equations (DAEs) integrate ordinary differential equations (ODEs) with algebraic constraints, providing a fundamental framework for developing models of dynamical systems characterized by time-scale separation, conservation laws and physical constraints. While sparse optimization has revolutionized model development by allowing data-driven discovery of parsimonious models from a library of possible equations, existing approaches for dynamical systems assume DAEs can be reduced to ODEs by eliminating variables before model discovery. This assumption limits the applicability of such methods for DAE systems with unknown constraints and time scales. We introduce sparse optimization for differential-algebraic systems (SODAs), a data-driven method for the identification of DAEs in their explicit form. By discovering the algebraic and dynamic components sequentially without prior identification of the algebraic variables, this approach leads to a sequence of convex optimization problems. It has the advantage of discovering interpretable models that preserve the structure of the underlying physical system. To this end, SODAs improves since SODAs is singular numerical stability when handling high correlations between library terms, caused by near-perfect algebraic relationships, by iteratively refining the conditioning of the candidate library. We demonstrate the performance of our method on biological, mechanical and electrical systems, showcasing its robustness to noise in both simulated time series and real-time experimental data.

DAE

A physically interpretable precursor framework for sub-seasonal prediction of Northern Hemisphere flash flourishing

Flash flourishing describes rapid vegetation increases that can quickly reshape land–atmosphere exchanges and impacts on ecosystem, yet its large-scale precursors, circulation context, and sub-seasonal predictability remain poorly understood. Here, we identified onset-stage circulation regimes across northern extratropical latitudes (NEL; >30°N) using 200 and 1000 hPa geopotential height, and examined their regional expressions over eastern Asia, western North America, and Europe. Flash flourishing onset in East Asian was associated with a baroclinic circulation regime and was preceded by a North Atlantic sea surface temperature (SST) precursor at a four-pentad lead. In contrast, onset in western North American and European preferentially occurred under barotropic regimes, preconditioned by Great Plains soil moisture at three-pentad lead and North Atlantic SST at a four-pentad lead, respectively. Ridge regression forecasts revealed regime-dependent sub-seasonal predictability, with mean out-of-sample R 2 exceeding 0.3 up to lead times of two pentads in East Asia, three pentads in western North America, and four pentads in Europe. Together, these findings established a mechanistic and regionally specific framework for anticipating rapid vegetation greening at sub-seasonal timescales.

Kong, Xiangxu [Nanjing Univ. of Information Scienc

Effect of Oxidizing Impurities on the Corrosion Behavior of Structural Alloys Exposed to MgCl 2 Molten Salt Vapor

The corrosion behavior of structural alloys SS304 and IN617 when exposed to vapors of MgCl 2 molten salt at 764 °C for 100 h in an argon atmosphere was investigated as a function of the level of oxidizing impurities (e.g., NaOH) present in the salt. Increasing the concentration of oxidizing impurities in the salt caused an increase in corrosion in both structural alloys, as measured by Cr depletion layers. A poorly adhered oxide scale layer containing Mg, Al, and Cr was observed on the surface of the metals, and local attack was observed as well. The Ni-based IN617 showed lower Cr depletion and less scale formation than the Fe-based SS304. These results demonstrate that oxidizing impurities in the molten salt will impact the vapor phase corrosion for structural members exposed to the molten salt headspace. This finding is highly relevant for understanding the effects of the vapor phase in molten salt applications, including molten salt thermal energy storage, molten salt nuclear reactors, and molten salt electrochemistry.

25 ENERGY STORAGE

Topographic effects on reflected acoustic waves from the OSIRIS-REx reentry observed from stratospheric balloons

During long-distance sound propagation in planetary atmospheres, acoustic waves may reflect off the air/surface interface one or more times. For low sound frequencies and flat interfaces, the incident and reflected wave tend to be nearly identical. However, this may not be the case when the downgoing acoustic wave encounters topography. Here, we describe a set of direct and reflected acoustic signals recorded on free-flying balloons during the hypersonic entry of the OSIRIS-REx sample return capsule (SRC). In two of the three cases presented here, an impulsive reflected arrival similar in form to the direct sonic boom of the SRC is observed, followed by a diffuse coda. In contrast, one of the floating stations lacked an impulsive reflection entirely, with only the coda present. We use full waveform modeling to show how reflection in the presence of complex topography can explain coda seen in all three examples as well as the lack of impulsive arrival on the third. Our results indicate that the complex signals often observed in long range acoustic propagation could be due, in part, to interactions with topography during transmission.

Lees, Jonathan M. [University of North Carolina, C

Carbon in GaN as a nonradiative recombination center

Trap-assisted nonradiative recombination has been shown to limit the efficiency of optoelectronic devices. While substitutional carbon (⁠C N ) has been suggested to be a nonradiative recombination center in GaN devices, a complete recombination cycle including the two charge-state transition levels has not been previously described. In this work, we investigate the trap-assisted recombination process due to C N in GaN, including multiphonon emission, radiative recombination, trap-assisted Auger–Meitner (TAAM) recombination as well as thermal emission of holes. Our study shows the key role of TAAM processes at the high carrier densities relevant for devices. We also reveal the carrier-density regimes where thermal emission and radiative recombination are expected to play an observable role. Our results highlight that carbon concentrations exceeding ~10 17 cm −3 can have a noticeable impact on device efficiency, not just in GaN active layers but also in InGaN and AlGaN. Furthermore, our comprehensive formalism not only offers detailed results for carbon but also provides a general framework for assessing the multiple processes that participate in trap-assisted recombination in semiconductors.

36 MATERIALS SCIENCE

Temperature-Dependent Transport Characteristics of 2D MoS2 Channel FETs Grown Using Salt-Based Precursors

D Transition Metal Dichalcogenides (TMDs), particularly MoS2, are promising candidates for sub-10 nm Gate All Around (GAA) CMOS FETs. Salt-assisted Chemical Vapor Deposition (CVD) enable lateral MoS2 growth at atmospheric pressure and low temperatures. This work analyzes salt-based precursor-driven CVD-grown MoS2 FETs at various temperatures. MoS2 was grown using Ammonium Molybdate salt, sulfurized at 750∘C, and transferred onto p−Si3/SiO2 substrates. At room temperature, threshold voltage (VT) ranged from -35 V to -25 V, with a peak drain current of 1.2μA/μm. As temperature increased above 325K, VT shifted exponentially, and carrier mobility dropped significantly. At 400 K, the gate lost channel control, though gate leakage current remained low. These results are compared with non-salt-based MoS2 growth to assess salt precursor effects.Notice: This manuscript has been authored by UT-Battelle, LLC, under contract DE-AC05-00OR22725 with the US Department of Energy (DOE). The US government retains and the publisher, by accepting the article for publication, acknowledges that the US government retains a nonexclusive, paid-up, irrevocable, worldwide license to publish or reproduce the published form of this manuscript, or allow others to do so, for US government purposes. DOE will provide public access to these results of federally sponsored research in accordance with the DOE Public Access Plan (https://www.energy.gov/doe-public-access-plan).

Jones, Andrew [ORNL] (ORCID:0009000233849687)

Multi-sensor Observations of a High-velocity Fireball over the South Atlantic on 2026 April 1

A high-velocity fireball was detected over the South Atlantic (41.9°S, 54.7°W) on 2026 April 1 at 02:13:14 UTC by U.S. Government (USG) sensors, with peak brightness at 90.5 km altitude. The event was well-observed from geostationary orbit by two civilian lightning imagers with near-orthogonal geometry, the Geostationary Operational Environmental Satellite-East Geostationary Lightning Mapper (GLM) and the Exploitation of Meteorological Satellites Meteosat Third Generation Imager 1 Lightning Imager, and it produced low-frequency acoustic signatures. The GLM measured a total radiated energy of 2.4 × 10 10 J, corresponding to a calculated impact energy of 0.086 kt TNT equivalent. Stereoscopic triangulation of the imager tracks yields an independent pre-atmospheric velocity of ~57 km s −1 , some 18% below the USG-reported value. Because orbital provenance is acutely sensitive to the entry velocity, whose reported uncertainty could be substantial, this discrepancy is notable. We document the multi-sensor record and identify the analysis required to assess provenance, which a subsequent study will present.

Silber, Elizabeth Allaryce [Sandia National Lab. (

Influence of controlled disorder on the dipolar spin-ice state of Ho-based pyrochlores

Pyrochlore magnets of the form 𝑅 2 ⁢𝐵 2 ⁢O 7 , in which rare-earth ions on the 𝑅 site form a three-dimensional network of corner-sharing tetrahedra, provide a canonical setting for geometrical frustration. Ho-based pyrochlores host a dipolar spin-ice ground state, characterized by Ising moments constrained by the ice rules and elementary excitations analogous to magnetic monopoles. Here, in this work, we examine how controlled chemical disorder influences this state by introducing site mixing on the nonmagnetic 𝐵 site in two compounds. Ho 2⁢ GaSbO 7 contains only Ga 3+ /Sb 5+ charge disorder, whereas Ho 2 ⁢ScSbO 7 exhibits both charge and substantial size disorder arising from the large ionic-radius mismatch between Sc 3+ and Sb 5+ . Although both materials retain the pyrochlore structure, neutron-scattering measurements reveal a reduced correlation length for the 𝑅/𝐵-site cation ordering and enhanced local structural distortions in Ho 2 ⁢ScSbO 7 . Despite these structural differences, bulk thermodynamic measurements and magnetic diffuse scattering demonstrate that both systems exhibit the defining signatures of a dipolar spin-ice state. Low-energy inelastic neutron spectroscopy further uncovers broad magnetic excitations that develop within the dipolar spin-ice regime, a feature absent in pristine Ho pyrochlores and indicative of disorder-induced splitting of the non-Kramers ground-state doublet. Together, these results show that controlled disorder generates tunable transverse-field-driven quantum fluctuations in Ho-based pyrochlores, although the dipolar spin-ice state is remarkably robust to this disorder.

magnetic anisotropy

Sparse non-Markovian Noise Modeling of Transmon-Based Multi-Qubit Operations

The influence of noise on quantum dynamics is one of the main factors preventing current quantum processors from performing accurate quantum computations. Sufficient noise characterization and modeling can provide key insights into the effect of noise on quantum algorithms and inform the design of targeted error protection protocols. However, constructing effective noise models that are sparse in model parameters, yet predictive can be challenging. In this work, we present an approach for effective noise modeling of multi-qubit operations on transmon-based devices. Through a comprehensive characterization of seven devices offered by the IBM Quantum Platform, we show that the model can capture and predict a wide range of single- and two-qubit behaviors, including non-Markovian effects resulting from spatiotemporally correlated noise sources. The model’s predictive power is further highlighted through multi-qubit dynamical decoupling demonstrations and an implementation of the variational quantum eigensolver. As a training proxy for the hardware, we show that the model can predict expectation values within a relative error of 0.5%; this is a sevenfold improvement over default hardware noise models. Through these demonstrations, we highlight key error sources in superconducting qubits and illustrate the utility of reduced noise models for predicting hardware dynamics.

open quantum systems & decoherence

Graphitization of banana peel – synergistic effect of Ca and other elements of the biomass on carbon structure transformation

Recently, growing interest has been observed in using calcium as an effective catalyst for graphitizing disordered carbon materials derived from biomass and bioprecursors, such as cellulose and lignin, at relatively low tem peratures (<2000 °C). Herein, it is demonstrated that in the presence of elements such as calcium, silicon, sulfur, and potassium, which are abundant in banana peel composed of cellulose and lignin base, graphitization occurs much more effectively. The addition of external calcium source, in the form of calcium carbonate, to the banana peel-derived carbon and heat-treatment of the mixture under protective atmosphere up to 1750 °C results in a significant increase in the degree of the graphitization order and local formation of graphite crystals, whereas the same preparation and heat-treatment procedure applied to carbon material derived from pristine cellulose does not lead to such effective graphitization. It is indicated that the graphitization of banana peel carbon occurs due to the synergistic effect of the external Ca-based catalyst and the internal elements present in the raw biomass, supporting the transformation of the disordered graphene-like layers into graphitic structures. These findings are important for the future production of green graphite from biomass

09 BIOMASS FUELS

Enhanced Catalytic Dechlorination of Polyvinyl Chloride (PVC) and H2 Production Enabled by Synergistic Gaδ+/Ga0 Active Sites in Liquid Metal Particles

Polyvinyl chloride (PVC) is ubiquitous yet challenging to recycle due to its tendency to thermally decompose above 250 °C, releasing toxic, corrosive chlorinated compounds, and its inability to melt. Here, we report a catalytic strategy for PVC upcycling at 160 °C using gallium liquid metal particles (Ga-LMP) featuring a dynamic Ga-GaOOH core–shell architecture. These catalysts enable concurrent dechlorination and hydrogen evolution, yielding up to 7% H2 (based on initial hydrogen atoms in PVC) along with a highly dechlorinated (>95%) carbonaceous solid and aqueous HCl. Mechanistic investigations combining X-ray photoelectron spectroscopy, infrared spectroscopy, solid-state NMR, inelastic neutron scattering, and ab initio molecular dynamics reveal a synergistic interplay between Gaδ+ sites in the GaOOH shell and metallic Ga0 in the core. Cationic Ga initiates C–Cl bond activation and HCl formation, while progressive reduction of the shell exposes Ga0 sites that promote C–H activation and H2 evolution. Control experiments with a Ga salt and bulk Ga liquid metal confirmed that neither oxidation state alone can achieve both transformations efficiently. This work establishes a dynamic dual-site paradigm for liquid metal catalysis, in which the in situ evolution and coexistence of oxidized and metallic species enable sequential and cooperative bond activation pathways. These findings provide a general design principle for novel liquid metal catalysts that target challenging polymer transformations under mild conditions.

Zingg, Benjamin [ORNL] (ORCID:0009000914530153)

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