Search NASA⌕ Search

SEARCH · Search NASA

Results for “torsional 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.

17 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

MSTor 2023 : A new version of the computer code for multistructural torsional anharmonicity, now with automatic torsional identification using redundant internal coordinates

The MSTor program is a computer program for calculating partition functions and thermodynamic functions of complex gas-phase molecules with multiple torsions; the multi-structural approximation with torsional anharmonicity (MS-T) can be used based on either a coupled torsional potential or an uncoupled torsional potential. The program can also carry out calculations in the multiple-structure local harmonic (MS-LH) approximation or multi-structural local quasiharmonic (MS-LQ) approximation or by the dual-level MS-T method. Furthermore, the program package includes eight utility codes that can be used as stand-alone programs to help the user to generate the input files for the code or to generate comparison results. The 2023 version of MSTor includes a new capability, namely automatic identification of torsional modes with redundant internal coordinates.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Unimolecular Reactions of E -Glycolaldehyde Oxide and Its Reactions with One and Two Water Molecules

The kinetics of Criegee intermediates are important for atmospheric modeling. However, the quantitative kinetics of Criegee intermediates are still very limited, especially for those with hydroxy groups. Here, we calculate rate constants for the unimolecular reaction of E-glycolaldehyde oxide [E-hydroxyethanal oxide, E-(CH 2 OH)CHOO], for its reactions with H 2 O and (H 2 O) 2 , and for the reaction of the E-(CH 2 OH)CHOO…H 2 O complex with H 2 O. For the highest level of electronic structure, we use W3X-L//CCSD(T)-F12a/cc-pVDZ-F12 for the unimolecular reaction and the reaction with water and W3X-L//DF-CCSD(T)-F12b/jun-cc-pVDZ for the reaction with 2 water molecules. For the dynamics, we use a dual-level strategy that combines conventional transition state theory with the highest level of electronic structure and multistructural canonical variational transition state theory with small-curvature tunneling with a validated density functional for the electronic structure. This dynamical treatment includes high-frequency anharmonicity, torsional anharmonicity, recrossing effects, and tunneling. We find that the unimolecular reaction of E-(CH 2 OH)CHOO depends on both temperature and pressure. The calculated results show that E-(CH 2 OH)CHOO…H 2 O + H 2 O is the dominant entrance channel, while previous investigations only considered Criegee intermediates + (H 2 O) 2 . In addition, we find that the atmospheric lifetime of E-(CH 2 OH)CHOO with respect to 2 water molecules is particularly short with a value of 1.71 × 10 –6 s at 0 km, which is about 2 orders of magnitude shorter than those usually assumed for Criegee intermediate reactions with water dimer. We also find that the OH group in E-(CH 2 OH)CHOO enhances its reactivity.

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

Experimental and computational studies of the kinetics of the reaction of hydrogen peroxide with the amidogen radical

In this work, the pulsed-laser photolysis/laser-induced fluorescence method is used to study the kinetics of the reaction of NH 2 with H 2 O 2 to yield a second-order rate constant of (2.42 ± 0.55) × 10 –14 cm 3 molecule –1 s –1 at 412 K in 10–22 mbar of Ar bath gas. There are no prior measurements for comparison. To check this value and enable reliable extrapolation to other temperatures, we also compute thermal rate constants for this process over the temperature range 298–3000 K via multi-structural canonical variational transition-state theory with small-curvature multidimensional tunneling (MS-CVT/SCT). The CVT/SCT rate constants are derived using a dual-level direct dynamics approach utilizing single-point CCSD(T)-F12b/cc-pVQZ-F12 energies—corrected for core-valence and scalar relativistic effects—and M06-2X/MG3S geometries, gradients, and Hessians—for all stationary and non-stationary points along the reaction path. The multistructural method with torsional anharmonicity, based on a coupled torsional potential, is then employed to calculate correction factors for the rate constants, accounting for the comprehensive effects of torsional anharmonicity on the kinetics of this reaction system. The final MS-CVT/SCT rate constants are found to be in good agreement with our measurements and can be expressed in modified Arrhenius form as 2.13 × 10 –15 (T/298 K) 4.02 exp(–513 K/T) cm 3 molecule –1 s –1 over the temperature range of 298–3000 K.

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

Pilgrim : A thermal rate constant calculator and a chemical kinetics simulator

Pilgrim is a program written in Python and designed to use direct dynamics in the calculation of thermal rate constants of chemical reactions by the variational transition state theory (VTST), based on electronic structure calculations for the potential energy surface. Pilgrim can also simulate reaction mechanisms using kinetic Monte Carlo (KMC). For reaction processes with many elementary steps, the rate constant of each of these steps can be calculated by means of conventional transition state theory (TST) or by using VTST. In the current version, Pilgrim can evaluate thermal rates using the canonical version of reaction-path VTST, which requires the calculation of the minimum energy path (MEP) associated with each elementary step or transition structure. Multi-dimensional quantum effects can be incorporated through the small-curvature tunneling (SCT) approximation. These methodologies are available both for reactions involving a single structure of the reactants and the transition state and also for reactions involving flexible molecules with multiple conformations of the reactant and/or of the transition state. For systems with many conformers, the program can evaluate each of the elementary reaction rate constants by multipath canonical VTST or multi-structural VTST. Moreover, the reactant can be unimolecular or bimolecular. Torsional anharmonicity can be incorporated through either the MSTor or the Q2DTor programs. Dual-level calculations are also available in Pilgrim: automatic high-level single-point energies can be used to correct the energy of reactants, transition states, products, and MEP points using the interpolated single-point energies (ISPE) algorithm. When the rate constants of all the chemical processes of interest are known, by means of their calculation using Pilgrim or alternatively through analytical fits to the rate constants as functions of temperature, it is possible to simulate a multistep mechanism under specified laboratory conditions using KMC. Finally, this algorithm allows performing a kinetic simulation to monitor the evolution of each chemical species with time and obtain the product yields.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Systematic exploration of the thermochemistry for a set of peroxy hydroperoxy-alkyl radicals

Here, the thermochemistry of peroxy hydroperoxy-alkyl ($\overline{O}$OQOOH) radicals has a significant influence on the reactivity of fuels and on the formation of highly oxygenated molecules (HOMs) in the atmosphere. Theoretical characterization of these radicals can be arduous due to their molecular size and complex fundamental interactions, such as hydrogen-bonding and torsional anharmonicity, and difficult to validate in the absence of any direct experimental thermochemical data. In this work, we systematically explore the thermochemistry of a set of $\overline{O}$OQOOH radicals with five increasingly affordable approaches with considerations for these interactions. In doing so, we present a novel conformer selection approach suited to predict properties at combustion temperatures. As a corollary, we also highlight the shortcomings in the standard choice of the ground conformer as the reference. The set of molecules is comprised of 149 C 2 -C 8 $\overline{O}$OQOOH radicals, selected to encompass a wide variety of branching and substitution patterns. Comparisons amongst the approaches help quantify the errors arising from various simplifying assumptions. For the largest of these radicals, and with the most affordable approach, final predictions of Gibbs energies are assigned a 2 σ uncertainty of 4 kcal mol -1 in the negative temperature coefficient region.

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

Kinetics of the benzyl + HO 2 and benzoxyl + OH barrierless association reactions: fate of the benzyl hydroperoxide adduct under combustion and atmospheric conditions

Radical–radical association reactions are challenging to address theoretically due to difficulties finding the bottleneck that variationally minimizes the reactive flux. For this purpose, the variable reaction coordinate (VRC) formulation of the variational transition state theory (VTST) represents an appropriate tool. Here, we revisited the kinetics of two radical–radical association reactions of importance in combustion modelling and poly-aromatic hydrocarbon (PAH) chemistry by performing VRC calculations: benzyl + HO 2 and benzoxyl + OH, both forming the adduct benzyl hydroperoxide. Our calculated rate constants are significantly lower than those previously reported based on VTST calculations, which results from a more efficient minimization of the reactive flux through the bottleneck achieved by the VRC formulation. Both reactions show different trends in the variation of their rate constants with temperature. We observed that if the pair of single occupied molecular orbitals (SOMOs) of the associating radicals show a similar nature, i.e. similar character, and thereby a small energy gap, a highly stabilized transition state structure is formed as the result of a very efficient SOMO–SOMO overlap, which may cancel out the free energy bottleneck for the formation of the adduct and result in large rate constants with a negative temperature dependence. This is the case of the benzoxyl and OH radical pair, whose SOMOs show O 2p nature with an energy gap of 20.2 kcal mol –1 . On the other hand, the benzyl and HO 2 radical pair shows lower rate constants with a positive temperature dependence due to the larger difference between both SOMOs (a 28.9 kcal mol –1 energy gap) as a consequence of the contribution of the multiple resonance structures of the benzyl radical. The reverse dissociation rate constants were also calculated using multi-structural torsional anharmonicity partition functions, which were not included in previous work, and the results show a much slower dissociation of benzyl hydroperoxide. Our work may help to improve kinetic models of interest in combustion and PAH formation, as well as to gain further understanding of radical–radical association reactions, which are ubiquitous in different environments.

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

TUMME: Tsinghua University Minnesota Master Equation program

We report that TUMME is a program for assembling and solving master equations for gas-phase chemical kinetics based on chemically significant eigenmodes. TUMME has interfaces to the Gaussian, Polyrate, and/or MSTor output files that allow the master equation code to obtain the microcanonical flux coefficients needed for the coefficient matrix of the master equation. The flux coefficients for reactions with barriers can be calculated by multi-structural variational transition state theory with small-curvature tunneling (MS-VTST/SCT) or by simpler approximations to this such as conventional transition state theory without tunneling (also called RRKM theory). The flux coefficients for barrierless reactions are provided by a hard-sphere model. TUMME is written in double precision with Python 3; quadruple and octuple precision are also available for some subtasks in C++. The Python code can run in serial or parallel (MP or MPI), and the C++ code can run on a single processor or on multiple processors with OpenMP.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Systematically derived thermodynamic properties for alkane oxidation

Key combustion properties, such as ignition delay time, show strong sensitivity to the thermochemistry of the main species in the standard radical oxidation pathway at low temperatures (600-1000K). Significant uncertainties persist in current estimates of thermodynamic properties, particularly for the larger species of relevance to common practical fuels. In this article, we use advanced computational schemes to evaluate thermodynamic properties for the fuel, fuel radical, peroxy, hydroperoxy-alkyl radical, and hydroperoxide species for a set of 17 fuels containing up to 9 carbon atoms, and with various degrees of branching in the alkane backbone. The procedure, termed STAR-1D, combines conformer sampling to find the minimum geometry, B2PLYP-D3/cc-pVTZ harmonic frequency evaluations, and omega B97X-D/cc-pVTZ one-dimensional torsional mapping. It includes two physically-based scaling routines: a frequency-dependent scaling to B2PLYP-D3/cc-pVTZ anharmonic frequencies and a scaling of one-dimensional omega B97X-D/ccpVTZ torsional profiles to reproduce the product of the B2PLYP-D3/cc-pVTZ frequencies in the harmonic limit. Substantive comparisons with existing experimental databases, together with careful examinations of key theoretical assumptions, are used to explore the accuracy of the predictions. These computationally intensive explorations of 195 species were facilitated by automated thermochemistry software. Further, high accuracy 0 K heats of formation from a separate study are used with the STAR-1D computations to generate NASA polynomial representations. In addition to their intrinsic value, the present results also provide a reliable database for the optimization of group additivity or machine learning schemes for scaling to larger combustion systems. Towards this end a complementary extensive conformational analysis is carried out for the medium sized species and the thermodynamic properties of the lowest energy hydrogen-bonded and non-hydrogen bonded conformers are contrasted for larger hydroperoxy-alkyl radical species.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1993

Molecular-beam spectrum of the 970/cm Fermi triad of CF3CH3

A tunable MW-sideband CO2 laser is used with an electric-resonance optothermal spectrometer to investigate the IR spectrum of CF3CH3 near 970/cm. A Fermi-coupled triad of states is observed, resulting from the interactions of 2nu6+nu11 and nu5+nu12 with the fundamental vibration, nu10, which is assumed to carry the oscillator strength in this region. The high resolution (about 3 MHz) of the spectrometer allows the observation of tunneling splittings associated with the nu6 torsional vibration. These splittings are used to identify the torsional character of the states observed. At the normal-mode level the nu10 and nu5+nu12 states are found to be nearly degenerate and interacting by an anharmonic matrix element of about 3/cm. The present investigation illustrates the utility of using resolved torsional splittings to unravel complex vibrational couplings in molecules.

Fraser, G. T.↗
DOE OSTI2026

Transient hydroperoxyalkyl intermediates (•QOOH) in isopentane oxidation. I. Conformer- and isomer-resolved infrared spectra

Transient carbon-centered hydroperoxyalkyl intermediates (•QOOH) in isopentane oxidation are characterized by infrared (IR) action spectroscopy under jet-cooled conditions with selective detection of hydroxyl (OH) radical products. Two distinct •QOOH isomers with radical sites at a primary carbon of one of the methyl groups (β-Me) or a secondary carbon (β-Et) of the ethyl group are identified by comparison with theoretically calculated IR absorption features for multiple conformers, including C–O, ethyl, and O–O torsions of the two isomers. Here, a master-equation analysis is developed and utilized to explore the conformational cooling process and the resultant conformer population distributions of the two isomers. Most of the IR features observed, including those in the strong fundamental and first overtone OH stretch regions as well as weaker combination bands involving OH stretch with torsion or OOH bend, are in good accord with computed anharmonic frequencies for the most populated conformers of the •QOOH_Et isomer. Relatively weak IR features most evident in the fundamental and first overtone OH stretch regions are ascribed to multiple conformers of the less stable •QOOH_Me isomer, along with a weak feature that is uniquely attributed to a combination band involving asymmetric CH 2 stretch and HCH bend of •QOOH_Me.

Chemical reaction jet spectroscopy↗
DOE OSTI2024

Infrared Fingerprint and Unimolecular Decay Dynamics of the Hydroperoxyalkyl Intermediate (•QOOH) in Cyclopentane Oxidation

A transient carbon-centered hydroperoxyalkyl intermediate (•QOOH) in the oxidation of cyclopentane is identified by IR action spectroscopy with time-resolved unimolecular decay to hydroxyl (OH) radical products that are detected by UV laser-induced fluorescence. Two nearly degenerate •QOOH isomers, β- and γ-QOOH, are generated by H atom abstraction of the cyclopentyl hydroperoxide precursor. Fundamental and first overtone OH stretch transitions and combination bands of •QOOH are observed and compared with anharmonic frequencies computed by second-order vibrational perturbation theory. An OH stretch transition is also observed for a conformer arising from torsion about a low-energy CCOO barrier. Definitive identification of the β-QOOH isomer relies on its significantly lower transition state (TS) barrier to OH products, which results in rapid unimolecular decay and near unity branching to OH products. A benchmarking approach is utilized to compute high-accuracy stationary point energies, most importantly TS barriers, for cyclopentane oxidation (C 5 H 9 O 2 ), building on higher level reference calculations for ethane oxidation (C 2 H 5 O 2 ). The experimental OH product appearance rates are compared with computed statistical microcanonical rates using RRKM theory, including heavy-atom tunneling, thereby validating the computed TS barrier. The results are extended to thermal unimolecular decay rate constants at temperatures and pressures relevant to cyclopentane combustion via master-equation modeling. Furthermore, the various torsional and ring puckering states of the wells and transition states are explicitly considered in these calculations.

Chemical calculations↗
DOE OSTI2022

Infrared spectroscopic signature of a hydroperoxyalkyl radical (•QOOH)

Infrared (IR) action spectroscopy is utilized to characterize a prototypical carbon-centered hydroperoxyalkyl radical (•QOOH) transiently formed in the oxidation of volatile organic compounds. The •QOOH radical formed in isobutane oxidation, 2-hydroperoxy-2-methylprop-1-yl, •CH 2 (CH 3 ) 2 COOH, is generated in the laboratory by H-atom abstraction from tert-butyl hydroperoxide (TBHP). IR spectral features of jet-cooled and stabilized •QOOH radicals are observed from 2950 to 7050 cm –1 at energies that lie below and above the transition state barrier leading to OH radical and cyclic ether products. The observed •QOOH features include overtone OH and CH stretch transitions, combination bands involving OH or CH stretch and a lower frequency mode, and fundamental OH and CH stretch transitions. Most features arise from a single vibrational transition with band contours well simulated at a rotational temperature of 10 K. In each case, the OH products resulting from unimolecular decay of vibrationally activated •QOOH are detected by UV laser-induced fluorescence. Assignments of observed •QOOH IR transitions are guided by anharmonic frequencies computed using second order vibrational perturbation theory, a 2 + 1 model that focuses on the coupling of the OH stretch with two low-frequency torsions, as well as recently predicted statistical •QOOH unimolecular decay rates that include heavy-atom tunneling. Most of the observed vibrational transitions of •QOOH are readily distinguished from those of the TBHP precursor. Furthermore, the distinctive IR transitions of •QOOH, including the strong fundamental OH stretch, provide a general means for detection of •QOOH under controlled laboratory and real-world conditions.

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

Experimental and theoretical studies of the doubly substituted methyl-ethyl Criegee intermediate: Infrared action spectroscopy and unimolecular decay to OH radical products

In this work the infrared (IR) action spectrum of the doubly substituted methyl-ethyl Criegee intermediate (MECI) is observed in the CH stretch overtone region with detection of OH products. The MECI exhibits four conformers, all of which undergo unimolecular decay via a 1,4 H-atom transfer mechanism, followed by the rapid release of OH products. Conformers with different orientations of the carbonyl oxide group with respect to the methyl and ethyl substituents (i.e., anti and syn ) decay via distinct transition state barriers (16.1 kcal mol –1 and 15.4 kcal mol –1 , respectively). The observed IR action spectrum is in good agreement with the predicted anharmonic IR absorption spectrum, but exhibits significant congestion, which is attributed to couplings between spectroscopic bright states and nearby dark states. Energy-dependent OH appearance rates are measured upon IR excitation of the strongest features in the IR action spectrum and are found to be on the order of 10 6 –10 7 s –1 . The experimental rates are in good agreement with computed Rice–Ramsperger–Kassel–Marcus rates for the unimolecular decay of MECI at these energies, which incorporate quantum mechanical tunneling and sophisticated hindered rotor treatments, as well as high-level theoretical calculations of the TS barrier heights, rovibrational properties, and torsional barriers associated with the MECI conformers. Master equation modeling is used to predict thermal rates for the unimolecular decay of anti - and syn -MECI of 473 s –1 and 660 s –1 , respectively. Comparison with other previously studied Criegee intermediate systems provides insights into substituent effects on unimolecular decay under both energy-dependent and thermal conditions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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 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↗