Search NASA⌕ Search

SEARCH · Search NASA

Results for “rate constants”

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
DOE OSTI2021

Probing the thermodynamics and kinetics of ethylene carbonate reduction at the electrode–electrolyte interface with molecular simulations

Understanding the formation of the solid-electrolyte interphase (SEI) in lithium-ion batteries (LIBs) is an ongoing area of research due to its high degree of complexity and the difficulties encountered by experimental studies. Herein, we investigate the initial stage of SEI growth, the reduction reaction of ethylene carbonate (EC), from both a thermodynamic and kinetic approach with theory and molecular simulation. We employed both the potential distribution theorem (PDT) and the SMD implicit solvent model to EC solvation for the estimation of reduction potentials of Li$^+$, EC, and Li$^+$-solvating EC (s-EC), as well as reduction rate constants of EC and s-EC. We find that solvation effects greatly influence these quantities of interest, particularly the Li$^+$/Li reference electrode potential in EC solvent. Further, we also compute the inner- and outer-sphere reorganization energies for both EC and s-EC at the interface of liquid EC and a hydroxyl-terminated graphite surface, where total reorganization energies are predicted to be 76.6 and 88.9 kcal/mol, respectively. With the computed reorganization energies, we estimate reduction rate constants across a range of overpotentials and show that EC has a larger electron transfer rate constant than s-EC at equilibrium, despite s-EC being more thermodynamically favorable. Furthermore, this manuscript demonstrates how ion solvation effects largely govern the prediction of reduction potentials and electron transfer rate constants at the electrode-electrolyte interface.

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

Revised estimates for ozone reduction by shuttle operation

Previous calculations by five different modeling groups of the effect of space shuttle operations on the ozone layer yielded an estimate of 0.2 percent ozone reduction for the Northern Hemisphere at 60 launches per year. Since these calculations were made, the accepted rate constant for the reaction between hydroperoxyl and nitric oxide to yield hydroxyl and nitrogen dioxide, HO2 + NO yields OH + NO2, was revised upward by more than an order of magnitude, with a resultant increase in the predicted ozone reduction for chlorofluoromethanes by a factor of approximately 2. New calculations of the shuttle effect were made with use of the new rate constant data, again by five different modeling groups. The new value of the shuttle effect on the ozone layer was found to be 0.25 percent. The increase resulting from the revised rate constant is considerably less for space shuttle operations than for chlorofluoromethane production, because the new rate constant also increases the calculated rate of downward transport of shuttle exhaust products out of the stratosphere.

Potter, A. E.↗
NASA NTRS1988

Kinetics of the reaction OH + HO2 yields H2O + O2 from 254 to 382 K

The discharge-flow resonance fluorescence technique has been used to determine the absolute rate constant for the reaction OH + HO2 yields H2O + O2 from 254 to 382 K at a total pressure of 1 Torr. Pseudo-=first-order conditions were used with HO2 in large excess over OH. The rate constant was obtained directly from observed decays of OH and measured concentrations of HO2. Since the observed rate constant was found to be very sensitive to small background concentrations of O and H atoms, NO2 was used to remove both atom species from the system. With added NO2 the result at 299 K is (1.1 + or - 0.3) x 10 to the -10th cu cm/molecule where the error limits are one standard deviation and include an estimate of overall experimental uncertainty. The temperature dependence expressed in Arrhenius form is (4.8 + or - 0.8) x 10 to the -11th exp (250 + or - 50)/T. The results are independent of the type of reactor surface and the precursor used to produce OH and HO2. The present results agree well with earlier measurements near 1-atm total pressure and suggest that this rate constant exhibits little or no pressure dependence between 1 and 1000 Torr.

Keyser, Leon F.↗
NASA NTRS1977

Absolute rate of the reaction of Cl/2P/ with methane from 200-500 K

Rate constants for the reaction Cl + CH4 yields HCl + CH3 were measured in the range 200-500 K by the plot flash photolysis-resonance fluorescence technique. An Arrhenius plot of the mean rate constant observed at 14 equally spaced temperatures is best fit by the expression k = 5.44 times 10 to the minus 19th times T to the exponent 2.50 times exp(-608/T) cu cm per molecule per sec. Appropriate expressions are also given for the rate constant in the two ranges 200-299 K and 299-500 K. These results are compared with a number of other existing determinations of the rate constant by other methods and in different temperature ranges. Most discrepancies are attributed to the nonlinear Arrhenius behavior.

Whytock, D. A.↗
DOE OSTI2023

Theory-Based Mechanism for Fluoromethane Combustion I: Thermochemistry and Abstraction Reactions

A new detailed chemical kinetic mechanism is presented for small fluorinated hydrocarbons. Ab initio electronic structure theory is used to provide heats of formation with subchemical accuracy. The ANL0 method is extended to include fluorine. The resulting heats of formation at 0 K are in excellent agreement with 36 benchmark species in the Active Thermochemical Tables, with a mean error of μ = –0.02 kJ/mol and a standard deviation of σ = 0.91 kJ/mol. The thermophysical properties for 92 small-molecule H/C/O/F species are computed. The rate coefficients for 40+ H-abstraction reactions involving H, O, F, OH, OF, HO 2 , and various methyl radicals with CH 4 , CH 3 F, CH 2 F 2 , CHF 3 , CH 2 O, and CHFO are discussed. Here, the computed rate constants are in excellent agreement with the available literature. Additionally, 30+ rate constants are provided for F abstraction, which are several orders of magnitude smaller than H abstraction. The thermophysical properties and rate constants are provided in a mechanism. This mechanism is the first in a series of theory-based investigations into the thermal destruction of per- and polyfluorinated species.

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

Temperature dependence of the reaction O/3P/ + OH/2 Pi/ yields O2 + H

Measurements of the absolute rate and temperature dependence of the rate constant for the reaction of ground state atomic oxygen with hydroxyl radicals to yield molecular oxygen and hydrogen atoms over the temperature range 211-499 K are presented. The reaction was monitored in a low-pressure discharge flow resonance fluorescence apparatus under pseudo-first-order conditions in which the oxygen concentration was greater than the hydroxyl concentration. An expression for the temperature dependence of the rate constant is derived from a least-squares fit to the observed bimolecular rate constants, which were found to range from 3.21 to 2.77 x 10 to the -11th cu cm/molecule per sec for temperatures from 211 to 499 K, indicative of a slight negative temperature dependence. Results are compared with those of previous workers, and the significance of the lower value of the rate constant obtained in the present study to models of HO(x) in the upper stratosphere is indicated.

Lewis, R. S.↗
DOE OSTI2026

Investigating shock-induced chemical reactions in Ni+Al multilayers: A continuum-based mesoscale approach with Arrhenius kinetics and artificial thermal conduction

A new continuum-based mesoscale modeling approach for shock-induced chemical reactions (SICRs) in Ni+Al multilayers is demonstrated in Sandia’s shock physics hydrocode, CTH. The approach utilizes Arrhenius-type kinetics and artificial thermal conduction. Our work builds upon previous efforts to parameterize equations of state for Ni x Al y [J. Appl. Phys. 137, 075102 (2025)], as well as simulations of inert shocks in realistic 2D microstructures [J. Appl. Phys. 137, 225301 (2025)]. To calibrate the reaction kinetics, pairs of the reaction coordinate, R′, vs time are extracted from the molecular dynamics (MD) literature. Here, the MD-informed kinetics are used to simulate the dynamic evolution of pressure and temperature in 2D mesoscale simulations. Overall, the MD-informed kinetics obtained for planar interfaces are too slow, as initial reaction is not observed on a nanosecond time scale. Even with quasi-periodic shock focusing leading to the formation of so-called “hot-spots,” the hot spots are unable to grow and coalesce using the fitted Arrhenius rate constants. However, by increasing the rate constants by two orders of magnitude, SICRs are observed at a shock pressure near 30 GPa, which is supported by experiments. Consequently, these mesoscale simulations suggest that unresolved shear-based mechanical mixing might possibly account for the discrepancies in kinetic rates, with shock-generated intense perturbations, interfacial vortical flows, and elevated temperatures serving as favorable reaction conditions. Future work will calibrate a shear-dependent reaction rate from the MD simulations with realistic interfaces that are wavy, diffuse, and disordered.

Kittell, David E. [Sandia National Laboratories (S↗
DOE OSTI2022

Bay capping via acetylene addition to polycyclic aromatic hydrocarbons: Mechanism and kinetics

The bay-capping mechanism on PAH armchair edges and the kinetics of acetylene addition to 6–6–5 and 5–6–5 bays have been explored by ab initio/RRKM-ME calculations. The bays on the edges were modeled by C 21 H 11 and C 20 H 9 radicals produced by H abstractions from 7H-benzo[c]cyclopenta[e]pyrene and dicyclopenta[cf]pyrene. The C 20 H 9 + C 2 H 2 reaction is shown to have a low entrance barrier and to rapidly form the capped product, indaceno[2,1,8,7-cdefg]pyrene, along with ethynyl substituted dicyclopenta[cf]pyrene at temperatures above 1400 K. The reactivity of C 21 H 11 is shown to be governed by the location of the unpaired electron; the π radical R1 formed by H abstraction from the CH 2 group in 7H-benzo[c]cyclopenta[e]pyrene reacts with C 2 H 2 very slowly owing to a high entrance barrier, with the bay-capping rate constant approaching 10 -16 cm 3 molecule -1 s -1 only at temperatures above 2000K. Further, this result reaffirms that the growth of π aryl radicals via acetylene addition is inefficient and reflects the generally low reactivity of such radicals where the spin density is highly delocalized over the entire polyaromatic system. Alternatively, the σ C 21 H 11 radical R2 produced by H abstraction from the five-membered ring at the bay rapidly reacts with C 2 H 2 forming the bay-capped product, with the rate constant on the order of 10 -12 cm 3 molecule -1 s -1 at T≥1500K. Rate constants for the capping reactions at the 6–6–5 and 5–6–5 bays are compared with those at the 6–0–6, 6–6–6, and 6–5–6 bays. The site-specific bay-capping rate constants have been utilized in kMC simulations of the PAH growth and the results showed measurable differences when the 6–6–5 and 5–6–5 bay-capping reactions are taken into account, including an increase of the growth rate and the formation of closed-shell PAH and a rise of the number of embedded five-membered rings accompanied with a slight decrease of their overall amount.

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

Modeling nanoparticle charge distribution in the afterglow of non-thermal plasmas and comparison with measurements

Particle charging in the afterglows of non-thermal plasmas typically take place in a non-neutral space charge environment. We model the same by incorporating particle-ion collision rate constant models, developed in prior work by analyzing particle-ion trajectories calculated using Langevin Dynamics simulations, into species transport equations for ions, electrons and charged particles in the afterglow. A scaling analysis of particle charging and additional Langevin Dynamics calculations of the particle-ion collision rate constant are presented to extend the range of applicability to ion electrostatic to thermal energy ratios of 300 and diffusive Knudsen number (that scales inversely with gas pressure) up to 2000. The developed collision rate constant models are first validated by comparing predictions of particle charge against measured values in a stationary, non-thermal DC plasma from past PK-4 campaigns published in Phys. Rev. Lett. 93(8): 085001 and Phys. Rev. E 72(1): 016406). The comparisons reveal excellent agreement within ±35% for particles of radius 0.6,1.0,1.3 μm in the gas pressure range of ~20-150 Pa. The experiments to probe particle charge distributions by Sharma et al. (J. Physics D: Appl. Phys. 53(24): 245204) are modeled using the validated particle-ion collision rate constant models and the calculated charge fractions are compared with measurements. The comparisons reveal that the ion/electron concentration and gas temperature in the afterglow critically influence the particle charge and the predictions are generally in qualitative agreement with the measurements. Along with critical assessment of the modeling assumptions, several recommendations are presented for future experimental design to probe charging in afterglows.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

N-Terminal Decarboxylation as a Probe for Intramolecular Contact Formation in $γ$-Glu-(Pro) n -Met Peptides

The kinetics of intramolecular-contact formation between remote functional groups in peptides with restricted conformational flexibility were examined using designed peptides with variable-length proline bridges. As probes for this motion, free radicals were produced using the • OH-induced oxidation at the C-terminal methionine residue of γ-Glu-(Pro) n -Met peptides ( n = 0–3). The progress of the radicals’ motion along the proline bridges was monitored as the radicals underwent reactions along the peptides’ backbones. Of particular interest was the reaction between the sulfur atom located in the side chain of the oxidized Met residue and the unprotonated amino group of the glutamic acid moiety. Interactions between them were probed by the radiation-chemical yields (expressed as G values) of the formation of C-centered, α-aminoalkyl radicals (αN) on the Glu residue. These radicals were monitored directly or via their reaction with p -nitroacetophenone (PNAP) to generate the optically detected PNAP •– radical anions. The yields of these αN radicals were found to be linearly dependent on the number of Pro residues. A constant decrease by 0.09 μM J –1 per spacing Pro residue of the radiation-chemical yields of G (αN) was observed. Previous reports support the conclusion that the αN radicals in these cases would have to result from (S$\therefore$N) + -bonded cyclic radical cations that arose as a result from direct contact between the ends of the peptides. Furthermore, by analogy with the rate constants for the formation of intramolecularly (S$\therefore$S) + -bonded radical cations in Met-(Pro) n -Met peptides ( J. Phys. Chem. B 2016 , 120 , 9732), the rate constants for the formation of intramolecularly (S$\therefore$N) + -bonded radical cations are activated to the same extent for all of the γ-Glu-(Pro) n -Met peptides. Thus, the continuous decrease of G (αN) with the number of Pro residues (from 0 to 3) suggests that the formation of a contact between the S-atom in the C-terminal Met residue and the N-atom of a deprotonated N-terminal amino group of Glu is controlled in peptides with 0 to 3 Pro residues by the relative diffusion of the S •+ and unoxidized N-atom. The overall rate constants of cyclization to form the (S$\therefore$N)-bonded radical cations were estimated to be 3.8 × 10 6 , 1.8 × 10 6 , and 8.1 × 10 5 s –1 for peptides with n = 0, 1, and 2 Pro residues, respectively. If activation is the same for all of the peptides, then these rate constants are a direct indication for the end-to-end dynamics along the chain.

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

Nonequilibrium boundary layer at a stagnation point for a hydrogen-helium stream over ablating graphite

The nonequilibrium axisymmetric stagnation point boundary layer over an ablating graphite surface is considered. The external stream is a high temperature mixture of hydrogen and helium. Variable thermodynamic and transport properties are assumed. Lennard-Jones potential model is used to calculate the transport coefficients of each species. Although the mixture rules for viscosity of the gas mixture are used, the weighting functions are more sophisticated than those commonly employed. For the conductivity of the mixture, generalized Wassiljewa coefficients are used. Seven species with 28 dissociation/recombination reactions are considered. Hansen's model for the dissociation rate constants is employed. The recombination rate constants are obtained by invoking detailed balance principles assisted by the JANAF thermodynamic data and the Hansen-Pearson thermodynamic data for C3.

Liu, T.-M.↗
DOE OSTI2022

Plutonium reactive transport in fractured granite: Multi-species experiments and simulations

Plutonium (Pu) in the subsurface environment can transport in different oxidation states as an aqueous solute or as colloidal particles. The transport behavior of Pu is affected by the relative abundances of these species and can be difficult to predict when they simultaneously exist. This study investigates the concurrent transport of Pu intrinsic colloids, Pu(IV) (aq) and Pu(V-VI) (aq) through a combination of controlled experiments and semi-analytical dual-porosity transport modeling. Pu transport experiments were conducted in a fractured granite at high and low flow rates to elucidate sorption processes and their scaling behavior. In the experiments, Pu(IV) (aq) was the least mobile of the Pu species, Pu(V-VI) (aq) had intermediate mobility, and the colloidal Pu, which consisted mainly of precipitated and/or hydrolyzed Pu(IV), was the most mobile. The semi-analytical modeling revealed that the sorption of each Pu species was rate-limited, as the sorption could not be described by assuming local equilibrium in the experiments. The model was able to describe the sorption of the different Pu species that occurring either on fracture surfaces, in the pores of the rock matrix, or simultaneously in both locations. While equally good fits to the data could be achieved using any of these assumptions, a fracture-dominated process was considered to be the most plausible because it provided the most reasonable estimates of sorption rate constants. Importantly, a key result of this work is that the sorption rate constant of all Pu species tends to decrease with increasing time scales, which implies that Pu will tend to be more mobile at longer time scales than observations at shorter time scales suggest. This result has important implications for predicting the environmental impacts of Pu in the safety assessments of geologic repositories for radioactive waste disposal, and we explore potential mechanistic bases for upscaling the sorption rate constants to time and distance scales that cannot be practically evaluated in experiments.

58 GEOSCIENCES↗
DOE OSTI2023

Constraining the Reaction Rate of Criegee Intermediates with Carboxylic Acids during the Multiphase Ozonolysis of Aerosolized Alkenes

Criegee intermediates (CIs) are central to the ozonolysis of unsaturated organic compounds, controlling the formation of small fragmentation reaction products and higher molecular weight oligomeric compounds through competing unimolecular and bimolecular pathways. In particular, the reaction of thermalized CI with carboxylic acids (RCOOH) to produce oligomeric α-acyloxyalkyl hydroperoxides (AAHPs) is of interest to the chemical transformation of organic compounds in the atmosphere. In the gas phase, the CI + RCOOH reaction proceeds at supercollisional rates, consistent with a barrierless reaction mechanism. However, the rate of this reaction in condensed phases remains uncertain. Here, we report the results of aerosol flow tube studies of the ozonolysis of submicron organic aerosol containing mixtures of an internal n-alkene (Z-9-tricosene, tri) and a saturated acid (2-hexyldecanoic acid, HDA). The decay of the reactants and formation of reaction products are monitored using atmospheric pressure chemical ionization mass spectrometry (APCI–MS). A 6-fold decrease in the reactive uptake of ozone is observed with increasing acid concentration, consistent with a previous study using alcohol additives. The decay of HDA with ozone exposure shows that HDA scavenges CIs at rates comparable to competing unimolecular pathways involving CIs. The decay kinetics of HDA are analyzed using a simple kinetic model in order to constrain the ratio of the CI + RCOOH and unimolecular CI rate constants. Empirical fitting to this model, combined with theoretical estimates of the unimolecular loss rate, yields a CI + RCOOH rate constant of 1.85 ± 0.27 × 10 –19 cm 3 molec –1 s –1 , which is six orders of magnitude slower than the expected diffusion limit in a liquid organic matrix. Stochastic kinetic simulations are conducted with different values of the CI + RCOOH rate constant to validate this result and show that the diffusion-limited rate is indeed much too rapid to reproduce the experimental results.

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

Non-Boltzmann Effects in Chain Branching and Pathway Branching for Diethyl Ether Oxidation

Low-temperature (LT) engine applications have several potential benefits, including reduced emissions and increased efficiency. Attaining these benefits requires accurate kinetic modeling of LT chain branching, which depends heavily on ketohydroperoxide (KHP) decomposition. For diethyl ether (DEE), a promising biofuel, current estimates of the KHP decomposition rate constant are largely based on empirical fits to data. In this study, we investigate the most important KHP isomer in DEE LT oxidation by applying variable reaction coordinate transition state theory to the main pathway for KHP decomposition: OO bond fission to produce •OH and a keto-alkoxy radical, •OQ'O. We also use ab initio kinetics methods to investigate the decomposition of •OQ'O, where we find dominant branching to acetic acid, with the remaining flux going to CH 3 C(O)OCHO. Additionally, new time-resolved measurements of DEE and acetic acid concentrations during LT (450–600 K) DEE oxidation are obtained in a laser photolysis flow reactor coupled with multiplexed photoionization mass spectrometry. These new experimental data, along with jet-stirred reactor data in the literature, are compared with the predictions of a recent DEE mechanism (Tran et al. Proc. Comb. Inst. 2019, 37, 511-519) that was modified with the newly calculated ab initio rate constants for KHP and •OQ'O decomposition. The predictions of the modified mechanism are quite poor when compared to the experimental data; this is primarily due to the new KHP ⇌ •OQ'O + •OH rate constant, which is 1–2 orders of magnitude slower than empirical values employed in recent mechanisms. To reconcile the new KHP rate constant and the experimental data, we explore and quantify the possible role of non-Boltzmann (nB) reaction sequences. The nB reactions have a substantial effect on both the overall mechanism reactivity and the •OQ'O branching to acetic acid. We also provide guidance on the proper implementation of nB reactions in kinetic mechanisms.

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

Laser flash photolysis of ozone - O/1D/ quantum yields in the fall-off region 297-325 nm

The wavelength dependence of the quantum yield for O(1D) production from ozone photolysis has been determined between 297.5 nm and 325 nm in order to resolve serious discrepancies among previous studies. The results of this investigation are compared to earlier work by calculating atmospheric production rate constants for O(1D). It is found that for the purpose of calculating this rate constant, there is now good agreement among three studies at 298 K. Furthermore, it appears that previous data on the temperature dependence of the O(1D) quantum yield fall-off is adequate for determining the vertical profile of the O(1D) production rate constant. Several experimental difficulties associated with using NO2(asterisk) chemiluminescence to monitor O(1D) have been identified.

Brock, J. C.↗
NASA NTRS1979

Heterogeneous atmospheric reactions - Sulfuric acid aerosols as tropospheric sinks

The reaction probabilities of various atmospheric species incident on a bulk sulfuric acid surface are measured in order to determine the role of sulfuric acid aerosols as pollutant sinks. Reaction products and unreacted starting materials leaving a Knudsen cell flow reactor after collision at 300 K with a H2SO4 surface or a soot surface were detected by mass spectrometry. Significant collision reaction probabilities are observed on a H2SO4 surface for H2O2, HNO3, HO2NO2, ClONO2, N2O5, H2O and NH3, and on soot for NH3. Estimates of the contribution of heterogeneous reactions to pollutant removal under atmospheric conditions indicate that while aerosol removal in the stratosphere is insignificant (loss rate constants approximately 10 to the -10th/sec), heterogeneous reactions may be the dominant loss process for several tropospheric species (loss rate constant approximately 10 to the -5th/sec, comparable to photolysis rate constants).

Baldwin, A. C.↗
DOE OSTI2024

Reactions of the ˙NO 3 radical with nuclear extraction ligands in alkane solution

Here, the reactivity of the nitrate radical (NO 3 ˙) with organophosphorus and amidic actinide and lanthanide complexing agents of interest to nuclear solvent extraction applications was measured, resulting in the first-ever reported bimolecular rate constants for this radicals’ reactions in dodecane solution. The order of reactivity for neutral organophosphorus compounds showed faster rate constants with increasing electron density on the phosphoryl phosphorus atom, indicating an increasing facility for electron abstraction reactions occurring in addition to H-atom abstraction from the ligand alkane chains. The only acidic organophosphorus compound investigated, HEH[EHP], showed low reactivity with the NO 3 ˙ radical, attributed to its dimerization in this non-polar solvent. Amide ligand reaction rates were faster than for organophosphorus molecules, suggesting more facile H-atom abstraction from carbonyl activated methylene and amyl groups. While all rate constants were slower than the diffusion-limited rate they were still rapid enough to result in significant oxidation of solvent extraction ligands in dodecane solution.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
NASA NTRS1981

Analysis of the tensile stress-strain behavior of elastomers at constant strain rates. I - Criteria for separability of the time and strain effects

A theoretical analysis of the tensile stress-strain relation of elastomers at constant strain rate is presented which shows that the time and the stress effect are separable if the experimental time scale coincides with a segment of the relaxation modulus that can be described by a single power law. It is also shown that time-strain separability is valid if the strain function is linearly proportional to the Cauchy strain, and that when time-strain separability holds, two strain-dependent quantities can be obtained experimentally. In the case where time and strain effect are not separable, superposition can be achieved only by using temperature and strain-dependent shift factors.

Hong, S. D.↗