The dissociation/recombination reaction CH4 (+M) ⇔ CH3 + H (+M): a case study for unimolecular rate theory

The dissociation/recombination reaction CH(4) (+M) ? CH(3) + H (+M) is modeled by statistical unimolecular rate theory completely based on dynamical information using ab initio potentials. The results are compared with experimental data. Minor discrepancies are removed by fine-tuning theoretical energy transfer data. The treatment accounts for transitional mode dynamics, adequate centrifugal barriers, anharmonicity of vibrational densities of states, weak collision and other effects, thus being "complete" from a theoretical point of view. Equilibrium constants between 300 and 5000 K are expressed as K(c) = k(rec)/k(dis) = exp(52,044 K/T) [10(-24.65) (T/300 K)(-1.76) + 10(-26.38) (T/300 K)(0.67)] cm(3) molecule(-1), high pressure recombination rate constants between 130 and 3000 K as k(rec,∞) = 3.34 × 10(-10) (T/300 K)(0.186) exp(-T/25,200 K) cm(3) molecule(-1) s(-1). Low pressure recombination rate constants for M = Ar are represented by k(rec,0) = [Ar] 10(-26.19) exp[-(T/21.22 K)(0.5)] cm(6) molecule(-2) s(-1), for M = N(2) by k(rec,0) = [N(2)] 10(-26.04) exp[-(T/21.91 K)(0.5)] cm(6) molecule(-2) s(-1) between 100 and 5000 K. Weak collision falloff curves are approximated by asymmetric broadening factors [J. Troe and V. G. Ushakov, J. Chem. Phys. 135, 054304 (2011)] with center broadening factors of F(c) ≈ 0.262 + [(T - 2950 K)/6100 K](2) for M = Ar. Expressions for other bath gases can also be obtained.     

Is there a unique melting temperature for two-state proteins?

Thermal unfolding (or folding) in many proteins occurs in an apparent two-state manner, suggesting that only two states, unfolded and folded, are populated. At the melting temperature, Tm, the two states coexist. Using lattice models with side chains we show that individual residues become structured at temperatures that deviate from Tm, which implies that partially folded conformations make substantial contribution to thermodynamic properties of two-state proteins. We also find that the folding cooperativity for a given residue is linked to its accessible surface area. These results are consistent with the experiments on GCN4-like zipper peptide, which showed that local melting temperatures differ from Tm. Analysis of thermal unfolding of six proteins shows that deltaT/Tm approximately N(-1), where deltaT is the transition width and N is the number of residues. This scaling allows us to conclude that, when corrected for finite size effects, folding cooperativity can be captured using coarse grained models.     

Solvent viscosity dependence of the folding rate of a small protein: distributed computing study

By using distributed computing techniques and a supercluster of more than 20,000 processors we simulated folding of a 20-residue Trp Cage miniprotein in atomistic detail with implicit GB/SA solvent at a variety of solvent viscosities (gamma). This allowed us to analyze the dependence of folding rates on viscosity. In particular, we focused on the low-viscosity regime (values below the viscosity of water). In accordance with Kramers' theory, we observe approximately linear dependence of the folding rate on 1/gamma for values from 1-10(-1)x that of water viscosity. However, for the regime between 10(-4)-10(-1)x that of water viscosity we observe power-law dependence of the form k approximately gamma(-1/5). These results suggest that estimating folding rates from molecular simulations run at low viscosity under the assumption of linear dependence of rate on inverse viscosity may lead to erroneous results.     

Refined analysis of the thermal dissociation of formaldehyde

New experimental results for the thermal dissociation of formaldehyde to radical and molecular products (Proc. Combust. Inst. 2007, 31, in press) form the basis of the present analysis of the respective low-pressure rate coefficients k(Rad,0) and k(Mol,0) of the reaction. The article supersedes an earlier analysis (J. Phys. Chem. A 2005, 109, 8320) which used less accurate and more preliminary input information. In addition, refined rotational factors F(rot) are determined and specific energy and angular momentum dependent branching ratios from a more detailed analysis of photolysis quantum yields (J. Phys. Chem. A 2007, 111, 3868) are employed as well. It is emphasized again that pyrolysis and photolysis are intimately related and should be analyzed in an internally consistent manner. The combination of the new with earlier experimental results for pyrolysis rates allows one to fit the height of the energy barrier for the molecular elimination channel with improved precision. A value of E0,1 = 81.7(+/-0.5) kcal mol(-1) is obtained. In addition, employing anharmonicity factors F(anh) from the earlier work, a total average energy transferred per collision of -DeltaE/hc = 100(+/-20) cm(-1) is fitted from the experiments in the bath gas Ar. This value is consistent with the value -DeltaE/hc = 80(+/-40) cm(-1) for the bath gas N(2) such as fitted from photolysis quenching experiments (using the same molecular parameters as for the pyrolysis). Rate coefficients for the temperature range 1200-3500 K are represented in the form k(Mol,0)/[Ar] = 7.3 x 10(14) T -6.1 exp(-47300 K/T) cm(3) molecule(-1) s(-1) and k(Rad,0)/[Ar] = 2.1 x 10(12) T -5.5 exp(-47300 K/T) cm(3) molecule(-1) s(-1) (accuracy +/-25%) and recommended for use in combustion chemistry.     

Statistical mechanics of helix bundles using a dynamic programming approach

Despite much study, biomolecule folding cooperativity is not well understood. There are quantitative models for helix-coil transitions and for coil-to-globule transitions, but no accurate models yet treat both chain collapse and secondary structure formation together. We develop here a dynamic programming approach to statistical mechanical partition functions of foldamer chain molecules. We call it the ascending levels model. We apply it to helix-coil and helix-bundle folding and cooperativity. For 14- to 50-mer Baldwin peptides, the model gives good predictions for the heat capacity and helicity versus temperature and urea. The model also gives good fits for the denaturation of Oas's three-helix bundle B domain of protein A (F13W*) and synthetic protein alpha3C by temperature and guanidine. The model predicts the conformational distributions. It shows that these proteins fold with transitions that are two-state, although the transitions in the Baldwin helices are nearly higher order. The model shows that the recently developed three-helix bundle polypeptoids of Lee et al. fold anti-cooperatively, with a predicted value of DeltaHvH/DeltaHcal = 0.72. The model also predicts that two-helix bundles are unstable in proteins but stable in peptoids. Our dynamic programming approach provides a general way to explore cooperativity in complex foldable polymers.     

Experimental and theoretical study of the carbon-13 and deuterium kinetic isotope effects in the Cl and OH reactions of CH3F

A laser flash photolysis-resonance fluorescence technique has been employed to determine absolute rate coefficients for the CH3F + Cl reaction in N2 bath gas in the temperature range of 200-700 K and pressure range of 33-133 hPa. The data were fitted to a modified Arrhenius expression k(T) = 1.14 x 10(-12) x (T/298)2.26 exp{-313/T}. The OH and Cl reaction rates of (13)CH3F and CD3F have been measured by long-path FTIR spectroscopy relative to CH3F at 298 +/- 2 K and 1013 +/- 10 hPa in purified air. The FTIR spectra were fitted using a nonlinear least-squares spectral fitting method including line data from the HITRAN database and measured infrared spectra as references. The relative reaction rates defined by alpha = k(light)/k(heavy) were determined to be k(OH+CH3F)/k(OH+CD3F) = 4.067 +/- 0.018, k(OH+CH3F)/k(OH+(13)CH3F) = 1.067 +/- 0.006, k(Cl+CH3F)/k(Cl+CD3F) = 5.11 +/- 0.07, and k(Cl+CH3F)/k(Cl+(13)CH3F) = 1.016 +/- 0.006. The carbon-13 and deuterium kinetic isotope effects in the OH and Cl reactions of CH3F have been further investigated by quantum chemistry methods and variational transition state theory.     

Kinetics of the unimolecular decomposition of the 2-chloroallyl radical

The thermal decomposition of the 2-chloroallyl radical, CH(2)CClCH(2) --> CH(2)CCH(2) + Cl (1), was studied using the laser photolysis/photoionization mass spectrometry technique. Rate constants were determined in time-resolved experiments as a function of temperature (720-840 K) and bath gas density ([He] = (3-12) x 10(16), [N(2)] = 6 x 10(16) molecule cm(-3)). C(3)H(4) was observed as a primary product of reaction 1. The rate constants of reaction 1 are in the falloff, close to the low-pressure limit, under the conditions of the experiments. The potential energy surface (PES) of reaction 1 was studied using a variety of quantum chemical methods. The results of the study indicate that the minimum energy path of the CH(2)CClCH(2) dissociation proceeds through a PES plateau corresponding to a weakly bound Cl-C(3)H(4) complex; a PES saddle point exists between the equilibrium CH(2)CClCH(2) structure and the Cl-C(3)H(4) complex. The results of quantum chemical calculations, the rate constant values obtained in the experimental study, and literature data on the reverse reaction of addition of Cl to allene were used to create a model of reactions 1 and -1. The experimental dependences of the rate constants on temperature and pressure were reproduced in RRKM/master equation calculations. The reaction model provides expressions for the temperature dependences of the high-pressure-limit and the low-pressure-limit rate constants and the falloff broadening factors (at T = 300-1600 K): k(infinity)(1) = 1.45 x 10(20)T(-1.75) exp(-19609 K/T) s(-1), k(infinity)(-)(1) = 8.94 x 10(-10)T(-0.40) exp(481 K/T) cm(3) molecule(-1) s(-1), k(1)(0)(He) = 5.01 x 10(-32)T(-12.02) exp(-22788 K/T) cm(3) molecule(-1) s(-1), k(1)(0)(N(2)) = 2.50 x 10(-32)T(-11.92) exp(-22756 K/T) cm(3) molecule(-1) s(-1), F(cent)(He) = 0.46 exp(-T/1001 K) + 0.54 exp(-T/996 K) + exp(-4008 K/T), and F(cent)(N(2)) = 0.37 exp(-T/2017 K) + 0.63 exp(-T/142 K) + exp(-4812 K/T). The experimental data are not sufficient to specify all the parameters of the model; consequently, some of the model parameters were obtained from quantum chemical calculations and from analogy with other reactions of radical decomposition. Thus, the parametrization is most reliable under conditions close to those used in the experiments.     

Robustness of atomistic Gō models in predicting native-like folding intermediates

Go? models are exceedingly popular tools in computer simulations of protein folding. These models are native-centric, i.e., they are directly constructed from the protein's native structure. Therefore, it is important to understand up to which extent the atomistic details of the native structure dictate the folding behavior exhibited by Go? models. Here we address this challenge by performing exhaustive discrete molecular dynamics simulations of a Go? potential combined with a full atomistic protein representation. In particular, we investigate the robustness of this particular type of Go? models in predicting the existence of intermediate states in protein folding. We focus on the N47G mutational form of the Spc-SH3 folding domain (x-ray structure) and compare its folding pathway with that of alternative native structures produced in silico. Our methodological strategy comprises equilibrium folding simulations, structural clustering, and principal component analysis.     

The spin-orbit and rotational constants for the N2 C" 5Pi(ui) (v = 3) state

The spin-orbit (A = -16.4 cm(-1)) and rotational (B = 1.017 cm(-1)) constants for the N2 C" 5Pi(ui)(v = 3) level are determined by a fit to rotational lines in the C" 5Pi(u)-A' 5Sigma(g)+(3-1) band that terminate in J'Omega' = 3(3), 4(3), 3(2), and 4(2) levels of the C" state. The C"-state spin-orbit constant is consistent with semi-empirical estimates, based on spin-orbit constants observed in several other electronic states of N2 and the atomic spin-orbit coupling constant, zeta(N 2p). The C"-A' bands exhibit the unusual feature of oppositely degraded sub-band heads, Omega' = 3 (red) and Omega' = 1, 0, and -1 (blue). The unusually wide range of B(Omega)eff values, from 0.85 cm(-1) (Omega = 3) to 1.28 cm(-1) (Omega = -1) for C" 5Pi(v = 3) should be diagnostically useful for Omega'-assignments. The C" 5Pi(v = 3) level lies 14257.17 and 90599 cm(-1) above A' 5Sigma(g)+(v = 1) and X 1Sigma(g)+(v = 0), respectively, and Re(C" 5Pi) = 1.50 A.     

Conformational entropy of a pseudoknot polymer

The thermodynamics and kinetics of ABAB pseudoknot formation owing to reversible intrachain reactions are investigated for a flexible polymer based on the off-lattice Monte Carlo simulations. The polymer is made of N hard spheres tethered by inextensible bonds and consists of two reactive pairs AA and BB with binding energies -epsilon1 and -epsilon2, respectively, and three loop lengths (l1, l2, and l3). Although two intermediate states, loops A and B, may be formed, the folding path goes mainly through the intermediate loop whose free energy reduction associated with coil-to-loop crossover is greater. The conformational entropy loss is found to follow DeltaS=alpha ln N+G, where alpha approximately 2.48 for coil-loop crossover and alpha approximately 2.43 for loop-pseudoknot crossover. The constant G depends on the three loop lengths and the two end-to-reactive site lengths (L1 and L2). For a given total loop length, G is maximum when the three loop lengths are equal (l1=l2=l3). When l1=l3, the entropy loss is minimum if l2=0. However, the condition l1 not equal l3 makes G even smaller. This consequence indicates that asymmetry in loop lengths is thermodynamically favorable and this fact is consistent with observations of pseudoknotted RNA structures.     

Surface patch binding induced intermolecular complexation and phase separation in aqueous solutions of similarly charged gelatin-chitosan molecules

The formation of selective surface patch binding induced complex coacervates between polyions, chitosan (cationic polyelectrolyte), and alkali-processed gelatin (polyampholyte), both carrying similar net charge, was investigated for two volumetric mixing ratios: r = [chitosan]/[gelatin] = 1:5 and 1:10. Formation of soluble intermolecular complexes between gelatin and chitosan molecules was observed in a narrow range of pH, though these biopolymers had the same kind of net charge, which was evidenced from electrophoretic measurement. This clearly established the role played by selective surface patch binding driven interactions. The temperature sweep measurements conducted on these coacervate samples through rheology and differential scanning calorimetry (DSC) studies yielded two characteristic melting temperatures located at approximately 68 +/- 3 degrees C and 82 +/- 3 degrees C. In the flow mode, the shear viscosity (eta) of the coacervate samples was found to scale with (power-law model) applied shear rate (gamma*) as eta(gamma*) approximately (gamma*)(-k); this yielded k = 0.76 +/- 0.2 (1 s(-1) < gamma* < 100 s(-1)), indicating non-Newtonian behavior. The static structure factor (I(q)) deduced from small angle neutron scattering (SANS) data in the low q (q is the scattering wavevector) (0.018 A(-1) < q < 0.072 A(-1)) region was fitted to the Debye-Bueche regime, I(q) approximately 1/(1 + zeta(2)q(2))2 that yielded a size of zeta approximately 215 +/- 20 A (for r = 1:10) and zeta approximately 260 +/- 20 A (for r = 1:5) samples, implying change in the size of inhomogeneities present with mixing ratio. In the intermediate q region, called the Ornstein-Zernike regime, I(q) approximately 1/(1 + xi(2)q(2)) gave a correlation length of xi approximately 10.0 +/- 2.0 A independent of the mixing ratio. The results taken together imply the existence of a weakly interconnected and heterogeneous network structure inside the coacervate phase separated by domains of polymer-poor regions.     

A kinetic study of the reactions of Ca+ ions with O3, O2, N2, CO2 and H2O

The reactions between Ca(+)(4(2)S(1/2)) and O(3), O(2), N(2), CO(2) and H(2)O were studied using two techniques: the pulsed laser photo-dissociation at 193 nm of an organo-calcium vapour, followed by time-resolved laser-induced fluorescence spectroscopy of Ca(+) at 393.37 nm (Ca(+)(4(2)P(3/2)-4(2)S(1/2))); and the pulsed laser ablation at 532 nm of a calcite target in a fast flow tube, followed by mass spectrometric detection of Ca(+). The rate coefficient for the reaction with O(3) is essentially independent of temperature, k(189-312 K) = (3.9 +/- 1.2) x 10(-10) cm(3) molecule(-1) s(-1), and is about 35% of the Langevin capture frequency. One reason for this is that there is a lack of correlation between the reactant and product potential energy surfaces for near coplanar collisions. The recombination reactions of Ca(+) with O(2), CO(2) and H(2)O were found to be in the fall-off region over the experimental pressure range (1-80 Torr). The data were fitted by RRKM theory combined with quantum calculations on CaO(2)(+), Ca(+).CO(2) and Ca(+).H(2)O, yielding the following results with He as third body when extrapolated from 10(-3)-10(3) Torr and a temperature range of 100-1500 K. For Ca(+) + O(2): log(10)(k(rec,0)/cm(6) molecule(-2) s(-1)) = -26.16 - 1.113log(10)T- 0.056log(10)(2)T, k(rec,infinity) = 1.4 x 10(-10) cm(3) molecule(-1) s(-1), F(c) = 0.56. For Ca(+) + CO(2): log(10)(k(rec,0)/ cm(6) molecule(-2) s(-1)) = -27.94 + 2.204log(10)T- 1.124log(10)(2)T, k(rec,infinity) = 3.5 x 10(-11) cm(3) molecule(-1) s(-1), F(c) = 0.60. For Ca(+) + H(2)O: log(10)(k(rec,0)/ cm(6) molecule(-2) s(-1)) = -23.88 - 1.823log(10)T- 0.063log(10)(2)T, k(rec,infinity) = 7.3 x 10(-11)exp(830 J mol(-1)/RT) cm(3) molecule(-1) s(-1), F(c) = 0.50 (F(c) is the broadening factor). A classical trajectory analysis of the Ca(+) + CO(2) reaction is then used to investigate the small high pressure limiting rate coefficient, which is significantly below the Langevin capture frequency. Finally, the implications of these results for calcium chemistry in the mesosphere are discussed.     

Model calculations of the spontaneous curvature, mean and Gaussian bending constants for a thermodynamically open surfactant film

The spontaneous curvature (H(0)), mean and Gaussian bending constants (k(c) and k (c)), as defined in the well-known Helfrich expression, have been calculated from a detailed model for a thermodynamically open surfactant layer. The effect of head group cross-section area, surfactant tail length and electrolyte concentration for monovalent ionic surfactants have been investigated. Geometrical packing constraints subjected to the aggregated hydrocarbon tails and electrostatics are found to be the dominant contributions to H(0), k(c) and k (c). In addition, the transition from spherocylindrical micelles to vesicles were investigated in terms of the three parameters and the following simple expressions were derived as criteria for coexistence between micelles and vesicles H(0)=1/4 xi and N(ves)/N(mic)=exp[4 pi(k(c)+k (c))/kT], where xi is the thickness of the hydrocarbon part of the film and N(mic) and N(ves) the average aggregation numbers of micelles and vesicles, respectively. However, it is found that the ratio N(ves)/N(mic) is order of magnitudes too large for vesicles to form at all in charged single-surfactant systems where the surfactant head is of moderate size.     

Falloff curves for the Reaction CH3 + O2 (+ M) --> CH3O2 (+ M) in the pressure range 2-1000 bar and the temperature range 300-700 K

The reaction CH(3) + O(2) (+M) --> CH(3)O(2) (+M) was studied in the bath gases Ar and N(2) in a high-temperature/high-pressure flow cell at pressures ranging from 2 to 1000 bar and at temperatures between 300 and 700 K. Methyl radicals were generated by laser flash photolysis of azomethane or acetone. Methylperoxy radicals were monitored by UV absorption at 240 nm. The falloff curves of the rate constants are represented by the simplified expression k/k(infinity) approximately [x/(1 + x)]F(cent)(1/{1+[(log)(x)/)(N)(]2}) with x = k(0)/k(infinity) F(cent) approximately 0.33, and N approximately 1.47, where k(0) and k(infinity) denote the limiting low and high-pressure rate constants, respectively. At low temperatures, 300-400 K, and pressures >300 bar, a fairly abrupt increase of the rate constants beyond the values given by the falloff expressions was observed. This effect is attributed to a contribution from the radical complex mechanism as was also observed in other recombination reactions of larger radicals. Equal limiting low-pressure rate constants k(0) = [M]7 x 10(-31)(T/300 K)(-3.0) cm(6) molecule(-2) s(-1) were fitted for M = Ar and N(2) whereas limiting high-pressure rate constants k(infinity) = 2.2 x 10(-12)(T/300 K)(0.9) cm(3) molecule(-1) s(-1) were approached. These values are discussed in terms of unimolecular rate theory. It is concluded that a theoretical interpretation of the derived rate constants has to be postponed until better information of the potential energy surface is available. Preliminary theoretical evaluation suggests that there is an "anisotropy bottleneck" in the otherwise barrierless interaction potential between CH(3) and O(2).     

Kinetics, mechanism, and thermochemistry of the gas phase reaction of atomic chlorine with dimethyl sulfoxide

A laser flash photolysis-resonance fluorescence technique has been employed to study the kinetics of the reaction of chlorine atoms with dimethyl sulfoxide (CH3S(O)CH3; DMSO) as a function of temperature (270-571 K) and pressure (5-500 Torr) in nitrogen bath gas. At T = 296 K and P > or = 5 Torr, measured rate coefficients increase with increasing pressure. Combining our data with literature values for low-pressure rate coefficients (0.5-3 Torr He) leads to a rate coefficient for the pressure independent H-transfer channel of k1a = 1.45 x 10(-11) cm3 molecule(-1) s(-1) and the following falloff parameters for the pressure-dependent addition channel in N2 bath gas: k(1b,0) = 2.53 x 10(-28) cm6 molecule(-2) s(-1); k(1b,infinity) = 1.17 x 10(-10) cm3 molecule(-1) s(-1), F(c) = 0.503. At the 95% confidence level, both k1a and k1b(P) have estimated accuracies of +/-30%. At T > 430 K, where adduct decomposition is fast enough that only the H-transfer pathway is important, measured rate coefficients are independent of pressure (30-100 Torr N2) and increase with increasing temperature. The following Arrhenius expression adequately describes the temperature dependence of the rate coefficients measured at over the range 438-571 K: k1a = (4.6 +/- 0.4) x 10(-11) exp[-(472 +/- 40)/T) cm3 molecule(-1) s(-1) (uncertainties are 2sigma, precision only). When our data at T > 430 K are combined with values for k1a at temperatures of 273-335 K that are obtained by correcting reported low-pressure rate coefficients from discharge flow studies to remove the contribution from the pressure-dependent channel, the following modified Arrhenius expression best describes the derived temperature dependence: k1a = 1.34 x 10(-15)T(1.40) exp(+383/T) cm3 molecule(-1) s(-1) (273 K < or = T < or = 571 K). At temperatures around 330 K, reversible addition is observed, thus allowing equilibrium constants for Cl-DMSO formation and dissociation to be determined. A third-law analysis of the equilibrium data using structural information obtained from electronic structure calculations leads to the following thermochemical parameters for the association reaction: delta(r)H(o)298 = -72.8 +/- 2.9 kJ mol(-1), deltaH(o)0 = -71.5 +/- 3.3 kJ mol(-1), and delta(r)S(o)298 = -110.6 +/- 4.0 J K(-1) mol(-1). In conjunction with standard enthalpies of formation of Cl and DMSO taken from the literature, the above values for delta(r)H(o) lead to the following values for the standard enthalpy of formation of Cl-DMSO: delta(f)H(o)298 = -102.7 +/- 4.9 kJ mol(-1) and delta(r)H(o)0 = -84.4 +/- 5.8 kJ mol(-1). Uncertainties in the above thermochemical parameters represent estimated accuracy at the 95% confidence level. In agreement with one published theoretical study, electronic structure calculations using density functional theory and G3B3 theory reproduce the experimental adduct bond strength quite well.     

Predicting protein folding rates using the concept of Chou's pseudo amino acid composition

One of the most important challenges in computational and molecular biology is to understand the relationship between amino acid sequences and the folding rates of proteins. Recent works suggest that topological parameters, amino acid properties, chain length and the composition index relate well with protein folding rates, however, sequence order information has seldom been considered as a property for predicting protein folding rates. In this study, amino acid sequence order was used to derive an effective method, based on an extended version of the pseudo-amino acid composition, for predicting protein folding rates without any explicit structural information. Using the jackknife cross validation test, the method was demonstrated on the largest dataset (99 proteins) reported. The method was found to provide a good correlation between the predicted and experimental folding rates. The correlation coefficient is 0.81 (with a highly significant level) and the standard error is 2.46. The reported algorithm was found to perform better than several representative sequence-based approaches using the same dataset. The results indicate that sequence order information is an important determinant of protein folding rates.     

Shock tube study of the reaction of CH with N2: overall rate and branching ratio

We have studied the reaction between CH and N2, (1) CH + N2 --> products, in shock tube experiments using CH and NCN laser absorption. CH was monitored by continuous-wave, narrow-line-width laser absorption at 431.1 nm. The overall rate coefficient of the CH + N2 reaction was measured between 1943 and 3543 K, in the 0.9-1.4 atm pressure range, using a CH perturbation approach. CH profiles recorded upon shock-heating dilute mixtures of ethane in argon and acetic anhydride in argon were perturbed by the addition of nitrogen. The perturbation in the CH concentration was principally due to the reaction between CH and N2. Rate coefficients for the overall reaction were inferred by kinetically modeling the perturbed CH profiles. A least-squares, two-parameter fit of the current overall rate coefficient measurements was k1 = 6.03 x 1012 exp(-11150/T [K]) (cm3 mol-1 s-1). The uncertainty in k1 was estimated to be approximately +/-25% and approximately +/-35% at approximately 3350 and approximately 2100 K, respectively. At high temperatures, there are two possible product channels for the reaction between CH and N2, (1a) CH + N2 --> HCN + N and (1b) CH + N2 --> H + NCN. The large difference in the rates of the reverse reactions enabled inference of the branching ratio of reaction 1, k1b/(k1b + k1a), in the 2228-2905 K temperature range by CH laser absorption in experiments in a nitrogen bath. The current CH measurements are consistent with a branching ratio of 1 and establish NCN and H as the primary products of the CH + N2 reaction. A detailed and systematic uncertainty analysis, taking into account experimental and mechanism-induced contributions, yields a conservative lower bound of 0.70 for the branching ratio. NCN was also detected by continuous-wave, narrow-line-width laser absorption at 329.13 nm. The measured NCN time histories were used to infer the rate coefficient of the reaction between H and NCN, H + NCN --> HCN + N, and to estimate an absorption coefficient for the NCN radical.     

Rate coefficients for the gas-phase reaction of OH radicals with dimethyl sulfide: temperature and O2 partial pressure dependence

Rate coefficients for the gas-phase reaction of hydroxyl (OH) radicals with dimethyl sulfide (CH(3)SCH(3), DMS) have been determined using a relative rate technique. The experiments were performed under different conditions of temperature (250-299 K) and O(2) partial pressure (approximately 0 Torr O(2)-380 Torr O(2)), at a total pressure of 760 Torr bath gas (N(2) + O(2)), in a 336 l reaction chamber, using long path in situ Fourier transform (FTIR) absorption spectroscopy to monitor the disappearance rates of DMS and the reference compounds (ethene, propene and 2-methylpropene). OH was produced by the photolysis of H(2)O(2). The following Arrhenius expressions adequately describe the rate coefficients as a function of temperature (units are cm(3) molecule(-1) s(-1)): k = (1.56 +/- 0.20) x 10(-12) exp[(369 +/- 27)/T], for approximately 0 Torr O(2); (1.31 +/- 0.08) x 10(-14) exp[(1910 +/- 69)/T], for 155 Torr O(2); (5.18 +/- 0.71) x 10(-14) exp[(1587 +/- 24)/T], for 380 Torr O(2). The results are compared with previous investigations.     

Convergence of fine-lattice discretization for near-critical fluids

In simulating continuum model fluids that undergo phase separation and criticality, significant gains in computational efficiency may be had by confining the particles to the sites of a lattice of sufficiently fine spacing, a(0) (relative to the particle size, say a). But a cardinal question, investigated here, then arises; namely, How does the choice of the lattice discretization parameter, zeta identical with a/a(0), affect the values of interesting parameters, specifically, critical temperature and density, T(c) and rho(c)? Indeed, for small zeta ( less, similar 4-8) the underlying lattice can strongly influence the thermodynamic properties. A heuristic argument, essentially exact in d = 1 and d = 2 dimensions, indicates that, for models with hard-core potentials, both T(c)(zeta) and rho(c)(zeta) should converge to their continuum limits as 1/zeta((d)(+1)/2) for d infinity; but the behavior of the error is highly erratic for d >/= 2. For smoother interaction potentials, the convergence is faster. Exact results for d = 1 models of van der Waals character confirm this; however, an optimal choice of zeta can improve the rate of convergence by a factor 1/zeta. For d >/= 2 models, the convergence of the second virial coefficients to their continuum limits likewise exhibits erratic behavior, which is seen to transfer similarly to T(c) and rho(c); but this can be used in various ways to enhance convergence and improve extrapolation to zeta = infinity as is illustrated using data for the restricted primitive model electrolyte.