Abstract

Entropy data are reported using different calculation methods for internal rotors on n-heptane, 2-methylhexane, and 2,3-dimethylpentane and on the different radical sites of each species corresponding to the loss of a hydrogen atom for temperatures between 298 and 1500 K. Structures, moments of inertia, vibration frequencies, and internal rotor potentials are calculated at the B3LYP/6-31G(d,p) level of theory. Comparisons with experimental literature data suggest limitations inuse of the rigid-rotor harmonic-oscillator (HO) approximation and advantages to the use of internal rotation contributions for entropy relative to torsion frequencies. The comparisons suggest the need to include contributions from all internal rotors where the barriers are at or below those of the above molecules. Calculation of entropy from the use of internal rotor contributions provides acceptable approximations to available literature values. Entropy values for radicals corresponding to carbon sites on these hydrocarbons are presented.

1. Introduction

Over the past decades advancements made in quantum chemistry, including the evolution of density functional theory (DFT) and compound methods, coupled with increasing computer processing power, have allowed for substantial advancements in estimating thermochemical and other chemical properties and systems [1].

Thermochemical properties including enthalpies, entropies, and heat capacities for stable molecules and for radicals corresponding to the loss of H atom through abstraction reactions are important and necessary for building and analyzing detailed chemical kinetic models. Several research studies have focused on linear and branched alkanes, including n-heptane, [26] because it has been considered as an important surrogate for major components in diesel and spark ignition transportation fuels. The parent hydrocarbons might be stable alkanes that have been studied experimentally, but oxidation under both atmospheric and combustion conditions creates hundreds of intermediate species that are fragmented, partially oxygenated, and/or unsaturated. The stability, thermochemical properties, and chemical kinetics of these intermediates is needed for the model used for optimization.

Enthalpy, entropy, and heat capacity are all important because of the influence they have on reactions’ paths. The TΔS entropy term in the Gibbs Energy influences the equilibrium constant and calculation of reverse reaction rate constants with a large lever because of the temperature multiplier. Entropy also has a large contribution to preexponential factors. The advances in high level computational methods have resulted in accurate determination of thermochemical properties for linear and branched alkanes [79].

When determining the standard enthalpy of formation for a species, for example, it is common to implement different theoretical chemical equations that incorporate mass balance, with equivalent hybridization, and bond type on each side of a work equation. This allows for significant error cancellation via the use of a due to similar environments, which result in the cancelation of error resulting, on both sides of the reaction equation. Hierarchies of the reaction classes have been created allowing for a better idea of the error cancelling provided [10]. Employing these techniques also allows lower level and less computationally demanding methods to be used without sacrificing the accuracy.

For entropy and heat capacity calculations, the simple rigid-rotor harmonic-oscillator (HO) is utilized, almost exclusively, to describe the 3N-6 vibrations for nonlinear species and the corresponding contributions to entropy and heat capacities. It is well known that there are accuracy issues in determining the lower frequency torsions corresponding to internal rotations using this approximation. This can be corrected by replacing these frequencies with methods to treat the internal rotations as hindered rotors to increase accuracy for thermodynamic properties and reaction kinetics [1217].

The initial research papers of Pitzer and Gwinn [1820] addressed contributions from symmetrical and unsymmetrical rotating groups on a rigid frame. Later work by Kilpatrick and Pitzer [21] expanded the method to include balanced and unbalanced linked rotating groups. These studies still serve as a basis for current research and development for new methods for treating internal rotations.

A number of studies utilizing different techniques for handling coupled and uncoupled internal rotor contributions have been reported [12, 13, 16, 2226]. Although more advanced methods addressing coupled internal rotators are available, basic treatment of internal rotations as uncoupled rotations provides improved accuracy over the HO approximation alone [22, 24, 27]. This has been shown through comparisons between experimental and computational entropies and heat capacities of n-alkanes [22, 2628].

The goal of the present study is to calculate entropy values for n-heptane, 2-methylhexane, and 2,3-dimethylpentane, compare them to available literature data, and then use the method(s) which appear to be optimal, to calculate entropy data for their carbon-centered radicals shown in Figure 1. String nomenclatures of the radical species are written with the letter j representing a radical site on the preceding carbon atom. Values from several calculation methods used to account for internal rotation contribution to entropy were compared to available literature values. We also show comparisons of the entropy when only the torsion frequency contributions from the methyl rotors are replaced by contributions from internal rotor analysis and torsion frequency contributions from other rotors is continued.

2. Computational Methods

The parent and radical species were optimized at the B3LYP [30, 31] level of theory with the 6-31G(d,p) basis set using the Gaussian 03 program suite [32]. Potential energy curves for all single bond rotations were determined for each species using scans at 10° intervals. During the rotation, the geometry was allowed to relax. The scans also serve to verify that our optimized structures were the lowest energy geometries. If there was a lower energy structure, species were reoptimized and scans were rerun to verify that we were working with the lowest energy structures. We note that the relaxed scan may not accurately represent the actual rotation motion for all rotors, as the vibration motion does not necessarily adjust to compensate for internal rotations.

3. Calculation Methods

The three methods that we will use in our entropy analysis are given below.

3.1. SMCPS

Translations, vibrations, and external rotation (denoted TVR) contributions to entropy are determined using the Statistical Mechanics for Heat Capacity and Entropy (SMCPS) program [33]. The SMCPS program applies the rigid-rotor harmonic-oscillator (HO) approximation using the geometry, frequencies, and moments of inertia from the optimized B3LYP/6-31G(d,p) structure along with the mass, electronic degeneracy, symmetry, and number of optical isomers for each species. The zero-point vibration energies (ZPVE) are scaled by 0.9806 as recommended by Scott and Radom [34]. Other than the frequency scaling factor this is used without corrections.

This code is also used for the vibrations with corresponding torsion frequencies removed in the VIBIR and ROTATOR calculation for the internal rotor contributions below.

3.2. VIBIR

This method will analyze internal rotation contributions to entropy with the Pitzer and Gwinn [1820] approximation method as applied in the VIBIR [35] code. The barriers and foldness of the potential energy graphs for each bond rotation are used; there is no adjustment for coupling of internal rotor motion with vibration. The reduced moments of inertia were calculated based on the optimized geometries using the mass and radius of rotation for the rotational groups.

VIBIR assumes that the rotational groups are symmetrical which is accurate for the primary and terminal methyl group rotation in these isomers. We also report results using VIBIR for nonsymmetrical internal rotors for comparison purposes.

3.3. ROTATOR

The second method for internal rotation treatment is with the ROTATOR [36] code. This uses the potential energy curves with expansion of the hindrance potential at discrete torsion angles in the truncated Fourier series: where , , and provide the minima and maxima of the torsion potentials with the allowance for a shift of the theoretical extreme angular positions. It calculates the Hamiltonian matrix on the basis of 1-dimensional wave functions of free internal rotor with subsequent calculation of energy levels by direct diagonalization of the Hamiltonian matrix. The code also calculates the moments of inertia by methods recommended by Pitzer and Gwinn. By calculating all the energy levels of the internal rotor within the potential barriers, the analysis also incorporates contributions from optical isomers such as those of peroxides.

This method allows for a more accurate description of the internal rotor potential in the calculation of energy levels. Both symmetrical and unsymmetrical barriers can be accurately described for calculated contributions to entropies and heat capacities.

4. Results and Discussion

4.1. Entropies of Parent Species

The calculation sets used to calculate entropy contributions in this study can be grouped into four methods. Abbreviations for these methods are denoted in italics below and used in Table 1 and Figure 2. Symmetry corrections are denoted in Table 1 for each method as well as in the Supplementary Materials.(i)SMCPS applies only the HO approximation. (There is no correction (subtraction) of for the three-fold symmetry of the primary methyl groups using the torsion frequencies, as there is in the use of Group Additivity.) (ii)SMCPS/Methyl applies the HO approximation for all vibrations except methyl torsions, which are substituted with values from VIBIR, which shows similar values to ROTATOR for the methyl rotors.(iii)SMCPS/VIBIR applies the HO approximation for vibrations with all torsion frequencies removed from SMCPS and replaced with contributions for internal rotors from use of VIBIR.(iv)SMCPS/ROTATOR applies the HO approximation for vibrations with all torsion frequencies removed from SMCPS and replaced with contributions for internal rotors from the use of ROTATOR.

Table 1 has a comparison between our calculated entropies with data from available reference literature values for the three stable hydrocarbons; this comparison provides a gauge on the methods. To determine the overall performance of the different methods, we calculated the relative difference in the entropies between 298–1500 K based on the literature values from Scott [11]. This is presented graphically in Figure 2.

The data in Table 1 for the parent molecules by the HO approximation method as implemented in the SMCPS code is shown to consistently underpredict entropies by 4–9 cal mol−1 K−1 in this 298–1500 K temperature range. It is noted that as the branching of these isomers increases, the HO approximation shows an increase in accuracy. Overall as the temperature increases we observe a growing deviation.

In a second analysis, each of the parent alkanes analyzed has the terminal methyl rotation contributions from the use of the Pitzer and Gwinn approximation method in the VIBIR code. These symmetrical terminal methyl rotations should be accurately calculated using this method. The contributions from the torsion frequencies associated with terminal methyl rotations are replaced with the VIBIR values and they are incorporated with the remaining vibration contributions at each temperature.

These SMCPS/Methyl values, in Figure 2, are observed to significantly underestimate the entropies by 7–12 cal mol−1 K−1 at 298 K and 17-18 cal mol−1 K−1 at 1500 K compared to the Scott values. In this case, it would be a better estimate to use the HO approximation values from SMCPS than to only consider the methyl rotations.

Calculation of entropy contributions using all of the internal rotors, which are considered to be uncoupled, was performed using both VIBIR and ROTATOR. Determining the corresponding rotation-vibration torsions using visualization in GaussView 5 was problematic due to the coupling of several internal rotors. To evaluate this, we tested several different vibration pair combinations where frequency evaluation in GaussView indicated torsions. It was observed that there was an insignificant change in the entropies using the various combinations. Internal rotations are typically characterized by lower vibrational frequencies as compared to other motions, so to maintain consistency the first six vibration frequencies for each species were replaced with the data from the respective rotor method. The data on these trials is included in the Supplementary Materials available online at http://dx.doi.org/10.1155/2013/673065.

At 298 K, the SMCPS/VIBIR values provide accuracy within 1 cal mol−1 K−1 to the reference values for all of the parent species; SMCPS/ROTATOR provides a slightly closer value for cc2c2cc. As the temperature increases, the SMCPS/VIBIR values trend to underestimate the Scott values by up to 11–13 cal mol−1 K−1 at 1500 K. In contrast, SMCPS/ROTATOR provides a much closer approximation throughout the temperature range. This suggests that modeling all of the rotors, where the rotational barriers are similar or lower to those determined for these hydrocarbon parents and radicals results in a better approximation to entropy.

The HO approximation with corrections for internal rotations by the ROTATOR code provided the best comparison to the available literature reference values. It is important to note that SMCPS/VIBIR provides very good values at 298 K.

4.2. Entropies of Radical Species

Entropies using the SMCPS/ROTATOR and SMCPS/VIBIR methods for the radical species are shown in Table 2. The SMCPS/ROTATOR values are given in bold, and it is seen that the SMCPS/VIBIR values are lower in comparison with an increasing deviation with increasing temperature.

SMCPS/VIBIR provides acceptable agreement with the SMCPS/ROTATOR values at 298 K as it did for the parent species. There is approximately a 2 cal mol−1 K−1 difference for cjc2cccc while the 2,3-dimethylpentane radicals all fall below a 1 cal mol−1 K−1 difference. This is not the case at higher temperatures, however, as the deviation increases with a minimum difference of 9 cal mol−1 K−1 at 1500 K between the two methods. VIBIR is based on the symmetrical barriers, while ROTATOR accurately models the complete potential energy curves for the internal rotation. In the SMCPS/VIBIR method all of the internal rotations which include asymmetrical nonmethyl rotations utilize the average barrier heights for the rotors. The error created by this approximation is not as pronounced at lower temperatures but is apparent in the larger deviation at increased temperatures. We attribute these differences to our use of the Pitzer and Gwinn method in our VIBIR code not being able to account for the nonsymmetrical internal rotor barriers, where the ROTATOR method does fully account for wells in nonsymmetrical barriers.

Since literature values are not available for the entropy of the radicals, our recommendation is to use the SMCPS/ROTATOR values in Table 2; this is based on the calculated results for the parent species.

5. Summary

Internal rotational barriers, moments of inertia, and vibrational frequencies were determined at the B3LYP/6-31G(d,p) level of theory for rotor contribution analysis. These potentials were used to determine entropies from 298 to 1500 K for n-heptane, 2-methylhexane, 2,3-dimethylpentane, and their radicals using torsion frequencies, Pitzer and Gwinn’s analysis, and a Schrodinger equation solution method. The HO approximation-torsion frequency calculation consistently underestimated the entropies compared to reference literature values. Values calculated which were only adjusted for terminal methyl rotations were substantially worse compared to the HO values, especially at increased temperatures. The HO approximation values which used internal rotational contributions from the ROTATOR code provided the best analysis over the 298–1500 K temperature range.

Entropies calculated using the SMCPS/ROTATOR method is recommended for the radical species. We plan on extending this analysis to include enthalpies, heat capacities, and bond dissociation energies in a future publication.

Supplementary Materials

The Supplementary Materials includes the structures, symmetry values, moments of inertia, and vibrational frequencies for the parent and radical heptane isomers in this study. A comparison of calculated values using different torsional frequencies which represented vibrations replaced by internal rotation contributions is also illustrated.

  1. Supplementary Materials
  2. Supplementary Materials