Search NASA⌕ Search

SEARCH · Search NASA

Results for “Range expansion”

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 217 records · Page 12

Role of effective mass and long-range interactions in the band-gap renormalization of photoexcited semiconductors

Understanding how to control changes in the electronic structure and related dynamical renormalizations by external driving fields is the key for understanding ultrafast spectroscopy and applications in electronics. Here, we focus on the band gap's modulation by external electric fields and uncover the effect of band dispersion on the gap renormalization. We employ the Green's function formalism using the real-time Dyson expansion to account for dynamical correlations induced by photodoping. The many-body formalism captures the dynamics of systems with long-range interactions, carrier mobility, and variable electron and hole effective mass. We also demonstrate that mean-field simulations based on the Hartree-Fock Hamiltonian, which lacks dynamical correlations, yields a qualitatively incorrect picture of band-gap renormalization. We find the trend that increasing effective mass, thus decreasing mobility, leads to as much as a 6% enhancement in band-gap renormalization. Further, the renormalization is strongly dependent on the degree of photodoping. As the screening induced by free electrons and holes effectively reduces any long-range and interband interactions for highly excited systems, we show that there is a specific turnover point with a minimal band gap. Here, we further demonstrate that the optical gap renormalization follows the same trend though its magnitude is altered by the Moss-Burstein effect.

Approximation methods for many-body systems↗

Nodal Capacity Expansion Modeling with ReEDS: A Case Study of the RTS-GMLC Test System

Test systems play an important role in power systems analysis and are used extensively for a wide range of purposes like reliability analysis, production cost modeling, studying the impacts of load shifting, and many more. These test systems have continued to evolve over the years to keep pace with changing system compositions due to technology and policy changes. In this work, we apply the Regional Energy Deployment System (ReEDS) model, a large-scale capacity expansion model (CEM), on the nodal RTS-GMLC system to perform capacity expansion of the bus-level test system. This work both demonstrates that a large-scale CEM can successfully be applied to nodal systems, and that the CEM can allow a test system to be evolved to meet desired criteria. The process allowed the coupling of CEM data not available in the test system (such as capital costs) with test system data to provide plausible system evolutions in line with scenario specifications.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Coverage-Dependent Adsorption of Hydrogen on Fe(100): Determining Catalytically Relevant Surface Structures via Lattice Gas Models

Hydrogen adatoms are a critical surface species for several reactions catalyzed by Fe surfaces such as Fischer–Tropsch synthesis, ammonia synthesis, and the hydrodeoxygenation of biomass-derived oxygenates. Parameterizing the energetics for H/Fe in terms of both coverage and configuration space can significantly aid in the development of multiscale models as well as provide atomic level insight into the dominant surface structures present under realistic reaction conditions. Here, we construct a lattice gas model for H/Fe(100), where the lateral interactions are determined from first-principles using density functional theory. Using 950 symmetrically unique H/Fe(100) configurations, we generate a cluster expansion with a predictive accuracy in terms of surface energy of 3.8 meV/site over a coverage range from 0 to 3 monolayers. Ten electronic ground state structures are identified from this thorough scan (including the structures at 0 and 3 monolayers), which were subsequently used to generate ab initio phase diagrams under a range of temperatures and pressures. Under reaction conditions typical of Fischer–Tropsch synthesis, ammonia synthesis, and biomass oxygenate hydrodeoxygenation, we find that the 1.0 monolayer structure is dominant. Furthermore, examination of the total H–H lateral interactions for the H/Fe(100) electronic ground state structures shows that H/Fe(100) can be accurately modeled via a mean-field ideal lattice gas model for coverages less than 1.0 monolayers. Altogether, this work enables the incorporation of H–H lateral interactions on Fe(100) into multiscale models, via either mean-field or site-dependent techniques, and provides atomic insight into the catalytically relevant H/Fe(100) structures for a range of heterogeneous reactions.

09 BIOMASS FUELS↗

One‐at‐a‐Time Parameter Perturbation Ensemble of the Community Land Model, Version 5.1

Comprehensive land models are subject to significant parametric uncertainty, which can be hard to quantify due to the large number of parameters and high model computational costs. We constructed a large parameter perturbation ensemble (PPE) for the Community Land Model version 5.1 with biogeochemistry configuration (CLM5.1-BGC). We performed more than 2,000 simulations perturbing 211 parameters across six forcing scenarios. This provides an expansive data set, which can be used to identify the most influential parameters on a wide range of output variables globally, by biome, or by plant functional type. We found that parameter effects can exceed scenario effects and that a small number of parameters explains a large fraction of variance across our ensemble. The most important parameters can differ regionally and also based on the forcing scenario. The software infrastructure developed for this experiment has greatly reduced the human and computer time needed for CLM PPEs, which can facilitate routine investigation of parameter sensitivity and uncertainty, as well as automated calibration.

Kennedy, Daniel [NSF National Center for Atmospher↗

An extensive assessment of the performance of pairwise and many-body interaction potentials in reproducing ab initio benchmark binding energies for water clusters n = 2–25

We assess the performance of 7 pairwise additive (TIP3P, TIP4P, TIP4P-ice, TIP5P, OPC, SPC, SPC/E) and 8 families of many-body potentials (q-AQUA, HIPPO, AMOEBA, EFP, TTM, WHBB, MB-pol, MB-UCB) in reproducing high-level ab initio benchmark values, CCSD(T) or MP2 at the complete basis set (CBS) limit for the binding energy and the many-body expansion (MBE) of water clusters n = 2–11, 16–17, 20, 25. By including a large range of cluster sizes having dissimilar hydrogen bonding networks, we obtain an understanding of how these potentials perform for different hydrogen bonding arrangements that are mostly outside of their parameterization range. While it is appropriate to compare the results of ab initio based many-body potentials directly to the electronic binding energies (D e 's), the pairwise additive ones are compared to the enthalpies at T = 298 K, ΔH(298 K), as the latter class of force fields are parametrized to reproduce enthalpies (implicitly accounting for zero-point energy corrections) rather than binding energies. We find that all pairwise additive potentials considered overestimate the reference ΔH values for the n = 2–25 clusters by >13%. For the water dimer (n = 2) in particular, the errors are in the range 83–119% for the pairwise additive potentials studied since these are based on an effective rather than the true 2-body interaction specifically designed as a means of partially accounting for the missing many-body terms. This stronger 2-body interaction is achieved by an enhanced monomer dipole moment that mimics its increase from the gas phase monomer to the condensed phase value. Indeed, for cluster sizes n ≥ 4 the percent deviations become slightly smaller (albeit all exceeding 13%). In contrast, we find that the many-body potentials perform more accurately in reproducing the electronic binding energies (D e 's) throughout the entire cluster range (n = 2–25), all reproducing the ab initio benchmark binding energies within ±7% of the respective CBS values. We further assess the ability of a subset of the many-body potentials (MB-UCB, q-AQUA, MB-pol, and TTM2.1-F) to also reproduce the magnitude of the ab initio many-body energy terms for water cluster sizes n = 7, 10, 16 and 17. The potentials show an overall good agreement with the available benchmark values. However, we identify characteristic differences upon comparing the many-body terms at both the ab initio-optimized geometries and the respective potential-optimized geometries to the reference ab initio values. Additionally, by applying this analysis to a wide range of cluster sizes, trends in the MBE of the potentials with increasing cluster size can be identified. Finally, in an attempt to draw a parallel between the pairwise additive and many-body potentials, we report the analysis of the individual molecular dipole moments for water clusters with 1 to ~4 solvation shells with the TTM2.1-F potential. We find that the internally solvated water molecules have in general a larger molecular dipole moment ranging from 2.6–3.0 D. This justifies the use of an enhanced, with respect to the gas-phase value, molecular dipole moment for the pairwise additive potentials, which is intended to fold in the many body terms into an effective (enhanced) pairwise interaction through the choice of the charges. These results have important implications for the development of future generations of efficient, transferable, and highly accurate classical interaction potentials in both the pairwise additive and many-body categories.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

AB-G0W0: A practical G0W0 method without frequency integration based on an auxiliary boson expansion

Common G0W0 implementations rely on numerical or analytical frequency integration to determine the G0W0 self-energy, which results in a variety of practical complications. Recently, we have demonstrated an exact connection between the G0W0 approximation and equation-of-motion quantum chemistry approaches [J. Tölle and G. Kin-Lic Chan, J. Chem. Phys. 158, 124123 (2023)]. Based on this connection, we propose a new method to determine G0W0 quasiparticle energies, which completely avoids frequency integration and its associated problems. To achieve this, we make use of an auxiliary boson (AB) expansion. We name the new approach AB-G0W0 and demonstrate its practical applicability in a range of molecular problems.

Chemistry↗

Scattering phase shifts from a quantum computer

We calculate two-body scattering phase shifts on a quantum computer using a leading order short-range effective field theory Hamiltonian. The algorithm combines the variational quantum eigensolver and the quantum subspace expansion. As an example, we consider scattering in the deuteron 3S1 partial wave. We calculate scattering phase shifts with a quantum simulator and on real hardware. Here, we also study how noise impacts these calculations and discuss noise mitigation required to extend our work to larger quantum processing units. With current hardware, up to five superconducting qubits can produce acceptable results, and larger calculations will require a significant noise reduction.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Nuclear effective field theory: Status and perspectives

The nuclear physics landscape has been redesigned as a sequence of effective field theories (EFTs) connected to the Standard Model through symmetries and lattice simulations of Quantum Chromodynamics (QCD). EFTs in this sequence are expansions around different low-energy limits of QCD, each with its own characteristics, scales, and ranges of applicability regarding energy and number of nucleons. We review each of the three main nuclear EFTs---Chiral, Pionless, Halo/Cluster---highlighting their similarities, differences, and connections. Here, in doing so, we survey the structural properties and reactions of nuclei that have been derived from the \abinitio solution of the few- and many-body problem built upon EFT input.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Climate Change Mitigation Potential of Wind Energy

Global wind resources greatly exceed current electricity demand and the levelized cost of energy from wind turbines has shown precipitous declines. Accordingly, the installed capacity of wind turbines grew at an annualized rate of about 14% during the last two decades and wind turbines now provide ~6–7% of the global electricity supply. This renewable electricity generation source is thus already playing a role in reducing greenhouse gas emissions from the energy sector. Here we document trends within the industry, examine projections of future installed capacity increases and compute the associated climate change mitigation potential at the global and regional levels. Key countries (the USA, UK and China) and regions (e.g., EU27) have developed ambitious plans to expand wind energy penetration as core aspects of their net-zero emissions strategies. The projected climate change mitigation from wind energy by 2100 ranges from 0.3–0.8 °C depending on the precise socio-economic pathway and wind energy expansion scenario followed. The rapid expansion of annual increments to wind energy installed capacity by approximately two times current rates can greatly delay the passing of the 2 °C warming threshold relative to pre-industrial levels. To achieve the required expansion of this cost-effective, low-carbon energy source, there is a need for electrification of the energy system and for expansion of manufacturing and installation capacity.

17 WIND ENERGY↗

Sample Size Effects on Petrophysical Characterization and Fluid-to-Pore Accessibility of Natural Rocks

Laboratory-scale analysis of natural rocks provides petrophysical properties such as density, porosity, pore diameter/pore-throat diameter distribution, and fluid accessibility, in addition to the size and shape of framework grains and their contact relationship with the rock matrix. Different types of laboratory approaches for petrophysical characterization involve the use of a range of sample sizes. While the sample sizes selected should aim to be representative of the rock body, there are inherent limitations imposed by the analytical principles and holding capacities of the different experimental apparatuses, with many instruments only able to accept samples at the μm–mm scale. Therefore, a total of nine (three limestones, three shales, two sandstones, and one dolomite) samples were collected from Texas to fill the knowledge gap of the sample size effect on the resultant petrophysical characteristics. The sample sizes ranged from 3 cm cubes to <75 μm particles. Using a combination of petrographic microscopy, helium expansion pycnometry, water immersion porosimetry, mercury intrusion porosimetry, and (ultra-) small-angle X-ray scattering, the impact of sample size on the petrophysical properties of these samples was systematically investigated here. The results suggest that the sample size effect is influenced by both pore structure changes during crushing and sample size-dependent fluid-to-pore connectivity.

36 MATERIALS SCIENCE↗

BACCHUS Analysis of Weak Lines in APOGEE Spectra (BAWLAS)

Elements with weak and blended spectral features in stellar spectra are challenging to measure and require specialized analysis methods to precisely measure their chemical abundances. In this work, we have created a catalog of approximately 120,000 giants with high signal-to-noise Apache Point Observatory Galactic Evolution Experiment (APOGEE) Data Release 17 (DR17) spectra, for which we explore weak and blended species to measure Na, P, S, V, Cu, Ce, and Nd abundances and 12 C/ 13 C isotopic ratios. We employ an updated version of the Brussels Automatic Code for Characterizing High-accuracy Spectra (BACCHUS) code to derive these abundances using the stellar parameters measured by APOGEE's DR17 Stellar Parameters and Chemical Abundances Pipeline, quality flagging to identify suspect spectral lines, and a prescription for upper limits. Combined, these allow us to provide our BACCHUS Analysis of Weak Lines in APOGEE Spectra catalog of precise chemical abundances for these weak and blended species, which agrees well with the literature and improves upon APOGEE abundances for these elements, some of which are unable to be measured with APOGEE's current, grid-based approach without computationally expensive expansions. This new catalog can be used alongside APOGEE and provides measurements for many scientific applications ranging from nuclear physics to Galactic chemical evolution and Milky Way population studies. To illustrate this we show some examples of uses for this catalog, such as showing that we observe stars with enhanced s-process abundances or that we can use the 12 C/ 13 C ratios to explore extra mixing along the red giant branch.

79 ASTRONOMY AND ASTROPHYSICS↗

Femtoscopic correlations of identical charged pions and kaons in $pp$ collisions at $\sqrt{s}$ = 13 TeV with event-shape selection

Collective behavior has been observed in high-energy heavy-ion collisions for several decades. Collectivity is driven by the high particle multiplicities that are produced in these collisions. At the CERN Large Hadron Collider (LHC), features of collectivity have also been seen in high-multiplicity proton-proton collisions that can attain particle multiplicities comparable to peripheral Pb-Pb collisions. One of the possible signatures of collective behavior is the decrease of femtoscopic radii extracted from pion and kaon pairs emitted from high-multiplicity collisions with increasing pair transverse momentum. This decrease can be described in terms of an approximate transverse mass scaling. In the present work, femtoscopic analyses are carried out by the ALICE Collaboration on charged pion and kaon pairs produced in pp collisions at $\sqrt{s}$=13TeV from the LHC to study possible collectivity in pp collisions. The event-shape analysis method based on transverse sphericity is used to select for spherical versus jetlike events, and the effects of this selection on the femtoscopic radii for both charged pion and kaon pairs are studied. This is the first time this selection method has been applied to charged kaon pairs. An approximate transverse-mass scaling of the radii is found in all multiplicity ranges studied when the difference in the Lorentz boost for pions and kaons is taken into account. This observation does not support the hypothesis of collective expansion of hot and dense matter that should only occur in high-multiplicity events. A possible alternate explanation of the present results is based on a scenario of common emission conditions for pions and kaons in pp collisions for the multiplicity ranges studied.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Structure of chalcogen overlayers on Au(111): Density functional theory and lattice-gas modeling

Ordering of different chalcogens, S, Se, and Te, on Au(111) exhibit broad similarities but also some distinct features, which must reflect subtle differences in relative values of the long-range pair and many-body lateral interactions between adatoms. We develop lattice-gas (LG) models within a cluster expansion framework, which includes about 50 interaction parameters. These LG models are developed based on density functional theory (DFT) analysis of the energetics of key adlayer configurations in combination with the Monte Carlo (MC) simulation of the LG models to identify statistically relevant adlayer motifs, i.e., model development is based entirely on theoretical considerations. The MC simulation guides additional DFT analysis and iterative model refinement. Given their complexity, development of optimal models is also aided by strategies from supervised machine learning. The model for S successfully captures ordering motifs over a broader range of coverage than achieved by previous models, and models for Se and Te capture the features of ordering, which are distinct from those for S. More specifically, the modeling for all three chalcogens successfully explains the linear adatom rows (also subtle differences between them) observed at low coverages of ~0.1 monolayer. The model for S also leads to a new possible explanation for the experimentally observed phase with a (5 × 5)-type low energy electron diffraction (LEED) pattern at 0.28 ML and to predictions for LEED patterns that would be observed with Se and Te at this coverage.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Revisiting proton–proton fusion in chiral effective field theory

We calculate the S-factor for proton–proton fusion using chiral effective field theory interactions and currents. By performing order-by-order calculations with a variety of chiral interactions that are regularized and calibrated in different ways, we assess the uncertainty in the S-factor from the truncation of the effective field theory expansion and from the sensitivity of the S-factor to the short-distance axial current determined from three- and four-nucleon observables. We find that S(0) = (4.100 ± 0.024(syst) ± 0.013(stat) ± 0.008(g$_{A}$)) × 10$^{−23}$ MeV fm$^{2}$, where the three uncertainties arise, respectively, from the truncation of the effective field theory expansion, use of the two-nucleon axial current fit to few-nucleon observables and variation of the axial coupling constant within the recommended range. The increased value of S(0) compared to previous calculations is mainly driven by an increase in the recommended value for the axial coupling constant and is in agreement with a recent analysis based on pionless effective field theory.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Climate Influences on Capacity Expansion Planning with Application to the Western U.S

Electric power system planners utilize a variety of planning tools to inform decisions concerning generation and transmission additions to the electric grid, the need for operational changes, and to evaluate potential stressors on the system. Numerous factors contribute to the planning process including projected fuel and technology costs, policy and load profiles. There is also a growing recognition of the interdependency of the electric grid with other natural and engineered systems. Here we explore how future climate change and hydropower operability might influence decisions related to electricity capacity expansion planning and operations. To do so we assemble a multi-model framework. Specifically, water resource modeling is used to simulate climate impacts on future water supply for thermoelectric and hydropower generation. Separately, temperature impacts on electricity load are evaluated. Together, these climate factors spatially constrain a capacity expansion model that projects generation and transmission additions to the grid. The projected new capacity-builds are then evaluated on their operations, reliability, and cost under average and extreme climate conditions using production cost modeling. This coupled framework is demonstrated on the electric grid in the Western U.S., supporting capacity expansion planning by WECC, the North American Electric Reliability Corporation (NERC) regional entity responsible for reliability assurance of the Western Interconnection. This region was selected in part because the West is unique in that it has high potential for renewable penetrations and is experiencing large retirements/displacements of baseload resources, primarily coal, leading to possible operational challenges in terms of changing resource mix and the need for resource flexibility. Toward this challenge, planning scenarios encompass a range of alternative energy, climate and drought futures. In this context we explore answers to two strategic questions: 1) How does changing climate influence electricity expansion planning (generation and transmission) and future operations, including type and capacity of new builds, system reliability, cost and environmental impacts? 2) How does the representation of hydropower in the modeling framework influence the evaluation of bulk power system operations? Results indicate that climate has a measurable influence on recommendations concerning the capacity, type and location of new generation and transmission additions, with up to 17 GW additional capacity needed by 2038 to meet peak loads (~6.6% increase over capacity-builds based on historical climate). The extent of additional infrastructure needs is strongly influenced by future water availability for hydropower and the potential deployment of demand response technologies. Systems designed for future climate conditions were found to maintain high system reliability under a range of electricity and water availability scenarios (including significant drought), with minimal system curtailments. Additional capacity needs due to higher load tend to increase cumulative 20-year investment and operating costs by $\$$5-$\$$17 billion and generation costs increase by 9 to 19%. Finally, changing the representation of hydropower flexibility has a relatively small influence on capacity expansion in the Western Interconnection through 2038, but hydropower flexibility impacts generation costs to a similar extent as climate.

13 HYDRO ENERGY↗

Temperature-Dependent Fracture Resistance of Silicon Nanopillars during Electrochemical Lithiation

During the lithation of silicon anodes, the solid-state diffusion of lithium into Li x Si follows the Arrhenius law, the resulting morphology and fracture behavior are determined by the silicon anode operation temperature. Here, we reveal the temperature dependence of the lithiation mechanics of crystalline silicon nanopillars (SiNPs) via microscopic observations of the anisotropic growth and fracture behavior. In this work, we fabricated 1D SiNP structures with various orientations ($\langle100\rangle, \langle110\rangle$, and $\langle111\rangle$) as working electrodes and operated them at temperatures ranging from –20 to 40 °C. The lithiation of crystalline silicon at low temperatures exhibited preferential volume expansion along $\langle110\rangle$ and decreased fracture resistance. Furthermore, low temperatures caused the catastrophic fracture of amorphous silicon after the second lithiation. Our findings demonstrate the importance of silicon anode temperature control to prevent mechanical fracture during the cycle of lithium-ion batteries in harsh environments (e.g., electric vehicles in winter).

25 ENERGY STORAGE↗

Autonomous Nanoparticle Synthesis Guided by In Situ Multiscale Structural Characterization

Autonomous synthesis platforms promise rapid exploration of vast parameter spaces; yet, integrating in situ structural characterization in closed-loop synthesis optimization remains challenging. We demonstrate a realization of such a closed-loop platform coupled with a droplet-flow microreactor, in situ X-ray scattering methods (SAXS/WAXS), and Gaussian process optimization to synthesize citrate-reduced Au nanoparticles with targeted characteristics. The system efficiently explored ∼19,000 synthesis recipes through 365 experiments, achieving precise control over size (4–60 nm) and polydispersity (σ < 0.11) across large citrate/gold ratios, exceeding traditional synthesis boundaries (1–10). Beyond confirming classical Turkevich–Frens trends, partial-dependence analysis revealed strong nonlinear coupling among precursor, citrate, and pH effects. Combining quantitative SAXS/WAXS analysis with electron microscopy characterization, we uncovered that crystallite size (d c ) and particle size (d) follow d c = 0.18d + β, where synthesis chemistry controls the intercept β while maintaining a universal slope. This parallel-band structure enables independent tuning of crystallite domain size at fixed particle diameter through a combination of chloride, gold precursor, citrate, and pH contributions (cross-validated Spearman ρ = 0.7 ± 0.1). High-resolution electron microscopy shows multiple lattice-fringe orientations within single particles, directly confirming polycrystalline domains and the ability to tune d c at the fixed d. The platform’s validation includes indistinguishable static versus flowing measurements, stable droplet transport at 100 °C, and <5% run-to-run variation, establishing a robust framework for mapping and controlling multiscale nanoparticle structure across expansive chemical spaces. In conclusion, the developed closed-loop platform can be applied to a borad range of nanosyntheis processes.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Infrared-active phonon modes in single-crystal thorium dioxide and uranium dioxide

The infrared-active phonon modes, in single-crystal samples of thorium dioxide (ThO 2 ) and uranium dioxide (UO 2 ), were investigated using spectroscopic ellipsometry and compared with density functional theory. Both ThO 2 and UO 2 are found to have one infrared-active phonon mode pair [consisting of one transverse optic (TO) and one associated longitudinal optic (LO) mode], which is responsible for the dominant features in the ellipsometric data. Furthermore, at room temperature, our results for the mode pair’s resonant frequencies and broadening parameters are comparable with previous reflectance spectroscopy characterizations and density functional theory predictions. For ThO 2 , our ellipsometry and density function theory results both show that the LO mode broadening parameter is larger than the TO mode broadening. This signifies mode anharmonicity, which can be attributed to the intrinsic phonon–phonon interaction. In addition to the main mode pair, a broad low-amplitude impurity-like vibrational mode pair is detected within the reststrahlen band for both ThO 2 and UO 2 . Elevated temperature measurements were performed for ThO 2 in order to study the mechanisms by which the phonon parameters evolve with increased heat. The observed change in the TO resonant frequency is in excellent agreement with previous density functional calculations, which only consider volume expansion of the crystal lattice. This suggests that the temperature-dependent change in the TO frequency is primarily due to volume expansion. The change in the main mode pair’s broadening parameters is nearly linear within the temperature range of this study, which indicates the intrinsic anharmonic scattering (via cubic anharmonicities) as the main decay mechanism.

36 MATERIALS SCIENCE↗