Search

The limitation of KoopmanEs theory like itEs application to some cases of a small donor - acceptor distance is discussed. Investigations shows that the two orbitals included must be properly selected in evaluating the electron-transfer (ET) matrix element.

The isodesmic reaction method is applied to calculate the potential energy surface (PES) along the reaction coordinates and the rate constants of the barrierless reactions for unimolecular dissociation reactions of alkanes to form two alkyl radicals and their reverse recombination reactions. The reaction class is divided into 10 subclasses depending upon the type of carbon atoms in the reaction centers. A correction scheme based on isodesmic reaction theory is proposed to correct the PESs at UB3LYP/6-31+G(d,p) level. To validate the accuracy of this scheme, a comparison of the PESs at B3LYP level and the corrected PESs with the PESs at CASPT2/aug-cc-pVTZ level is performed for 13 representative reactions, and it is found that the deviations of the PESs at B3LYP level are up to 35.18 kcal/mol and are reduced to within 2 kcal/mol after correction, indicating that the PESs for barrierless reactions in a subclass can be calculated meaningfully accurately at a low level of ab initio method using our correction scheme. High-pressure limit rate constants and pressure dependent rate constants of these reactions are calculated based on their corrected PESs and the results show the pressure dependence of the rate constants cannot be ignored, especially at high temperatures. Furthermore, the impact of molecular size on the pressure-dependent rate constants of decomposition reactions of alkanes and their reverse reactions has been studied. The present work provides an effective method to generate meaningfully accurate PESs for large molecular system.

Intramolecular H-migration reaction of hydroperoxyalkylperoxy radicals ( • O 2 QOOH) is one of the most important reaction families in the low-temperature oxidation of hydrocarbon fuels. This reaction family is first divided into classes depending upon H atom transfer from -OOH bonded carbon or non-OOH bonded carbon, and then the two classes are further divided depending upon the ring size of the transition states and the types of the carbons from which the H atom is transferred. High pressure limit rate rules and pressure-dependent rate rules for each class are derived from the rate constants of a representative set of reactions within each class using electronic structure calculations performed at the CBS-QB3 level of theory. For the intramolecular H-migration reactions of • O 2 QOOH radicals for abstraction from an -OOH substituted carbon atom (-OOH bonded case), the result shows that it is acceptable to derive the rate rules by taking the average of the rate constants from a representative set of reactions with different sizes of the substitutes. For the abstraction from a non-OOH substituted carbon atom (non-OOH bonded case), rate rules for each class are also derived and it is shown that the difference between the rate constants calculated by CBS-QB3 method and rate constants estimated from the rate rules may be large; therefore, to get more reliable results for the low-temperature combustion modeling of alkanes, it is better to assign each reaction its CBS-QB3 calculated rate constants, instead of assigning the same values for the same reaction class according to rate rules. The intramolecular H-migration reactions of • O 2 QOOH radicals (a thermally equilibrated system) are pressure-dependent, and the pressure-dependent rate constants of these reactions are calculated by using the Rice-Ramsberger-Kassel-Marcus/master-equation theory at pressures varying from 0.01 to 100 atm. The impact of molecular size on the pressure-dependent rate constants of the intramolecular H-migration reactions of • O 2 QOOH radicals has been studied, and it is shown that the pressure dependence of the rate constants of intramolecular H-migration reactions of • O 2 QOOH radicals decreases with the molecular size at low temperatures and the impact of molecular size on the pressure-dependent rate constants decreases as temperature increases. It is shown that it is acceptable to derive the pressure-dependent rate rules by taking the average of the rate constants from a representative set of reactions with different sizes of the substitutes. The barrier heights follow the Evans-Polanyi relationship for each type of intramolecular hydrogen-migration reaction studied. All calculated rate constants are fitted by a nonlinear least-squares method to the form of a modified Arrhenius rate expression at pressures varying from 0.01 to 100 atm and at the high-pressure limit. Furthermore, thermodynamic parameters for all species involved in these reactions are calculated by the composite CBS-QB3 method and are given in NASA format.

Hydrogen abstraction from toluene by OH, H, O, CH3, and HO2 radicals are important reactions in oxidation process of toluene. Geometries and corresponding harmonic frequencies of the reactants, transition states as well as products involved in these reactions are determined at the B3LYP/6-31G(2df,p) level. To achieve highly accurate thermochemical data for these stationary points on the potential energy surfaces, the Gaussian-4(G4) composite method was employed. Torsional motions are treated either as free rotors or hindered rotors in calculating partion functions to determine thermodynamic properties. The obtained standard enthalpies of formation for reactants and some prodcuts are shown to be in excellent agreement with experimental data with the largest error of 0.5 kcal mol(-1). The conventional transition state theory (TST) with tunneling effects was adopted to determine rate constants of these hydrogen abstraction reactions based on results from quantum chemistry calculations. To faciliate its application in kinetic modeling, the obtained rate constants are given in Arrhenius expression: k(T) = AT(n) exp(-EaR/T). The obtained reaction rate constants also agree reasonably well with available expermiental data and previous theoretical values. Branching ratios of these reactions have been determined. The present reaction rates for these reactions have been used in a toluene combustion mechanism, and their effects on some combustion properties are demonstrated.

Within the framework of constrained density functional theory (CDFT), the diabatic or charge localized states of electron transfer (ET) have been constructed. Based on the diabatic states, inner reorganization energy λin has been directly calculated. For solvent reorganization energy λs, a novel and reasonable nonequilibrium solvation model is established by introducing a constrained equilibrium manipulation, and a new expression of λs has been formulated. It is found that λs is actually the cost of maintaining the residual polarization, which equilibrates with the extra electric field. On the basis of diabatic states constructed by CDFT, a numerical algorithm using the new formulations with the dielectric polarizable continuum model (D-PCM) has been implemented. As typical test cases, self-exchange ET reactions between tetracyanoethylene (TCNE) and tetrathiafulvalene (TTF) and their corresponding ionic radicals in acetonitrile are investigated. The calculated reorganization energies λ are 7293 cm(-1) for TCNE/TCNE(-) and 5939 cm(-1) for TTF/TTF(+) reactions, agreeing well with available experimental results of 7250 cm(-1) and 5810 cm(-1), respectively.

We present a further interpretation of reaction class transition state theory (RC-TST) proposed by Truong et al. for the accurate calculation of rate coefficients for reactions in a class. It is found that the RC-TST can be interpreted through the isodesmic reaction method, which is usually used to calculate reaction enthalpy or enthalpy of formation for a species, and the theory can also be used for the calculation of the reaction barriers and reaction enthalpies for reactions in a class. A correction scheme based on this theory is proposed for the calculation of the reaction barriers and reaction enthalpies for reactions in a class. To validate the scheme, 16 combinations of various ab initio levels with various basis sets are used as the approximate methods and CCSD(T)/CBS method is used as the benchmarking method in this study to calculate the reaction energies and energy barriers for a representative set of five reactions from the reaction class: R(c)CH(R(b))CR(a)CH2 + OH(•) → R(c)C(•)(R(b))CR(a)CH2 + H2O (R(a), R(b), and R(c) in the reaction formula represent the alkyl or hydrogen). Then the results of the approximate methods are corrected by the theory. The maximum values of the average deviations of the energy barrier and the reaction enthalpy are 99.97 kJ/mol and 70.35 kJ/mol, respectively, before correction and are reduced to 4.02 kJ/mol and 8.19 kJ/mol, respectively, after correction, indicating that after correction the results are not sensitive to the level of the ab initio method and the size of the basis set, as they are in the case before correction. Therefore, reaction energies and energy barriers for reactions in a class can be calculated accurately at a relatively low level of ab initio method using our scheme. It is also shown that the rate coefficients for the five representative reactions calculated at the BHandHLYP/6-31G(d,p) level of theory via our scheme are very close to the values calculated at CCSD(T)/CBS level. Finally, reaction barriers and reaction enthalpies and rate coefficients of all the target reactions calculated at the BHandHLYP/6-31G(d,p) level of theory via the same scheme are provided.

Aromatic hydrocarbon fuels, such as toluene, are important components in real jet fuels. In this work, reactive molecular dynamics (MD) simulations employing the ReaxFF reactive force field have been performed to study the high-temperature oxidation mechanisms of toluene at different temperatures and densities with equivalence ratios ranging from 0.5 to 2.0. From the ReaxFF MD simulations, we have found that the initiation consumption of toluene is mainly through three ways, (1) the hydrogen abstraction reactions by oxygen molecules or other small radicals to form the benzyl radical, (2) the cleavage of the C-H bond to form benzyl and hydrogen radicals, and (3) the cleavage of the C-C bond to form phenyl and methyl radicals. These basic reaction mechanisms are in good agreement with available chemical kinetic models. The temperatures and densities have composite effects on toluene oxidation; concerning the effect of the equivalence ratio, the oxidation reaction rate is found to decrease with the increasing of equivalence ratio. The analysis of the initiation reaction of toluene shows that the hydrogen abstraction reaction dominates the initial reaction stage at low equivalence ratio (0.5-1.0), while the contribution from the pyrolysis reaction increases significantly as the equivalence ratio increases to 2.0. The apparent activation energies, E(a), for combustion of toluene extracted from ReaxFF MD simulations are consistent with experimental results.

Thermal cracking of n-decane and n-decane in the presence of several fuel additives are studied in order to improve the rate of thermal cracking by using reactive molecular dynamics (MD) simulations employing the ReaxFF reactive force field. From MD simulations, we find the initiation mechanisms of pyrolysis of n-decane are mainly through two pathways: (1) the cleavage of a C-C bond to form smaller hydrocarbon radicals, and (2) the dehydrogenation reaction to form an H radical and the corresponding decyl radical. Another pathway is the H-abstraction reactions by small radicals including H, CH(3), and C(2)H(5). The basic reaction mechanisms are in good agreement with existing chemical kinetic models of thermal decomposition of n-decane. Quantum mechanical calculations of reaction enthalpies demonstrate that the H-abstraction channel is easier compared with the direct C-C or C-H bond-breaking in n-decane. The thermal cracking of n-decane with several additives is further investigated. ReaxFF MD simulations lead to reasonable Arrhenius parameters compared with experimental results based on first-order kinetic analysis. The different chemical structures of the fuel additives greatly affect the apparent activation energy and pre-exponential factors. The presence of diethyl ether (DEE), methyl tert-butyl ether (MTBE), 1-nitropropane (NP), 3,6,9-triethyl-3,6,9-trimethyl-1,2,4,5,7,8-hexaoxonane (TEMPO), triethylamine (TEA), and diacetonediperodixe (DADP) exhibit remarkable promoting effect on the thermal cracking rates, compared with that of pure n-decane, in the following order: NP > TEMPO > DADP > DEE (∼MTBE) > TEA, which coincides with experimental results. These results demonstrate that reactive MD simulations can be used to screen for fuel additives and provide useful information for more comprehensive chemical kinetic model studies at the molecular level., (© 2012 American Chemical Society)

To provide insight on the reaction mechanism of the methyperoxy (CH(3)O(2)•) self-reaction, stationary points on both the spin-singlet and the spin-triplet potential energy surfaces of 2(CH(3)O(2)•) have been searched at the B3LYP/6-311++G(2df,2p) level. The relative energies, enthalpies, and free energies of these stationary points are calculated using CCSD(T)/cc-pVTZ. Our theoretical results indicate that reactions on a spin-triplet potential energy surface are kinetically unfavorable due to high free energy barriers, while they are more complicated on the spin-singlet surface. CH(3)OOCH(3) + O(2)(1) can be produced directly from 2(CH(3)O(2)•), while in other channels, three spin-singlet chain-structure intermediates are first formed and subsequently dissociated to produce different products. Besides the dominant channels producing 2CH(3)O• + O(2) and CH(3)OH + CH(2)O + O(2) as determined before, the channels leading to CH(3)OOOH + CH(2)O and CH(3)O• + CH(2)O + HO(2)• are also energetically favorable in the self-reaction of CH(3)O(2)• especially at low temperature according to our results.

In our recent work, a new form of the electrostatic solvation energy for the nonequilibrium polarization has been derived by introducing the method of constrained equilibrium state in the framework of continuous medium theory. Up until now, the idea of the constrained equilibrium state method has not been introduced into the explicit solvent model by others; therefore this nonequilibrium energy form was further equivalently extended to the explicit solvent model in this work based on the discrete representation of the solvent permanent charges and induced dipoles. Making use of this expression in explicit solvent model, we modified the nonequilibrium module in the averaged solvent electrostatic potential/molecular dynamics program to implement numerical calculations. Subsequently, the new codes were applied to study the solvatochromic shifts of the n → π* absorption spectra for acetone and trans-formic acid in aqueous solution. The calculation results show a good agreement with the experimental observations. When our results of spectral shift are compared with those achieved directly from the continuum model, it can be seen that both the explicit solvent model and continuum model derived based on the constrained equilibrium approach can give reasonable predictions. The hydrogen bond effect was also discussed and deemed to be a dominant contribution to the spectral shift by calculating the n → π* absorption spectra of acetone-water complexes.

Propenols have been found to be common intermediates in the hydrocarbon combustion and they are present in substantial concentrations in a wide range of flames. However, the kinetics properties of these species in combustion flames have not received much attention. In this work, the mechanism and kinetics of the OH hydrogen abstraction from propenols are investigated. Three stable conformations of propenols, (E)-1-propenol, (Z)-1-propenol, and syn-propen-2-ol, are taken into consideration. The potential energy profiles for the three reaction systems have been first investigated by the CCSD(T) method. The geometric parameters and relative energies of the reactants, reactant complexes, transition states, product complexes, and products have been investigated theoretically. The rate constants are calculated in the temperature range of 200-3000 K by the Variflex code based on the weak collision master equation/microcanonical variational RRKM theory. For all considered reactions, our results support a stepwise mechanism involving the formation of a reactant complex in the entrance channel and a product complex in the exit channel. In the reaction of OH with (E)-1-propenol, the hydrogen abstractions from the -CH(3) and -OH sites are dominant and competitive with each other in the temperature range from 500 to 2000 K. Above 2000 K, the hydrogen abstraction from the -CH group bonded to O atom becomes dominant with a relative yield of 51.1% at 3000 K. In the reaction of OH with (Z)-1-propenol, the hydrogen abstractions from -CH(3), -CH bonded to O atom, and -OH are preferable in the temperature range from 500 to 1800 K, with the first two channels being competitive with each other. Above 1800 K, the hydrogen abstraction reaction from the CH group bonded to the CH(3) group becomes dominant with the branching ratio of 90.3% at 3000 K. In the reaction of OH with syn-propen-2-ol, the abstractions from the -CH(3) and -OH sites are competitive with each other when the temperature is higher than 500 K, and they become dominant above 800 K with the relative yields of 70.5% and 29.5% at 3000 K, respectively. The predicted total rate constants at the pressure of 1 atm fitted by modified three-parameter Arrhenius expressions in two different temperature ranges are also provided.

Recently, enols have been found to be the common intermediates in hydrocarbon combustion flames (Taatjes et al. Science 2005, 308, 1887), but the knowledge of kinetic properties for such species in combustion flames is rare. Therefore in this work, particular attention is paid to the formation of enols in combustion flames. Starting with HO and propene (CH(3)CH=CH(2)), the reaction mechanism involving eight product channels has been investigated systematically. It is revealed that the electrophilic addition of OH to the double bond of CH(3)CH=CH(2) is unselective and the chemically activated adducts, CH(3)CHOH=CH(2) and CH(3)CH=CH(2)OH, may undergo dissociation in competition with H-abstractions. The kinetics and product branching ratios of the HO and propene reaction have been evaluated in the temperature range of 200-3000 K by Variflex code, based on the weak collision master equation/microcanonical variational RRKM theory. Available experimental kinetic data can be quantitatively reproduced by this study, with a minor adjustment (1.0 kcal/mol) of the OH central addition barrier. From the theoretical calculations with multiple reflection correction included, the total rate constant is fitted to k(t) = 6.07 x 10(-5)T(-2.54) exp(108/T) cm(3) x molecule(-1) x s(-1) in the range of 200-800 K and k(t) = 7.11 x 10(-23)T(3.38) exp(-1097/T) cm(3) x molecule(-1) x s(-1) in the range of 800-3000 K, which are in close agreement with experimental data. The branching ratios of enol channels are consistent with the observation in low-pressure flames and hence the reaction mechanisms presented here provide valuable descriptions of enol formations in hydrocarbon combustion chemistry.

Density functional theory was applied to the investigation of photoinduced electron transfer (ET) and the absorption spectrum for a bis-naphthalimide spermine conjugate. The multichannel feature of ET excitation in this system was focused on because four groups may act as electron donors and acceptors. The segment in this conjugate, N-(N-methylpropyl)-1,8-naphthalimide, which contains one donor and acceptor pair, was studied at first. Through theoretical calculation, the absorption band at 340 nm was assigned to the pi-->pi* transition. For the whole system involving four chromophores, this work suggested three types of ET. From the theoretical investigation, the naphthalimide radical anion turned out to be formed via intramolecular ET between the two terminal naphthalimide groups, rather than via the electron transfer between the dialkylamine moiety and the naphthalimide one. Furthermore, the electronic coupling matrix elements according to the generalized Mulliken-Hush theory were estimated and the detailed analyses showed that the strongest absorption was due to the local excitation of the naphthalimide chromophore.

The excited states of a donor-bridge-acceptor (DBA) model system have been investigated using time-dependent density-functional theory (TD-DFT) in vacuo and in solution. It is found that the MPW1PW91 functional always gives higher excitation energies than those with a B3LYP functional. Results from both TD-B3LYP and TD-MPW1PW91 are found consistent with the experimental observations. The two most intense absorptions of the DBA system, one resulting from the local excitation of catechol moiety and the other from that of dicyanoethylene, possess the pipi* transition feature. It seems that the solvent polarity does not remarkably influence the positions of absorption peaks. The spectroscopic properties of isolated donor, acceptor, and bridge and the donor-bridge compound have been investigated at the TD-B3LYP/6-31+G* and TD-MPW1PW91/6-31+G* levels. Results indicate that the donor and the acceptor are weakly coupled with the bridge. Therefore, it is more likely that the electron transfer takes place through a superexchange mechanism. In addition, we calculate the electronic coupling matrix elements according to the generalized Mulliken-Hush theory, and the detailed analyses also predict that the strong absorptions are due to the local excitation of the DBA system.