Speed dependent rotational angular momentum polarization of the O2 (a 1Deltag) fragment following ozone photolysis in the wavelength range 248-265 nm

The translational anisotropy and rotational angular momentum polarization of a selection of rotational states of the O2 (a 1Deltag; v=0) photofragment formed from ozone photolysis at 248, 260, and 265 nm have been determined using the technique of resonance enhanced multiphoton ionization in combination with time of flight mass spectrometry. At 248 nm, the dissociation is well described as impulsive in nature with all rotational states exhibiting similarly large, near-limiting values for the bipolar moments describing their angular momentum alignment and orientation. At 265 nm, however, the angular momentum polarization parameters determined for consecutive odd and even rotational states exhibit clear differences. Studies at the intermediate wavelength of 260 nm strongly suggest that such a difference in the angular momentum polarization is speed dependent and this proposal is consistent with the angular momentum polarization parameters extracted and reported previously for longer photolysis wavelengths [G. Hancock et al., Phys. Chem. Chem. Phys. 5, 5386 (2003); S. J. Horrocks et al., J. Chem. Phys. 126, 044308 (2007)]. The alternation of angular momentum polarization for successive odd and even J states may be a consequence of the different mechanisms leading to the formation of the two O2 (a 1Deltag) Lambda doublets. Specifically, the involvement of out of plane parent rotational motion is proposed as the origin for the observed depolarization for the Delta- relative to the Delta+ state.     

FTIR spectroscopic and theoretical study of the photochemistry of matrix-isolated coumarin

The infrared spectrum of monomeric unsubstituted coumarin (C9H6O2; 2H-1-benzopyran-2-one), isolated in solid argon at 10 K is presented and assigned. The UV-induced (lambda>200 nm) unimolecular photochemistry of the matrix-isolated compound was studied experimentally. Three main photoreactions were observed: (a) decarboxylation of the compound and formation of benzocyclobutadiene and CO2, with the Dewar form of coumarin as intermediate; (b) isomerization of the compound, leading to production of a conjugated ketene; and (c) decarbonylation, leading to formation of CO and benzofuran complex. Further decomposition of benzofuran to produce ethynol is suggested. Photochannels (a) and (b) correspond to those previously observed for matrix-isolated alpha-pyrone and its sulfur analogs (Phys. Chem. Chem. Phys. 2004, 6, 929; J. Phys. Chem. A 2006, 110, 6415), while route (c) is similar to the UV-induced photochemistry of coumarin in the gaseous phase (J. Phys. Chem. A 2000, 104, 1095). Interpretation of the experimental data is supported by extensive calculations performed at the B3LYP/6-311++G(d,p), MP2/6-31G(d,p) and MP2/6-311++G(d,p) levels.     

Complexes of type C6H7(+)·L (L = N2 and CO2) studied by explicitly correlated coupled cluster theory

Complexes of the benzenium ion (C(6)H(7)(+)) with N(2) or CO(2) have been studied by explicitly correlated coupled cluster theory at the CCSD(T)-F12x (x = a, b) level [T. B. Adler et al., J. Chem. Phys. 127, 221106 (2007)] and the double-hybrid density functional B2PLYP-D [T. Schwabe and S. Grimme, Phys. Chem. Chem. Phys. 9, 3397 (2007)]. Improved harmonic vibrational wavenumbers for C(6)H(7)(+) have been obtained by CCSD(T?)-F12a calculations with the VTZ-F12 basis set. Combining them with previous B2PLYP-D anharmonic contributions we arrive at anharmonic wavenumbers which are in excellent agreement with recent experimental data from p-H(2) matrix isolation IR spectroscopy [M. Bahou et al., J. Chem. Phys. 136, 154304 (2012)]. The energetically most favourable conformer of C(6)H(7)(+)·N(2) shows a π-bonded structure similar to C(6)H(7)(+)·Rg (Rg = Ne, Ar) [P. Botschwina and R. Oswald, J. Phys. Chem. A 115, 13664 (2011)] with D(e) ≈ 870 cm(-1). For C(6)H(7)(+)·CO(2), a slightly lower energy is calculated for a conformer with the CO(2) ligand lying in the ring-plane of the C(6)H(7)(+) moiety (D(e) ≈ 1508 cm(-1)). It may be discriminated from other conformers through a strong band predicted at 1218 cm(-1), red-shifted by 21 cm(-1) from the corresponding band of free C(6)H(7)(+).     

Experimental and modelling study of the recombination reaction H + O(2) (+M) --> HO(2) (+M) between 300 and 900 K, 1.5 and 950 bar, and in the bath gases M = He, Ar, and N(2)

The recombination reaction H + O(2) (+M) --> HO(2) (+M) was studied by laser flash photolysis in a high pressure flow cell, over the temperature range 300-900 K, the pressure range 1.5-950 bar and in the bath gases M = He and N(2). Earlier experiments by Hahn et al. (Phys. Chem. Chem. Phys. 2004, 6, 1997) in the bath gas M = Ar were also extended. The data were analyzed in terms of unimolecular rate theory employing new calculations of relevant molecular parameters. Improved energy transfer parameters for the bath gases M = He, Ar, N(2), and H(2)O could thus be obtained and complete falloff curves were constructed. In the case of water, the high pressure rates well connect with pulse radiolysis results obtained in supercritical water by Janik et al. (J. Phys. Chem. A 2007, 111, 79).     

Structure factor for hard hyperspheres in higher dimensions

The structure factor for hard hyperspheres in two to eight dimensions is computed by Fourier transforming the pair correlation function obtained by computer simulation at a variety of densities. The resulting structure factors are compared to the known Percus-Yevick equations for odd dimensions and to the model proposed by Leutheusser [J. Chem. Phys. 84, 1050 (1986)] and Rosenfeld [J. Chem. Phys. 87, 4865 (1987)] in even dimensions. It is found that there is fine agreement among all these approaches at low to moderate densities but that the accuracy of the analytical models breaks down as the freezing transition is approached. The structure factor gives another insight into the decrease in the ordering of the hyperspheres as the dimension is increased.     

Identity SN2 reactions X- + CH3X --> XCH3 + X- (X=F, Cl, Br, and I) in vacuum and in aqueous solution: a valence bond study

The recently developed (L. Song, W. Wu, Q. Zhang, S. Shaik, J. Phys. Chem. A 2004, 108, 6017) valence bond method coupled with a polarized continuum model (VBPCM) has been applied to the identity SN2 reaction of halides in the gas phase and in aqueous solution. The barriers computed at the level of the breathing orbital VB method (P. C. Hiberty, J. P. Flament, E. Noizet, Chem. Phys. Lett. 1992, 189, 259), BOVB and VBPCM//BOVB, are comparable to CCSD(T) and CCSD(T)//PCM results and to experimentally derived barriers in solution (W. J. Albery, M. M. Kreevoy, Adv. Phys. Org. Chem. 1978, 16, 85). The reactivity parameters needed to apply the valence bond state correlation diagram (VBSCD) method (S. Shaik, J. Am. Chem. Soc. 1984, 106, 1227), were also determined by VB calculations. It has been shown that the reactivity parameters along with their semiempirical derivations provide a satisfactory qualitative and quantitative account of the barriers.     

Side-chain radical losses from radical cations allows distinction of leucine and isoleucine residues in the isomeric peptides Gly-XXX-Arg

Sequencing of peptides via low-energy collision-induced dissociation of protonated peptides typically yields b(n) and y(n) sequence ions. The isomeric residues leucine and isoleucine rarely can be distinguished in these experiments since they give b(n) and y(n) sequence ions of the same m/z. Siu's pioneering work on electrospray ionization of copper complexes of peptides (Chu IK, Rodriquez CF, Lau TC, Hopkinson AC, Siu KWM. J. Phys. Chem. B 2000; 104: 3393) provides a way of forming radical cations of peptides in the gas phase. This method was used to generate M(+ small middle dot) ions of the two isomeric peptides Gly-Leu-Arg and Gly-Ile-Arg in order to compare their fragmentation reactions. Both radical cations fragment to give even electron y(2) and y(1) sequence ions as well as side-chain radical losses of CH(3) and CH(3)CH(2) for isoleucine and (CH(3))(2)CH for leucine. In contrast the [M + H](+) and [M + 2H](2+) ions do not allow distinction between the isomeric leucine and isoleucine peptides.     

Theoretical and experimental study of weakly bound CO2-(pH2)2 trimers

The infrared spectrum of CO(2)-(pH(2))(2) trimers is predicted by performing exact basis-set calculations on a global potential energy surface defined as the sum of accurately known two-body pH(2)-CO(2) (J. Chem. Phys. 2010, 132, 214309) and pH(2)-pH(2) potentials (J. Chem. Phys. 2008, 129, 094304). These results are compared with new spectroscopic measurements for this species, for which 13 transitions are now assigned. A reduced-dimension treatment of the pH(2) rotation has been employed by applying the hindered-rotor averaging technique of Li, Roy, and Le Roy (J. Chem. Phys. 2010, 133, 104305). Three-body effects and the quality of the potential are discussed. A new technique for displaying the three-dimensional pH(2) density in the body-fixed frame is used, and shows that in the ground state the two pH(2) molecules are localized much more closely together than is the case for the two He atoms in the analogous CO(2)-(He)(2) species. A clear tunneling splitting is evident for the torsional motion of the two pH(2) molecules on a ring about the CO(2) molecular axis, in contrast to the case of CO(2)-(He)(2) where a more regular progression of vibrational levels reflects the much lower torsional barrier.     

Theoretical study of pattern formation during the catalytic oxidation of CO on Pt{100} at low pressures

Theoretical studies have thus far been unable to model pattern formation during the reaction in this system on physically feasible length and time scales. In this paper, we derive a computational reaction-diffusion model for this system in which most of the input parameters have been determined experimentally. We model the surface on a mesoscopic scale intermediate between the microscopic size of CO islands and the macroscopic length scale of pattern formation. In agreement with experimental investigations [M. Eiswirth et al., Z. Phys. Chem., Neue Folge 144, 59 (1985)], the results from our model divide the CO and O(2) partial pressure parameter space into three regions defined by the level of CO coverage or the presence of sustained oscillations. We see CO fronts moving into oxygen-covered regions, with the 1 x 1 to hex phase change occurring at the leading edge. There are also traveling waves consisting of successive oxygen and CO fronts that move into areas of relatively high CO coverage, and in this case, the phase change is more gradual and of lower amplitude. The propagation speed of these reaction waves is similar to those observed experimentally for CO and oxygen fronts [H. H. Rotermund et al., J. Chem. Phys. 91, 4942 (1989); H. H. Rotermund et al., Nature (London) 343, 355 (1990); J. Lauterbach and H. H. Rotermund, Surf. Sci. 311, 231 (1994)]. In the two-dimensional version of our model, the traveling waves take the form of target patterns emitted from surface inhomogeneities.     

Quasiclassical trajectory calculations comparing the reactivity and dynamics of symmetric and asymmetric stretch and the role of the bending mode excitations of methane in the Cl + CH4 reaction

To analyze the effects of the symmetric (nu(1)) and asymmetric (nu(3)) stretch mode excitations and the role played by the "umbrella" bending (nu(4)) mode excitation in the reactivity and the dynamics of the gas-phase Cl+CH(4) reaction, an exhaustive dynamics study was performed. Quasiclassical trajectory (QCT) calculations, including corrections to avoid zero-point energy leakage along the trajectories, were used in this work on an analytical potential energy surface previously developed by Espinosa-Garcia et al. [J. Chem. Phys. (to be published)]. First, with respect to the reactivity, we found that the nu(1) mode excitation is more reactive than the nu(3) mode by a factor of 1.20, in agreement with the experimental tendency between these modes. The inclusion of the nu(4) bending mode practically does not affect this relative reactivity, (nu(1+)nu(4))(nu(3+)nu(4)) = 1.16. Second, with respect to the dynamics (rotovibrational and angular distributions of the products), the two stretch modes, nu(1) and nu(3), give very similar pictures, reproducing the experimental behavior, and the nu(4) "umbrella" mode does not affect the dynamics. The satisfactory reproduction (always qualitatively acceptable and sometimes even quantitatively) of a great variety of experimental data by the QCT study presented here lends confidence to the potential energy surface constructed by Espinosa-Garcia et al. [J. Chem. Phys. (to be published)].     

Parent-molecule rotational depolarization of photofragment angular momentum distributions: diatomic and polyatomic molecules

We extend the A(q)(k) polarization-parameter model, which describes product angular momentum polarization from one photon photodissociation of polyatomic molecules in the molecular frame [J. Chem. Phys., 2010, 132, 224310], to the case of rotating parent molecules. The depolarization of the A(q)(k) is described by a set of rotational depolarization factors that depend on the angle of rotation of the molecular axis γ. We evaluate these rotational depolarization factors for the case of dissociating diatomic molecules and demonstrate that they are in complete agreement with the results of Kuznetsov and Vasyutinskii [J. Chem. Phys., 2005, 123, 034307] obtained from a fully quantum mechanical approach of the same problem, showing the effective equivalence of the two approaches. We further evaluate the set of rotational depolarization factors for the case of dissociating polyatomic molecules that have three (near) equal moments of inertia, thus extending these calculations to polyatomic systems. This ideal case yields insights for the dissociation of polyatomic molecules of various symmetries when we compare the long lifetime limit with the results obtained for the diatomic case. In particular, in the long lifetime limit the depolarization factors of the A(0)(k) (odd k), Re(A(1)(k)) (even k) and Im(A(1)(k)) (odd k) for diatomic molecules vanish; in contrast, for polyatomic molecules the depolarization factors for the A(0)(k) (odd k) reduce to a value of 1/3, whereas for the Re(A(1)(k)) (even k) and Im(A(1)(k)) (odd k) they reduce to 1/5.     

Reactivity enhancement of ultracold O(3P)+H2 collisions by van der Waals interactions

The role of van der Waals forces in O((3)P)+H(2)(upsilon=1,j=0) collisions is investigated theoretically at low and ultralow temperatures. Quantum scattering calculations have been performed for zero total angular momentum using the lowest London-Eyring-Polanyi-Sato double-polynomial (3)A(") potential-energy surface reported by [Rogers et al., J. Phys. Chem. A 104, 2308 (2000)] and its recent BMS1 and BMS2 extensions developed by [Brandao et al., J. Chem. Phys. 121, 8861 (2004)] which provide a more accurate treatment of the van der Waals interaction. Our calculations show that van der Waals forces strongly influence chemical reactivity at ultracold translational energies. The presence of a zero-energy resonance for the BMS1 surface is found to enhance reactivity in the ultracold regime and shift the Wigner threshold to lower temperatures.     

The use of atomic intrinsic polarizabilities in the evaluation of the dispersion energy

The recent approach presented by Becke and Johnson [J. Chem. Phys. 122, 154104 (2005); 123, 024101 (2005); 123, 154101 (2005); 124, 174104 (2006); 124, 014104 (2006)] for the evaluation of dispersion interactions based on the properties of the exchange-hole dipole moment is combined with a Hirshfeld-type partitioning for the molecular polarizabilities into atomic contributions, recently presented by some of the present authors [A. Krishtal et al., J. Chem. Phys. 125, 034312 (2006)]. The results on a series of nine dimers, involving neon, methane, ethene, acetylene, benzene, and CO(2), taken at their equilibrium geometry, indicate that when the C(6), C(8), and C(10) terms are taken into account, the resulting dispersion energies can be obtained deviating 3% or 8% from high level literature data [E. R. Johnson and A. D. Becke, J. Chem. Phys. 124, 174104 (2006)], without the use of a damping function, the only outlier being the parallel face-to-face benzene dimer.     

Accurate prediction of heats of formation by a combined method of B3LYP and neural network correction

Recently, we proposed the X1 method which combines the B3LYP/6‐311+G(3df,2p)//B3LYP/6‐311+G(d,p) method with a neural network correction for an accurate yet efficient prediction of heats of formation (Wu and Xu, J Chem Phys 2007, 127, 214105). In this contribution, we discuss in detail how to set up the X1 neural network. We give examples, showing how to apply the X1 method and how the applicability of X1 can be extended. The overall mean absolute deviation of the X1 method from experiment for the 488 heats of formation is 1.52 kcal/mol compared with 9.44 kcal/mol for the original B3LYP results. © 2008 Wiley Periodicals, Inc. J Comput Chem 2009     

Using heat capacity and compressibility to choose among two-state models of liquid water

We extend to heat capacity Cp the model of Vedamuthu, Singh, and Robinson (J. Phys. Chem. 1994, 98, 2222). This model and that of Bartell (J. Phys. Chem. B 1997, 101, 7573) fit successfully, even in the supercooled region, the temperature dependence of Cp, volume, and isothermal compressibility kappa(T). The Robinson model is superior for kappa(T). Tanaka's model (J. Chem. Phys. 2000, 112, 799) fails for C(p) even after correction of a derivational error. All three models assume that the liquid consists of low-density component 1 and high-density component 2. We conclude that Robinson's tactics, ignoring the intercomponent equilibrium constant and determining compositions solely from volumes, yield the most reliable compositions and individual-component properties. Our fits of the Robinson model to C(p) yield at 0 degrees C H(2) - H(1) of (135 +/- 35) J/g, H(1) - H(ice) of 0.8DeltaH(fus), and C(2) - C(1) of (0.1 +/- 0.7) J/K.g. The enthalpy difference between the components is largely responsible for the rapid change of C(p) at the lowest supercooled temperatures. We propose an adjustment to Speedy and Angell's (J. Chem. Phys. 1976, 65, 851) experimental values of kappa(T) for supercooled water.     

Product distributions, rate constants, and mechanisms of LiH+H reactions

We present a quantum-mechanical investigation of the LiH depletion reaction LiH+H-->Li+H2 and of the H exchange reaction LiH+H'-->LiH'+H. We report product distributions, rate constant, and mechanism of the former, and rate constant and mechanism of the latter reaction. We use the potential-energy surface by Dunne et al. [Chem. Phys. Lett. 336, 1 (2001)], the real-wave-packet method by Gray and Balint-Kurti [J. Chem. Phys. 108, 950 (1998)], and the J-shifting approximation. The 1H2 nuclear-spin statistics and progressions of vib-rotational states (v',j') rule both initial-state-resolved and thermal product distributions, which have saw-toothed shapes with odd j' preferred with respect to even j'. At high collision energies and temperatures, we obtain a regular 3-to-1 intensity alternation of rotational states. At low collision energies and temperatures, the degeneracy and density of many H2 levels can, however, give more irregular distributions. During the collision, the energy flows from the reactant translational mode to the product vibration and recoil ones. The rate constants of both reactions are not Arrhenius type because the reactions are barrier-less. The low-temperature, LiH depletion rate constant is larger than the H exchange one, whereas the contrary holds at high temperature. The real-time mechanisms show the nuclear rearrangements of the nonreactive channel and of the reactive ones, and point out that the LiH depletion is preferred over the H exchange at short times. This confirms the rate-constant results.     

Monte Carlo simulations of critical cluster sizes and nucleation rates of water

We have calculated the critical cluster sizes and homogeneous nucleation rates of water at temperatures and vapor densities corresponding to experiments by Wolk and Strey [J. Phys. Chem B 105, 11683 (2001)]. The calculations have been done with an expanded version of a Monte Carlo method originally developed by Vehkamaki and Ford [J. Chem. Phys. 112, 4193 (2000)]. Their method calculates the statistical growth and decay probabilities of molecular clusters. We have derived a connection between these probabilities and kinetic condensation and evaporation rates, and introduce a new way for the calculation of the work of formation of clusters. Three different interaction potential models of water have been used in the simulations. These include the unpolarizable SPC/E [J. Phys. Chem. 91, 6269 (1987)] and TIP4P [J. Chem. Phys. 79, 926 (1983)] models and a polarizable model by Guillot and Guissani [J. Chem. Phys. 114, 6720 (2001)]. We show that TIP4P produces critical cluster sizes and a temperature and vapor density dependence for the nucleation rate that agree well with the experimental data, although the magnitude of nucleation rate is constantly overestimated by a factor of 2 x 10(4). Guissani and Guillot's model is somewhat less successful, but both the TIP4P and Guillot and Guissani models are able to reproduce a much better experimental temperature dependency of the nucleation rate than the classical nucleation theory. Using SPC/E results in dramatically too small critical clusters and high nucleation rates. The water models give different average binding energies for clusters. We show that stronger binding between cluster molecules suppresses the decay probability of a cluster, while the growth probability is not affected. This explains the differences in results from different water models.     

Characterization of molecularly imprinted polymers employing crosslinkers with nonsymmetric polymerizable groups

Crosslinking monomers have been developed with a combination of methacrylamide and methacrylate or vinyl ketone polymerizable groups that provide molecularly imprinted polymers (MIPs) with improved binding and selectivity. The differential reactivity rates of the polymerizable groups prompted an investigation into the time‐dependent behavior of the crosslinkers, which suggests a new mechanism for MIP formation. The mechanism involves the formation of long sections of linear poly(vinyl ketone) with pendant methacrylamide groups that form a highly crosslinked network in a subsequent step. This has implications for the sequence morphology of polymers, affecting the structure and improving the binding properties of MIPs. © 2004 Wiley Periodicals, Inc. J Polym Sci Part A: Polym Chem 42: 3668–3675, 2004     

Thermochemistry of halomethanes CF(n)Br(4-n) (n = 0-3) based on iPEPICO experiments and quantum chemical computations

Internal energy selected bromofluoromethane cations were prepared and their internal energy dependent fragmentation pathways were recorded by imaging photoelectron photoion coincidence spectroscopy (iPEPICO). The first dissociation reaction is bromine atom loss, which is followed by fluorine atom loss in CF(3)Br and CF(2)Br(2) at higher energies. Accurate 0 K appearance energies have been obtained for these processes, which are complemented by ab initio isodesmic reaction energy calculations. A thermochemical network is set up to obtain updated heats of formation of the samples and their dissociative photoionization products. Several computational methods have been benchmarked against the well-known interhalogen heats of formation. As a corollary, we stumbled upon an assignment issue for the ClF heat of formation leading to a 5.7 kJ mol(-1) error, resolved some time ago, but still lacking closure because of outdated compilations. Our CF(3)(+) appearance energy from CF(3)Br confirms the measurements of Asher and Ruscic (J. Chem. Phys. 1997, 106, 210) and Garcia et al. (J. Phys. Chem. A 2001, 105, 8296) as opposed to the most recent result of Clay et al. (J. Phys. Chem. A 2005, 109, 1541). The ionization energy of CF(3) is determined to be 9.02-9.08 eV on the basis of a previous CF(3)-Br neutral bond energy and the CF(3) heat of formation, respectively. We also show that the breakdown diagram of CFBr(3)(+), a weakly bound parent ion, can be used to obtain the accurate adiabatic ionization energy of the neutral of 10.625 ± 0.010 eV. The updated 298 K enthalpies of formation Δ(f)H(o)(g) for CF(3)Br, CF(2)Br(2), CFBr(3), and CBr(4) are reported to be -647.0 ± 3.5, -361.0 ± 7.4, -111.6 ± 7.7, and 113.7 ± 4 kJ mol(-1), respectively.     

Weak interactions in ion–ligand complexes of C3H3(+) isomers: competition between H-bound and C-bound structures in c-C3H3(+)·L and H2CCCH(+)·L (L = Ne, Ar, N2, CO2, and O2)

Explicitly correlated coupled cluster theory at the CCSD(T)-F12x level (T. B. Adler, G. Knizia, and H.-J. Werner, J. Chem. Phys.127, 221106, 2007) has been employed to study structures and vibrations of complexes of type c-C(3)H(3)(+)·L and H(2)C(3)H(+)·L (L = Ne, Ar, N(2), CO(2), and O(2)). Both cations have different binding sites, allowing for the formation of weak to moderately strong hydrogen bonds as well as "C-bound" or "π-bound" structures. In contrast to previous expectations, the energetically most favourable structures of all H(2)C(3)H(+)·L complexes investigated are "C-bound", with the ligand bound to the methylenic carbon atom. The theoretical predictions enable a more detailed interpretation of infrared photodissociation (IRPD) spectra than was possible hitherto. In particular, the bands observed in the range 3238-3245 cm(-1) (D. Roth and O. Dopfer, Phys. Chem. Chem. Phys.4, 4855, 2002) are assigned to essentially free acetylenic CH stretching vibrations of the propargyl cation in "C-bound" H(2)C(3)H(+)·L complexes.