Direct measurement of the energy thresholds to conformational isomerization in tryptamine: experiment and theory

The methods of stimulated emission pumping-hole filling spectroscopy (SEP-HFS) and stimulated emission pumping population transfer spectroscopy (SEP-PTS) were applied to the conformation-specific study of conformational isomerization in tryptamine [TRA, 3-(2-aminoethyl)indole]. These experimental methods employ stimulated emission pumping to selectively excite a fraction of the population of a single conformation of TRA to well-defined ground-state vibrational levels. This produces single conformations with well-defined internal energy, tunable over a range of energies from near the zero-point level to well above the lowest barriers to conformational isomerization. When the SEP step overcomes a barrier to isomerization, a fraction of the excited population isomerizes to form that product. By carrying out SEP excitation early in a supersonic expansion, these product molecules are subsequently cooled to their zero-point vibrational levels, where they can be detected downstream with a third tunable laser that probes the ground-state population of a particular product conformer via a unique ultraviolet transition using laser-induced fluorescence. The population transfer spectra (recorded by tuning the SEP dump laser while holding the pump and probe lasers fixed) exhibit sharp onsets that directly determine the energy thresholds for conformational isomerization in a given reactant-product conformer pair. In the absence of tunneling effects, the first observed transition in a given X-Y PTS constitutes an upper bound to the energy barrier to conformational isomerization, while the last transition not observed constitutes a lower bound. The bounds for isomerizing conformer A of tryptamine to B(688-748 cm(-1)), C(1)(860-1000 cm(-1)), C(2)(1219-1316 cm(-1)), D(1219-1282 cm(-1)), E(1219-1316 cm(-1)), and F(688-748 cm(-1)) are determined. In addition, thresholds for isomerizing from B to A(<1562 cm(-1)), B to F(562-688 cm(-1)), and out of C(2) to B(<747 cm(-1)) are also determined. The A-->B and B-->A transitions are used to place bounds on the relative energies of minima B relative to A, with B lying at least 126 cm(-1) above A. The corresponding barriers have been computed using both density functional and second-order many-body perturbation theory methods in order to establish the level of theory needed to reproduce experimental results. While most of the computed barriers match experiment well, the barriers for the A-F and B-F transitions are too high by almost a factor of 2. Possible reasons for this discrepancy are discussed.     

Laser probes of conformational isomerization in flexible molecules and complexes

Molecules with several flexible coordinates have potential energy surfaces with a large number of minima and many transition states separating them. A general experimental protocol is described that is capable of studying conformational isomerization in such circumstances, measuring the product quantum yields following conformation-specific infrared excitation, and measuring energy thresholds for isomerization of specific X --> Y reactant-product isomer pairs following excitation via stimulated emission pumping (SEP). These methods have been applied to a series of molecules of varying size and conformational complexity, including 3-indolepropionic acid (IPA), meta-ethynylstyrene, N-acetyltryptophan methyl amide (NATMA), N-acetyltryptophan amide (NATA), and melatonin. Studies of isomerization in solute-solvent complexes are also described, including a measurement of the barrier to isomerization in the IPA-H2O complex, and a unique isomerization reaction in which a single water molecule is shuttled between H-bonding sites on the trans-formanilide (TFA) molecule.     

Flexing the muscles of m-divinylbenzene: direct measurement of the barriers to conformational isomerization

The energy thresholds to isomerization of the three conformational isomers of m-divinylbenzene (cis-cis, cis-trans, and trans-trans) were directly measured by stimulated emission pumping-population transfer (SEP-PT) spectroscopy. The experimentally determined isomerization thresholds are Ethresh(cc --> ct, tt) = 1080-1232 cm(-1), Ethresh(tt --> ct, cc) = 1130-1175 cm(-1), Ethresh(ct --> cc) = 997-1175 cm(-1), and Ethresh(ct --> tt) = 997-1232 cm(-1). On the basis of the threshold values for X --> Y and Y --> X isomerization, the relative energies of the conformational isomers are -102 < or = E(ct) < or = +178 cm(-1) and -102 < or = E(cc) < or = +95 cm(-1) relative to E(tt) = 0. UV-hole-filling (UVHF) spectroscopy was also used to determine the effect of population returning to the ground state via fluorescence. A full set of governing equations for SEP-PT and UVHF spectroscopy is reported that will be generally useful for future studies using these methods. By comparison of these results with the computed stationary points on a calculated surface (DFT B3LYP/6-31+G*), the isomerization pathway was determined to involve sequential isomerization of each vinyl group rather than concerted motion. The energy thresholds were also combined with the ground state torsional vibrational energy levels to obtain a new fitted two-dimensional torsional potential for mDVB.     

Isomer-specific spectroscopy and conformational isomerization energetics of o-, m-, and p-ethynylstyrenes

The infrared and ultraviolet spectroscopy of o-, m-, and p-ethynylstyrene isomers (oES, mES, and pES) were studied by a combination of methods, including resonance-enhanced two-photon ionization (R2PI), UV-UV hole-burning spectroscopy (UVHB), resonant ion-dip infrared spectroscopy (RIDIRS), and rotationally resolved fluorescence excitation spectroscopy. In addition, the newly developed method of stimulated emission pumping-population transfer spectroscopy (SEP-PTS) was used to determine the energy threshold to conformational isomerization in m-ethynylstyrene. The S(1) <-- S(0) origin transitions of oES and pES occur at 32 369 and 33 407 cm(-1), respectively. In mES, the cis and trans conformations are calculated to be close in energy. In the R2PI spectrum of mES, the two most prominent peaks (32672 and 32926 cm(-1)) were confirmed by UVHB spectroscopy to be S(1) <-- S(0) origins of these two conformers. The red-shifted conformer was identified as the cis structure by least-squares fitting of the rotationally resolved fluorescence excitation spectrum of the origin band. There are also two possible conformations in oES, but transitions due to only one were observed experimentally, as confirmed by UVHB spectroscopy. Density functional theory calculations (B3LYP/6-31+G) predict that the cis-ortho conformer, in which the substituents point toward each other, is about 8 kJ/mol higher in energy than the trans-ortho isomer, and should only be about 5% of the room temperature population of oES. Ground-state infrared spectra in the C-H stretch region (3000-3300 cm(-1)) of each isomer were obtained with RIDIRS. In all three structural isomers, the acetylenic C-H stretch fundamental was split by Fermi resonance. Infrared spectra were also recorded in the excited electronic state, using a UV-IR-UV version of RIDIR spectroscopy. In all three isomers the acetylenic C-H stretch fundamental was unshifted from the ground state, but no Fermi resonance was seen. The first observed and last unobserved transitions in the SEP-PT spectrum were used to place lower and upper bounds on the barrier to cis --> trans isomerization in m-ethynylstyrene of 990-1070 cm(-1). Arguments are given for the lack of a kinetic shift in the measurement. The analogous trans --> cis barrier is in the same range (989-1065 cm(-1)), indicating that the relative energies of the zero-point levels of the two isomers are (E(ZPL)(cis) - E(ZPL)(trans))= -75 to +81 cm(-1). Both the barrier heights and relative energies of the minima are close to those determined by DFT (Becke3LYP/6-31+G) calculations.     

Methods for studying reaction kinetics in gas chromatography, exemplified by using the 1-chloro-2,2-dimethylaziridine interconversion reaction

In this paper, methods are described that are used for studying first-order reaction kinetics by gas chromatography. Basic theory is summarized and illustrated using the interconversion of 1-chloro-2,2-dimethylaziridine enantiomers as a representative example. For the determination of the kinetic and thermodynamic activation data of interconversion the following methods are reviewed: (i) classical kinetic methods where samples of batch-wise kinetic studies are analyzed by enantioselective gas chromatography, (ii) stopped-flow methods performed on one chiral column, (iii) stopped-flow methods performed on an achiral column or empty capillary coupled in series with two chiral columns, (iv) on-flow method performed on an achiral column coupled in series with two chiral columns, and (v) reaction gas chromatography, known as a dynamic gas chromatography, where the interconversion is performed on chiral column during the separation process. The determination of kinetic and thermodynamic activation data by methods (i) through (iv) is straightforward as the experimental data needed for the evaluation (particularly the concentration of reaction constituents) are accessible from the chromatograms. The evaluation of experiments from reaction chromatography method (v) is complex as the concentration bands of reaction constituents are overlapped. The following procedures have been developed to determination peak areas of reaction constituents in such complex chromatograms: (i) methods based on computer-assisted simulations of chromatograms where the kinetic activation parameters for the interconversion of enantiomers are obtained by iterative comparison of experimental and simulated chromatograms, (ii) stochastic methods based on the simulation of Gaussian distribution functions and using a time-dependent probability density function, (iii) approximation function and unified equation, (iv) computer-assisted peak deconvolution methods. Evaluation of the experimental data permits the calculation of apparent rate constants for both the interconversion of the first eluted (k (A-->B)(app)) as well as the second eluted (k(B-->A)(app)) enantiomer. The mean value for all the rate constants (from all the reviewed methods) was found for 1-chloro-2,2-dimethylaziridine A-->B enantiomer interconversion at 100 degrees C: k (A-->B)(app)=21.2 x 10(-4)s(-1) with a standard deviation sigma=10.7 x 10(-4). Evaluating data for reaction chromatography at 100 degrees C {k (app)=k(A-->B)(app)=k(B-->A)(app)=13.9 x 10(-4)s(-1), sigma=3.0 x 10(-4)s(-1)} shows that differences between k(A-->B)(app) and k(B-->A)(app) are the same within experimental error. It was shown both theoretically and experimentally that the Arrhenius activation energy (E(a)) calculated from Arrhenius plots (lnk(app) versus 1/T) is proportional to the enthalpy of activation {E(a)=DeltaH+RT}. Statistical treatment of Gibbs activation energy values gave: DeltaG (app)=110.5kJmol(-1), sigma=2.4kJmol(-1), DeltaG (A-->B)(app)=110.5kJmol(-1), sigma=2.2kJmol(-1), DeltaG (B-->A)(app)=110.3kJmol(-1), sigma=2.8kJmol(-1). This shows that the apparent Gibbs energy barriers for the interconversion of 1-chloro-2,2-dimethylaziridine enantiomers are equal DeltaG (app)=DeltaG(A-->B)(app)=DeltaG(B-->A)(app) and within the given precision of measurement independent of the experimental method used.     

Experimental determination of conformational isomerization energy thresholds in serotonin

Stimulated emission pumping-population transfer (SEP-PT) and hole-filling (SEP-HF) spectroscopies were used to determine the energy thresholds to isomerization between thirteen reactant-product conformer pairs in the biomolecule serotonin (SERO). Serotonin is a close structural analog of tryptamine (TRA), differing in having a hydroxyl group in the 5 position of the indole ring. A previous spectroscopic study (LeGreve; et al. J. Am. Chem. Soc. 2007, 129 (13), 4028) identified eight conformational isomers of SERO, whose interconversion involves motion of the 3-ethylamine side chain, the 5-OH group, or both. In the cases in which only an ethylamine side chain reorientation occurred, the barriers were found to be similar to, but systematically somewhat smaller than, those in TRA, which has been studied by similar methods (Dian; et al. Science 2004, 303 (5661), 1169; Clarkson; et al. J. Chem. Phys. 2005, 122 (21), Art. No. 214311). In most cases, the experimental thresholds are well reproduced by calculated transition states separating the conformational wells; however, tunneling effects may artificially reduce the thresholds observed for isomerization of SERO(A,Gpy(out)) and SERO(B,Gpy(up)) into SERO(C,Gph(out)). The A --> A' isomerization involving only the OH rotation from anti to syn was found to be 721-761 cm-1, in accordance with the calculated classical barrier. For isomerizations in which the ethylamine side chain reorients as does the OH group, the barriers to isomerization were consistent with sequential rather than concerted motion of both groups. Finally, some evidence for mode-specific effects in the product quantum yields near threshold is presented.     

The singlet electronic ground state isomers of dialuminum monoxide: AlOAl, AlAlO, and the transition state connecting them

The singlet electronic ground state isomers, X (1)Sigma(g) (+) (AlOAl D(infinityh)) and X (1)Sigma(+) (AlAlO C(infinitynu)), of dialuminum monoxide have been systematically investigated using ab initio electronic structure theory. The equilibrium structures and physical properties for the two molecules have been predicted employing self-consistent field (SCF) configuration interaction with single and double excitations (CISD), multireference CISD (MRCISD), coupled cluster with single and double excitations (CCSD), CCSD with perturbative triples [CCSD(T)], CCSD with iterative partial triple excitations (CCSDT-3 and CC3), and full triples (CCSDT) coupled cluster methods. Four correlation consistent polarized valence (cc-pVXZ) type basis sets were used. The AlAlO system is rather challenging theoretically. The two isomers are confirmed to have linear structures at all levels of theory. The symmetric isomer AlOAl is predicted to lie 81.9 kcal mol(-1) below the asymmetric isomer AlAlO at the cc-pV(Q+d)Z CCSD(T) level of theory. The predicted harmonic vibrational frequencies for the X (1)Sigma(g) (+) AlOAl molecule, omega(1)=517 cm(-1), omega(2)=95 cm(-1), and omega(3)=1014 cm(-1), are in good agreement with experimental values. The harmonic vibrational frequencies for the X (1)Sigma(+) AlAlO structure, omega(1)=1042 cm(-1), omega(2)=73 cm(-1), and omega(3)=253 cm(-1), presently have no experimental values with which to be compared. With the same methods the barrier heights for the isomerization AlOAl-->AlAlO and AlAlO-->AlOAl reactions were predicted to be 84.3 and 2.4 kcal mol(-1), respectively. The dissociation energies D(0) for AlOAl (X (1)Sigma(g) (+)) and AlAlO (X (1)Sigma(+))-->AlO (X (2)Sigma(+))+Al ((2)P) were determined to be 130.8 and 48.9 kcal mol(-1), respectively. Thus, both symmetric AlOAl (X (1)Sigma(g) (+)) and asymmetric AlAlO (X (1)Sigma(+)) isomers are expected to be thermodynamically stable with respect to the dissociation into AlO (X (2)Sigma(+)) + Al ((2)P) and kinetically stable for the isomerization reaction (AlAlO-->AlOAl) at sufficiently low temperatures.     

Prediction of FCC gasoline octane numbers using FT-MIR and PLS

A method for predicting "octane numbers" (RON and MON) in fluid catalytic cracking (FCC) gasolines is proposed. Using FT-MIR and PLS, improvements have been obtained in sample throughput, reduced delay times, accuracy (repeatability and reproducibility), amounts of samples and reagents and environmental working conditions when compared with current standard methods. A total number of 140 daily production samples were taken; and from there, a learning group was prepared (44 samples); a validation set (96 samples) was prepared, as well. Sample spectra were recorded from 4000 to 600 cm(-1) at 4 cm(-1) intervals (traditional sealed NaCl cells). The PLS technique was used in its two variants (1 and 2-block). Both provided similar results. Their predictive characteristics are very good: SEP(RON)=0.38; SEP(MON)=0.40; repeatability <0.1 O.N.; reproducibility <0.3 O.N. (SEP=Standard Error of Prediction).     

Quantum energy flow and the kinetics of water shuttling between hydrogen bonding sites on trans-formanilide

A potential energy surface for trans-formanilide (TFA)-H2O is calculated and applied to study energy flow in the complex as well as the kinetics of water shuttling between hydrogen bonding sites on TFA. In addition to the previously identified H2O-TFA(C[Double Bond]O) and H2O-TFA(NH) minima, with the water monomer bound to the C[Double Bond]O and NH groups, respectively, the new surface reveals a second local minimum with the water bound to the C[Double Bond]O group, and which lies energetically 310 cm(-1) above the previously identified H2O-TFA(C[Double Bond]O) global minimum. On this surface, the energy barrier for water shuttling from H2O-TFA(C[Double Bond]O) global minimum to H2O-TFA(N-H) is 984 cm(-1), consistent with the experimental bounds of 796 and 988 cm(-1) [J. R. Clarkson et al. Science 307, 1443 (2005)]. The ergodicity threshold of TFA is calculated to be 1450 cm(-1); for the TFA-H2O complex, the coupling to the water molecule is found to lower the ergodicity threshold to below the isomerization barrier. Energy transfer between the activated complex and the vibrational modes of TFA is calculated to be sufficiently rapid that the Rice-Ramsperger-Kassel-Marcus (RRKM) theory does not overestimate the rate of water shuttling. The possibility that the rate constant for water shuttling is higher than the RRKM theory estimate is discussed in light of the relatively high energy of the ergodicity threshold calculated for TFA.     

B4O分子异构化稳定性的理论研究

采用CCSD(T)/6-311+G(2df)//B3LYP/6-311+G(d)方法, 系统研究了B₄O体系各个异构体的结构和能量, 以及重要异构体的解离和异构化稳定性. 结果表明, 单态平面三角形含-BO单元的异构体cB₃-BO ¹1能量最低, 其次是带状的异构体B₄O ¹2(10.9 kJ/mol), 并且¹1和¹2结构都具有良好的动力学稳定性. 因此对于B₄O体系, ¹1和¹2都是有可能存在的. 而文献报道的三态直线型结构BBBBO(146.0 kJ/mol)的能量比异构体¹1和¹2高得多.     

Phase-space overlap measures. I. Fail-safe bias detection in free energies calculated by molecular simulation

We consider ways to quantify the overlap of the parts of phase space important to two systems, labeled A and B. Of interest is how much of the A-important phase space lies in that important to B, and how much of B lies in A. Two measures are proposed. The first considers four total-energy distributions, formed from all combinations made by tabulating either the A-system or the B-system energy when sampling either the A or B system. Measures for A in B and B in A are given by two overlap integrals defined on pairs of these distributions. The second measure is based on information theory, and defines two relative entropies which are conveniently expressed in terms of the dissipated work for free-energy perturbation (FEP) calculations in the A-->B and B-->A directions, respectively. Phase-space overlap is an important consideration in the performance of free-energy calculations. To demonstrate this connection, we examine bias in FEP calculations applied to a system of independent particles in a harmonic potential. Systems are selected to represent a range of overlap situations, including extreme subset, subset, partial overlap, and nonoverlap. The magnitude and symmetry of the bias (A-->B vs B-->A) are shown to correlate well with the overlap, and consequently with the overlap measures. The relative entropies are used to scale the amount of sampling to obtain a universal bias curve. This result leads to develop a simple heuristic that can be applied to determine whether a work-based free-energy measurement is free of bias. The heuristic is based in part on the measured free energy, but we argue that it is fail-safe inasmuch as any bias in the measurement will not promote a false indication of accuracy.     

Electronic spectroscopy and photoisomerization of trans-urocanic acid in a supersonic jet.

trans-Urocanic acid (trans-UA), a component of the epidermal layer of skin, exhibits wavelength-dependent photochemistry. The quantum efficiency of isomerization to cis-UA is greatest when the molecule is excited on the long wavelength tail of its absorption profile in solution (300-320 nm). However, exciting the molecule where it absorbs UV light most efficiently (260-285 nm) causes almost no isomerization. We have used fluorescence excitation and dispersed emission methods in a supersonic jet to investigate the electronic states involved in this complex and interesting photochemistry. Three distinct regions are present in the excitation spectrum. Region I, which is below the isomerization barrier, contains sharp, well-resolved peaks that upon excitation emit from the S(1) state of trans-UA. Region II exhibits peaks that increase in broadness and decrease in intensity with increasing excitation energy. Upon excitation these peaks produce dual emission from the S(1) states of both trans- and cis-UA. The trans to cis isomerization barrier is estimated to be 1400 cm(-1). Region III exhibits excitation to the S(2) electronic state and has a broad structure that spans 3000 cm(-1) and occurs 4000 cm(-1) above S(1). S(2) excitation results in essentially no trans to cis isomerization.     

Vibrations in the B4 rhombic structure

A double minimum six-dimensional potential energy surface (PES) is determined in symmetry coordinates for the most stable rhombic (D2h) B4 isomer in its 1Ag electronic ground state by fitting to energies calculated ab initio. The PES exhibits a barrier to the D4h square structure of 255 cm(-1). The vibrational levels (J=0) are calculated variationally using an approach which involves the Watson kinetic energy operator expressed in normal coordinates. The pattern of about 65 vibrational levels up to 1600 cm(-1) for all stable isotopomers is analyzed. Analogous to the inversion in ammonia-like molecules, the rhombus rearrangements lead to splittings of the vibrational levels. In B4 it is the B1g (D4h) mode which distorts the square molecule to its planar rhombic form. The anharmonic fundamental vibrational transitions of 11B4 are calculated to be (splittings in parentheses): G(0)=2352(22) cm(-1), nu1(A1g)=1136(24) cm(-1), nu2(B1g)=209(144) cm(-1), nu3(B2g)=1198(19) cm(-1), nu4(B2u)=271(24) cm(-1), and nu5(Eu)=1030(166) cm(-1) (D4h notation). Their variations in all stable isotopomers were investigated. Due to the presence of strong anharmonic resonances between the B1g in-plane distortion and the B2u out-of-plane bending modes, the higher overtones and combination levels are difficult to assign unequivocally.     

Laser spectroscopic study on the conformations and the hydrated structures of benzo-18-crown-6-ether and dibenzo-18-crown-6-ether in supersonic jets

The laser-induced fluorescence spectra of jet-cooled benzo-18-crown-6 (B18C6) and dibenzo-18-crown-6 (DB18C6) exhibit a number of vibronic bands in the 35 000-37 000 cm(-1) region. We attribute these bands to monomers and hydrated clusters by fluorescence-detected IR-UV and UV-UV double resonance spectroscopy. We found four and two conformers for bare B18C6 and DB18C6, and the hydration of one water molecule reduces the number of isomers to three and one for B18C6-(H(2)O)(1) and DB18C6-(H(2)O)(1), respectively. The IR-UV spectra of B18C6-(H(2)O)(1) and DB18C6-(H(2)O)(1) suggest that all isomers of the monohydrated clusters have a double proton-donor type (bidentate) hydration. That is, the water molecule is bonded to B18C6 or DB18C6 via two O-H[dot dot dot]O hydrogen bonds. The blue shift of the electronic origin of the monohydrated clusters and the quantum chemical calculation suggest that the water molecule in B18C6-(H(2)O)(1) and DB18C6-(H(2)O)(1) prefers to be bonded to the ether oxygen atoms near the benzene ring.     

Photophysical studies of the trans to cis isomerization of the push-pull molecule: 1-(pyridin-4-yl)-2-(N-methylpyrrol-2-yl)ethene (mepepy)

Organic molecules possessing intramolecular charge-transfer properties (D-pi-A type molecules) are of key interest particularly in the development of new optoelectronic materials as well as photoinduced magnetism. One such class of D-pi-A molecules that is of particular interest contains photoswitchable intramolecular charge-transfer states via a photoisomerizable pi-system linking the donor and acceptor groups. Here we report the photophysical and electronic properties of the trans to cis isomerization of 1-(pyridin-4-yl)-2-(N-methylpyrrol-2-yl)ethene ligand (mepepy) in aqueous solution using photoacoustic calorimetry (PAC) and theoretical methods. Density functional theory (DFT) calculations demonstrate a global energy difference between cis and trans isomers of mepepy to be 8 kcal mol(-1), while a slightly lower energy is observed between the local minima for the trans and cis isomers (7 kcal mol(-1)). Interestingly, the trans isomer appears to exhibit two ground-state minima separated by an energy barrier of approximately 9 kcal mol(-1). Results from the PAC studies indicate that the trans to cis isomerization results in a negligible volume change (0.9 +/- 0.4 mL mol(-1)) and an enthalpy change of 18 +/- 3 kcal mol(-1). The fact that the acoustic waves associated with the trans to cis transition of mepepy overlap in frequency with those of a calorimetric reference implies that the conformational transition occurs faster than the approximately 50 ns response time of the acoustic detector. Comparison of the experimental results with theoretical studies provide evidence for a mechanism in which the trans to cis isomerization of mepepy results in the loss of a hydrogen bond between a water molecule and the pyridine ring of mepepy.     

A theoretical survey on the structures, energetics, and isomerization pathways of the B(5)O radical

Boron-centered radicals have received growing interest. Recently, two groups reported density functional theory investigations (GGA-PW91 and B3LYP) on a hexa-atomic boron-oxide radical, B(5)O, which has led to great discrepancies on the type of low-lying structures. In this article, we not only explore the energetics of doublet and quartet B(5)O isomers at high electron-correlated levels (CCSD(T)/6-311+G(2df), CCSD(T)/aug-cc-pVTZ, and G3B3) but also investigate the isomerization and fragmentation stability of the low-lying B(5)O isomers. All the high-level studies consistently show that the B(5)O radical possesses a belt-like ground structure (2)01 in doublet electronic state followed by isomer (2)02 with an exocyclic - BO moiety at around 3.0 kcal/mol. Kinetically, (2)01 and (2)02 are separated by a considerable barrier of about 20 kcal/mol. Thus, the two isomeric forms of B(5)O radical should be very promising for isolation in laboratory. However, the other four isomers reported recently are all kinetically unstable toward conversion to (2)01 and (2)02. The high thermodynamic and kinetic stability of (2)01 and (2)02 might make them as important building cores in the growth of boron-oxide clusters. This results would also help deeply understand the oxidation and doping mechanism of pure boron clusters.     

Nonradiative decay dynamics of methyl-4-hydroxycinnamate and its hydrated complex revealed by picosecond pump-probe spectroscopy

The lifetimes of methyl 4-hydroxycinnamate (OMpCA) and its mono-hydrated complex (OMpCA-H(2)O) in the S(1) state have been measured by picosecond pump-probe spectroscopy in a supersonic beam. For OMpCA, the lifetime of the S(1)-S(0) origin is 8-9 ps. On the other hand, the lifetime of the OMpCA-H(2)O complex at the origin is 930 ps, which is ～100 times longer than that of OMpCA. Furthermore, in the complex the S(1) lifetime shows rapid decrease at an energy of ～200 cm(-1) above the origin and finally becomes as short as 9 ps at ～500 cm(-1). Theoretical calculations with a symmetry-adapted cluster-configuration interaction (SAC-CI) method suggest that the observed lifetime behavior of the two species is described by nonradiative decay dynamics involving trans → cis isomerization. That is both OMpCA and OMpCA-H(2)O in the S(1) state decay due to the trans → cis isomerization, and the large difference of the lifetimes between them is due to the difference of the isomerization potential energy curve. In OMpCA, the trans → cis isomerization occurs smoothly without a barrier on the S(1) surface, while in the OMpCA-H(2)O complex, there exists a barrier along the isomerization coordinate. The calculated barrier height of OMpCA-H(2)O is in good agreement with that observed experimentally.     

First spectroscopic studies of the B2Sigma+ --> X2Sigma+ system in the 12C17O+ molecule spectrum

The emission spectrum of the B2Sigma+ --> X2Sigma+ system in the 12C17O+ has been recorded at high resolution with conventional spectroscopic techniques. A graphite hollow cathode lamp filled with oxygen, enriched approximately 45% with the 17O2 isotope, has given strong spectra of hitherto unreported bands system of 12C17O+. Four bands, which form the 0-v(') (v(')=1 -4) progression, have been rotationally analysed with effective Hamiltonians of Brown [J. Mol. Spectrosc. 74 (1979) 294]. Next the rovibronic structure parameters, for the B2Sigma+ and X2Sigma+ states, have been determined from merge calculations. Finally, principal equilibrium molecular constants for the ground state: omegae=2186.032(26) cm(-1), omegae xe=14.7848(43) cm(-1), Be=1.927271(73) cm(-1) , alphae=1.8432(24) x 10(-2) cm(-1), De=6.018(26) x 10(-6) cm(-1), betae=1.11(92) x 10(-8) cm(-1) and the sigma(e)(B-->X)=45876.563(34) cm(-1) value of the 12CO17+ molecule have been derived for the first time in this investigation.     

Gas-phase infrared and ab initio study of the unstable CF3CNO molecule and its stable furoxan ring dimer

The unstable trifluoroacetonitrile N-oxide molecule, CF3CNO, has been generated in high yield in the gas phase from CF3BrC=NOH and studied for the first time by gas-phase mid-infrared spectroscopy. Cold trapping of this molecule followed by slow warming forms the stable ring dimer, bis(trifluoromethyl)furoxan, also investigated by gas-phase infrared spectroscopy. The spectroscopy provides an investigation into the vibrational character of the two molecules, the assignments supported by calculations of the harmonic vibrational frequencies using in the case of CF3CNO both ab initio (CCSD(T)) and density functional theory (B3LYP) and B3LYP for the ring dimer. The ground-state structures of both molecules were investigated at the B3LYP level of theory, with CF3CNO further investigated using coupled-cluster. The CCSD(T) method suggests a slightly bent (C(s)) structure for CF3CNO, while the B3LYP method (with basis sets ranging from 6-311G(d) to cc-pVTZ) suggests a close-to-linear or linear CCNO chain. The CCN bending potential in CF3CNO was explored at the CCSD(T)(fc)/cc-pVTZ level, with the results suggesting that CF3CNO exhibits strong quasi-symmetric top behavior with a barrier to linearity of 174 cm(-1). Since both isomerization and dimerization are feasible loss processes for this unstable molecule, the relative stability of CF3CNO with respect to the known cyanate (CF3OCN), isocyanate (CF3NCO), and fulminate (CF3ONC) isomers and the mechanism of the dimerization process to the ring furoxan and other isomers were studied with density functional theory.     

TD-DFT calculations of the potential energy curves for the trans-cis photo-isomerization of protonated Schiff base of retinal

One-dimensional potential energy curves for the isomerization of protonated Schiff base of retinal (PSBR) in bacteriorhodopsin (bR), i.e., isomerization from all-trans- to 13-cis-forms, have been calculated by means of time-dependent density functional theory (TD-DFT) calculations, in order to elucidate the mechanism of initial step in photo-absorption. The transition state of the isomerization in the first excited state is located at theta(13-14)=58 degrees , where theta(13-14) means twist angle around the C(13)=C(14) double bond of PSBR The potential barrier is formed by the avoided crossing between S(1) (B(u)-like) and S(2) (A(g)-like) states. The mechanism of the isomerization was discussed on the basis of theoretical results.