Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “anharmonicity”

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
NASA NTRS1995

An Accurate ab initio Quartic Force Field and Vibrational Frequencies for CH4 and Isotopomers

A very accurate ab initio quartic force field for CH4 and its isotopomers is presented. The quartic force field was determined with the singles and doubles coupled-cluster procedure that includes a quasiperturbative estimate of the effects of connected triple excitations, CCSD(T), using the correlation consistent polarized valence triple zeta, cc-pVTZ, basis set. Improved quadratic force constants were evaluated with the correlation consistent polarized valence quadruple zeta, cc-pVQZ, basis set. Fundamental vibrational frequencies are determined using second-order perturbation theory anharmonic analyses. All fundamentals of CH4 and isotopomers for which accurate experimental values exist and for which there is not a large Fermi resonance, are predicted to within +/- 6 cm(exp -1). It is thus concluded that our predictions for the harmonic frequencies and the anharmonic constants are the most accurate estimates available. It is also shown that using cubic and quartic force constants determined with the correlation consistent polarized double zeta, cc-pVDZ, basis set in conjunction with the cc-pVQZ quadratic force constants and equilibrium geometry leads to accurate predictions for the fundamental vibrational frequencies of methane, suggesting that this approach may be a viable alternative for larger molecules. Using CCSD(T), core correlation is found to reduce the CH4 r(e), by 0.0015 A. Our best estimate for r, is 1.0862 +/- 0.0005 A.

Lee, Timothy J.↗
NASA NTRS2002

Theoretical Study of Infrared and Raman Spectra of Hydrated Magnesium Sulfate Salts

Harmonic and anharmonic vibrational frequencies, as well as infrared and Raman intensities, are calculated for MgSO4.nH20 (n=1-3). Electronic structure theory at the second order Moller-Plesset perturbation theory (MP2) level with a triple-zeta + polarization (TZP) basis set is used to determine the geometry, properties, and vibrational spectra of pure and hydrated MgSO4 salts. The direct vibrational self-consistent field (VSCF) method and its correlation corrected (CC-VSCF) extension are used to determine anharmonic corrections to vibrational frequencies and intensities for the pure MgSO4 and its complex with one water molecule. Very significant differences are found between vibrational of water molecules in complexes with MgSO4 and pure water. Some of the O-H stretching frequencies are shifted to the red very significantly (by up to 1500-2000/cm) upon complexation with magnesium sulfate. They should be observed between 1700 and 3000/cm in a region very different from the corresponding O-H stretch frequency region of pure water (3700-3800/cm). In addition, the SO2 stretching vibrations are found at lower frequency regions than the water vibrations. They can serve as unique identifiers for the presence of sulfate salts. The predicted infrared and Raman spectra should be of valuable help in the design of future missions and analysis of observed data from the ice surface of Jupiter's moon Europa that possibly contains hydrated MgSO4 salts.

Chaban, Galina M.↗
NASA NTRS2005

Hydrogen Bonding and Vibrational Spectroscopy: A Theoretical Study

Effects of hydrogen bonding on vibrational spectra are studied for several hydrogen-bonded complexes, in which hydrogen bonding ranges from weak (<5 kcal/mol) to very strong (>25 kcal/mol). The systems studied include complexes of inorganic acids and salts with water and ammonia, as well as complexes of several organic molecules (nitriles and amino acids) with water. Since anharmonic effects are very strong in hydrogen-bonded systems, anharmonic vibrational frequencies and infrared intensities are computed using the correlation-corrected vibrational self-consistent field (CC-VSCF) method with ab initio potential surfaces at the MP2 and CCSD(T) levels. The most common spectral effects induced by hydrogen bonding are red shifts of stretching vibrational frequencies ranging from approx.200/cm to over 2000/cm and significant increases of infrared intensities for those bonds that participate in hydrogen bonding. However, some systems (e.g. nitrile-water complexes) exhibit shifts in the opposite direction (to the blue) upon formation of hydrogen bonds.

Chaban, Galina M.↗
NASA NTRS2019

Quantifying the Influence of Redox Conditions on the Seismic Properties of Olivine

A series of eight olivine specimens were fabricated by hot-pressing at 1200 degrees Celsius and 300 megapascals (MPa). Each hot-pressed specimen was then wrapped in Pt, Ni or NiFe foil to vary oxygen fugacity (fO2), and interrogated via forced torsional oscillation. Mechanical testing was conducted at 10 oscillation periods between 1 and 1000 seconds, at a confining pressure of 200 MPa, during a slow staged-cooling from a maximum temperature of 1200 degrees Celsius to room temperature. After mechanical testing, each specimen was axially sectioned and EBSD (Electron BackScatter Diffraction) was used for the determination of the representative grain size, and grain size distribution of each sample. In addition, each longitudinal section was mapped via FTIR (Fourier Transform InfraRed spectroscopy) to determine the spatial distribution and concentrations of chemically bound and molecular water. Amongst these specimens, chemically bound ‘water’ contents were observed to vary between 0 and 1150 atom ppm (parts per million) H/Si, and molecular water concentrations varied between 0 and 245 atom ppm H/Si. Our forced-oscillation results demonstrate that the measured magnitude of anelastic relaxation within the experimental ‘window’ of oscillation periods is unrelated to the water content. Rather, a relationship was observed between the magnitude of anelastic relaxation and the prevailing redox conditions, which is influenced by the choice of metal sleeving used during the mechanical test. Further, regardless of water content or metal sleeving, each specimen exhibits coupled variations in shear modulus and dissipation within the observational window, indicative of ‘high-temperature background’ behavior, that can be described by a Burgers-type model. During initial fitting of the Burgers models, the unrelaxed shear modulus at a reference temperature of 900 degrees Celsius (G (sub UR)) and the temperature derivative of the unrelaxed shear modulus (dG (sub U)/dT), were treated as adjustable parameters. For all Fe-bearing olivine samples (but not a hydrous and oxidized Fe-free sample) we observe deficits of G (sub UR), and increased values of dG (sub U)/dT, relative to the expected elastic (anharmonic) behavior of Fo90 olivine. This behavior is indicative of anelastic relaxation occurring at shorter periods than observable within the ‘window’ of oscillation periods used in the mechanical test. Moving towards a comprehensive seismologically applicable Burgers model, which includes this newly observed effect of redox conditions on anelastic relaxation, we will present our progress in reconciling truly anharmonic and elastic behavior of Fo (sub 90) olivine with our observed forced-oscillation data.

Cline, Christopher, II↗
NASA NTRS2019

The Influence of Redox Conditions on the Seismic Properties of Polycrystalline Olivine

Eight olivine specimens were fabricated by use of a solgel method and hot-pressing at 1200 degrees Centigrade and 300 megapascals (MPa) inside of welded Pt capsules. Each hot-pressed specimen was then recovered, precision ground, and wrapped in Pt, Ni or NiFe foil to vary oxygen fugacity (fO2) during the subsequent forced torsional oscillation measurements. Mechanical testing was conducted at 10 oscillation periods between 1 and 1000 seconds, at a confining pressure of 200 MPa, during a slow staged-cooling from a maximum temperature of 1200 degrees Centigrade down to room temperature. After mechanical testing, each specimen was axially sectioned and EBSD (Electron BackScatter Diffraction) was used for the determination of the representative grain size, and grain size distribution of each sample. In addition, each longitudinal section was mapped via FTIR (Fourier Transform Infrared Spectroscopy) to determine the spatial distribution and concentrations of chemically bound and molecular water. Amongst these eight specimens, chemically bound ‘water’ contents were observed to vary between 0 and 1150 atom parts per million (ppm) H/Si, and molecular water concentrations varied between 0 and 245 atom ppm H/Si. Our forced-oscillation results demonstrate that the measured magnitude of anelastic relaxation within the experimental ‘window’ of oscillation periods is unrelated to the water content. Rather, a relationship was observed between the magnitude of anelastic relaxation and the prevailing redox conditions, which is influenced by the choice of metal sleeving used during the mechanical test. Further, regardless of water content or metal sleeving, each specimen exhibits coupled variations in shear modulus and dissipation within the observational window, indicative of ‘high-temperature background’ behavior, that can be described by a Burgers-type model. During initial fitting of the Burgers models, the unrelaxed shear modulus at a reference temperature of 900°C (elastic unloading/reloading shear modulus G (sub UR)) and the temperature derivative of the unrelaxed shear modulus (dGU/dT), were treated as adjustable parameters. For all Fe-bearing olivine samples we observe deficits of G (sub UR) and increased values of dGU/dT, relative to the expected elastic (anharmonic) behavior for Fo (Forsterite content percentage) (sub 90) olivine. This behavior is indicative of anelastic relaxation occurring at shorter periods than observable within the ‘window’ of oscillation periods used in the mechanical test. Moving towards a comprehensive and seismologically applicable Burgers model, which includes this newly observed effect of redox conditions on anelastic relaxation, we will present our progress on reconciling truly anharmonic and elastic behavior of Fo (sub 90) olivine with our observed forced-oscillation data.

Cline, C. J., II↗
NASA NTRSDate not supplied

Metal Hydroxide and Oxyhydroxide Gas Database

In a high temperature combustion situation such as a gas turbine or rocket engine, combustion products include water vapor. This water vapor reacts with oxides in the hot-gas-path to form volatile hydroxides and oxyhydroxides. The oxides may be structural oxide components or protective oxides on alloys or non-oxide ceramics. In either case, this enhanced volatilization leads to degradation of components. Reliable thermochemical data for many metal hydroxides and oxyhydroxides has not been available. We review a database put together at NASA over many years based on experimental and ab initio computations for thermochemical data of Al, Si, Ti, Cr, Mn, Y, Zr, La, Ta, Gd, and Yb hydroxides and oxyhydroxides. Experimental methods for measuring thermodynamic quantities of metal hydroxides and oxyhydroxides are limited since the required high temperatures and high pressures preclude many conventional methods. Nonetheless, transpiration has proven a valuable tool in our laboratories for the Si-O-H, Cr-O-H, and Ti-O-H systems. Data are analyzed with second and third law methods. These data can be leveraged with quantum chemistry to provide reliable free energy functions for the third law calculation and comparison to the derived enthalpies of formation. Both wave function and DFT methods can be utilized for deriving the thermodynamic functions for the various metal hydroxides and oxyhydroxides. We have found that, in general, reasonable molecular geometries and spectroscopic data can be derived with a relatively low level of theory, i.e. DFT with the B3LYP functional. Accurate calculations of the enthalpy of formation require higher levels of theory, such as the Coupled Cluster, Singles, Doubles, and perturbative Triples CCSD(T) approach with complete basis sets (CBS). Also, corrections for hindered rotations and anharmonicity are necessary. For some symmetries, the anharmonic correction has proven to be a challenge and we have found that this is only treated properly by certain codes. The database lists the fHo(298), So(298), and Cp in the form Cp = A + BT + CT2 + D/T + E/T2, where A, B, C, D, and E are the fitting constants. Thus, it can be easily used in any of the common computational thermochemistry codes. To date, the database contains 51 gaseous species of the formula MOx(OH)y.

Metal Hydroxides↗
DOE OSTI2022

Dynamically preferred state with strong electronic fluctuations from electrochemical synthesis of sodium manganate

Electrochemical (de)intercalation is a delicate method to precisely control the alkaline ion composition in alkaline transition metal oxides. Because of complicated interactions, metal charge ordering patterns can form spontaneously at special fractional alkaline compositions and orderings. In this study, we show that this elegant electrochemical process can create dynamically preferred structures in an anharmonic energy landscape that conventional syntheses and computations can rarely visit. Specifically, electrochemically prepared Na 1/2 MnO 2 ordering exhibits abnormal structure distortions, charge orderings, and dynamical activities. Strong magnetic fluctuations and lattice dynamics are observed in an unusually wide temperature range in Na 1/2 MnO 2 , which distinguishes it from all other Na x MnO 2 at higher or lower Na compositions. The results emphasize the unique opportunity of using electrochemical processes to design and create novel quantum states with strongly coupled and mutually enhanced electronic and lattice fluctuations, likely through a special dynamic charge flux functional, as suggested by our computational investigations.

25 ENERGY STORAGE↗
DOE OSTI2023

Isotope Effects in the Zundel–Eigen Isomerization of H + (H 2 O) 6

The isomerization pathway between the energetically low-lying Zundel and Eigen isomers of the protonated water hexamer was investigated using high-level ab initio calculations including a treatment of zero-point corrections. On the basis of these calculations, the Zundel–Eigen isomerization was found to proceed through a stable intermediate isomer, which consists of a four-membered ring with two single acceptor water molecules. Further, the inclusion of vibrational zero-point energy is shown to be important for accurately establishing the relative energies of the three relevant isomers involved in the Zundel–Eigen isomerization. Diffusion Monte Carlo calculations including anharmonic vibrational effects show that all three isomers of H + (H 2 O) 6 and D + (D 2 O) 6 have well-defined structures. The energetic ordering of the three isomers changes upon deuteration. The implications of these results for the vibrational spectra of H + (H 2 O) 6 and D + (D 2 O) 6 are also discussed.

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

Multimode two-dimensional vibronic spectroscopy. I. Orientational response and polarization-selectivity

Two-dimensional Electronic–Vibrational (2D EV) spectroscopy and two-dimensional Vibrational–Electronic (2D VE) spectroscopy are among the newest additions to the coherent multidimensional spectroscopy toolbox, and they are directly sensitive to vibronic couplings. In this first of two papers, the complete orientational response functions are developed for a model system consisting of two coupled anharmonic oscillators and two electronic states in order to simulate polarization-selective 2D EV and 2D VE spectra with arbitrary combinations of linearly polarized electric fields. Here, we propose analytical methods to isolate desired signals within complicated spectra and to extract the relative orientation between vibrational and vibronic dipole moments of the model system using combinations of polarization-selective 2D EV and 2D VE spectral features. Time-dependent peak amplitudes of coherence peaks are also discussed as means for isolating desired signals within the time-domain. This paper serves as a field guide for using polarization-selective 2D EV and 2D VE spectroscopies to map coupled vibronic coordinates on the molecular frame.

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

PBE-GGA predicts the B8↔B2 phase boundary of FeO at Earth’s core conditions

FeO is a crucial component of the Earth’s core, and its thermodynamic properties are essential to developing more accurate core models. It is also a notorious correlated insulator in the NaCl-type (B1) phase at ambient conditions. It undergoes two polymorphic transitions at 300 K before it becomes metallic in the NiAs-type (B8) structure at ~100 GPa. Although its phase diagram is not fully mapped, it is well established that the B8 phase transforms to the CsCl-type (B2) phase at core pressures and temperatures. Here, we report a successful ab initio calculation of the B8↔B2 phase boundary in FeO at Earth’s core pressures. We show that fully anharmonic free energies computed with the Perdew–Burke–Ernzerhof-generalized gradient approximation coupled with thermal electronic excitations reproduce the experimental phase boundary within uncertainties at P > 255 GPa, including the largely negative Clapeyron slope of –52 MPa/K. This study validates the applicability of a standard density functional theory functional to FeO under Earth’s core conditions and demonstrates the theoretical framework that enables complex predictive studies of this region.

58 GEOSCIENCES↗
DOE OSTI2022

Stability and distortion of fcc LaH 10 with path-integral molecular dynamics

The synthesis of the high-temperature superconductor LaH 10 requires pressures in excess of 100 GPa, wherein it adopts a face-centered cubic structure. Upon decompression, this structure undergoes a distortion, which still supports superconductivity, but with a lower critical temperature. Previous calculations have shown that quantum and anharmonic effects are necessary to stabilize the cubic structure, but have not resolved the low pressure distortion. Using large scale path-integral molecular dynamics enabled by a machine learned potential, we show that a rhombohedral distortion appears at sufficiently low pressures. Here, we also highlight the importance of quantum zero-point motion in stabilizing the cubic structure.

08 HYDROGEN↗
DOE OSTI2025

Thermal disorder and phonon softening in the ferroelectric phase transition of lead titanate

We report a molecular dynamics study of ab initio quality of the ferroelectric phase transition in crystalline PbTi⁢O3. We model anharmonicity accurately in terms of potential energy and polarization surfaces trained on density functional theory data with modern machine learning techniques. Our simulations demonstrate that the transition has a strong order-disorder character, in agreement with diffraction experiments, and provide fresh insight into the approach to equilibrium across the phase transition. We find that the emergence and disappearance of the macroscopic polarization is driven by dipolar switching at the nanometer scale. We also computed the infrared optical absorption spectra in both the ferroelectric and the paraelectric phases, finding good agreement with the experimental Raman frequencies. Often, the almost ideal displacive character of the soft mode detected by Raman scattering in the paraelectric phase has been contrasted with the order-disorder character of the transition suggested by diffraction experiments. We settle this issue by showing that the soft mode coexists with a strong Debye relaxation associated with thermal disordering of the dipoles. The Debye relaxation feature is centered at zero frequency and appears near the transition temperature in both the ferroelectric and the paraelectric phases.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Phonon second harmonic generation in NaBr studied by inelastic neutron scattering and computer simulation

The phenomenon of second harmonic generation (SHG) was found for phonons in anharmonic NaBr by inelastic neutron scattering. The temperature dependence of this phonon SHG was measured from 300 K to 650 K. At 300 K the second harmonic (SH) is seen as a high-energy branch around 33 meV, nearly independent of $\overrightarrow{Q}$. The temperature effective potential (TDEP) method and classical molecular dynamics (MD) simulation with machine learning interatomic potential were able to reproduce the SH, and showed that SHG occurs with the flat transverse optical (TO) phonon branch. A classical model of a nonlinear medium explains the intensity and lifetime of the SH, compared to those of the TO modes. Also successful was a quantum model based on the Heisenberg-Langevin equation for interacting phonons coupled to a thermal bath, which also predicts a spectral distribution of the SH. In conclusion, the measured temperature dependence of the intensity of the second harmonic showed that it follows the Planck distribution of a one-phonon quasiparticle, and not two TO phonons.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Ultralow Thermal Conductivity, Multiband Electronic Structure and High Thermoelectric Figure of Merit in TlCuSe

The entanglement of lattice thermal conductivity, electrical conductivity, and Seebeck coefficient complicates the process of optimizing thermoelectric performance in most thermoelectric materials. Semiconductors with ultralow lattice thermal conductivities and high power factors at the same time are scarce but fundamentally interesting and practically important for energy conversion. In this work, an intrinsic p-type semiconductor TlCuSe that has an intrinsically ultralow thermal conductivity (0.25 W m -1 K -1 ), a high power factor (11.6 μ W cm -1 K -2 ), and a high figure of merit, ZT (1.9) at 643 K is described. The weak chemical bonds, originating from the filled antibonding orbitals p-d* within the edge-sharing CuSe 4 tetrahedra and long Tl-Se bonds in the PbClF-type structure, in conjunction with the large atomic mass of Tl lead to an ultralow sound velocity. Strong anharmonicity, coming from Tl + lone-pair electrons, boosts phonon-phonon scattering rates and further suppresses lattice thermal conductivity. The multiband character of the valence band structure contributing to power factor enhancement benefits from the lone-pair electrons of Tl + as well, which modify the orbital character of the valence bands, and pushes the valence band maximum off the Gamma-point, increasing the band degeneracy. The results provide new insight on the rational design of thermoelectric materials.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Sensitive Thermochromic Behavior of InSeI, a Highly Anisotropic and Tubular 1D van der Waals Crystal

Thermochromism, the change in color of a material with temperature, is the fundamental basis of optical thermometry. A longstanding challenge in realizing sensitive optical thermometers for widespread use is identifying materials with pronounced thermometric optical performance in the visible range. Herein, it is demonstrated that single crystals of indium selenium iodide (InSeI), a 1D van der Waals (vdW) solid consisting of weakly bound helical chains, exhibit considerable visible range thermochromism. A strong temperature-dependent optical band edge absorption shift ranging from 450 to 530 nm (2.8 to 2.3 eV) over a 380 K temperature range with an experimental (dE g /dT) max value extracted to be 1.26 × 10 −3 eV K −1 is shown. This value lies appreciably above most dense conventional semiconductors in the visible range and is comparable to soft lattice solids. The authors further seek to understand the origin of this unusually sensitive thermochromic behavior and find that it arises from strong electron–phonon interactions and anharmonic phonons that significantly broaden band edges and lower the E g with increasing temperature. The identification of structural signatures resulting in sensitive thermochromism in 1D vdW crystals opens avenues in discovering low-dimensional solids with strong temperature-dependent optical responses across broad spectral windows, dimensionalities, and size regimes.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Thickness-Dependent Thermal Conductivity and Phonon Mean Free Path Distribution in Single-Crystalline Barium Titanate

Nanosized perovskite ferroelectrics are widely employed in several electromechanical, photonics, and thermoelectric applications. Scaling of ferroelectric materials entails a severe reduction in the lattice (phonon) thermal conductivity, particularly at sub-100 nm length scales. Such thermal conductivity reduction can be accurately predicted using the information of phonon mean free path (MFP) distribution. The current understanding of phonon MFP distribution in perovskite ferroelectrics is still inconclusive despite the critical thermal management implications. Here, high-quality single-crystalline barium titanate (BTO) thin films, a representative perovskite ferroelectric material, are grown at several thicknesses. Using experimental thermal conductivity measurements and first-principles based modeling (including four-phonon scattering), the phonon MFP distribution is determined in BTO. The simulation results agree with the measured thickness-dependent thermal conductivity. The results show that the phonons with sub-100 nm MFP dominate the thermal transport in BTO, and phonons with MFP exceeding 10 nm contribute ≈35% to the total thermal conductivity, in significant contrast to previously published experimental results. The experimentally validated phonon MFP distribution is consistent with the theoretical predictions of other complex crystals with strong anharmonicity. This work paves the way for thermal management in nanostructured and ferroelectric-domain-engineered systems for oxide perovskite-based functional materials.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Suppressed Lone Pair Electrons Explain Unconventional Rise of Lattice Thermal Conductivity in Defective Crystalline Solids

Manipulating thermal properties of materials can be interpreted as the control of how vibrations of atoms (known as phonons) scatter in a crystal lattice. Compared to a perfect crystal, crystalline solids with defects are expected to have shorter phonon mean free paths caused by point defect scattering, leading to lower lattice thermal conductivities than those without defects. While this is true in many cases, alloying can increase the phonon mean free path in the Cd-doped AgSnSbSe 3 system to increase the lattice thermal conductivity from 0.65 to 1.05 W m –1 K –1 by replacing 18% of the Sb sites with Cd. It is found that the presence of lone pair electrons leads to the off-centering of cations from the centrosymmetric position of a cubic lattice. X-ray pair distribution function analysis reveals that this structural distortion is relieved when the electronic configuration of the dopant element cannot produce lone pair electrons. Furthermore, a decrease in the Grüneisen parameter with doping is experimentally confirmed, establishing a relationship between the stereochemical activity of lone pair electrons and the lattice anharmonicity. The observed “harmonic” behavior with doping suggests that lone pair electrons must be preserved to effectively suppress phonon transport in these systems.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Electrochemical Control of the Ultrafast Lattice Response of a Layered Semimetal

The unique layer-stacking in two-dimensional (2D) van der Waals materials facilitates the formation of nearly degenerate phases of matter and opens novel routes for the design of low-power, reconfigurable functional materials. Electrochemical ion intercalation between stacked layers offers a promising approach to stabilize bulk metastable phases and to explore the effects of extreme carrier doping and strain. However, in situ characterization methods to study the structural evolution and dynamical functional properties of these intercalated materials remains limited. Here a novel experimental platform is presented capable of simultaneously performing electrochemical lithium-ion intercalation and multimodal ultrafast characterization of the lattice using both electron diffraction and nonlinear optical techniques. Using the layered semimetal WTe 2 as a model system, the interlayer shear phonon mode that modulates stacking between 2Dlayers is probed, showing that small amounts of lithiation enhance the amplitude and lifetime of the phonon, contrary to expectations. This results from the dynamically fluctuating and anharmonic structure between nearly degenerate phases at room temperature, which can be stabilized by electronic carriers accompanying the inserted lithium ions. At high lithiation, the T d ’ structure emerges and quenches the phonon response. This work defines new approaches for using electrochemistry to engineer the dynamic structure of 2D materials.

36 MATERIALS SCIENCE↗