A variational method for the calculation of spin-rovibronic energy levels of any triatomic molecule in an electronic triplet state

The triatomic spin-rovibronic variational code RVIB3 has been extended to include the effect of two uncoupled electrons, for both (3)Sigma(-) and (3)Pi (Renner-Teller) electronic states. The spin-orbital-rotational kinetic energy is included in the usual way, via terms (J+L+S). The phenomenological terms AL.S and lambda 2/3(3S(z)(2)) are introduced to reproduce the 3 spin-orbit and spin-spin splittings, respectively. Calculations are performed to evaluate the spin-rovibronic energy levels of CCO (X) over tilde (3) Sigma(-) and CCO (A) over tilde (3) Pi for which the Born-Oppenheimer potentials are derived from high-accuracy ab initio calculations.

Calculations of rovibrational energies and dipole transition intensities for polyatomic molecules using MULTIMODE

We report rigorous calculations of rovibrational energies and dipole transition intensities for three molecules using a new version of the code MULTIMODE. The key features of this code which permit, for the first time, such calculations for moderately sized but otherwise general polyatomic molecules are briefly described. Calculations for the triatomic molecule BF(2) are done to validate the code. New calculations for H(2)CO and H(2)CS are reported; these make use of semiempirical potentials but ab initio dipole moment surfaces. The new dipole surface for H(2)CO is a full-dimensional fit to the dipole moment obtained with the coupled-cluster with single and double excitations and a perturbative treatment of triple excitations method with the augmented correlation consistent triple zeta basis set. Detailed comparisons are made with experimental results from a fit to relative data for H(2)CS and absolute intensities from the HITRAN database for H(2)CO.

Intramolecular proton transfer in [Ph2B-OH2](+). Participation of a phenyl group

Ab initio calculations at the HF/6-31+G* level on [Ph2B-OH2](+) show that in the gas phase the structure with the proton attached to an ipso C is lower in energy than the one with the proton on the oxygen atom by 8.40 kcal mol(-1). The transition states and reaction paths for intramolecular proton transfer in [Ph2B-OH2](+) have also been studied.

The gas-phase reactions Of SiCl4 and Si2Cl6 with CH3OH and C2H5OH: an investigation by mass spectrometry and matrix-isolation infrared spectroscopy

The gas phase reactions Of SiCl4 and Si2Cl6 With CH3OH and C2H5OH have been investigated using both mass spectrometry and matrix isolation techniques. SiCl4 reacts with both CH3OH and C2H5OH upon mixing of the vapours for times in excess of 3 h to generate the HCl-elimination products SiCl3OR (R = CH3 or C2H5). The identity of these products is confirmed by deuteration experiments and by ab initio calculations at the HF/6-31G(d) level. Further products are generated when the mixture is passed through a tube heated to 750degreesC. Si2Cl6 reacts with CH3OH and C2H5OH via a different mechanism in which the Si-Si bond is cleaved to yield SiCl3OR and HCl. Other products of the type SiCl4-n(OCH3)(n) are tentatively identified by a combination of mass spectrometric and matrix isolation measurements. These latter products indicate further replacement of Cl atoms by OR groups as a result of reaction of CH3OH or C2H5OH with the initial product.

Reactions of SiCl2 with N2O, NO and O-2

The thermal route to dichlorosilylene by pyrolysis Of Si2Cl6 has been investigated using both mass spectrometry and matrix isolation techniques. The formation Of SiCl2 in the gas phase was confirmed by employing a known "trapping" agent, namely buta-1,3-diene, which gave the product 1, 2-dichloro-1-silacyclopent-3 -ene. Dichlorosilylene was then reacted with N2O and NO. The observed products in the mass spectrum from the N2O reaction were SiCl2O and its polymers and N-2. On reacting SiCl2 with NO, SiCl2O and its polymers, Cl-2 and N2O were all observed. Infrared spectra of argon matrices supported these findings from mass spectrometry. A mechanism is proposed for this reaction based on these observations involving the intermediacy of cyclo-Cl2SiO2 and is supported by ab initio calculations at the MP2 and G3 levels. The reaction between SiCl2 and O-2 has also been investigated. The products seen in this case were SiCl2O and Cl-2. Ab initio calculations again suggest that cyclo-Cl2SiO2 is involved, and a chain mechanism seems the most likely route to Cl-2 formation. The calculations lead to DeltaH(f)degrees (SiO2,g) = -276 +/- 4- 6 kJ mol(-1).

Rovibrational energy levels for the electronic ground state of A1OH

The vibrational-rotational energy levels of aluminum monohydroxide in its electronic ground state, (A) over tilde (1)A' AlOH, have been predicted using the variational method. The potential energy surface of the (X) over tilde (1)A' ground state of AlOH was determined employing the ab initio coupled cluster method with single, double, and perturbative triple excitations [CCSD(T)] and the correlation-consistent polarized valence quadruple zeta (cc-pVQZ) basis set. Low-lying J= 0 and J= 1 vibrational levels are reported. These are analyzed in terms of the quasilinearity of the molecule. Coriolis effects are shown to be significant. We hope that our predictions will be of value in the future when assigning rovibrational transitions in spectroscopic studies. (c) 2006 Elsevier B.V. All rights reserved.

The chemistry of intrinsically chiral surfaces

Intrinsically chiral metal and mineral surfaces show enantioselective behaviour without modifiers. Examples are artificial high-Miller-index surfaces of metal single crystals with cubic bulk lattice symmetry, which have no mirror planes and are therefore chiral, or surfaces of naturally occurring crystallites of some common minerals, such as alpha-quartz or calcite. Recent findings with regards to the surface geometry, reactivity and thermal stability of intrinsically chiral surfaces are discussed. A number of enantioselective effects have been reported in connection with the adsorption of small chiral molecules (e.g. alanine, cysteine) on intrinsically chiral surfaces under well-defined conditions. From a combination of experimental surface science techniques and theoretical ab initio model calculations it emerges that these effects are due to a combination of attractive and repulsive adsorbate-substrate and inter-adsorbate interactions.

Vibrations in the B-4 rhombic structure

A double minimum six-dimensional Potential energy surface (PES) is determined in symmetry coordinates for the most stable rhombic (D-2h) B-4 isomer in its (1)A(g) electronic ground state by fitting to energies calculated ab initio. The PES exhibits a barrier to the D-4h 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 B-4 it is the B-1g (D-4h mode which distorts the square molecule to its planar rhombic form. The anharmonic fundamental vibrational transitions of B-11(4) are calculated to be (splittings in parentheses): G(O) = 2352(22) cm(-1), v(1)(A(1g)) - 1136(24) cm(-1,) v(2)(B-1g)=209(144) cm(-1) v(3)(B-2g)=1198(19)cm(-1), v(4)(B-2u) = 271(24) cm(-1), and v(5) (E-u) = 1030( 166) cm(-1) (D-4h notation). Their variations in all stable isotoporners were investigated. Due to the presence of strong anharmonic resonances between the B-1g in-plane distortion and the B-2u, out-of-plane bending modes. the hiaher overtones and combination levels are difficult to assign unequivocally. (C) 2005 American Institute of Physics.

Potential energy surface and MULTIMODE vibrational analysis of C2H3+

A full dimensional, ab initio-based semiglobal potential energy surface for C2H3+ is reported. The ab initio electronic energies for this molecule are calculated using the spin-restricted, coupled cluster method restricted to single and double excitations with triples corrections [RCCSD(T)]. The RCCSD(T) method is used with the correlation-consistent polarized valence triple-zeta basis augmented with diffuse functions (aug-cc-pVTZ). The ab initio potential energy surface is represented by a many-body (cluster) expansion, each term of which uses functions that are fully invariant under permutations of like nuclei. The fitted potential energy surface is validated by comparing normal mode frequencies at the global minimum and secondary minimum with previous and new direct ab initio frequencies. The potential surface is used in vibrational analysis using the "single-reference" and "reaction-path" versions of the code MULTIMODE. (c) 2006 American Institute of Physics.

C-H...O hydrogen bonding in 4-phenyl-benzaldehyde: a comprehensive crystallographic, spectroscopic and computational study

The crystal structure of 4-phenyl-benzaldehyde reveals the presence of a dimer linked by the C=O and C( 9)-H groups of adjacent molecules. In the liquid phase, the presence of C-(HO)-O-... bonded forms is revealed by both vibrational and NMR spectroscopy. A Delta H value of - 8.2 +/- 0.5 kJ mol(-1) for the dimerisation equilibrium is established from the temperature-dependent intensities of the bands assigned to the carbonyl-stretching modes. The NMR data suggest the preferential engagement of the C(2,6)-H and C(10/12)/C(11)-H groups as hydrogen bond donors, instead of the C(9)-H group. While ab initio calculations for the isolated dimers are unable to corroborate these NMR results, the radial distribution functions obtained from molecular dynamics simulations show a preference for C(2,6)-H and C(10/12)/C(11)-(HO)-O-... contacts relative to the C(9)-(HO)-O-... ones.

Tests of MULTIMODE calculations of rovibrational energies of CH4

We report variational calculations of rovibrational energies of CH4 using the code MULTIMODE and an ab initio force field of Schwenke and Partridge. The systematic convergence of the energies with respect to the level of mode coupling is presented. Converged vibrational energies calculated using the five-mode representation of the potential for zero total angular momentum are compared with previous, benchmark calculations based on Radau coordinates using this force field for zero total angular momentum and for J = 1. Very good agreement with the previous benchmark calculations is found. (c) 2006 Elsevier B.V. All rights reserved.

Studying reactive processes with classical dynamics: rebinding dynamics in MbNO

A new surface-crossing algorithm suitable for describing bond-breaking and bond-forming processes in molecular dynamics simulations is presented. The method is formulated for two intersecting potential energy manifolds which dissociate to different adiabatic states. During simulations, crossings are detected by monitoring an energy criterion. If fulfilled, the two manifolds are mixed over a finite number of time steps, after which the system is propagated on the second adiabat and the crossing is carried out with probability one. The algorithm is extensively tested (almost 0.5 mu s of total simulation time) for the rebinding of NO to myoglobin. The unbound surface ((FeNO)-N-...) is represented using a standard force field, whereas the bound surface (Fe-NO) is described by an ab initio potential energy surface. The rebinding is found to be nonexponential in time, in agreement with experimental studies, and can be described using two time constants. Depending on the asymptotic energy separation between the manifolds, the short rebinding timescale is between 1 and 9 ps, whereas the longer timescale is about an order of magnitude larger. NO molecules which do not rebind within 1 ns are typically found in the Xenon-4 pocket, indicating the high affinity of NO to this region in the protein.

Potential energy surface and moleculardynamics of MbNO: existence of an unsuspected FeON minimum

Ligands such as CO, O2, or NO are involved in the biological function of myoglobin. Here we investigate the energetics and dynamics of NO interacting with the Fe(II) heme group in native myoglobin using ab initio and molecular dynamics simulations. At the global minimum of the ab initio potential energy surface (PES), the binding energy of 23.4 kcal/mol and the Fe-NO structure compare well with the experimental results. Interestingly, the PES is found to exhibit two minima: There exists a metastable, linear Fe-O-N minimum in addition to the known, bent Fe-N-O global minimum conformation. Moreover, the T-shaped configuration is found to be a saddle point, in contrast to the corresponding minimum for NO interacting with Fe(III). To use the ab initio results for finite temperature molecular dynamics simulations, an analytical function was fitted to represent the Fe-NO interaction. The simulations show that the secondary minimum is dynamically stable up to 250 K and has a lifetime of several hundred picoseconds at 300 K. The difference in the topology of the heme-NO PES from that assumed previously (one deep, single Fe-NO minimum) suggests that it is important to use the full PES for a quantitative understanding of this system. Why the metastable state has not been observed in the many spectroscopic studies of myoglobin interacting with NO is discussed, and possible approaches to finding it are outlined.

Computational chemistry for elucidating protein function: energetics and dynamics of myoglobin-ligand systems

State-of-the-art computational methodologies are used to investigate the energetics and dynamics of photodissociated CO and NO in myoglobin (Mb···CO and Mb···NO). This includes the combination of molecular dynamics, ab initio MD, free energy sampling, and effective dynamics methods to compare the results with studies using X-ray crystallography and ultrafast spectroscopy metho ds. It is shown that modern simulation techniques along with careful description of the intermolecular interactions can give quantitative agreement with experiments on complex molecular systems. Based on this agreement predictions for as yet uncharacterized species can be made.

Theoretical investigation of infrared spectra and pocket dynamics of photodissociated carbonmonoxy myoglobin

Molecular dynamics simulations of the photodissociated state of carbonmonoxy myoglobin (MbCO) are presented using a fluctuating charge model for CO. A new three-point charge model is fitted to high-level ab initio calculations of the dipole and quadrupole moment functions taken from the literature. The infrared spectrum of the CO molecule in the heme pocket is calculated using the dipole moment time autocorrelation function and shows good agreement with experiment. In particular, the new model reproduces the experimentally observed splitting of the CO absorption spectrum. The splitting of 3–7 cm−1 (compared to the experimental value of 10 cm−1) can be directly attributed to the two possible orientations of CO within the docking site at the edge of the distal heme pocket (the B states), as previously suggested on the basis of experimental femtosecond time-resolved infrared studies. Further information on the time evolution of the position and orientation of the CO molecule is obtained and analyzed. The calculated difference in the free energy between the two possible orientations (Fe···CO and Fe···OC) is 0.3 kcal mol−1 and agrees well with the experimentally estimated value of 0.29 kcal mol−1. A comparison of the new fluctuating charge model with an established fixed charge model reveals some differences that may be critical for the correct prediction of the infrared spectrum and energy barriers. The photodissociation of CO from the myoglobin mutant L29F using the new model shows rapid escape of CO from the distal heme pocket, in good agreement with recent experimental data. The effect of the protein environment on the multipole moments of the CO ligand is investigated and taken into account in a refined model. Molecular dynamics simulations with this refined model are in agreement with the calculations based on the gas-phase model. However, it is demonstrated that even small changes in the electrostatics of CO alter the details of the dynamics.