Evaluation of antioxidants using oxidation reaction rate constants

An evaluation method for the capacity of antioxidants to protect drugs against oxidation is presented. As a new viewpoint, to determine the priority of the competitive oxidations between the antioxidant and the protected drug, and to compare the drug-protection capacity of antioxidants, it is important to determine their oxidation rate constants using chemical kinetics instead of standard oxidation (or reduction) potentials. Sodium sulfite, sodium bisulfite and sodium pyrosulfite were used as models for the determination of oxidation reaction rate constants in aqueous solutions. In the experiments, sufficient air was continually infused into the solution to keep the concentration of dissolved oxygen constant. The residual concentrations of the antioxidants were determined by iodimetry, and the concentration of dissolved oxygen by oxygen electrode. The data were fitted by linear regressions to obtain the reaction rate constants. It was found that the degradation of sodium sulfite, sodium bisulfite or sodium pyrosulfite obeyed pseudo zero-order kinetics in the buffer solutions. Because of the ionization equilibrium, these three antioxidants have the same ion form in solutions at a definite pH value and therefore their apparent rate constants were essentially the same. The average apparent rate constants of the three antioxidants at 25°C are (1.34 ± 0.03) × 10⁻³ at pH 6.8, (1.20 ± 0.02) × 10⁻³ at pH 4.0 and (6.58 ± 0.02) × 10⁻³ mol·L⁻¹·h⁻¹ at pH 9.2, respectively. Translated from Acta Chimica Sinica, 2006, 64(6): 496–500 (in Chinese)     

氟代乙炔氢迁移反应微正则系综速率常数的计算

以氟代乙炔的氢迁移反应为例, 计算了隧道效应、曲率因子和转动选态对微正则系综单分子反应速率常数的影响。结果表明, 在低能区域, 隧道效应和曲率因子对反应速率有较大影响, 均使反应速率增加, 而在高能区域作用不明显。与此同时还研究了不同对称态下的反应速率常数, 结果表明, 本例中的两对称态(A'态和A"态)的反应速率常数在能量较低时差异较大, 随能量升高, 两者趋于相等。     

Comparison between hydrogen-transfer and fluoro-transfer reactions on the microcanonical unimolecular rate constants

The microcanonical rate constants for the hydrogen-transfer process of HCCF (reaction 7) and the fluoro-transfer process of FCCF (reaction 8) are carried out with tunneling correction and curvature correction. The results show that the tunneling effects and curvature effects on the rate constant of reaction 7 is quite different from that of reaction 8. The rate constants for different rotational states are also studied for these reactions.     

Theoretical calculations of the unimolecular canonical rate constants

Microcanonical rate constants k(E) and canonical rate constants k(T) for unimolecular reactions have been obtained through the calculations of cumulative reaction probabilities N(E) with the unsymmetrical Eckart potential tunneling correction. By way of example, the reactions HCN→CNH (I) and FNC→NCF (II) have been employed. For reaction (I), the calculated rate constants are in agreement with the experimental data; for reaction (II), the results are in accordance with the rate constants kCVT/MEPSAG(T) calculated by the common program POLYRATE.     

An analytic hindered rotor model for calculating microcanonical variational unimolecular rate constants from reaction path properties

A model is proposed for performing microcanonical variational transition state theory calculations which incorporates ideas from vibrator and flexible variational transition state models. Vibrational frequencies, moments of inertia, and potential energy for the variational transition state are found by reaction path following as for the vibrator model. However, the transitional modes are treated as hindered rotors using an analytic potential and an analytic density of states, which are fit to barriers for hindered rotation determined from reaction path following. The model proposed here differs from the flexible transition state model in that the density of states for the transitional modes is analytic and transitional modes and external rotational angular momenta are uncoupled. For the H + CH₃ ? CH₄ system, rate constants calculated with this new model are only 6–23% smaller on average from those of the flexible transition state model for values of total angular momentum which correspond to average rotational temperatures of 0–2000 K. Harmonic frequencies calculated for the transitional modes from the hindered rotor Hamiltonian are in good agreement with the exact values found by a reaction path analysis. © 1994 John Wiley & Sons, Inc.     

A metropolis Monte Carlo method for computing microcanonical statistical rate constants

A new method of computing microcanonical statistical rate constants is presented. The method utilizes the Metropolis Monte Carlo algorithm in a manner which circumvents some of the numerical inefficiency associated with other Metropolis and “shot-gun” Monte Carlo based schemes. It is therefore expected to be useful in studies of many degree of freedom systems where numerical efficiency is crucial. Optimization of the method efficiency with respect to its adjustable parameters is examined in detail, both theoretically and in a numerical study of the T-shaped Ar₃ “inversion” process. The energy dependence of the T-shaped Ar₃ inversion rate is studied in a sample application of the method. An application to full three dimensional Ar₃ will be presented in a future study.     

On reaction rate constants and rate kernels

A simple model for a reversible isomerization reaction A ? B is used to show how rate constants and their time-dependent generalizations can be related to microscopic reactive steps and nonreactive internal dynamics (e.g., vibrational relaxation).     

Implementation of a microcanonical variational transition state theory for direct dynamics calculations of rate constants

Calculation of microcanonical rate constants has been an important field in chemical dy-namic studies for many years because it can be used not only to give good prediction of rate con-stants in microcanonical assembly, but also to calculate rate constants with certain conserved quantum numbers such as the total angular momentum, and in turn, can be easily converted into thermal rate constants[1—3]. The widely used method for calculating microcanonical rate constants of unimolecular reac-tions…     

Semiclassical calculation of nonadiabatic thermal rate constants: application to condensed phase reactions

The nonadiabatic transition state theory proposed recently by Zhao et al. [J. Chem. Phys. 121, 8854 (2004)] is extended to calculate rate constants of complex systems by using the Monte Carlo and umbrella sampling methods. Surface hopping molecular dynamics technique is incorporated to take into account the dynamic recrossing effect. A nontrivial benchmark model of the nonadiabatic reaction in the condensed phase is used for the numerical test. It is found that our semiclassical results agree well with those produced by the rigorous quantum mechanical method. Comparing with available analytical approaches, we find that the simple statistical theory proposed by Straub and Berne [J. Chem. Phys. 87, 6111 (1987)] is applicable for a wide friction region although their formula is obtained using Landau-Zener [Phys. Z. Sowjetunion 2, 46 (1932); Proc. R. Soc. London, Ser. A 137, 696 (1932)] nonadiabatic transition probability along a one-dimensional diffusive coordinate. We also investigate how the nuclear tunneling events affect the dependence of the rate constant on the friction.     

Model system-bath Hamiltonian and nonadiabatic rate constants for proton-coupled electron transfer at electrode-solution interfaces

An extension of the Anderson-Newns-Schmickler model for electrochemical proton-coupled electron transfer (PCET) is presented. This model describes reactions in which electron transfer between a solute complex in solution and an electrode is coupled to proton transfer within the solute complex. The model Hamiltonian is derived in a basis of electron-proton vibronic states defined within a double adiabatic approximation for the electrons, transferring proton, and bath modes. The interaction term responsible for electronic transitions between the solute complex and the electrode depends on the proton donor-acceptor vibrational mode within the solute complex. This model Hamiltonian is used to derive the anodic and cathodic rate constants for nonadiabatic electrochemical PCET. The derivation is based on the master equations for the reduced density matrix of the electron-proton subsystem, which includes the electrons of the solute complex and the electrode, as well as the transferring proton. The rate constant expressions differ from analogous expressions for electrochemical electron transfer because of the summation over electron-proton vibronic states and the dependence of the couplings on the proton donor-acceptor vibrational motion. These differences lead to additional contributions to the total reorganization energy, an additional exponential temperature-dependent prefactor, and a temperature-dependent term in the effective activation energy that has different signs for the anodic and cathodic processes. This model can be generalized to describe both nonadiabatic and adiabatic electrochemical PCET reactions and provides the framework for the inclusion of additional effects, such as the breaking and forming of other chemical bonds.     

Determination of the elementary reaction rate constants on the basis of the summary values of the rate constants

Conclusions The rate constants of reactions of atomic oxygen with methane, ethane, ethylene, and propylene along different pathways were determined on the basis of the ratios between the concentrations of the products formed and the summary rate constants measured earlier.Translated from Izvestiya Akademii Nauk SSSR, Seriya Khimicheskaya, No. 12, pp. 2693–2699, December, 1971.     

Relation between lyoluminescence and radical reaction rate constants

Light emission during the dissolving of irradiated sugars (lyoluminescence, LL) allows the estimation of absorbed dose. The use of 1-mannose as LL substance and the correlation between the concentration of paramagnetic centres and LL yield in the presence of [Fe(CN)₆]^4? and CNS^? anions demonstrated the possibility to measure relative rate constants of the reaction of mannose peroxy-radicals with different substances soluble in water.     

The three-dimensional nonadiabatic dynamics calculation of DH(2)(+) and HD(2)(+) systems by using the trajectory surface hopping method based on the Zhu-Nakamura theory

A theoretical investigation on the nonadiabatic processes of the full three-dimensional D(+)+H(2) and H(+)+D(2) reaction systems has been performed by using trajectory surface hopping (TSH) method based on the Zhu-Nakamura (ZN) theory. This ZN-TSH method refers to not only classically allowed hops but also classically forbidden hops. The potential energy surface constructed by Kamisaka et al. is employed in the calculation. A new iterative method is proposed to yield the two-dimensional seam surface from the topography of the adiabatic potential surfaces, in which the inconvenience of directly solving the first-order partial differential equation is avoided. The cross sections of these two systems are calculated for three competing channels of the reactive charge transfer, the nonreactive charge transfer, and the reactive noncharge transfer, for ground rovibrational state of H(2) or D(2). Also, this study provides reaction probabilities of these three processes for the total angular momentum J=0 and ground initial vibrational state of H(2) or D(2). The calculated results from ZN-TSH method are in good agreement with the exact quantum calculations and the experimental measurements.     

Evaluation of the rate constants in chemical reactions

This article is concerned with the application of a new method to recover the rate constants in chemical reactions. The method is based on treating the unknown parameters as time dependent. With appropriate experimental data the unknown rate constants are guided from an arbitrary initial condition to their true value at a final time. An explicit equation describing the time evolution of the parameters is obtained by minimizing the error along the trajectory. The method leads to an iterative algorithm which is described in detail. Numerical results with the method indicate that accurate estimates of the rate constants can be obtained directly from experimental data. © 1998 John Wiley & Sons, Inc. Int J Chem Kinet 30: 151–159, 1998.     

Semiclassical theory of electronically nonadiabatic chemical dynamics: incorporation of the Zhu-Nakamura theory into the frozen Gaussian propagation method

The title theory is developed by combining the Herman-Kluk semiclassical theory for adiabatic propagation on single potential-energy surface and the semiclassical Zhu-Nakamura theory for nonadiabatic transition. The formulation with use of natural mathematical principles leads to a quite simple expression for the propagator based on classical trajectories and simple formulas are derived for overall adiabatic and nonadiabatic processes. The theory is applied to electronically nonadiabatic photodissociation processes: a one-dimensional problem of H2+ in a cw (continuous wave) laser field and a two-dimensional model problem of H2O in a cw laser field. The theory is found to work well for the propagation duration of several molecular vibrational periods and wide energy range. Although the formulation is made for the case of laser induced nonadiabatic processes, it is straightforwardly applicable to ordinary electronically nonadiabatic chemical dynamics.     

The hydroxyl radical reaction rate constants and atmospheric reaction products of three siloxanes

The relative rate technique has been used to measure the hydroxyl radical (OH) reaction rate constant of hexamethyldisiloxane (MM, (CH₃)₃Si-O-Si(CH₃)₃), octamethyltrisiloxane (MDM, (CH₃)₃Si-O-Si(CH₃)₂-O-Si(CH₃)₃), and decamethyltetrasiloxane (MD₂M, (CH₃)₃Si-O-Si(CH₃)₂-O-Si(CH₃)₂-O-Si(Ch₃)₃). Hexamethyldisiloxane, octamethyltrisiloxane, and decamethyltetrasiloxane react with OH with bimolecular rate constants of 1.32 ± 0.05 × 10⁻¹² cm³molecule⁻¹s⁻¹, 1.83 ± 0.09 × 10⁻¹² cm³ molecule⁻¹s⁻¹, and 2.66 ± 0.13 × 10⁻¹² cm³molecule⁻¹s⁻¹, respectively. Investigation of the OH + siloxane reaction products yielded trimethylsilanol, pentamethyldisiloxanol, heptamethyltetrasiloxanol, hexamethylcyclotrisiloxane, octamethylcyclotetrasiloxane, and other compounds. Several of these products have not been reported before because these siloxanes and the proposed reaction mechanisms yielding these products are complicated. Some unusual cyclic siloxane products were observed and their formation pathways are discussed in light of current understanding of siloxane atmospheric chemistry. © 1997 John Wiley & Sons, Inc. Int J Chem Kinet 29: 445–451, 1997.     

Measurement of the reaction rate constants of oxygen atoms with chlorine and iodomethane using a resonance fluorescence technique

The following reaction rate constants of oxygen atoms with iodomethane and chlorine were measured using resonance fluorescence under jet conditions at 298 K: k ₁ = (2.4 ± 0.5) × 10–15 and k ₂ = (6.9 ± 0.2) × 10⁻¹⁴ cm³/s, respectively.     

Differential method for the determination of the order of reaction and the rate constants

A method is proposed whereby the orders and rate constants for processes obeying the rate law ?dA/dt = kAⁿ may be determined. The method is illustrated in two ways. First, simulated data for processes of various orders are treated, and the treatment is shown to be capable of reproducing orders and rate constants to a high degree of accuracy. The factors affecting the accuracy with which n and k can be determined are considered. These are inaccuracy in the determination of concentration values, irregularity of the time intevals between concentration determinations, and the length of those time intervals. It is shown that if concentrations are determined at times that are close together, the effect of the other two factors is small, but if the time intervals are made longer, the errors due to the other two factors affect the calculated values of n and k much more seriously. Second, the method was applied to two homogeneous reactions, of which one was first-order and one was second order, and three heterogeneous reactions, of which one was found by the original workers to be first order, one to be zero order, and one to vary between zero and first order, depending on the initial pressure. The present method gives results in agreement with these conclusions and reproduces the rate constants to within ±5% in all cases.     

Reaction rate constants and mechanistic detail of the Ni+ + butanone reaction

The unimolecular decomposition kinetics of the jet-cooled Ni(+)-butanone cluster ion has been monitored over a range of internal energies (16000-18800 cm?1). First-order rate constants are acquired for the precursor ion dissociation into three product channels. The temporal growth of each fragment ion is selectively monitored in a custom instrument and yields similar valued rate constants at a common ion internal energy. The decomposition reaction is proposed to proceed along two parallel reaction coordinates. Each dissociative pathway is rate-limited by the initial Ni(+) oxidative addition into either the C-CH? or C-C?H? σ-bond in the butanone molecule. Ratios of integrated product ion intensities as well as the measured rate constants are used to determine values for each σ-bond activation rate constant. The lowest energy measurement presented in this study occurs when the binary complex ion possesses an internal energy of 16000 cm?1. Under this condition, the Ni(+) assisted decomposition of the butanone molecule is rate limited by k(act)(C-C?H?) = (0.92 ± 0.08) × 10? s?1 and k(act)(C-CH?) = (0.37 ± 0.03) × 10? s?1. The relative magnitudes of the two rate constants reflect the greater probability for reaction to occur along the C-C?H? σ-bond insertion pathway, consistent with thermodynamic arguments. DFT calculations at the B3LYP/6-311++G(d,p) level of theory suggest the most likely geometries and relative energies of the reactants, intermediates, and products.