Search NASA⌕ Search

SEARCH · Search NASA

Results for “Kinetics Parameters”

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 1,153 records · Page 64

Tuning the many-body interactions in a helical Luttinger liquid

In one-dimensional (1D) systems, electronic interactions lead to a breakdown of Fermi liquid theory and the formation of a Tomonaga-Luttinger Liquid (TLL). The strength of its many-body correlations can be quantified by a single dimensionless parameter, the Luttinger parameter K, characterising the competition between the electrons’ kinetic and electrostatic energies. Recently, signatures of a TLL have been reported for the topological edge states of quantum spin Hall (QSH) insulators, strictly 1D electronic structures with linear (Dirac) dispersion and spin-momentum locking. Here we show that the many-body interactions in such helical Luttinger Liquid can be effectively controlled by the edge state’s dielectric environment. This is reflected in a tunability of the Luttinger parameter K, distinct on different edges of the crystal, and extracted to high accuracy from the statistics of tunnelling spectra at tens of tunnelling points. The interplay of topology and many-body correlations in 1D helical systems has been suggested as a potential avenue towards realising non-Abelian parafermions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Probabilistic effects of porosity and chemical kinetics on the shock initiation of an octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine (HMX) based explosive

In this work, we investigate the effects of porosity and chemical reactivity on the shock initiation of an HMX (octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine)-based polymer bonded explosive due to thin flyer impact using two-dimensional mesoscale simulations. Microstructural samples are computationally generated, and the calculated shock response is compared to experiments. The sample dimensions are 0.1×1.05 mm 2 , and each sample accounts for grain scale heterogeneities, including randomized grain morphology and pore spacing. Simulations are performed using the Lawrence Livermore National Laboratory (LLNL) arbitrary Lagrangian–Eulerian multi-physics code ALE3D. Chemistry and equation of states are handled with LLNL's thermochemical code Cheetah. We employ a modified Arrhenius-based chemical kinetics model that scales the temperature exponent fitting parameter with flyer shock conditions. This new model is calibrated to the mean response of experimental shock initiation thresholds for three separate Kapton flyer thicknesses. A simplified shock initiation criterion is proposed, which classifies the macroscale initiation response based on cross-sectional temperature profiles of the shocked sample. This criterion is used to further investigate the probabilistic response of the sample to variations in the mean porosity and chemical kinetics rate, and we find an inverse correlation between the shock pressure and the confidence interval of the experimental threshold velocity. Furthermore, the methodology presented in this work for quantifying the probabilistic nature of shock initiation is material agnostic. These mesoscale studies may be used to inform continuum scale models and expedite the development of future, tailored energetic formulations.

36 MATERIALS SCIENCE↗

Excitation of sodium D-line radiation in collisions of sodium atoms with internally excited H2, D2 and N2.

Excitation of D-line radiation in collisions of Na atoms with vibrationally excited N2, H2, and D2 was studied in two modulated crossed-beam experiments. In both experiments, the molecular excitation was provided by heating the molecular beam to temperatures which were assumed to give populations corresponding to the Boltzmann distribution. In the first experiment, a total rate coefficient was measured as a function of molecular beam temperature. The second experiment achieved partial separation of internal vs kinetic energy transfer effects by using a velocity-selected molecular beam. The data from both experiments were used to determine parameters for the dependence of the transfer cross section upon the kinetic energy for a given change in the vibrational quantum number. Results indicate that most of the Na excitation energy comes from internal vibrational energy, with the remainder coming from kinetic energy.

Krause, H. F.↗

Event-by-Event Evaluation of the Prompt Fission Neutron Spectrum From 239 Pu(n,f)

We have developed an improved evaluation method for the spectrum of neutrons emitted in fission of 239 Pu induced by incident neutrons with energies up to 20 MeV. The <v> covariance data, including incident energy correlations introduced by the <v> evaluation method, were used to fix the input parameters in our event-by-event model of fission, FREYA, by applying formal statistical methods. Formal estimates of uncertainties in the evaluation were developed by randomly sampling model inputs and calculating likelihood functions based on agreement with the evaluated <v>. Our approach is able to employ a greater variety of fission measurements than the relatively coarse spectral data alone. It also allows the study of numerous fission observables for more accurate model validation. The combination of an event-by-event Monte Carlo fission model with a statistical-likelihood analysis is thus a powerful tool for evaluation of fission-neutron data. Our empirical model FREYA follows the complete fission event from birth of the excited fragments through their decay via neutron emission until the fragment excitation energy is below the neutron separation energy when neutron emission can no longer occur. The most recent version of FREYA incorporates pre-equilibrium neutron emission, the emission of the first neutron before equilibrium is reached in the compound nucleus, and multi-chance fission, neutron evaporation prior to fission when the incident neutron energy is above the neutron separation energy. Energy, momentum, charge and mass number are conserved throughout the fission process. The best available values of fragment masses and total kinetic energies are used as inputs to FREYA. We fit three parameters that are not well under control from previous measurements: the shift in the total fragment kinetic energy; the energy scale of the asymptotic level density parameter, controlling the fragment 'temperature' for neutron evaporation; and the relative excitation of the light and heavy fragments, governing the number and energy of neutrons emitted from each fragment. The latter two parameters are assumed to be independent of the incident neutron energy while the first varies with incident energy. We describe our method and the subsequent spectral evaluation and present the results of several standard validation calculations that test our new evaluation. These benchmarks include critical assemblies, sensitive to criticality in fast systems; pulsed sphere measurements testing the spectra at incident neutron energies of 14 MeV; and other tests.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Carbon Dioxide Insertion into Rhenium Hydrides as a Probe for the Impact of Solvent on Linear Free Energy Relationships between Thermodynamic and Kinetic Hydricity

The kinetics of CO 2 insertion into electronically different Re( R bpy)(CO) 3 H ( R bpy = 4,4'-R-2,2'- bipyridine; R = OMe, t Bu, Me, H, Br, COOMe, CF 3 ) complexes to form Re( R bpy)(CO) 3 {OC(H)O} compounds were determined in acetone, dimethylacetamide (DMAc), dimethylformamide (DMF), dimethylsulfoxide (DMSO), and 3-methoxypropionitrile (3-MPN) and compared with previous data in acetonitrile (MeCN). The rates of CO 2 insertion for any one complex of the type Re(Rbpy)(CO) 3 H in different solvents correlates with the Dimroth-Reichardt (E T (30)) solvent parameter. Hammett plots in each solvent indicate that insertion reactions are faster for bpy ligands with electron-donating groups. There is, however, no correlation between the slope of the Hammett plot in different solvents and any common solvent parameter. Similarly, the enthalpies and entropies of activation and kinetic isotope effects associated with CO 2 insertion into Re(bpy)(CO) 3 H in different solvents do not correlate with any common solvent parameters. Theoretical calculations were used to determine the relative thermodynamic hydricities of Re( R bpy)(CO) 3 H type complexes in MeCN, acetone, DMF, and DMSO and in each solvent complexes with more electron-donating substituents on the bpy ligand are stronger hydride donors. Linear Free Energy Relationships (LFERs) between calculated thermodynamic and experimental kinetic hydricity, as measured through CO 2 insertion reactions, were observed in MeCN, acetone, DMF, and DMSO. Furthermore, the slopes of the LFERs correlate with the dielectric constant of the solvent. Overall, this work provides fundamental information about the thermodynamics and kinetics of hydride transfer reactions in different solvents, which is valuable for catalyst design.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Embedded Error Bayesian Calibration of Thermal Decomposition of Organic Materials

Organic materials are an attractive choice for structural components due to their light weight and versatility. However, because they decompose at low temperatures relative to tradiational materials they pose a safety risk due to fire and loss of structural integrity. To quantify this risk, analysts use chemical kinetics models to describe the material pyrolysis and oxidation using thermogravimetric analysis. This process requires the calibration of many model parameters to closely match experimental data. Previous efforts in this field have largely been limited to finding a single best-fit set of parameters even though the experimental data may be very noisy. Furthermore the chemical kinetics models are often simplified representations of the true de- composition process. The simplification induces model-form errors that the fitting process cannot capture. In this work we propose a methodology for calibrating decomposition models to thermogravimetric analysis data that accounts for uncertainty in the model-form and experimental data simultaneously. The methodology is applied to the decomposition of a carbon fiber epoxy composite with a three-stage reaction network and Arrhenius kinetics. The results show a good overlap between the model predictions and thermogravimetric analysis data. Uncertainty bounds capture devia- tions of the model from the data. The calibrated parameter distributions are also presented. In conclusion, the distributions may be used in forward propagation of uncertainty in models that leverage this material.

36 MATERIALS SCIENCE↗

Study of nitrogen flowing afterglow with mercury vapor injection

The reaction kinetics in nitrogen flowing afterglow with mercury vapor addition was studied by optical emission spectroscopy. The DC flowing post-discharge in pure nitrogen was created in a quartz tube at the total gas pressure of 1000 Pa and discharge power of 130 W. The mercury vapors were added into the afterglow at the distance of 30 cm behind the active discharge. The optical emission spectra were measured along the flow tube. Three nitrogen spectral systems – the first positive, the second positive, and the first negative, and after the mercury vapor addition also the mercury resonance line at 254 nm in the spectrum of the second order were identified. The measurement of the spatial dependence of mercury line intensity showed very slow decay of its intensity and the decay rate did not depend on the mercury concentration. In order to explain this behavior, a kinetic model for the reaction in afterglow was developed. This model showed that the state Hg(6 {sup 3}P{sub 1}), which is the upper state of mercury UV resonance line at 254 nm, is produced by the excitation transfer from nitrogen N{sub 2}(A{sup 3}Σ{sup +}{sub u}) metastables to mercury atoms. However, the N{sub 2}(A{sup 3}Σ{sup +}{sub u}) metastables are also produced by the reactions following the N atom recombination, and this limits the decay of N{sub 2}(A{sup 3}Σ{sup +}{sub u}) metastable concentration and results in very slow decay of mercury resonance line intensity. It was found that N atoms are the most important particles in this late nitrogen afterglow, their volume recombination starts a chain of reactions which produce excited states of molecular nitrogen. In order to explain the decrease of N atom concentration, it was also necessary to include the surface recombination of N atoms to the model. The surface recombination was considered as a first order reaction and wall recombination probability γ = (1.35 ± 0.04) × 10{sup −6} was determined from the experimental data. Also sensitivity analysis was applied for the analysis of kinetic model in order to reveal the main control parameters in the model.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Effects of alteration product precipitation on glass dissolution

Understanding the mechanisms that control the durability of nuclear waste glass is paramount if reliable models are to be constructed so that the glass dissolution rate in a given geological repository can be calculated. Presently, it is agreed that (boro)silicate glasses dissolve in water at a rate dependent on the solution concentration of orthosilicic acid (H 4 SiO 4 ) with higher [H 4 SiO 4 ] leading to lower dissolution rates. Once the reaction has slowed as a result of the buildup of H 4 SiO 4 , another increase in the rate has been observed that corresponds to the precipitation of certain silica-bearing alteration products. However, it has also been observed that the concentration of silica-bearing solution species does not significantly decrease, indicating saturation, while other glass tracer elements concentrations continue to increase, indicating that the glass is still dissolving. In this study, we have used the Geochemist’s Workbench code to investigate the relationship between glass dissolution rates and the precipitation rate of a representative zeolitic silica-bearing alteration product, analcime [Na(AlSi 2 O 6 )∙H 2 O]. To simplify the calculations, we suppressed all alteration products except analcime, gibbsite (Al(OH) 3 ), and amorphous silica. The pseudo-equilibrium-constant matrix for amorphous silica was substituted for the glass pseudo-equilibrium-constant matrix because it has been shown that silicate glasses act as a silica-only solid with respect to kinetic considerations. In this article, we present the results of our calculations of the glass dissolution rate at different values for the analcime precipitation rate constant and the effects of varying the glass dissolution rate constant at a constant analcime precipitation rate constant. From the simulations we conclude, firstly, that the rate of glass dissolution is dependent on the kinetics of formation of the zeolitic phase. Therefore, the kinetics of secondary phase formation is an important parameter that should be taken into account in future glass dissolution modeling efforts. Secondly, the results indicate that, in the absence of a gel layer, the glass dissolution rate controls the rate of analcime precipitation in the long term. Finally, the meaning of these results pertinent to long-term glass durability is discussed.

analcime↗

Directing Charge Transfer in Quantum Dot Assemblies

The optical and electronic properties of semiconductor quantum dots (QDs) make them attractive candidates for applications in photovoltaics, spintronics, photocatalysis, and optoelectronics. Understanding how to control the flow of charge in QD assemblies is essential for realizing novel applications. This Account explores some unique characteristics of charge transport in QD dyads, triads, and their assemblies. The emerging features of these assemblies that provide new opportunities to manipulate charge flow at the nanoscale are (1) cascading energy landscapes and band offsets to inhibit charge recombination, (2) electrostatic fields that direct charge flow through QD–QD and QD-conjugated polymer junctions, and (3) QD chirality and chiral imprinting that promotes vectorial electron and spin selective transport. Charge flow kinetics is determined by a combination of familiar electron transfer parameters (reaction free energy, reorganization energy, and electronic coupling), donor and acceptor electronic densities of states, and internal electric fields. Electron transfer and electronic structure theory, combined with kinetic modeling, place the measured kinetics of QD electron transfer donor–acceptor assemblies into a unified conceptual context. The experimental transfer rates measured in these systems depend upon structure and the internal electric fields that are present in the assemblies. A negatively charged donor and positively charged acceptor, for example, facilitates (inhibits) electron (hole) transfer, while an electric field of opposite orientation (reversal of charges) inhibits (promotes) electron (hole) transfer. These and other emerging rules that govern charge flow in NP assemblies provide a strategy to design the directionality and yield of interfacial charge transport. Chirality at the nanoscale can induce spin selective charge transport, providing new ways to direct charge (and spin) flow in QD assemblies. Magnetoresistance and magnetic conductive probe atomic force microscopy experiments show spin selective electron transport for chirally imprinted QD assemblies. Photoinduced electron transfer from achiral donor-QDs to chiral acceptor-QDs depends on the electron spin and chiroptical properties of the acceptor-QDs. These assemblies show transport characteristics that correlate with features of the QDs’ circular dichroism spectra, presenting intriguing challenges to theory, and indicating that spectroscopic signatures may assist in the design and diagnosis of functional molecular assemblies. Theoretical and experimental studies of charge transport in well-defined QD assemblies are establishing design principles for vectorial charge transport and are also refining questions surrounding the mechanism and control of these processes. Finally, these intensified efforts are forging links between fundamental discoveries regarding mechanism and practical applications for these novel assembled nanostructures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

QUANTIFYING OBSERVATIONAL PROJECTION EFFECTS USING MOLECULAR CLOUD SIMULATIONS

The physical properties of molecular clouds are often measured using spectral-line observations, which provide the only probes of the clouds' velocity structure. It is hard, though, to assess whether and to what extent intensity features in position-position-velocity (PPV) space correspond to 'real' density structures in position-position-position (PPP) space. In this paper, we create synthetic molecular cloud spectral-line maps of simulated molecular clouds, and present a new technique for measuring the reality of individual PPV structures. Using a dendrogram algorithm, we identify hierarchical structures in both PPP and PPV space. Our procedure projects density structures identified in PPP space into corresponding intensity structures in PPV space and then measures the geometric overlap of the projected structures with structures identified from the synthetic observation. The fractional overlap between a PPP and PPV structure quantifies how well the synthetic observation recovers information about the three-dimensional structure. Applying this machinery to a set of synthetic observations of CO isotopes, we measure how well spectral-line measurements recover mass, size, velocity dispersion, and virial parameter for a simulated star-forming region. By disabling various steps of our analysis, we investigate how much opacity, chemistry, and gravity affect measurements of physical properties extracted from PPV cubes. For the simulations used here, which offer a decent, but not perfect, match to the properties of a star-forming region like Perseus, our results suggest that superposition induces a ∼40% uncertainty in masses, sizes, and velocity dispersions derived from {sup 13}CO (J = 1-0). As would be expected, superposition and confusion is worst in regions where the filling factor of emitting material is large. The virial parameter is most affected by superposition, such that estimates of the virial parameter derived from PPV and PPP information typically disagree by a factor of ∼2. This uncertainty makes it particularly difficult to judge whether gravitational or kinetic energy dominate a given region, since the majority of virial parameter measurements fall within a factor of two of the equipartition level α ∼ 2.

79 ASTRONOMY AND ASTROPHYSICS↗

Heat- and light-induced transformations of Yb trapping sites in an Ar matrix

The low-lying electronic states of Yb isolated in a solid Ar matrix grown at 4.2 K are characterized through absorption and emission spectroscopy. Yb atoms are found to occupy three distinct thermally stable trapping sites labeled “red,” “blue,” and “violet” according to the relative positions of the absorption features they produce. Classical simulations of the site structure and relative stability broadly reproduced the experimentally observed matrix-induced frequency shifts and thus identified the red, blue, and violet sites as due to respective single substitutional (SS), tetravacancy (TV), and hexavacancy (HV) occupation. Prolonged excitation of the {sup 1}S → {sup 1}P transition was found to transfer the Yb population from HV sites into TV and SS sites. The process showed reversibility in that annealing to 24 K predominantly transferred the TV population back into HV sites. Population kinetics were used to deduce the effective rate parameters for the site transformation processes. Experimental observations indicate that the blue and violet sites lie close in energy, whereas the red one is much less stable. Classical simulations identify the blue site as the most stable one.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Sensitivity Analysis of FEAST-Metal Fuel Performance Code: Initial Results

This memo documents the completion of the LANL milestone, M3FT-12LA0202041, describing methodologies and initial results using FEAST-Metal. The FEAST-Metal code calculations for this work are being conducted at LANL in support of on-going activities related to sensitivity analysis of fuel performance codes. The objective is to identify important macroscopic parameters of interest to modeling and simulation of metallic fuel performance. This report summarizes our preliminary results for the sensitivity analysis using 6 calibration datasets for metallic fuel developed at ANL for EBR-II experiments. Sensitivity ranking methodology was deployed to narrow down the selected parameters for the current study. There are approximately 84 calibration parameters in the FEAST-Metal code, of which 32 were ultimately used in Phase II of this study. Preliminary results of this sensitivity analysis led to the following ranking of FEAST models for future calibration and improvements: fuel conductivity, fission gas transport/release, fuel creep, and precipitation kinetics. More validation data is needed to validate calibrated parameter distributions for future uncertainty quantification studies with FEAST-Metal. Results of this study also served to point out some code deficiencies and possible errors, and these are being investigated in order to determine root causes and to improve upon the existing code models.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Data for KETCHUP: Parameterizing of Large-Scale Kinetic Models Using Multiple Datasets with Different Reference States

Repository for Kinetic Estimation Tool Capturing Heterogeneous Datasets Using Pyomo (KETCHUP), a flexible parameter estimation tool that leverages a primal-dual interior-point algorithm to solve a nonlinear programming (NLP) problem that identifies a set of parameters capable of recapitulating the steady-state fluxes and concentrations in wild-type and perturbed metabolic networks. KETCHUP can use K-FIT [2] input files. Example K-FIT input files are located in the K-FIT repository at https://github.com/maranasgroup/K-FIT.

Metabolomics↗

Kinklike structures in an arcsin real scalar dynamics

Highlights: • Kinklike analyses in a real scalar theory with an arcsin dynamic. • First-order formalism implementation and deformation procedure in an arcsin scalar dynamics. In this paper, we analyze kink-like analytical solutions in a real scalar theory with an arcsin dynamics inspired by the arcsin electrodynamics presented in Kruglov (2015). This analysis is done by means of the first-order formalism. This formalism provides a framework where the equations of motion can be simplified by preserving the linear stability of the theory. In this work, the deformation procedure is implemented with the aim of finding exact solutions in systems with generalized dynamics. Along the paper, we explore how the first-order formalism is implemented in the arcsin kinetics and how such a term influences the kink-like solutions. As a part of the result of our paper, we show that the kink-like solutions are similar to the ones obtained in the standard scalar kinetic theory. We also show that the extra parameter, that controls the non-linearities of the model, plays an essential role in the energy densities and stability potentials. These quantities vary according to this parameter. The goals here are to show how the first-order framework is implemented in this arcsin scenario and to present the analytical kink-like solutions that can be found by means of the first-order framework and deformation method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Molecular Dynamics Study of Combustion Reactions in a Supercritical Environment. Part 2: Boxed MD Study of CO + OH → CO 2 + H Reaction Kinetics

Oxy-fuel combustion technology holds a great promise in both increasing the efficiency of the energy conversion and reducing environmental impact. Yet, effects of the higher pressures and replacement of the nitrogen with carbon dioxide diluent are not well understood at present. The title reaction is one of the most important processes in combustion. Despite numerous studies, the effects of supercritical carbon dioxide environment did not receive much attention in the past. In this work, we report the results of boxed molecular dynamics simulations of these effects at QM/MM theory level with periodical boundary conditions. The free energy barriers for HOCO intermediate formation and decomposition were tabulated in a wide range of pressures (1–1000 atm) and temperatures (400–1600 K). Pressure dependence of calculated rate constants for these reaction steps and overall reaction were analyzed. We discovered that the CO 2 environment may increase these rate constants up to a factor of 25, at near critical conditions. At higher temperatures, this effect weakens significantly. Numerical values for parameters of extended Arrhenius equation, suitable for combustion kinetic modeling are reported.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Cellular Automaton Modeling of Dynamic Recrystallization of Nimonic 80A Superalloy Based on Inhomogeneous Distribution of Dislocations Inside Grains

A modified cellular automaton (CA) model in the time and space scale is developed to simulate the dynamic recrystallization (DRX) mechanism of Nimonic 80A nickel-base superalloy during hot deformation. Some necessary thermo-mechanical parameters used as input data in CA approach such as work hardening and softening coefficient are fitted according to the previous hot compression experiments. Based on the experimental results, nucleation kinetics is modified as a function of Zener-Holloman parameter. A dislocation density model at transient recrystallization boundary is proposed, in which the dislocation density is a function of temperature, strain rate and annealed dislocation density. The simulated results are in good agreement with experimental results, proving that the modified CA model is capable of predicting both flow behavior and microstructural evolution of Nimonic 80A during hot deformation. According to the experimental and simulated results, the DRX behavior is discussed in detail. It demonstrates that the flow stress and microstructural evolution of DRX have a relationship with the deformation temperature and strain rate.

36 MATERIALS SCIENCE↗

Lutetium(III) aqua ion: On the dynamical structure of the heaviest lanthanoid hydration complex

The structure and dynamics of the lutetium(III) ion in aqueous solution have been investigated by means of a polarizable force field molecular dynamics (MD). An 8-fold square antiprism (SAP) geometry has been found to be the dominant configuration of the lutetium(III) aqua ion. Nevertheless, a low percentage of 9-fold complexes arranged in a tricapped trigonal prism (TTP) geometry has been also detected. Dynamic properties have been explored by carrying out six independent MD simulations for each of four different temperatures: 277 K, 298 K, 423 K, 632 K. The mean residence time of water molecules in the first hydration shell at room temperature has been found to increase as compared to the central elements of the lanthanoid series in agreement with previous experimental findings. Water exchange kinetic rate constants at each temperature and activation parameters of the process have been determined from the MD simulations. The obtained structural and dynamical results suggest that the water exchange process for the lutetium(III) aqua ion proceeds with an associative mechanism, in which the SAP hydration complex undergoes temporary structural changes passing through a 9-fold TTP intermediate. Such results are consistent with the water exchange mechanism proposed for heavy lanthanoid atoms.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Bulk Attenuation in the Earth’s Mantle Due to Phase Transitions

The calculated bulk attenuation due to phase transition from mineral physics data is reported here. With relaxation time less than 1 s, the calculated value for a pyrolite mantle is consistent with the inverted bulk attenuation of the upper mantle from seismic observations. The two important mechanisms of phase transitions, diffusion-controlled and kinetics-controlled, have different relaxation time as indicated by the models here. The diffusion controlled is more likely to contribute to the observed seismic bulk attenuation than the kinetic-controlled process based on the available diffusivity and kinetics data. The correlation between the bulk attenuation and relaxation time emphasizes the importance of a number of parameters in the mineral physics database such as Fe-Mg diffusivity and kinetics in olivine-wadsleyite-ringwoodite-perovskite, Mg-Ca-Al-Si diffusivity and kinetics in pyroxene-garnet-Ca perovskite, some of which are still unknown.

58 GEOSCIENCES↗