Rotationally resolved C2 symmetric conformers of bis-(4-hydroxyphenyl)methane: prototypical examples of excitonic coupling in the S1 and S2 electronic states

Rotationally resolved microwave and ultraviolet spectra of jet-cooled bis-(4-hydroxyphenyl)methane (b4HPM) have been obtained using Fourier-transform microwave and UV laser/molecular beam spectrometers. A recent vibronic level study of b4HPM [Rodrigo, C. P.; Mu?ller, C. W.; Pillsbury, N. R.; James, W. H., III; Plusquellic, D. F.; Zwier, T. S. J. Chem. Phys. 2011, 134, 164312] has assigned two conformers distinguished by the orientation of the in-plane OH groups and has identified two excitonic origins in each conformer. In the present study, the rotationally resolved bands of all four states have been well-fit to asymmetric rotor Hamiltonians. For the lower exciton (S(1)) levels, the transition dipole moment (TDM) orientations are perpendicular to the C(2) symmetry axes and consist of 41(2):59(2) and 34(2):66(2)% a:c hybrid-type character. The S(1) levels are therefore delocalized states of B symmetry and represent the antisymmetric combinations of the zero-order locally excited states of the p-cresol-like chromophores. The TDM polarizations of bands located at ≈132 cm(-1) above the S(1) origins are exclusively b-type and identify them as the upper exciton S(2) origin levels of A symmetry. The TDM orientations and the relative band strengths from the vibronic study have been analyzed within a dipole-dipole coupling model in terms of the localized TDM orientations, μ(loc), on the two chromophores. The out-of-the-ring plane angles of μ(loc) are both near 20° and are similar to results for diphenylmethane [Stearns, J. A.; Pillsbury, N. R.; Douglass, K. O.; Mu?ller, C. W.; Zwier, T. S.; Plusquellic, D. F. J. Chem. Phys. 2008, 129, 224305]. The in-plane angles are, however, rotated by 14 and 18° relative to DPM and, in part, explain the smaller than expected exciton splittings of these two conformers.     

Kinetics and mechanism of the oxidation of acetylacetone by permanganate ion

The kinetics and mechanism of permanganate ion oxidation of acetylacetone (Acac) was studied in acidic and alkaline media. The rate constants for keto, enol, and enolate anions were determined and discussed. Delocalization of the π‐electrons of the double bond by conjugation results in a slower oxidation rate of enol than can be usually observed for unsaturated compounds. In the case of the keto form, the acid‐catalyzed nucleophilic attack of permanganate ion occurs on the carbonyl‐C atom. For enolate anion a mechanism with basis‐catalyzed electron abstraction is suggested. © 2006 Wiley Periodicals, Inc. Int J Chem Kinet 38: 444–450, 2006     

Detection of vibrational bending mode ν8 and overtone bands of the propargyl radical, HCCCH2 X̃ 2B1

Infrared (IR) absorption spectra of matrix-isolated HCCCH(2) have been measured. Propargyl radicals were generated in a supersonic pyrolysis nozzle, using a method similar to that described in a previous study (Jochnowitz, E. B.; Zhang, X.; Nimlos, M. R.; Varner, M. E.; Stanton, J. F.; Ellison, G. B. J. Phys. Chem. A 2005, 109, 3812-3821). Besides the nine vibrational modes observed in the previous study, this investigation detected the HCCCH(2) X? (2)B(1) out-of-plane bending mode (ν(8)) at 378.0 (±1.9) cm(-1) in a cryogenic argon matrix. This is the first experimental observation of ν(8) for the propargyl radical. In addition, seven overtone and combination bands have also been detected and assigned. Ab initio coupled-cluster anharmonic force field calculations were used to guide the analysis. Furthermore, ν(12), the HCCCH(2) in-plane bending mode, has been assigned to 333 (±10) cm(-1) based on the detection of its overtone (2ν(12), 667.7 ± 1.0 cm(-1)) and a possible combination band (ν(10) + ν(12), 1339.0 ± 0.8 cm(-1)). This is the first experimental estimation of ν(12) for the propargyl radical.     

High-resolution absorption spectrum of jet-cooled CH3Cl between 70,000 and 85,000 cm(-1): new assignments

The absorption spectrum of jet-cooled CH(3)Cl was photographed from 165 to 117 nm (or 60,000 - 85,000 cm(-1), 7.5-10.5 eV) at a resolution limit of 0.0008 nm (0.3-0.6 cm(-1) or 0.04-0.08 meV). Even in the best structured region of the spectrum, from 70,000 to 85,000 cm(-1) (8.7-10.5 eV), observed bandwidths (full width at half maximum) are large, from 50 to 150 cm(-1). No rotational feature could be resolved. The spectrum is dominated by two strong bands near 9 eV, 140 nm, the D and E bands of Mulliken [J. Chem. Phys. 8, 382 (1940)] or the spectral region D of Price [J. Chem. Phys.4, 539 (1936)]. Their relative intensity is incompatible with previous assignments, namely, to a triplet and a singlet state belonging to the same configuration. On the basis of the present ab initio calculations, those bands are now assigned to two singlet states, the (1)A(1) and (1)E excited states resulting from the 2e(3)4pe Rydberg configuration. The present calculations also reveal that the two (1)E states issued from 2e(3)4sa(1) and 2e(3)4pa(1) are quasidegenerate and strongly mixed. They should be assigned to the two broad bands near 8 eV, 160 nm, the B and C bands of Mulliken and Price. Three vibrational modes are observed to be active: the CCl bond stretch nu(3)(a(1)), and the CH(3) umbrella and rocking vibrations, respectively, nu(2)(a(1)) and nu(6)(e). The fundamental frequencies deduced are well within the ranges defined by the corresponding values in the neutral and ion ground states. The possibility of a dynamical Jahn-Teller effect induced by the nu(6)(e) vibrational mode in the (1)E Rydberg states is discussed.     

Theoretical study of [Si,O,C,O] species: Prediction of new species on triplet potential energy surface

The intermediates [Si,O,C,O] of the Si + CO₂ reaction have been studied in detail using high level ab iniitio methods. Both singlet and triplet [Si,O,C,O] species are characterized structurally and energetically. On the singlet potential energy surface (PES), the vdw‐OSi–CO isomer and in the triplet PES, the bent‐SiOCO isomer is found to be thermodynamically as well as kinetically most stable species. All possible isomerization transition states (TS) are located on both singlet and triplet potential surfaces. On the triplet surface, the stability of the trans‐OSiCO isomer is comparable with that of the bent‐SiOCO isomer. A non‐planar cis‐SiOCO isomer is located on the triplet PES, which is predicted for the first time. Heats of formation at 0 K (Δ_fH°, 0 K) for all singlet and triplet species are computed using G3B3, G3MP2, and CBS‐Q theories. The discrepancy between G3B3 and the other two methods for the heat of formation value for triplet trans‐OSiCO is discussed. The PESs for singlet as well as triplet species with their dissociation asymptotes are explored at the CCSD(T)/6‐311G(d,p)//MP2/6‐311G(d,p) level of theory. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

Rate constants and H atom branching ratios of the gas-phase reactions of methylidyne CH(X2Pi) radical with a series of alkanes

The reactions of the CH radical with several alkanes were studied, at room temperature, in a low-pressure fast-flow reactor. CH(X2Pi, v = 0) radicals were obtained from the reaction of CHBr(3) with potassium atoms. The overall rate constants at 300 K are (0.76 +/- 0.20) x 10(-10) [Fleurat-Lessard, P.; Rayez, J. C.; Bergeat, A.; Loison, J. C. Chem. Phys. 2002, 279, 87],1 (1.60 +/- 0.60) x 10(-10)[Galland, N.; Caralp, F.; Hannachi, Y.; Bergeat, A.; Loison, J.-C. J. Phys. Chem. A 2003, 107, 5419],2 (2.20 +/- 0.80) x 10(-10), (2.80 +/- 0.80) x 10(-10), (3.20 +/- 0.80) x 10(-10), (3.30 +/- 0.60) x 10(-10), and (3.60 +/- 0.80) x 10(-10) cm3 molecule(-1) s(-1), (errors refer to +/-2sigma) for methane, ethane, propane, n-butane, n-pentane, neo-pentane, and n-hexane respectively. The experimental overall rate constants correspond to those obtained using a simple classical capture theory. Absolute atomic hydrogen production was determined by V.U.V. resonance fluorescence, with H production from the CH + CH4 reaction being used as a reference. Observed H branching ratios were for CH4, 1.00[Fleurat-Lessard, P.; Rayez, J. C.; Bergeat, A.; Loison, J. C. Chem. Phys. 2002, 279, 87];1 C(2)H(6), 0.22 +/- 0.08 [Galland, N.; Caralp, F.; Hannachi, Y.; Bergeat, A.; Loison, J.-C. J. Phys. Chem. A 2003, 107, 5419];2 C(3)H(8), 0.19 +/- 0.07; C(4)H(10) (n-butane), 0.14 +/- 0.06; C(5)H(12) (n-pentane), 0.52 +/- 0.08; C(5)H(12) (neo-pentane), 0.51 +/- 0.08; C(5)H(12) (iso-pentane), 0.12 +/- 0.06; C(6)H(14) (n-hexane), 0.06 +/- 0.04.     

What is the best DFT functional for vibronic calculations? A comparison of the calculated vibronic structure of the S1-S0 transition of phenylacetylene with cavity ringdown band intensities

The sensitivity of vibronic calculations to electronic structure methods and basis sets is explored and compared to accurate relative intensities of the vibrational bands of phenylacetylene in the S(1)(A(1)B(2)) ← S(0)(X(1)A(1)) transition. To provide a better measure of vibrational band intensities, the spectrum was recorded by cavity ringdown absorption spectroscopy up to energies of 2000 cm(-1) above the band origin in a slit jet sample. The sample rotational temperature was estimated to be about 30 K, but the vibrational temperature was higher, permitting the assignment of many vibrational hot bands. The vibronic structure of the electronic transition was simulated using a combination of time-dependent density functional theory (TD-DFT) electronic structure codes, Franck-Condon integral calculations, and a second-order vibronic model developed previously [Johnson, P. M.; Xu, H. F.; Sears, T. J. J. Chem. Phys. 2006, 125, 164331]. The density functional theory (DFT) functionals B3LYP, CAM-B3LYP, and LC-BLYP were explored. The long-range-corrected functionals, CAM-B3LYP and LC-BLYP, produced better values for the equilibrium geometry transition moment, but overemphasized the vibronic coupling for some normal modes, while B3LYP provided better-balanced vibronic coupling but a poor equilibrium transition moment. Enlarging the basis set made very little difference. The cavity ringdown measurements show that earlier intensities derived from resonance-enhanced multiphoton ionization (REMPI) spectra have relative intensity errors.     

Doppler-free two-photon excitation spectroscopy and the Zeeman effects. Perturbations in the 14(0)1 and 1(0)(1)14(0)1 bands of the S1 <-- S0 transition of C6D6

Doppler-free two-photon excitation spectra and the Zeeman effects for the 1 band of the S1 1B2u <-- S0 1A1g transition in gaseous benzene-d6 were measured. Although the spectral lines were strongly perturbed, almost all of the lines near the band origin could be assigned. From a deperturbation analysis, the perturbation near the band origin was identified as originating from an anharmonic resonance interaction. Perturbation centered at K = 28-29 in the 14(0)1 band was analyzed, and it was identified as originating from a perpendicular Coriolis interaction. The symmetry and the assignment of the perturbing state proposed by Schubert et al. (Schubert, U.; Riedle, E.; Neusser, H. J. J. Chem. Phys. 1989, 90, 5994.) were confirmed. No perturbation originating from an interaction with a triplet state was observed in both bands. From the Zeeman spectra and the analysis, it is demonstrated that rotationally resolved levels are not mixed with a triplet state. The intersystem mixing is not likely to occur at levels of low excess energy in the S1 state of an isolated benzene. Nonradiative decay of an isolated benzene in the low vibronic levels of the S1 state will occur through the internal mixing followed by the rotational and vibrational relaxation in the S0 state.     

Extension of the renormalized coupled-cluster methods exploiting left eigenstates of the similarity-transformed hamiltonian to open-shell systems: a benchmark study

The recently formulated completely renormalized coupled-cluster method with singles, doubles, and noniterative triples, exploiting the biorthogonal form of the method of moments of coupled-cluster equations (Piecuch, P.; W?och, M. J. Chem. Phys. 2005, 123, 224105; Piecuch, P.; W?och, M.; Gour, J. R.; Kinal, A. Chem. Phys. Lett. 2006, 418, 467), termed CR-CC(2,3), is extended to open-shell systems. Test calculations for bond breaking in the OH radical and the F2+ ion and singlet-triplet gaps in the CH2, HHeH, and (HFH)- biradical systems indicate that the CR-CC(2,3) approach employing the restricted open-shell Hartree--Fock (ROHF) reference is significantly more accurate than the widely used CCSD(T) method and other noniterative triples coupled-cluster approximations without making the calculations substantially more expensive. A few molecular examples, including the activation energies of the C2H4 + H --> C2H5 forward and reverse reactions and the triplet states of the CH2 and H2Si2O2 biradicals, are used to show that the dependence of the ROHF-based CR-CC(2,3) energies on the method of canonicalization of the ROHF orbitals is, for all practical purposes, negligible.     

Ab initio equation of state of an organic molecular crystal: 1,1-diamino-2,2-dinitroethylene

A complete equation of state for the molecular crystal 1,1-diamino-2,2-dinitroethylene has been calculated from first principles for temperatures between 0 and 400 K, and for specific volumes from 61 to 83 cm3/mol, corresponding to relative volumes from 0.78 to 1.06. The calculated 300 K isotherm agrees very well with the experimentally measured pressure-volume relation reported by Peiris et al. (Peiris, S. M.; Wong, C. P.; Zerilli, F. J. J. Chem. Phys. 2004, 120, 8060). The volumetric thermal expansion coefficient is calculated to be 140 ppm/K at 300 K and atmospheric pressure and varies considerably with specific volume as well as temperature. The Grüneisen parameter varies significantly with temperature, but its variation with specific volume is small. The calculated specific heat (160 J/mol/K at 300 K and atmospheric pressure) has only a very small dependence on specific volume.     

Do spin state changes matter in organometallic chemistry? A computational study

Spin changes occur often in organometallic chemistry, and their effect on kinetics is not well understood. We report computations on the singlet and triplet potential energy surfaces of several processes of this type and show that the topology of the individual surfaces, as well as of the crossing regions between them, can be used to rationalize the observed reactivity in all cases. In particular, the slow addition of dihydrogen to W[N(CH(2)CH(2)NSiMe(3))(3)]H (Schrock, R. R.; Shih, K. Y.; Dobbs, D. A.; Davis, W. M. J. Am. Chem. Soc. 1995, 117, 6609) is shown to be a "spin-blocked" reaction with a high barrier due to the crossing between reactant triplet and product singlet surfaces. In contrast, addition of CO to TpCo(CO) (Detrich, J. L.; Reinaud, O. M.; Rheingold, A. L.; Theopold, K. H. J. Am. Chem. Soc. 1995, 117, 11745) is fast because the triplet and singlet surfaces cross at low energy. Particular care is taken to use DFT methods which are in adequate agreement with experimental and high-level computational energetics for these systems.     

Two-state reactivity in alkane hydroxylation by non-heme iron-oxo complexes

Density functional theory is used to explore the mechanisms of alkane hydroxylation for four synthetic non-heme iron(IV)-oxo complexes with three target substrates (Kaizer, J.; Klinker, E. J.; Oh, N. Y.; Rohde; J.-U.; Song, W. J.; Stubna, A.; Kim, J.; Münck, E.; Nam, W.; Que, L., Jr. J. Am. Chem. Soc. 2004, 126, 472-473; Rohde, J.-U.; Que, L., Jr. Angew. Chem. Int. Ed. 2005, 44, 2255-2258.). The iron-oxo reagents possess triplet ground states and low-lying quintet excited states. The set of experimental and theoretical reactivity trends can be understood if the reactions proceed on the two spin states, namely two-state reactivity (TSR); an appropriate new model is presented. The TSR model makes testable predictions: (a) If crossing to the quintet state occurs, the hydroxylation will be effectively concerted; however, if the process transpires only on the triplet surface, stepwise hydroxylation will occur, and side products derived from radical intermediates would be observed (e.g., loss of stereochemistry). (b) In cases of crossing en route to the quintet transition state, one expects kinetic isotope effects (KIEs) typical of tunneling. (c) In situations where the two surfaces contribute to the rate, one expects intermediate KIEs and radical scrambling patterns that reflect the two processes. (d) Solvent effects on these reactions are expected to be very large.     

Internal electric field in cytochrome C explored by visible electronic circular dichroism spectroscopy

Electronic circular dichroism (ECD) is a valuable tool to explore the secondary and tertiary structure of proteins. With respect to heme proteins, the corresponding visible ECD spectra, which probe the chirality of the heme environment, have been used to explore functionally relevant structural changes in the heme vicinity. While the physical basis of the obtained ECD signal has been analyzed by Woody and co-workers in terms of multiple electronic coupling mechanism between the electronic transitions of the heme chromophore and of the protein (Hsu, M.C.; Woody, R.W. J. Am. Chem. Soc. 1971, 93, 3515), a theory for a detailed quantitative analysis of ECD profiles has only recently been developed (Schweitzer-Stenner, R.; Gorden, J. P.; Hagarman, A. J. Chem. Phys. 2007, 127, 135103). In the present study this theory is applied to analyze the visible ECD-spectra of both oxidation states of three cytochromes c from horse, cow and yeast. The results reveal that both B- and Q-bands are subject to band splitting, which is caused by a combination of electronic and vibronic perturbations. The B-band splittings are substantially larger than the corresponding Q-band splittings in both oxidation states. For the B-bands, the electronic contribution to the band splitting can be assigned to the internal electric field in the heme pocket, whereas the corresponding Q-band splitting is likely to reflect its gradient (Manas, E. S.; Vanderkooi, J. M.; Sharp, K. A. J. Phys. Chem. B 1999, 103, 6344). We found that the electronic and vibronic splitting is substantially larger in the oxidized than in the reduced state. Moreover, these states exhibit different signs of electronic splitting. These findings suggest that the oxidation process increases the internal electric field and changes its orientation with respect to the molecular coordinate system associated with the N-Fe-N lines of the heme group. For the reduced state, we used our data to calculate electric field strengths between 27 and 31 MV/cm for the investigated cytochrome c species. The field of the oxidized state is more difficult to estimate, owing to the lack of information about its orientation in the heme plane. Based on band splitting and the wavenumber of the band position we estimated a field-strength of ca. 40 MV/cm for oxidized horse heart cytochrome c. The thus derived difference between the field strengths of the oxidized and reduced state would contribute at least -55 kJ/mol to the enthalpic stabilization of the oxidized state. Our data indicate that the corresponding stabilization energy of yeast cytochrome c is smaller.     

Metal ion binding by a bicyclic diamide: deep UV Raman spectroscopic characterization

Deep UV resonance Raman spectroscopy was used for characterizing ligand-metal ion complexes. The obtained results demonstrated a strong intrinsic sensitivity and selectivity of a Raman spectroscopic signature of a bicyclic diamide, a novel chelating agent for lanthanides and actinides (Lumetta, G. J.; Rapko, B. M.; Garza, P. A.; Hay, B. P.; Gilbertson, R. D.; Weakley, T. J. R.; Hutchison, J. E. J. Am. Chem. Soc. 2002, 124, 5644). Molecular modeling, which included structure optimization and calculation of Raman frequencies and resonance intensities, allowed for assigning all strong Raman bands of the bicyclic diamide as well as predicting the band shifts observed because of complex formation with metal ions. A comparative analysis of Raman spectra and the results of the molecular modeling could be used for elucidating the structure of complexes in solution.     

Reaction of OH + NO2: high pressure experiments and falloff analysis

High pressure experiments on the OH + NO2 reaction are presented for 3 different temperatures. At 300 K, experiments in He (p = 2-500 bar) as well as in Ar (p = 2-4 bar) were performed. The rate constants obtained in Ar agree well with values which have been reported earlier by our group (Forster, R.; Frost, M.; Fulle, D.; Hamann, H. F.; Hippler, H.; Schlepegrell, A.; Troe, J. J. Chem. Phys. 1995, 103, 2949. Fulle, D.; Hamann, H. F.; Hippler, H.; Troe, J. J. Chem. Phys. 1998, 108, 5391). In contrast, the rate coefficients determined in He were found to be 15-25% lower than the values given in our earlier publications. Additionally, results for He as bath gas at elevated temperatures (T = 400 K, p = 3-150 bar; T = 600 K, p = 3-150 bar) are reported. The results obtained at elevated pressures are found to be in good agreement with existing literature data. The observed falloff behavior is analyzed in terms of the Troe formalism taking into account two reaction channels: one yielding HNO3 and one yielding HOONO. It is found that the extracted parameters are in agreement with rate constants for vibrational relaxation and isotopic scrambling as well as with experimentally determined branching ratios. Based on our analysis we determine falloff parameters to calculate the rate constant for atmospheric conditions.     

The Richest Collection of Tautomeric Polymorphs: The Case of 2‐Thiobarbituric Acid

Five new polymorphs and one hydrated form of 2‐thiobarbituric acid have been isolated and characterised by solid‐state methods. In both the crystalline form II and in the hydrate form, the 2‐thiobarbituric molecules are present in the enol form, whereas only the keto isomer is present in crystalline forms I (reported in 1967 by Calas and Martinex), III , V and VI . In form IV , on the other hand, a 50:50 ordered mixture of enol/keto molecules is present. All new forms have been characterised by single‐crystal X‐ray diffraction, 1D and 2D (¹H, ¹³C, and ¹⁵N) solid‐state NMR spectroscopy, Raman spectroscopy and X‐ray powder diffraction at variable temperature. It has been possible to induce keto–enol conversion between the forms by mechanical methods. The role of hydrogen‐bond interactions in determining the relative stability of the polymorphs and as a driving force in the conversions has been ascertained. To the best of the authors’ knowledge, the 2‐thiobarbituric family of crystal forms represents the richest collection of examples of tautomeric polymorphism so far reported in the literature.     

Fewest-switches surface hopping and decoherence in multiple dimensions

In a recent article (Subotnik, J. E.; Shenvi, N. J. Chem. Phys.2011, 134, 24105), we introduced a new approach for incorporating decoherence into the fewest-switches surface-hopping (FSSH) algorithm, titled augmented FSSH (A-FSSH). The A-FSSH algorithm was designed to correct the well-known overcoherence problem in traditional FSSH, and thus allow wave packets on different surfaces to separate naturally subject to different forces. As presented earlier, however, the A-FSSH algorithm was restricted to two electronic states. We now extend the method to more than two electronic states and present several new model problems with multiple electronic and nuclear dimensions. Lastly, starting with the quantum Liouville equation, we rederive and implement the new phase correction suggested by Shenvi (Shenvi, N.; Subotnik, J. E.; Yang, W. J. Chem. Phys.2011, 135, 24101) and co-workers for propagating the electronic amplitude along a specified nuclear trajectory and find much improved results in certain cases.