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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

The Reaction Rates of Amidogen and Ammonia with Nitrous Oxide: Implications for Combustion Mechanisms

Pulsed laser photolysis experiments with laser-induced fluorescence detection of NH 2 set an upper limit to the rate constant for reaction with N 2 O of k < 1 × 10 −15 cm 3 molecule −1 s −1 at 513 K. Computations were based on geometries and anharmonic frequency analysis (B2PLYP-D3/cc-pVTZ) followed by coupled cluster calculations extrapolated to the infinite basis set limit, with corrections for core−valence electron correlation, scalar relativistic effects, and correlation up to CCSDT(Q). Species that showed multireference character were quantified with MRCI(7,7)+Q/cc-pVTZ theory. Rate constants were obtained for the dominant product channel H 2 NN + NO, along with HN 3 + OH, H 2 NO + N 2 , NNH + HNO, ON(NH) 2 and HNNH + NO. The last channel is slow even at 2500 K, contrary to an early empirical estimate and confirming recent suggestions. Modeling of literature experiments on oxidation of NH 3 by N 2 O shows that all channels are too slow to make a significant impact on the loss of N 2 O in ammonia flames. Similarly, the direct NH 3 + N 2 O reaction is found to be negligibly slow.

Ammonia↗
DOE OSTI2025

Ab Initio Molecular Dynamics Study of the Reduction of Acetone by the Hydrated Electron

We have investigated the reaction dynamics of the reduction of acetone by the solvated electron in water using ab initio molecular dynamics simulations at 298 and 373 K. The rate constants derived from the simulations are consistent with experimental observations that the reaction has a low activation energy. Detailed analyses carried out to shed light on the mechanism of the reaction show that solvent reorganization plays a key role as a reaction coordinate, as expected from the Marcus electron transfer theory. Furthermore, constrained density functional theory calculations indicate that the electronic coupling is large, placing the reaction in the adiabatic limit. Indeed, the activation energies and rate constants of the simulations are in accordance with the predictions of adiabatic Marcus theory.

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

Data analytics approach to predict high-temperature cyclic oxidation kinetics of NiCr-based Alloys

Abstract Although of practical importance, there is no established modeling framework to accurately predict high-temperature cyclic oxidation kinetics of multi-component alloys due to the inherent complexity. We present a data analytics approach to predict the oxidation rate constant of NiCr-based alloys as a function of composition and temperature with a highly consistent and well-curated experimental dataset. Two characteristic oxidation models, i.e., a simple parabolic law and a statistical cyclic oxidation model, have been chosen to numerically represent the high-temperature oxidation kinetics of commercial and model NiCr-based alloys. We have successfully trained machine learning (ML) models using highly ranked key input features identified by correlation analysis to accurately predict experimental parabolic rate constants ( k p ). This study demonstrates the potential of ML approaches to predict oxidation kinetics of alloys over wide composition and temperature ranges. This approach can also serve as a basis for introducing more physically meaningful ML input features to predict the comprehensive cyclic oxidation behavior of multi-component high-temperature alloys with proper constraints based on the known underlying mechanisms.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Hierarchical memories: Simulating quantum LDPC codes with local gates

Constant-rate low-density parity-check (LDPC) codes are promising candidates for constructing efficient fault-tolerant quantum memories. However, if physical gates are subject to geometric-locality constraints, it becomes challenging to realize these codes. In this paper, we construct a new family of [[N,K,D]] codes, referred to as hierarchical codes, that encode a number of logical qubits K=Ω(N/log(N) 2 ). The N th element of this code family is obtained by concatenating a constant-rate quantum LDPC code with a surface code; nearest-neighbor gates in two dimensions are sufficient to implement the corresponding syndrome-extraction circuit and achieve a threshold. Below threshold the logical failure rate vanishes superpolynomially as a function of the distance D(N). We present a bilayer architecture for implementing the syndrome-extraction circuit, and estimate the logical failure rate for this architecture. Under conservative assumptions, we find that the hierarchical code outperforms the basic encoding where all logical qubits are encoded in the surface code.

Pattison, Christopher A. [California Institute of ↗
DOE OSTI2021

Validation of the TOPAS Monte Carlo toolkit for HDR brachytherapy simulations

The goal of this work is to validate the user-friendly Geant4-based Monte Carlo toolkit TOol for PArticle Simulation (TOPAS) for brachytherapy applications. Brachytherapy simulations performed with TOPAS were systematically compared with published TG-186 reference data. The photon emission energy spectrum, the air-kerma strength, and the dose-rate constant of the model-based dose calculation algorithm (MBDCA)-WG generic Ir-192 source were extracted. For dose calculations, a track-length estimator was implemented. Additionally, the four Joint AAPM/ESTRO/ABG MBDCA-WG test cases were evaluated through histograms of the local and global dose difference volumes. A prostate, a palliative lung, and a breast case were simulated. For each case, the dose ratio map, the histogram of the global dose difference volume, and cumulative dose-volume histograms were calculated. The air-kerma strength was (9.772 ± 0.001) × 10 -8 U Bq -1 (within 0.3% of the reference value). The dose-rate constant was 1.1107 ± 0.0005 cGy h -1 U -1 (within 0.01% of the reference value). For all cases, at least 96.9% of voxels had a local dose difference within [-1%, 1%] and at least 99.9% of voxels had a global dose difference within [-0.1%, 0.1%]. The implemented track-length estimator scorer was more efficient than the default analog dose scorer by a factor of 237. For all clinical cases, at least 97.5% of voxels had a global dose difference within [-1%, 1%]. Dose-volume histograms were consistent with the reference data. TOPAS was validated for high-dose-rate brachytherapy simulations following the TG-186 recommended approach for MBDCAs. Built on top of Geant4, TOPAS provides broad access to a state-of-the-art Monte Carlo code for brachytherapy simulations.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

The kinetics of the reactions of ground-state NH with N 2 O and implications for ammonia combustion

Quantum chemistry and canonical transition state theory have been applied to derive rate constants k for reactions of 3 Σ − imidogen with nitrous oxide. 3 channels were quantified. The fastest leads to NNH+NO and may be summarized as k = 2.8 × 10 7 T 1.68 exp.(−23.1 kcal mol −1 /RT) cm 3 mole −1 s −1 . Two much slower pathways lead to N 2 + 3 HNO and spin-forbidden 1 HNO. Modeling of literature data, from an NH 3 /O 2 flame and jet-stirred reactor measurements on dilute NH 3 /N 2 O mixtures, with the new kinetic information suggests that these channels are too slow to consume N 2 O significantly, and confirms that a widely used proposed rate constant for NH+N 2 O overestimates the reactivity.

Marshall, Paul [University of North Texas, Denton,↗
DOE OSTI2021

Substitution Reactions in the Pyrolysis of Acetone Revealed through a Modeling, Experiment, Theory Paradigm

The development of high-fidelity mechanisms for chemically reactive systems is a challenging process that requires the compilation of rate descriptions for a large and somewhat ill-defined set of reactions. The present unified combination of modeling, experiment, and theory provides a paradigm for improving such mechanism development efforts. Here we combine broadband rotational spectroscopy with detailed chemical modeling based on rate constants obtained from automated ab initio transition state theory-based master equation calculations and high-level thermochemical parametrizations. Broadband rotational spectroscopy offers quantitative and isomer-specific detection by which branching ratios of polar reaction products may be obtained. Using this technique, we observe and characterize products arising from H atom substitution reactions in the flash pyrolysis of acetone (CH 3 C(O)CH 3 ) at a nominal temperature of 1800 K. The major product observed is ketene (CH 2 CO). Minor products identified include acetaldehyde (CH 3 CHO), propyne (CH 3 CCH), propene (CH 2 CHCH 3 ), and water (HDO). Literature mechanisms for the pyrolysis of acetone do not adequately describe the minor products. The inclusion of a variety of substitution reactions, with rate constants and thermochemistry obtained from automated ab initio kinetics predictions and Active Thermochemical Tables analyses, demonstrates an important role for such processes. Furthermore, the pathway to acetaldehyde is shown to be a direct result of substitution of acetone's methyl group by a free H atom, while propene formation arises from OH substitution in the enol form of acetone by a free H atom.

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

Deposition Behavior of Supersaturated Silicic Acid on the Surfaces of Chlorite and Sericite around Geological Repository - 20138

Around the repository for HLW underground, cementitious materials used for the repository may complicate the migration of radionuclides in the host rock by the alteration of groundwater to a high alkaline condition. Specially, the alkaline groundwater dissolves silicate minerals of host rocks but the plenty of silicic acid dissolved from host rocks in groundwater will become supersaturated and deposit on the surface of flow-paths after decreasing the pH value to about 8 with dispersion and dilution. The deposition of silicic acid is presumed to recover and improve the barrier performance of host rocks for the radionuclide migration because of the narrowing of flow paths and the decrease of the permeability. In this study, the deposition behavior of supersaturated silicic acid was examined by the deposition experiments and the deposition rate constants, k [m/s], were evaluated from the experimental data. As the minerals of the surface of flow-paths, this study focused on chlorite and sericite, hydrothermally-altered minerals distributed to fill fractures of granite. In the experiments, deposition behavior of silicic acid on the powder of chlorite or sericite was observed after the pH of the solution was adjusted to 8 in order to make the silicic acid supersaturated. The solution of CaCl{sub 2} and NaCl were added into a silicic acid to simulate Ca{sup 2+} from the cementitious materials of the repository, Na{sup +} of the saline groundwater in Japan. From the results, the value of k was estimated to be about the order of 10{sup -10}∼10{sup -9} m/s for both chlorite and sericite and the value was almost the same as those for amorphous silica in the previous study. Therefore, the deposition of silicic acid may be assessed by similar apparent deposition rate constants even in actual underground containing various silicate minerals. Moreover, the Damkoehler number which is the dimensionless number expressed by the ratio of characteristic time of mass transfer and chemical reaction was calculated to be 10{sup 4} or more from the value of k. This value shows that the deposition reaction of silicic acid was estimated as much more predominant than the groundwater flow under the situation simulated in this study. These results suggest that the deposition of supersaturated silicic acid on the flow paths composed with various silicate minerals will be able to retard the migration of radionuclides because of narrowing micro flow-paths. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2021

Electrochemical Study on the Kinetic Properties of Fe 2+ /Fe, Ni 2+ /Ni, Cr 2+ /Cr and Cr 3+ /Cr 2+ in Molten MgCl 2 -KCl-NaCl Salts

Here, the fundamental kinetic parameters of corrosion product ions (Fe, Ni and Cr ions) are of significant importance to understand the corrosion mechanisms of the structural materials in the molten salts. In the present study, the fundamental data of Fe 2+ , Ni 2+ , Cr 2+ and Cr 3+ (diffusion coefficient D, exchange current density i 0 , charge transfer coefficient α, limiting current density i L and standard rate constant k0) were measured at different concentrations (1.53 × 10 –6 –7.48 × 10 –4 mol cm –3 ) and temperatures (600 °C–800 °C). The values of D are independent of concentrations and follow Arrhenius law with temperature, which descend in the order of D Cr 2+ > D Fe 2+ > D Cr 3+ > D Ni 2+ with values ranging from 0.94 × 10 –5 to 3.31 × 10 –5 cm 2 s –1 . Both i 0 and k 0 depend on the temperatures, which also follow Arrhenius law.

14 SOLAR ENERGY↗
DOE OSTI2020

Formation of Phenanthrene via Recombination of Indenyl and Cyclopentadienyl Radicals: A Theoretical Study

This work presents quantum chemical G3(MP2,CC)//B2PLYPD3/6-311G(d,p) calculations of the potential energy surface for the indenyl (C 9 H 7 ) + cyclopentadienyl (C 5 H 5 ) reaction followed by unimolecular decomposition of the C 14 H 11 radicals formed as the primary products, as well as the Rice–Ramsperger–Kassel–Marcus master equation (RRKM-ME) calculations to predict temperature- and pressure-dependent reaction rate constants and product branching ratios. The reaction begins with the barrierless recombination of indenyl and cyclopentadienyl forming a C 14 H 12 molecule with a new C–C bond connecting two five-membered rings, which subsequently dissociates to C 14 H 11 radicals by H losses. The primary products of the C 9 H 7 + C 5 H 5 → C 14 H 11 + H reaction can directly decompose by another H loss to benzofulvalene, and this pathway is most favorable in terms of the entropy factor and hence is preferable at higher temperatures. Otherwise, the initial C 14 H 11 isomers can undergo significant structural rearrangements before eliminating an H atom and producing phenanthrene, anthracene, or benzoazulenes, among which the formation of phenanthrene via the “spiran” pathway is clearly preferred. The calculated barriers along the computed favorable dissociation pathways are relatively low, in the ~30–40 kcal/mol range, making the C 14 H 11 radicals unstable at temperatures above 1000–1250 K at 1 atm. The results of RRKM-ME calculations show that, under typical combustion conditions, the decomposition of the C 14 H 11 radicals predominantly leads to benzofulvalene. However, the latter can be rapidly converted to phenanthrene via H-assisted isomerization with the rate constant for the benzofulvalene + H → phenanthrene + H reaction being close to 10 –11 cm 3 molecule –1 s –1 at 1000–1500 K and 1 atm. Here, the results provide further support for the hypothesis that recombination of two π radicals containing five-membered rings can lead to a growth of PAH with the formation of two fused six-membered rings, but the reaction mechanism may not be direct and is likely to involve two consecutive H atom losses leading to a fulvalene-like product, with subsequent H-assisted isomerization of the latter to a benzenoid PAH.

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

Reductive Dynamic and Static Excited State Quenching of a Homoleptic Ruthenium Complex Bearing Aldehyde Groups

A new homoleptic Ru polypyridyl complex bearing two aldehyde groups on each bipyridine ligand, [Ru(dab) 3 ](PF 6 ) 2 , where dab is 4,4′-dicarbaldehyde-2,2′-bipyridine, was synthesized, characterized, and utilized for iodide photo-oxidation studies. In acetonitrile (CH 3 CN) solution, the complex displayed an intense metal-to-ligand charge transfer (MLCT) absorbance maximum at 475 nm (ε = 22,000 M –1 cm –1 ) and an infrared (IR) band at 1712 cm –1 assigned to the pendent aldehyde groups. Visible light excitation in air-saturated solution resulted in room temperature photoluminescence (PL) with a maximum at 675 nm, a quantum yield, ϕ PL = 0.048, and an excited state lifetime, τ ο = 440 ns, from which radiative and nonradiative relaxation rate constants were extracted, k r = 9.1 × 10 4 s –1 and knr = 1.8 × 10 6 s –1 . Pulsed visible light excitation yielded transient UV–vis and IR absorption spectra consistent with an MLCT excited state; relaxation occurred with the maintenance of two isosbestic points in the visible region, and a lifetime that agreed with that measured by time-resolved PL. Cyclic voltammetry studies in a CH 3 CN solution with 0.1 M TBAPF 6 electrolyte revealed a quasi-reversible oxidation, E°(Ru III/II ) = +1.25 V vs. Fc +/0 , and three sequential one-electron reductions at −1.10, −1.25, and −1.54 V vs. Fc +/0 . Here, an excited state reduction potential of E°(Ru *2+/+ ) = +0.89 V vs. Fc +/0 was estimated with the Rehm–Weller expression. Titration of tetrabutylammonium iodide, TBAI, into a CD 3 CN solution of [Ru(dab) 3 ](PF 6 ) 2 resulted in significant shifts in the aldehyde H atom and 3,3′-biypridyl resonances that were analyzed with a 1:1 equilibrium model, from which K eq = 460 M –1 was extracted, increasing to 5800 M –1 when the solvent was changed to acetone-d 6 . Iodide titrations resulted in a significant quenching of the [Ru(dab) 3 ] *2+ lifetime and quantum yield in both CH 3 CN and acetone solvents. In CH 3 CN, the quenching was mainly dynamic and well described by the Stern–Volmer model, from which a quenching rate constant, k q , of 4.5 × 10 10 M –1 s –1 and an equilibrium constant, K eq , of 8.3 × 10 3 M –1 were obtained. In acetone, the static quenching pathway by iodide was greatly enhanced, with a K eq of 1.2 × 10 4 M –1 and a higher k q of 9.2 × 10 10 M –1 s –1 .

Aldehydes↗
DOE OSTI2023

Highly Coupled Heterobicycle‐Fused Porphyrin Dimers: Excitonic Coupling and Charge Separation with Coordinated Fullerene, C 60

Abstract Porphyrin dimers have been widely explored and studied owing to their importance in photosynthetic systems. A vast variety of dimers linked by different groups and at different angles have been synthesized and studied; however, the means by which to synthesize rigidly fused porphyrins with direct conjugation of the chromophores remains limited. Such a class of porphyrins may possess interesting properties that unconjugated or stacked dimers may not exhibit. In this study, bisbenzimidazole‐fused porphyrin dimers and their mono‐ and bis−zinc derivatives are synthesized and characterized. As a consequence of excitonic coupling, these dimers exhibit a split Soret band irrespective of the metal ion in the porphyrin cavity. Steady‐state fluorescence and excitation spectra followed by femtosecond transient absorption spectral studies of the heterometallated dimer, (free‐base and zinc porphyrin) reveals the occurrence of efficient singlet‐singlet energy transfer (>95 % efficiency and rate constant >10 12 s −1 ) within the dyad. Further, donor‐acceptor conjugates were formed by metal‐ligand axial coordination of phenyl imidazole functionalized C 60 and were characterized by a variety of physicochemical techniques. Excited state charge separation from both singlet and triplet excited states of ZnP in the conjugates has been established. The lifetime of the final charge‐separated state was in the 30–40 μs range revealing charge stabilization. Interestingly, no charge separation in the conjugate derived from the heterometallated dimer was observed wherein excitation transfer dominated the process. The present study brings out the importance of the rigid π‐spacer connecting porphyrin dimers in governing the energy and electron transfer events when coupled with an electron acceptor.

Chemistry↗
DOE OSTI2025

Catalytic Resonance Theory: Experimental and Kinetic Interpretation of Programmable Catalysis

The complexity of programmable catalysts that change in the physical or electronic state on the time scale of a catalytic turnover results in uncertainty in the kinetic outcome of common measurements with dynamic catalyst experiments. To provide an interpretation of utility to experimentalists, a general programmable catalytic reaction exhibiting high turnover efficiency was simulated to understand the kinetic behavior of catalysts with varying frequency, binding energy, oscillation amplitude, and minimum oscillation binding energy at different temperatures and reactant gas pressures. The time-averaged catalytic rate for varying temperature in the form of an Arrhenius plot exhibited three distinct kinetic regimes, with the intermediate temperatures or applied frequencies exhibiting zero slope, indicative of barrier-free kinetic control. Transitions in experimentally measurable kinetic regimes with both temperature and applied catalyst oscillation frequency were associated with transitions in the sensitivity of reaction rate constants and degrees of rate control. Furthermore, these distinct kinetic macroscopic features specific to programmable catalysts were identified for observation in experimental dynamic kinetic measurements.

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

Response of 304L stainless steel and 6061-T651 aluminum alloy at -40 °C

The objective of this work is to extend the thermal-mechanical, elastic-plastic calibrations for 304L stainless steel [1] and and 6061-T651 aluminum alloy [2] to the regime between room temperature and -40 °C. The basis to extend the calibration consisted of new uniaxial tension tests conducted at -40 °C using the same plate material stocks, circular cylindrical specimen geometries and testing apparatus as previously, followed by attempts to fit power-law hardening functions to replicate the response observed in the specimens and then extend the yield, hardening constant, hardening exponent and rate constant functions in the calibrations to cover the new temperature regime.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Formation of late-generation atmospheric compounds inhibited by rapid deposition

Reactive organic carbon species are important fuel for atmospheric chemical reactions, including the formation of secondary organic aerosol. However, in parallel to atmospheric oxidation processes, deposition can remove compounds from the atmosphere and impact downstream environments. To understand the impact of deposition on atmospheric oxidation, we present a framework for predicting and visualizing the fate of a molecule on the basis of the physicochemical properties of compounds (Henry’s law constant, vapour pressure and reaction rate constants), which are used to estimate timescales for oxidation and deposition. Further, by implementing our deposition rates in chemical models, we show that deposition substantially suppresses atmospheric reactivity and aerosol formation by removing early-generation products and preventing the formation of large fractions (up to 90%) of downstream, late-generation compounds. Deposition is frequently missing in the laboratory experiments and detailed chemical modelling, which probably biases our understanding of atmospheric composition.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Simple Approximation for the Ideal Reference State of Gases Adsorbed on Solid-State Surfaces

Reference states are useful as models for facilitating calculations of equilibrium constants, and they may also serve as standard states that are convenient for organizing and tabulating thermodynamic data; however, standard state conventions and appropriate reference states for adsorbed species have received less attention than those for pure substances and solutes. Here, we compare seven choices of reference states for calculations of equilibrium constants and transition state theory rate constants for flat surfaces, in particular (1) an ideal 2D harmonic oscillator, (2) an ideal rigid-molecule harmonic oscillator, (3) an ideal 2D harmonic oscillator with separable surface modes, (4) a 2D ideal gas, (5) an ideal 2D hindered translator, (6) an ideal 2D hindered translator with lowest-order barriers, and (7) a simple ideal 2D hindered translator proposed in this work. The advantage of models 5–7 is that they can treat both mobile and localized adsorbates in a consistent way, whereas models 1–3 are only appropriate for localized adsorbates, and model 4 is only appropriate for a freely translating adsorbate. Furthermore, models 6 and 7 reduce the computational cost without the user having to calculate barrier heights for diffusion. An advantage of the simple ideal 2D hindered translator is that it has a physical high-temperature limit. We also propose a reference state for nonflat surfaces. Here, the user is encouraged to choose a reference state based on the appropriateness of the model and the practicality of the calculations.

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

Modeling Predictive Maintenance for NuScale’s Condensate and Feedwater System Using EMRALD

Event Modeling Risk Assessment using Linked Diagrams (EMRALD) is used to model of NuScale’s feedwater and condenser system to capture component degradation and repairing process. To model the reliability of the components, a three-stage failure rate is used contrary to a constant rate leading to component failure. Results from EMRALD will be used to compare different maintenance strategies to reduce system downtime. The outcome of the project will assist in maximizing remaining useful life of the overall system and increasing the availability and revenue of the plant.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2026

Novel insight into the kinetics of amide bond glycolysis for nylon-6 depolymerization

Chemical recycling of nylon-6 to short-chain oligomers and monomer ε-caprolactam via catalytic glycolysis is a potential solution for plastic waste remediation. Here, in this work, the kinetics of amide bond glycolysis (with ethylene glycol) in nylon-6 and the model compound N-phenethyl-3-phenylpropanamide (M1) were evaluated at 473 K in the presence of the cyclic amidine catalyst 1,5,7-triazabicyclo[4.4.0]dec-5-ene. Rates of polymer glycolysis were determined by the time-dependent shift in molecular weight distribution, whereas rates of M1 glycolysis were determined using liquid chromatography. The similarity of the first-order rate constants for glycolysis of nylon-6 and M1 at 473 K with 0.03 M amidine catalyst (5 mol% relative to amide bonds), 1.22 × 10 −5 s −1 and 2.18 × 10 −5 s −1 , respectively, confirmed the suitability of M1 as a model compound for nylon-6 glycolysis. Similar rates of glycolysis in the presence of other cyclic amidine catalysts as well as sodium methoxide revealed little influence of base strength. Glycolysis rates were unexpectedly non-linear in catalyst loading and deactivation occurred with long reaction times, presumably by non-selective decomposition of products as detected by liquid chromatography.

Depolymerization rate and rate constant↗