The Al+ -H(2) cation complex: rotationally resolved infrared spectrum, potential energy surface, and rovibrational calculations

The infrared spectrum of the Al(+)-H(2) complex is recorded in the H-H stretch region (4075-4110 cm(-1)) by monitoring Al(+) photofragments. The H-H stretch band is centered at 4095.2 cm(-1), a shift of -66.0 cm(-1) from the Q(1)(0) transition of the free H(2) molecule. Altogether, 47 rovibrational transitions belonging to the parallel K(a)=0-0 and 1-1 subbands were identified and fitted using a Watson A-reduced Hamiltonian, yielding effective spectroscopic constants. The results suggest that Al(+)-H(2) has a T-shaped equilibrium configuration with the Al(+) ion attached to a slightly perturbed H(2) molecule, but that large-amplitude intermolecular vibrational motions significantly influence the rotational constants derived from an asymmetric rotor analysis. The vibrationally averaged intermolecular separation in the ground vibrational state is estimated as 3.03 A, decreasing by 0.03 A when the H(2) subunit is vibrationally excited. A three-dimensional potential energy surface for Al(+)-H(2) is calculated ab initio using the coupled cluster CCSD(T) method and employed for variational calculations of the rovibrational energy levels and wave functions. Effective dissociation energies for Al(+)-H(2)(para) and Al(+)-H(2)(ortho) are predicted, respectively, to be 469.4 and 506.4 cm(-1), in good agreement with previous measurements. The calculations reproduce the experimental H-H stretch frequency to within 3.75 cm(-1), and the calculated B and C rotational constants to within approximately 2%. Agreement between experiment and theory supports both the accuracy of the ab initio potential energy surface and the interpretation of the measured spectrum.     

Can ortho-para transitions for water be observed?

The spectrum of water can be considered as the juxtaposition of the spectra of two molecules, with different total nuclear spin: ortho-H2O, and para-H2O. No transitions have ever been observed between the two different nuclear-spin isotopomers. The interconversion time is unknown and it is widely assumed that interconversion is forbidden without some other intervention. However, weak nuclear spin-rotation interaction occurs and can drive ortho to para transitions. Ab initio calculations show that the hyperfine nuclear spin-rotational coupling constants are about 30 kHz. These constants are used to explore the whole vibration-rotation spectrum with special emphasis on the coupling between nearby levels. Predictions are made for different spectral regions where the strongest transitions between ortho and para levels of water could be experimentally observed.     

Spectroscopy and potential energy surface of the H2-CO2 van der Waals complex: experimental and theoretical studies

A 4-D ab initio potential energy surface is calculated for the intermolecular interaction of hydrogen and carbon dioxide, using the CCSD(T) method with a large basis set. The surface has a global minimum with a well depth of 212 cm(-1) and an intermolecular distance of 2.98 A for a planar configuration with both the O-C-O and H-H axes perpendicular to the intermolecular axis. Bound state calculations are performed for the H(2)-CO(2) van der Waals complex with H(2) in both the para and ortho spin states, and the binding energy of paraH(2)-CO(2)(50.4 cm(-1)) is found to be significantly less than that of orthoH(2)-CO(2)(71.7 cm(-1)). The surface supports 7 bound intermolecular vibrational states for paraH(2)-CO(2) and 19 for orthoH(2)-CO(2), and the lower rotational levels with J< or = 4 follow an asymmetric rotor pattern. The calculated infrared spectrum of paraH(2)-CO(2) agrees well with experiment. For orthoH(2)-CO(2), the ground state rotational levels allowed by symmetry are found to have (K(a), K(c))=(even, odd) or (odd, even). This somewhat unexpected fact enables the previously observed experimental spectrum to be assigned for the first time, in good agreement with theory, and indicates that the orientation of hydrogen is perpendicular to the intermolecular axis in the ground state of the orthoH(2)-CO(2) complex.     

Quantum mechanical calculations for the H2O + hnu --> O(1D) + H2 photodissociation process

Quantum mechanical wavepacket calculations for the photodissociation of water in the second absorption band are presented. Using O + H2 Jacobi coordinates, partial cross sections for the O(1D) + H2 channel are calculated for different initial rotational states. Conical intersection and Renner-Teller effects are included. The branching ratios for the four accessible dissociation channels at 121.6 nm are in good agreement with experiment (J. Chem. Phys. 1982, 77, 2432). The calculations predict significant rotational and vibrational excitation of the H2 fragments. Photodissociation of ortho and para water produces predominantly, but not exclusively, ortho and para H2 fragments, respectively.     

Rovibrational states of the H2O-H2 complex: an ab initio calculation

All bound rovibrational levels of the H(2)O-H(2) dimer are calculated for total angular momentum J = 0-5 on two recent intermolecular potential surfaces reported by Valiron et al. [J. Chem. Phys. 129, 134306 (2008)] and Hodges et al. [J. Chem. Phys. 120, 710 (2004)] obtained through ab initio calculations. The method used handles correctly the large amplitude internal motions in this complex; it involves a discrete variable representation of the intermolecular distance coordinate R and a basis of coupled free rotor wave functions for the hindered internal rotations and the overall rotation of the dimer. The basis is adapted to the permutation symmetry associated with the para/ortho (p/o) nature of both H(2)O and H(2) as well as to inversion symmetry. Dimers containing oH(2) are more strongly bound than dimers with pH(2), as expected, with dissociation energies D(0) of 33.57, 36.63, 53.60, and 59.04 cm(-1)for pH(2)O-pH(2), oH(2)O-pH(2), pH(2)O-oH(2), and oH(2)O-oH(2), respectively, on the potential of Valiron et al. that corresponds to a binding energy D(e) of 235.14 cm(-1). Rovibrational wave functions are computed as well and the nature of the bound states in the four different dimer species is discussed. Converged rovibrational levels on both potentials agree well with the high-resolution spectrum reported by Weida and Nesbitt [J. Chem. Phys. 110, 156 (1999)]; the hindered internal rotor model that was used to interpret this spectrum is qualitatively correct.     

Infrared spectrum of the I- -D2 anion complex

The infrared spectrum of the I(-)-D(2) anion complex is measured in the D(2) stretch region by monitoring production of I(-) photofragments. The rotationally resolved spectrum consists of two overlapping Sigma-Sigma subbands, redshifted by approximately 58 cm(-1) from the free D(2) fundamental vibrational band. These subbands are associated with absorptions by I(-)-D(2) complexes containing ortho and para forms of D(2). The measured rotational constants are consistent with a 3.79 A separation between I(-) and the D(2) center of mass, contracting by 0.08 A when the D(2) subunit is vibrationally excited. Spectroscopic data are used to generate effective radial potential energy curves describing the interaction of ortho and para D(2) with I(-) from which the dissociation energies of I(-)-D(2)(ortho) and I(-)-D(2)(para) are estimated as D(0)=236 and 297 cm(-1), respectively.     

Calculation of some thermodynamic properties of air plasmas: Internal partition functions,plasma composition,and thermodynamic functions

The thermodynamic properties of air plasmas are calculated in the framework of the C.L.T.E. hypothesis, for pressures varying from 1 to 200 atm, in a temperature range from 1000 to 30,000 K. In these calculations, the neutral mon--, di-, and triatomic species, as well as their positive and negative ions, are taken into account. From the internal partition function values, the plasma compositions and the thermodynamic functions (specific enthalpy, Gibbs and Helmholtz free energies and entropy) are calculated.     

State-to-state differential and relative integral cross sections for rotationally inelastic scattering of H2O by hydrogen

State-to-state differential cross sections (DCSs) for rotationally inelastic scattering of H(2)O by H(2) have been measured at 71.2 meV (574 cm(-1)) and 44.8 meV (361 cm(-1)) collision energy using crossed molecular beams combined with velocity map imaging. A molecular beam containing variable compositions of the (J = 0, 1, 2) rotational states of hydrogen collides with a molecular beam of argon seeded with water vapor that is cooled by supersonic expansion to its lowest para or ortho rotational levels (J(KaKc) = 0(00) and 1(01), respectively). Angular speed distributions of fully specified rotationally excited final states are obtained using velocity map imaging. Relative integral cross sections are obtained by integrating the DCSs taken with the same experimental conditions. Experimental state-specific DCSs are compared with predictions from fully quantum scattering calculations on the most complete H(2)O-H(2) potential energy surface. Comparison of relative total cross sections and state-specific DCSs show excellent agreement with theory in almost all details.     

Ternary alkali metal transition metal acetylides A2MC2 (A = Na, K; M = Pd, Pt)

Ternary transition metal acetylides A2MC2 (A = Na, K; M = Pd, Pt) can be synthesised by reaction of the respective alkali metal acetylide A2C2 with palladium or platinum in an inert atmosphere at about 350 degrees C. The crystal structures are characterised by (infinity)1[M(C2)(2/2)2-] chains, which are separated by the alkali metals (P3m1, Z = 1). The refinement of neutron powder diffraction data gave C-C = 1.263(3) A for Na2PdC2 (Na2PtC2: 1.289(4) A), which is distinctively longer than the expected value for a C-C triple bond (1.20 A). On the basis of band-structure calculations this can be attributed to a strong back-bonding from the metal into the anti-bonding orbitals of the C2 unit. This was further confirmed by Raman spectroscopic investigations, which showed that the wavenumbers of the C-C stretching vibrations in Na2PdC2 and Na2PtC2 are about 100 cm(-1) smaller than in acetylene. 13C MAS-NMR spectra demonstrated that the acetylenic C2 units in the title compounds are very different from those in acetylene. Electrical conductivity measurements and band-structure calculations showed that the black title compounds are semiconductors with a small indirect band gap (approximately 0.2 eV).     

Effects of fluorination on iridium(III) complex phosphorescence: magnetic circular dichroism and relativistic time-dependent density functional theory

We use a combination of low temperature, high field magnetic circular dichroism, absorption, and emission spectroscopy with relativistic time-dependent density functional calculations to reveal a subtle interplay between the effects of chemical substitution and spin-orbit coupling (SOC) in a family of iridium(III) complexes. Fluorination at the ortho and para positions of the phenyl group of fac-tris(1-methyl-5-phenyl-3-n-propyl-[1,2,4]triazolyl)iridium(III) cause changes that are independent of whether the other position is fluorinated or protonated. This is demonstrated by a simple linear relationship found for a range of measured and calculated properties of these complexes. Further, we show that the phosphorescent radiative rate, k(r), is determined by the degree to which SOC is able to hybridize T(1) to S(3) and that k(r) is proportional to the inverse fourth power of the energy gap between these excitations. We show that fluorination in the para position leads to a much larger increase of the energy gap than fluorination at the ortho position. Theory is used to trace this back to the fact that fluorination at the para position increases the difference in electron density between the phenyl and triazolyl groups, which distorts the complex further from octahedral symmetry, and increases the energy separation between the highest occupied molecular orbital (HOMO) and the HOMO-1. This provides a new design criterion for phosphorescent iridium(III) complexes for organic optoelectronic applications. In contrast, the nonradiative rate is greatly enhanced by fluorination at the ortho position. This may be connected to a significant redistribution of spectral weight. We also show that the lowest energy excitation, 1A, has almost no oscillator strength; therefore, the second lowest excitation, 2E, is the dominant emissive state at room temperature. Nevertheless the mirror image rule between absorption and emission is obeyed, as 2E is responsible for both absorption and emission at all but very low (<10 K) temperatures.     

Chlorination of 2-phenoxypropanoic acid with NCP in aqueous acetic acid: using a novel ortho-para relationship and the para/meta ratio of substituent effects for mechanism elucidation

Rate constants were measured for the oxidative chlorodehydrogenation of (R,S)-2-phenoxypropanoic acid and nine ortho-, ten para- and five meta-substituted derivatives using (R,S)-1-chloro-3-methyl-2,6-diphenylpiperidin-4-one (NCP) as chlorinating agent. The kinetics was run in 50% (v/v) aqueous acetic acid acidified with perchloric acid under pseudo-first-order conditions with respect to NCP at temperature intervals of 5 K between 298 and 318 K, except at the highest temperature for the meta derivatives. The dependence of rate constants on temperature was analyzed in terms of the isokinetic relationship (IKR). For the 20 reactions studied at five different temperatures, the isokinetic temperature was estimated to be 382 K, which suggests the preferential involvement of water molecules in the rate-determining step. The dependence of rate constants on meta and para substitution was analyzed using the tetralinear extension of the Hammett equation. The parameter lambda for the para/meta ratio of polar substituent effects was estimated to be 0.926, and its electrostatic modeling suggests the formation of an activated complex bearing an electric charge near the oxygen atom belonging to the phenoxy group. A new approach is introduced for examining the effect of ortho substituents on reaction rates. Using IKR-determined values of activation enthalpies for a set of nine pairs of substrates with a given substituent, a linear correlation is found between activation enthalpies of ortho and para derivatives. The correlation is interpreted in terms of the selectivity of the reactant toward para- or ortho-monosubstituted substrates, the slope of which being related to the ortho effect. This slope is thought to be approximated by the ratio of polar substituent effects from ortho and para positions in benzene derivatives. Using the electrostatic theory of through-space interactions and a dipole length of 0.153 nm, this ratio was calculated at various positions of a charged reaction center along the benzene C1-C4 axis, being about 2.5 near the ring and decreasing steeply with increasing distance until reaching a minimum value of -0.565 at 1.3 nm beyond the aromatic ring. Activation enthalpies and entropies were estimated for substrates bearing the isoselective substituent in either ortho and para positions, being demonstrated that they are much different from the values for the parent substrate. The electrophilic attack on the phenolic oxygen atom by the protonated chlorinating agent is proposed as the rate-determining step, this step being followed by the fast rearrangement of the intermediate thus formed, leading to products containing chlorine in the aromatic ring.     

Water line lists close to experimental accuracy using a spectroscopically determined potential energy surface for H2(16)O, H2(17)O, and H2(18)O

Line lists of vibration-rotation transitions for the H(2) (16)O, H(2) (17)O, and H(2) (18)O isotopologues of the water molecule are calculated, which cover the frequency region of 0-20 000 cm(-1) and with rotational states up to J=20 (J=30 for H(2) (16)O). These variational calculations are based on a new semitheoretical potential energy surface obtained by morphing a high accuracy ab initio potential using experimental energy levels. This potential reproduces the energy levels with J=0, 2, and 5 used in the fit with a standard deviation of 0.025 cm(-1). Linestrengths are obtained using an ab initio dipole moment surface. That these line lists make an excellent starting point for spectroscopic modeling and analysis of rotation-vibration spectra is demonstrated by comparison with recent measurements of Lisak and Hodges [J. Mol. Spectrosc. (unpublished)]: assignments are given for the seven unassigned transitions and the intensity of the strong lines are reproduced to with 3%. It is suggested that the present procedure may be a better route to reliable line intensities than laboratory measurements.     

Millimeter-wave spectroscopy of H(2)C=CD: Tunneling splitting and ortho-para mixing interaction

The H(2)C=CD isotopic species of vinyl radical produced in a supersonic jet expansion by ultraviolet laser photolysis was studied by millimeter-wave spectroscopy. Due to the tunneling motion of the α deuteron, the ground state is split into two components, 0(+) and 0(-). Tunneling-rotation transitions connecting the lower (0(+)) and upper (0(-)) components of the tunneling doublet were observed in the frequency region of 184-334 GHz, including three R- and two Q-branch transitions. Three and two pure rotational transitions in the K(a)=0 and 1 stacks, respectively, were also observed for each of the 0(+) and 0(-) states in the frequency region of 52-159 GHz. Least-squares analysis of the observed frequencies for the tunneling-rotation and pure rotational transitions with well resolved hyperfine structures yielded a set of precise molecular constants, among which the tunneling splitting in the ground state was determined to be ΔE(0)=1187.234(17)?MHz, which is 1/14 that for H(2)C=CH. The potential barrier height derived from the observed tunneling splitting by an analysis of the tunneling dynamics using a one-dimensional model is 1545?cm(-1), consistent with the value 1568?cm(-1) obtained for the normal vinyl. The observed spectrum was found to be perturbed by a hyperfine interaction connecting ortho and para levels. The constant for the interaction, which we call the ortho-para mixing Fermi contact interaction, has been determined to be δa(F) ((β))=68.06(53)?MHz. This is believed to be the first definite detection of such an interaction. By this interaction the ortho and para states of H(2)C=CD are mixed up to about 0.1%. The constant is more than 1000 times larger than spin-rotation interaction constants that cause ortho-para mixing in closed shell molecules and suggests extremely rapid conversion between the ortho and para nuclear spin isomers of H(2)C=CD.     

Near-resonant energy transfer from highly vibrationally excited OH to N2

The probability per collision P(T) of near-resonant vibration-to-vibration energy transfer (ET) of one quantum of vibrational energy from vibrational levels nu=8 and nu=9 of OH to N(2)(nu=0), OH(nu)+N(2)(0)-->OH(nu-1)+N(2)(1), is calculated in the 100-350 K temperature range. These processes represent important steps in a model that explains the enhanced 4.3 microm emission from CO(2) in the nocturnal mesosphere. The calculated energy transfer is mediated by weak long-range dipole-quadrupole interaction. The results of this calculation are very sensitive to the strength of the two transition moments. Because of the long range of the intermolecular potential, the resonance function, a measure of energy that can be efficiently exchanged between translation and vibration-rotation degrees of freedom, is rather narrow. A narrow resonance function coupled with the large rotational constant of OH is shown to render the results of the calculation very sensitive to the rotational distribution, or the rotational temperature if one exists, of this molecule. The calculations are carried out in the first and second orders of perturbation theory with the latter shown to give ET probabilities that are an order of magnitude larger than the former. The reasons for the difference in magnitude and temperature dependence of the first- and second-order calculations are discussed. The results of the calculations are compared with room temperature measurements as well as with an earlier calculation. Our calculated results are in good agreement with the room temperature measurements for the transfer of vibrational energy for the exothermic OH(nu=9) ET process but are about an order lower than the room temperature measurements for the exothermic OH(nu=8) ET process. The cause of this discrepancy is explored. This calculation does not give the large values of the rate coefficients needed by the model that explains the enhanced 4.3 microm emission from CO(2) in the nocturnal mesosphere.     

Gating of Quantum Interference in Molecular Junctions by Heteroatom Substitution

To guide the choice of future synthetic targets for single‐molecule electronics, qualitative design rules are needed, which describe the effect of modifying chemical structure. Here the effect of heteroatom substitution on destructive quantum interference (QI) in single‐molecule junctions is, for the first time experimentally addressed by investigating the conductance change when a “parent” meta‐phenylene ethylene‐type oligomer (m‐OPE) is modified to yield a “daughter” by inserting one nitrogen atom into the m‐OPE core. We find that if the substituted nitrogen is in a meta position relative to both acetylene linkers, the daughter conductance remains as low as the parent. However, if the substituted nitrogen is in an ortho position relative to one acetylene linker and a para position relative to the other, destructive QI is alleviated and the daughter conductance is high. This behavior contrasts with that of a para‐connected parent, whose conductance is unaffected by heteroatom substitution. These experimental findings are rationalized by transport calculations and also agree with recent “magic ratio rules”, which capture the role of connectivity in determining the electrical conductance of such parents and daughters.     

Structural and energetic aspects of the protonation of phenol,catechol, resorcinol,and hydroquinone

The various protonated forms of phenol (1), catechol (2), resorcinol (3), and hydroquinone (4) were explored by ab initio quantum chemical calculations at the MP2/6-31G(d) and B3LYP/6-31G(d) levels. Proton affinities (PA) of 1-4 were calculated by the combined G2(MP2,SVP) method, and their gas-phase basicities were estimated after calculation of the change in entropy on protonation. These theoretical data were compared with the corresponding experimental values determined in a high-pressure mass spectrometer. This comparison confirmed that phenols are essentially carbon bases and that protonation generally occurs in a position para to the hydroxyl group. Resorcinol is the most effective base (PA = 856 kJ mol-1) due to the participation of both oxygen atoms in the stabilization of the protonated form. Since protonation is accompanied by a freezing of the two internal rotations, a significant decrease in entropy is observed. The basicity of catechol (PA = 823 kJ mol-1) is due to the existence of an intramolecular hydrogen bond, which is strengthened upon protonation. The lower basicity of hydroquinone (PA = 808 kJ mol-1) is a consequence of the fact that protonation necessarily occurs in a position ortho to the hydroxyl group. When the previously published data are reconsidered and a corrected protonation entropy is used, a proton affinity value of 820 kJ mol-1 is obtained for phenol.     

Theoretical study on the molecular structures of toluene,para‐fluorotoluene,para‐chlorotoluene,and 4‐methylpyridine and their sixfold internal rotational barriers

Two kinds of sixfold internal rotational configurations of toluene, para‐fluorotoluene, para‐chlorotoluene, and 4‐methylpyridine were calculated using Hartree–Fock (HF), second‐order Møller–Plesset (MP2), and Beck's three parameter hybrid functional using the LYP correlation functional (B3LYP) theory methods with various high‐level basis sets. Structures and energies were compared for different configurations. Calculations indicate that the orthogonal configuration has a local minimum while the planar configuration is a transition structure. Furthermore, geometries of the orthogonal and the planar configurations are quite similar, except for a methyl CH bond. Sixfold internal rotational barriers were calculated from the energy difference of two different configurations. For the calculated results, HF methods underestimated the rotational barriers, but MP2 calculations overestimated them. However, the density functional theory (DFT) method is a reliable method since the calculated internal rotational barriers are similar to the experimental ones. © 2000 John Wiley & Sons, Inc. Int J Quant Chem 77: 772–778, 2000     

H atom yields from the reactions of CN radicals with C2H2, C2H4, C3H6, trans-2-C4H8, and iso-C4H8

The kinetics and H atom channel yield at both 298 and 195 K have been determined for reactions of CN radicals with C2H2 (1.00+/-0.21, 0.97+/-0.20), C2H4 (0.96+/-0.032, 1.04+/-0.042), C3H6 (pressure dependent), iso-C4H8 (pressure dependent), and trans-2-C4H8 (0.039+/-0.019, 0.029+/-0.047) where the first figure in each bracket is the H atom yield at 298 K and the second is that at 195 K. The kinetics of all reactions were studied by monitoring both CN decay and H atom growth by laser-induced fluorescence at 357.7 and 121.6 nm, respectively. The results are in good agreement with previous studies where available. The rate coefficients for the reaction of CN with trans-2-butene and iso-butene have been measured at 298 and 195 K for the first time, and the rate coefficients are as follows: k298K=(2.93+/-0.23)x10(-10) cm3 molecule(-1) s(-1), k195K=(3.58+/-0.43)x10(-10) cm3 molecule(-1) s(-1) and k298K=(3.17+/-0.10)x10(-10) cm3 molecule(-1) s(-1), k195K=(4.32+/-0.35)x10(-10) cm3 molecule(-1) s(-1), respectively, where the errors represent a combination of statistical uncertainty (2sigma) and an estimate of possible systematic errors. A potential energy surface for the CN+C3H6 reaction has been constructed using G3X//UB3LYP electronic structure calculations identifying a number of reaction channels leading to either H, CH3, or HCN elimination following the formation of initial addition complexes. Results from the potential energy surface calculations have been used to run master equation calculations with the ratio of primary:secondary addition, the average amount of downward energy transferred in a collision DeltaEd, and the difference in barrier heights between H atom elimination and an H atom 1, 2 migration as variable parameters. Excellent agreement is obtained with the experimental 298 K H atom yields with the following parameter values: secondary addition complex formation equal to 80%, DeltaEd=145 cm(-1), and the barrier height for H atom elimination set 5 kJ mol(-1) lower than the barrier for migration. Finally, very low temperature master equation simulations using the best fit parameters have been carried out in an increased precision environment utilizing quad-double and double-double arithmetic to predict H and CH3 yields for the CN+C3H6 reaction at temperatures and pressures relevant to Titan. The H and CH3 yields predicted by the master equation have been parametrized in a simple equation for use in modeling.     

Structural characteristics of phenol formaldehyde novolak resin depending on polycondensation type using molecular simulation

Phenol formaldehyde novolak resins have various structures depending on the polycondensation types. Their structures were characterized using molecular mechanics and molecular dynamics. Dimer, tetramer, hexamer, octamer, and decamer of the resins with the ortho–ortho, ortho–para, and para–para sequences were calculated. The ortho–ortho resins have the structural characteristics of intramolecular hydrogen bonds between hydroxyl groups of the adjacent phenolic units. For the ortho–para and para–para resins, the intramolecular hydrogen bonds are formed mainly between hydroxyl groups of the backbone phenolic units. The para–para resins also have intramolecular hydrogen bonds between hydroxyl groups of the branched phenolic units. A factor determining the structural characteristics of the resins was found to be the geometry of the basic unit (dimer). The order of the end‐to‐end distances between hydrogen atoms on the para‐position of the basic units of the resins is ortho–ortho resin < ortho–para resin < para–para resin. The calculational results were found to be consistent with the gel permeation chromatography (GPC) analysis. Copyright © 2001 John Wiley & Sons, Ltd.     

Nuclear spin dependence of the reaction of H(3)+ with H2. I. Kinetics and modeling

The chemical reaction H(3)(+) + H(2) → H(2) + H(3)(+) is the simplest bimolecular reaction involving a polyatomic, yet is complex enough that exact quantum mechanical calculations to adequately model its dynamics are still unfeasible. In particular, the branching fractions for the "identity," "proton hop," and "hydrogen exchange" reaction pathways are unknown, and to date, experimental measurements of this process have been limited. In this work, the nuclear-spin-dependent steady-state kinetics of the H(3)(+) + H(2) reaction is examined in detail, and employed to generate models of the ortho:para ratio of H(3)(+) formed in plasmas of varying ortho:para H(2) ratios. One model is based entirely on nuclear spin statistics, and is appropriate for temperatures high enough to populate a large number of H(3)(+) rotational states. Efforts are made to include the influence of three-body collisions in this model by deriving nuclear spin product branching fractions for the H(5)(+) + H(2) reaction. Another model, based on rate coefficients calculated using a microcanonical statistical approach, is appropriate for lower-temperature plasmas in which energetic considerations begin to compete with the nuclear spin branching fractions. These models serve as a theoretical framework for interpreting the results of laboratory studies on the reaction of H(3)(+) with H(2).