Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “acoustical properties”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Shivalik Plasma Device-I, a glow discharge device to study the collective dynamics of dusty plasma

Using a uniquely configured glow discharge-based Shivalik Plasma Device-I, we demonstrate a variety of collective phenomena in dusty plasma away from the glow discharge region. The cylindrical glass device produces plasma using parallel disc-shaped electrodes with a smaller anode size than the cathode. The dust microparticles are initially sprinkled over the grounded cathode. These particles acquire a significant negative charge upon plasma formation, resulting in their levitation due to the balance between the Coulomb force and gravity. The new device supports the levitation of a big-sized (10 × 8 × 5 cm 3 ) three-dimensional dust cloud over the glass surface. It contrasts the dusty plasma formations in-between electrodes reported earlier. As the discharge voltage varies from high to low, the dust cloud travels from over the glass surface to between the electrodes. A complex interplay of dust void over the cathode, a sharp density gradient, and gravity lead to self-excitation of collective dust phenomena. It includes dust density waves (phase velocity, v ph ~4 cm/s), dust cloud oscillation (frequency, f = 5 Hz), sheared flow (flow velocity, v f ~ 1 cm/s), and multiple-sized dust vortices. These dust vortices provided an excellent platform for studying turbulent mixing phenomena. The power spectrum analysis agreed with two-dimensional Kolmogorov power-law scaling. This is an ideal dusty plasma apparatus where we can create or move the dust cloud to a location of choice from the glass surface to in-between the electrodes and excite one among many collective dust dynamics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Operando Acoustic Monitoring of SEI Formation and Long-Term Cycling in NMC/SiGr Composite Pouch Cells

Stable long-term cycling and solid-electrolyte-interphase (SEI) formation are key challenges in the design of Si/graphite composites as Li-ion battery (LIB) anode materials. Typically, these long-term cycling properties are examined in flooded half-cell settings making use of a Li-metal counter electrode and a Si/graphite working electrode. This form factor has the advantage of offering an unlimited supply of Li-ions and electrolyte, thus isolating performance degradation to the passivation of the working electrode. However, half-cell studies are ineffective in revealing performance and degradation mechanisms of the Si/graphite composite in a more commercially realistic full cell setting. This paper outlines an operando acoustic technique that can offer insights on SEI formation and capacity degradation of Si/graphite composites in a full cell setting. Through a combination of electrochemical and chemical analyses, we show that increasing passivation of the silicon particles in the Si/graphite composite anode is correlated with an increase in the acoustic time-of-flight shift. We further show that temporary loss of the acoustic signal during the first cycle is associated with significant gassing of the cell. The operando acoustic technique outlined here is low-cost, simple to setup and has the potential for localized resolution, indicating usefulness in commercial-scale Si/graphite cell quality control and diagnosis.

42 ENGINEERING↗
DOE OSTI2021

Using Dense Array Waveform Correlations to Build a Velocity Model with Stochastic Variability

In an ongoing effort to improve 3D seismic-wave propagation modeling for frequencies up to 10 Hz, we used cross correlations between vertical-component waveforms from an underground chemical explosion to estimate the statistical properties of small-scale velocity heterogeneities. The waveforms were recorded by a dense 2D seismic array deployed during the Source Physics Experiments for event number 5 (SPE-5) in a series of six underground chemical explosions, conducted at the Nevada National Security Site. The array consisted of 996 geophones with a 50–100 m grid spacing, deployed at the SPE site at the north end of the Yucca Flat basin. The SPE were conducted to investigate the generation and propagation of seismic and acoustic waves from underground explosions. Comparisons of decay rates of waveform cross correlations as function of interstation distance, computed for observed and synthetic seismograms from the SPE-5 chemical explosion, were used to constrain statistical properties of correlated stochastic velocity perturbations representing small-scale heterogeneities added to a geology-based velocity model of the Yucca Flat basin. Using comparisons between recorded and simulated waveform cross correlations, we were able to recover sets of statistical properties of small-scale velocity perturbations in the velocity model that produce the best-fit between the recorded and simulated ground motion. The stochastic velocity fluctuations in the velocity model that produced the smallest misfits have a horizontal correlation distance of between 400 and 800 m, a vertical correlation distance between 100 and 200 m, and a standard deviation of 10% from the nominal model velocity in the alluvium basin layers. They also have a horizontal correlation distance of 1000 m, a vertical correlation distance of 250 m, and a standard deviation of 6% in the underlying and consolidated sedimentary layers, up to a depth of 4 km. Comparisons between observed and simulated wavefields were used to assess the proposed small-scale heterogeneity enhancements to the Yucca Flat basin model. We found that adding a depth-resolved stochastic variability to the geology-based velocity model improves the overall performance of ground-motion simulations of an SPE-5 explosion in the modeled frequency range up to 10 Hz. The results may be applicable to other similar basins.

58 GEOSCIENCES↗
DOE OSTI2022

Persistent hot carrier diffusion in boron arsenide single crystals imaged by ultrafast electron microscopy

Cubic boron arsenide (BAs) is promising for microelectronics thermal management because of its high thermal conductivity. Recently, its potential as an optoelectronic material is also being explored. However, it remains challenging to measure its photocarrier transport properties because of small sizes of available high-quality crystals. Here, we use scanning ultrafast electron microscopy (SUEM) to directly visualize the diffusion of photocarriers in BAs single crystals. Surprisingly, we observed ambipolar diffusion at low optical fluence with persistent hot carrier dynamics for above 200 ps, which can likely be attributed to the large frequency gap between acoustic and optical phonons, the same feature that is responsible for the high thermal conductivity. At higher optical fluence, we observed spontaneous electron-hole separation. Our results show BAs is an attractive optoelectronic material combining high thermal conductivity and excellent photocarrier transport properties. Finally, our study also demonstrates the capability of SUEM to probe photocarrier transport in emerging materials.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Extended linear-in- T resistivity due to electron-phason scattering in moiré superlattices

Due to its incommensurate nature, moiré superlattices host not only acoustic phonons but also another type of soft collective modes called phasons. Here, we investigate the impact of electron-phason scattering on the transport properties of moiré systems. We show that the resistivity can scale linearly with temperature down to temperatures much lower than the Bloch-Grüneisen scale defined by electron kinematics on the Fermi surface. This result stems from the friction between layers, which transfers phason spectral weight to a broad diffusive low-energy peak in the mechanical response of the system. As a result, phason scattering becomes a very efficient channel for entropy production at low temperatures. We also consider the contributions of phasons to thermodynamic properties at low temperatures and find a “metallic-like” linear-in-T behavior for the specific heat, despite the fact that this behavior is due to mechanical and not electronic degrees of freedom. Finally, we discuss the implications of this finding to reports of linear-in-T resistivity in the phase diagram of twisted bilayer graphene.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Mitigating the noise of DESI mocks using analytic control variates

In order to address fundamental questions related to the expansion history of the Universe and its primordial nature with the next generation of galaxy experiments, we need to model reliably large-scale structure observables such as the correlation function and the power spectrum. Cosmological N-body simulations provide a reference through which we can test our models, but their output suffers from sample variance on large scales. Fortunately, this is the regime where accurate analytic approximations exist. To reduce the variance, which is key to making optimal use of these simulations, we can leverage the accuracy and precision of such analytic descriptions using Control Variates (CV). The power of control variates stems from utilizing inexpensive but highly correlated surrogates of the statistics one wishes to measure. The stronger the correlation between the surrogate and the statistic of interest, the larger the variance reduction delivered by the method. We apply two control variate formulations to mock catalogs generated in anticipation of upcoming data from the Dark Energy Spectroscopic Instrument (DESI) to test the robustness of its analysis pipeline. Our CV-reduced measurements offer a factor of 5-10 improvement in the measurement error compared with the raw measurements. We explore the relevant properties of the galaxy samples that dictate this reduction and comment on the improvements we find on some of the derived quantities relevant to Baryon Acoustic Oscillation (BAO) analysis.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Measurements of sound propagation in Mars' lower atmosphere

Acoustics has become extraterrestrial and Mars provides a new natural laboratory for testing sound propagation models compared to those ones on Earth. Owing to the unique combination of a microphone and two sound sources, the Ingenuity helicopter and the SuperCam laser-induced sparks, the Mars 2020 Perseverance rover payload enables the in situ characterization of unique sound propagation properties of the low-pressure CO 2 -dominated Mars atmosphere. In this study, we show that atmospheric turbulence is responsible for a large variability in the sound amplitudes from laser-induced sparks. This variability follows the diurnal pattern of turbulence. In addition, acoustic measurements acquired over one Martian year reveal a variation of the sound intensity by a factor of 1.8 from a constant source due to the seasonal cycle of pressure and temperature that significantly modifies the acoustic impedance and shock-wave formation. Finally, we show that the evolution of the Ingenuity tones and laser spark amplitudes with distance is consistent with one of the existing sound absorption models, which is a key parameter for numerical simulations applied to geophysical experiments on CO 2 -rich atmospheres. Overall, these results demonstrate the potential of sound propagation to interrogate the Mars environment and will therefore help in the design of future acoustic-based experiments for Mars or other planetary atmospheres such as Venus and Titan.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Electrical Conductivity and Sound Velocities of Talc Under High Pressure and High Temperature Conditions and Application to the Subducting Cocos Plate

Talc is expected to be an important water carrier in Earth's upper mantle, and understanding its electrical and seismic properties under high pressure and temperature conditions is required to detect possible talc-rich regions in subduction zones imaged using geophysical observations. Here, we conducted acoustic and electrical experiments on natural talc aggregates at relevant pressure-temperature conditions. Compressional wave velocity (V p ) was measured using ultrasonic interferometry in a Paris-Edinburgh press at pressures up to 3.4 GPa and temperatures up to 873 K. Similar Vp values are obtained regardless of the initial crystallographic preferred orientation of the samples, which can be explained by talc grain reorientation during the experiment, with the (001) plane becoming perpendicular to the uniaxial compression axis. Electrical conductivity of the same starting material was determined using impedance spectroscopy in a multi-anvil press up to 6 GPa and 1263 K. Two conductivity jumps are observed, at ∼860–1025 K and ∼940–1080 K, depending on pressure, and interpreted as talc dehydroxylation and decomposition, respectively. Electrical anisotropy is observed at low temperature and decreases with increasing pressure (∼10 at 1.5 GPa and ∼2 at 3.5 GPa). Comparison of acoustic and electrical results with geophysical observations in central Mexico supports the presence of a talc-bearing layer atop the subducted Cocos plate.

Cocos plate↗
DOE OSTI2021

Dispersion and the speed-limited particle-in-cell algorithm

This paper discusses temporally continuous and discrete forms of the speed-limited particle-in-cell (SLPIC) method first treated by Werner et al. [Phys. Plasmas 25, 123512 (2018)]. The dispersion relation for a 1D1V electrostatic plasma whose fast particles are speed-limited is derived and analyzed. By examining the normal modes of this dispersion relation, we show that the imposed speed-limiting substantially reduces the frequency of fast electron plasma oscillations while preserving the correct physics of lower-frequency plasma dynamics (e.g. ion acoustic wave dispersion and damping). We then demonstrate how the timestep constraints of conventional electrostatic particle-in-cell methods are relaxed by the speed-limiting approach, thus enabling larger timesteps and faster simulations. Here, these results indicate that the SLPIC method is a fast, accurate, and powerful technique for modeling plasmas wherein electron kinetic behavior is nontrivial (such that a fluid/Boltzmann representation for electrons is inadequate) but evolution is on ion timescales.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Vibrational properties and thermal transport in quaternary chalcogenides: The case of Te-based compositions

Vibrational thermal properties of CuZn 2 InTe 4 , AgZn 2 InTe 4 , and Cu 2 CdSnTe 4 , derived from binary II-VI zinc-blendes, are reported based on first-principles calculations. While the chalcogenide atoms in these materials have the same lattice positions, the cation atom arrangements vary, resulting in different crystal symmetries and subsequent properties. The compositional differences have important effects on the vibrational thermal characteristics of the studied materials, which demonstrate that low-frequency optical phonons hybridize with acoustic phonons and lead to enhanced phonon-phonon scattering and low lattice thermal conductivities. The phonon density of states, mode Grüneisen parameters, and phonon scattering rates are also calculated, enabling deeper insight into the microscopic thermal conduction processes in these materials. Compositional variations drive differences among the three materials considered here; nonetheless, their structural similarities and generally low thermal conductivities (0.5–4 W/m K at room temperature) suggest that other similar II-VI zinc-blende derived materials will also exhibit similarly low values, as also corroborated by experimental data. Finally, this, combined with the versatility in designing a variety of motifs on the overall structure, makes quaternary chalcogenides interesting for thermal management and energy conversion applications that require low thermal conductivity.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Magnetostrictive Ultrasonic Waveguide Transducer for In-Pile Thermometry

Real-time and reliable temperature measurement on nuclear fuels is crucial for the safe operation of existing pressurized-water reactors and future advanced nuclear reactors. Magnetostrictive materials deform when subjected to a magnetic field or exhibit magnetization variation when stressed. Here, based on these properties, this study prototyped an ultrasonic thermometer (UT) consisting of a magnetostrictive waveguide, a dc coil providing appropriate magnetic biasing, and an ac coil generating an acoustic impulse and detecting the resulting acoustic echoes. By tracking the time of flight between the excitation pulse and the echoes, the UT can potentially detect nuclear fuel cladding temperature from a long distance. In this study, magnetostrictive iron–gallium alloys, or Galfenol, were selected as the waveguide due to their large magnetostriction, superior temperature survivability, excellent radiation resilience, and high mechanical robustness. A multiphysics finite-element model considering electrical, magnetic, and mechanical dynamics in the magnetostrictive UT was then developed. The model exhibited an error of 0.24% in time-of-flight simulation and, therefore, enabled computer-aided design and guided signal processing. Between room temperature and 120 °C, the new Galfenol-based UT exhibits a linear sensitivity of 162.8×10 ₋6 °C ₋1 , which is 51.7% higher than a previous magnetostrictive UT based on iron–cobalt–vanadium alloys.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2023

Single crystal growth and thermoelectric properties of Nowotny chimney-ladder compound Fe 2 Ge 3

Fe 2 Ge 3 with an incommensurate Nowotny chimney-ladder (NCL) structure is a promising material for thermoelectric applications due to its low thermal conductivity. Previous experimental studies on Fe 2 Ge 3 have mainly focused on polycrystalline samples, resulting in a limited understanding of the material's intrinsic thermoelectric properties and the underlying causes of its low thermal conductivity. Here we report the synthesis and thermoelectric properties of single crystalline Fe 2 Ge 3 . Millimeter-sized Fe 2 Ge 3 single crystals grown by the chemical vapor transport method enable the study of the intrinsic thermoelectric properties. The Seebeck coefficient of Fe 2 Ge 3 is negative and its magnitude increases linearly with temperature, showing a degenerate n-type semiconductor behavior. Analysis of the electrical resistivity and specific heat data indicates the existence of an Einstein mode with a characteristic temperature of about 60 K, suggesting the presence of low-energy optical phonons. Further, the thermal conductivity of Fe 2 Ge 3 along the c axis is as low as 1.9Wm -1 K -1 at 300 K and exhibits a nearly temperature-independent characteristic, which is distinct from the previous theoretical calculations with a stronger temperature dependence. The low thermal conductivity may be attributed to the scattering of acoustic phonons by low-energy optical modes and the presence of non-extended diffuson modes, as reported in another NCL compound, MnSi 1.74 . This study provides valuable insights into the electrical and thermal properties of Fe 2 Ge 3 , which can open up possibilities for future advances in thermoelectric applications.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Theoretical Optimization of Trapped-Bubble-Based Acoustic Metamaterial Performance

Acoustic metamaterials have proven to be a versatile tool for the precise control and manipulation of sound waves. One of the promising designs of acoustic metamaterials employ the arrays of bubbles and find applications for soundproofing, blast mitigation, and many others. An obvious advantage of bubble-based metamaterials is their ability to be relatively thin while absorbing low-frequency sound waves. The vast majority of theories developed to predict resonant behavior of bubble-based metamaterials capitalize on Minnaert frequency. Here, we propose a novel theoretical approach to characterize bubble-based metamaterials that are based on our previous findings for a single bubble trapped in circular cavity modeled as a thin clamped plate. We obtain analytical expressions for resonant frequencies of bubble metascreens using self-consistent approximation. Two geometry factors, distance between bubble centers and distance between bubble center and interface of acoustic impedance change, are taken into account. We demonstrate the existence of multiple bandgaps and possibility of switching between them via adjustment of geometry parameters and reflector properties.

Gritsenko, Dmitry (ORCID:0000000248181512)↗
DOE OSTI2024

Atomic-scale visualization of defect-induced localized vibrations in GaN

Phonon engineering is crucial for thermal management in GaN-based power devices, where phonon-defect interactions limit performance. However, detecting nanoscale phonon transport constrained by III-nitride defects is challenging due to limited spatial resolution. Here, we used advanced scanning transmission electron microscopy and electron energy loss spectroscopy to examine vibrational modes in a prismatic stacking fault in GaN. By comparing experimental results with ab initio calculations, we identified three types of defect-derived modes: localized defect modes, a confined bulk mode, and a fully extended mode. Additionally, the PSF exhibits a smaller phonon energy gap and lower acoustic sound speeds than defect-free GaN, suggesting reduced thermal conductivity. Our study elucidates the vibrational behavior of a GaN defect via advanced characterization methods and highlights properties that may affect thermal behavior.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Stability and pulsation of the first dark stars

ABSTRACT The first bright objects to form in the Universe might not have been ‘ordinary’ fusion-powered stars, but ‘dark stars’ (DSs) powered by the annihilation of dark matter (DM) in the form of weakly interacting massive particles (WIMPs). If discovered, DSs can provide a unique laboratory to test DM models. DSs are born with a mass of the order of M⊙ and may grow to a few million solar masses; in this work we investigate the properties of early DSs with masses up to $\sim \! 1000 \, \mathrm{ M}_\odot$, fueled by WIMPS weighing 100 GeV. We improve the previous implementation of the DM energy source into the stellar evolution code mesa. We show that the growth of DSs is not limited by astrophysical effects: DSs up to $\sim \!1000 \, \mathrm{ M}_{\odot }$ exhibit no dynamical instabilities; DSs are not subject to mass-loss driven by super-Eddington winds. We test the assumption of previous work that the injected energy per WIMP annihilation is constant throughout the star; relaxing this assumption does not change the properties of the DSs. Furthermore, we study DS pulsations, for the first time investigating non-adiabatic pulsation modes, using the linear pulsation code gyre. We find that acoustic modes in DSs of masses smaller than $\sim \! 200 \, \mathrm{ M}_\odot$ are excited by the κ − γ and γ mechanism in layers where hydrogen or helium is (partially) ionized. Moreover, we show that the mass-loss rates potentially induced by pulsations are negligible compared to the accretion rates.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Method to account for natural fracture induced elastic anisotropy in geomechanical characterization of shale gas reservoirs

Shale has been usually recognized as a transverse isotropic (TI) medium in conventional geomechanical log interpretation due to its laminated nature. However, when natural fractures exist in the shale rock, additional elastic anisotropy is introduced, converting laminated Shale to an orthorhombic (OB) medium. Previous studies illustrate that neglecting the natural fracture induced anisotropy in shale geomechanical log interpretation could lead to inaccurate evaluations of elastic moduli and in-situ stresses. In this paper, a new method is developed to account for the natural fracture induced anisotropy in geomechanical log interpretation based upon the TI acoustic model developed by the author and a characterization technique of elastic wave anisotropy (Sayers, 1991). The new OB model incorporates the four acoustic log data inputs and five modeling constraints in a nonlinear optimization algorithm to solve for the nine independent stiffness coefficients of an OB rock, and further to solve for the geomechanical properties and in-situ stress profiles in an OB formation. The new method was validated with a Marcellus Gas Shale field case. Both the new OB model and the conventional TI model were applied to interpret the minimum horizontal stress profile for the same formation. By comparing the results, the OB model is more robust than the TI from two aspects. First, the average stress magnitude predicted by the OB model is closer to the one measured by the Diagnostic Fracture Injection Test (DFIT). Second, the OB model predicts a more obvious stress barrier between the lower Marcellus and upper Onondaga Limestone than the TI model does. Finally, the predicted stress barrier is consistent with the observation of the microseismic events of a horizontal well drilled and completed nearby, which reveals that no hydraulic fracture propagates downward through the bottom boundary of Marcellus Shale into the underlying Onondaga Limestone.

03 NATURAL GAS↗
DOE OSTI2021

A surface acoustic wave hydrogen sensor with tin doped indium oxide layers for intermediate temperatures

We report high temperature surface acoustic wave (SAW) gas sensors with conducting sensing layers require tuning of the sheet conductivity for optimal response. Conducting metal oxides are attractive sensing materials for their tunable electronic properties and high thermal stability, amongst others. Here, we have investigated the application of indium oxide (IO) and indium tin oxide (ITO) films on langasite (LGS)-based SAW reflective delay line sensor devices for monitoring hydrogen at 350 °C. Specifically, we modeled the effect of the IO and ITO sensing layer thickness on the wave velocity, attenuation, and effective electromechanical coefficient. This was followed by an experimental demonstration of tuning of the ITO film sheet conductivity by controlling the dopant concentration and yielding an improvement in the sensor sensitivity. The current study provides a pathway towards the development of conductivity-based sensing layers for high temperature SAW gas sensors with improved sensitivity.

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

Higher-order phonon scattering: advancing the quantum theory of phonon linewidth, thermal conductivity and thermal radiative properties

Phonon scattering plays a central role in the quantum theory of phonon linewidth, which in turn governs important properties including infrared spectra, Raman spectra, lattice thermal conductivity, thermal radiative properties, and also signifi- cantly affects other important processes such as hot electron relaxation. Since Maradudin and Fein’s classic work in 1962, three-phonon scattering had been considered as the dominant intrinsic phonon scattering mechanism and has seen tremendous advances. However, the role of the higher-order four-phonon scattering had been persistently unclear and so was ignored. The tremendous complexity of the formalism and computational challenges stood in the way, prohibiting the direct and quantitative treatment of four-phonon scattering. In 2016, a rigorous four-phonon scattering formalism was developed, and the prediction was realized using empirical potentials. In 2017, the method was extended using first-principles calculated force constants, and the thermal conductivities of boron arsenides (BAs), Si and diamond were predicted. The predictions for BAs were later confirmed by several independent experiments. Four-phonon scattering has since been investigated in a range of materials and established as an important intrinsic scattering mechanism for thermal transport and radiative properties. Specifically, four-phonon scattering is important when the fourth-order scattering potential or phase space becomes relatively large. The former scenario includes: (i) nearly all materials when the temperature is high; (ii) strongly anharmonic (low thermal conductivity) materials, including most rocksalt compounds, halides, hydrides, chalcogenides and oxides. The latter scenario includes: (iii) materials with large acoustic–optical phonon band gaps, such as XY compounds with a large atomic mass ratio between X and Y; (iv) two- dimensional materials with reflection symmetry, such as single-layer graphene, single-layer boron nitride and carbon nanotubes; and (v) phonons with a large density of states, such as optical phonons, which are important for Raman, infrared and thermal radiative properties. Four-phonon scattering is expected to gain broad interest in various technologically important materials for thermoelectrics, thermal barrier coatings, thermal energy storage, phase change, nuclear power, ultra-high temperature ceramics, infrared spectra, Raman spectra, radiative transport, hot electron relaxation and radiative cooling. Four-phonon scattering has been, and will continue to be, established as an important intrinsic phonon scattering mechanism beyond three-phonon scattering. The prediction of four-phonon scattering will transition from a breakthrough to a new routine in the next decade.

Feng, Tianli↗