Linkage disequilibrium assessment via log-linear modeling of SNP haplotype frequencies

Analyses of high-density single-nucleotide polymorphism (SNP) data, such as genetic mapping and linkage disequilibrium (LD) studies, require phase-known haplotypes to allow for the correlation between tightly linked loci. However, current SNP genotyping technology cannot determine phase, which must be inferred statistically. In this paper, we present a new Bayesian Markov chain Monte Carlo (MCMC) algorithm for population haplotype frequency estimation, particulary in the context of LD assessment. The novel feature of the method is the incorporation of a log-linear prior model for population haplotype frequencies. We present simulations to suggest that 1) the log-linear prior model is more appropriate than the standard coalescent process in the presence of recombination (>0.02cM between adjacent loci), and 2) there is substantial inflation in measures of LD obtained by a "two-stage" approach to the analysis by treating the "best" haplotype configuration as correct, without regard to uncertainty in the recombination process. Genet Epidemiol 25:106-114, 2003. (C) 2003 Wiley-Liss, Inc.

Molybdenum eta(3)-allyl dicarbonyl complexes as a new class of precursors for highly reactive epoxidation catalysts with tert-butyl hydroperoxide

New Mo(II) diimine derivatives of [Mo(q (3)allyl)X(CO)(2)(CH3CN)(2)] (allyl = C3H5 and C5H5O; X = Cl, Br) were prepared, and [MO(eta(3)-C3H5)Cl(CO)(2)(BIAN)] (BIAN = 1,4-(4-chloro)phenyl-2,3-naphthalene-diazabutadiene) (7) was structurally characterized by single-crystal X-ray diffraction. This complex adopted an equatorial-axial arrangement of the bidentate ligand (axial isomer), in contrast with the precursors, found as the equatorial isomer in the solid and fluxional in solution. The new complexes of the type [Mo(eta(3)-allyl)X(CO)(2)(N-N)l (N-N is a bidentate chelating dinitrogen ligand) were tested for the catalytic epoxidation of cyclooctene using tert-butyl hydroperoxide as oxidant. All catalytic systems were 100% selective toward epoxide formation. While their turnover frequencies paralleled those of related Mo(eta) carbonyl compounds or Mo(VI) compounds bearing similar N-donor ligands, they exhibited similar olefin conversions in consecutive catalytic runs. The acetonitrile precursors were generally more active than the diimine complexes, and the chloro derivatives more active than the bromo ones. Combined vibrational and NMR spectroscopy and computational studies (DFT) were used to investigate the nature of the molybdenum species formed in the catalytic system with [Mo(eta(3)-C3H5)Cl(CO)(2){1,4-(2,6-dimethyl)phenyl-2.3-dimethyldiazabuta diene}] (4) and to propose that the resulting species may be dimeric bearing oxide bridges.

The photodimerisation of the alpha- and beta-forms of trans-cinnamic acid: a study of single crystals by vibrational microspectroscopy

Infrared and Raman microspectroscopy have been used to follow the photodimerisation reactions of single crystals, the alpha- and beta-forms of trans-cinnamic acid. This approach allows the starting materials and products -alpha-truxillic acid that has C-i symmetry and beta-truxinic acid, which has C-s symmetry-to be identified. It also allows the topotactic nature of the reaction to be confirmed. Attempts to produce the poorly-defined unreactive gamma-form of trans-cinnamic acid resulted only in a mixture of the alpha- and beta-forms. The findings suggest a wide role for these spectroscopic methods in monitoring solid-state organic reactions. (C) 2002 Elsevier Science B.V. All rights reserved.

Calculation of converged rovibrational energies and partition function for methane using vibrational-rotational configuration interaction

The rovibration partition function of CH4 was calculated in the temperature range of 100-1000 K using well-converged energy levels that were calculated by vibrational-rotational configuration interaction using the Watson Hamiltonian for total angular momenta J=0-50 and the MULTIMODE computer program. The configuration state functions are products of ground-state occupied and virtual modals obtained using the vibrational self-consistent field method. The Gilbert and Jordan potential energy surface was used for the calculations. The resulting partition function was used to test the harmonic oscillator approximation and the separable-rotation approximation. The harmonic oscillator, rigid-rotator approximation is in error by a factor of 2.3 at 300 K, but we also propose a separable-rotation approximation that is accurate within 2% from 100 to 1000 K. (C) 2004 American Institute of Physics.

Large vibrational variational calculations using 'multimode' and an iterative diagonalization technique

We have favoured the variational (secular equation) method for the determination of the (ro-) vibrational energy levels of polyatomic molecules. We use predominantly the Watson Hamiltonian in normal coordinates and an associated given potential in the variational code 'Multimode'. The dominant cost is the construction and diagonalization of matrices of ever-increasing size. Here we address this problem, using pertubation theory to select dominant expansion terms within the Davidson-Liu iterative diagonalization method. Our chosen example is the twelve-mode molecule methanol, for which we have an ab initio representation of the potential which includes the internal rotational motion of the OH group relative to CH3. Our new algorithm allows us to obtain converged energy levels for matrices of dimensions in excess of 100 000.

Anharmonic vibrational levels of the two cyclic isomers of SiC3

Using coupled-cluster approach full six-dimensional analytic potential energy surfaces for two cyclic SiC3 isomers [C-C transannular bond (I) and Si-C transannular bond (II)] have been generated and used to calculate anharmonic vibrational wave functions. Several strong low-lying anharmonic resonances have been found. In both isomers already some of the fundamental transitions cannot be described within the harmonic approximation. Adiabatic electron affinities and ionization energies have been calculated as well. The Franck-Condon factors for the photodetachment processes c-SiC3-(I)-> c-SiC3(I) and c-SiC3-(II)-> c-SiC3(II) are reported. (c) 2006 American Institute of Physics.

Spectroscopy of highly excited vibrational states of HCN in its ground electronic state

An experimental technique based on a scheme of vibrationally mediated photodissociation has been developed and applied to the spectroscopic study of highly excited vibrational states in HCN, with energies between 29 000 and 30 000 cm(-1). The technique consists of four sequential steps: in the first one, a high power laser is used to vibrationally excite the sample to an intermediate state, typically (0,0,4), the nu(3) mode being approximately equivalent to the C-H stretching vibration. Then a second laser is used to search for transitions between this intermediate state and highly vibrationally excited states. When one of these transitions is found, HCN molecules are transferred to a highly excited vibrational state. Third, a ultraviolet laser photodissociates the highly excited molecules to produce H and CN radicals in its A (2)Pi electronic state. Finally, a fourth laser (probe) detects the presence of the CN(A) photofragments by means of an A-->B-->X laser induced fluorescence scheme. The spectra obtained with this technique, consisting of several rotationally resolved vibrational bands, have been analyzed. The positions and rotational parameters of the states observed are presented and compared with the results of a state-of-the-art variational calculation. (C) 2004 American Institute of Physics.

Highly excited vibrational states of HCN around 30 000 cm(-1)

This article describes the analysis and interpretation of rovibrational spectra involving highly excited vibrational states in the molecule of HCN. The spectra were obtained by means of the vibrationally mediated photodissociation technique. Analysis of the spectra revealed four bands with Sigma-Sigma structures that, once fitted, provided the energies and rotational constants of four new, highly excited vibrational states in the region of the potential energy surface near and above 30 000 cm(-1). All the states could be identified with the help of a state-of-the-art variational calculation. Together with the states already observed in previous works, eight highly excited states have so far been identified in this region. (c) 2006 American Institute of Physics.

Quantum studies of the vibrations in H3O2- and D3O2-

The vibrations of H3O2- and D3O2- are investigated using diffusion Monte Carlo (DMC) and vibrational configuration-interaction approaches, as implemented in the program MULTIMODE. These studies use the potential surface recently developed by Huang [ J. Am. Chem. Soc. 126, 5042 (2004)]. The focus of this work is on the vibrational ground state and fundamentals which occur between 100 and 3700 cm(-1). In most cases, excellent agreement is obtained between the fundamental frequencies calculated by the two approaches. This serves to demonstrate the power of both methods for treating this very anharmonic system. Based on the results of the MULTIMODE and DMC treatments, the extent and nature of the couplings in H3O2- and D3O2- are investigated. (C) 2005 American Institute of Physics.

Full-dimensional vibrational calculations for H5O2+ using an ab initio potential energy surface

We report quantum diffusion Monte Carlo (DMC) and variational calculations in full dimensionality for selected vibrational states of H5O2+ using a new ab initio potential energy surface [X. Huang, B. Braams, and J. M. Bowman, J. Chem. Phys. 122, 044308 (2005)]. The energy and properties of the zero-point state are focused on in the rigorous DMC calculations. OH-stretch fundamentals are also calculated using "fixed-node" DMC calculations and variationally using two versions of the code MULTIMODE. These results are compared with infrared multiphoton dissociation measurements of Yeh [L. I. Yeh, M. Okumura, J. D. Myers, J. M. Price, and Y. T. Lee, J. Chem. Phys. 91, 7319 (1989)]. Some preliminary results for the energies of several modes of the shared hydrogen are also reported.

Hydrogen spillover on carbon-supported metal catalysts studied by inelastic neutron scattering. Surface vibrational states and hydrogen riding modes

Hydrogen spillover on carbon-supported precious metal catalysts has been investigated with inelastic neutron scattering (INS) spectroscopy. The aim, which was fully realized, was to identify spillover hydrogen on the carbon support. The inelastic neutron scattering spectra of Pt/C, Ru/C, and PtRu/C fuel cell catalysts dosed with hydrogen were determined in two sets of experiments: with the catalyst in the neutron beam and, using an annular cell, with carbon in the beam and catalyst pellets at the edge of the cell excluded from the beam. The vibrational modes observed in the INS spectra were assigned with reference to the INS of a polycyclic aromatic hydrocarbon, coronene, taken as a molecular model of a graphite layer, and with the aid of computational modeling. Two forms of spillover hydrogen were identified: H at edge sites of a graphite layer (formed after ambient dissociative chemisorption of H-2), and a weakly bound layer of mobile H atoms (formed by surface diffusion of H atoms after dissociative chernisorption of H-2 at 500 K). The INS spectra exhibited characteristic riding modes of H on carbon and on Pt or Ru. In these riding modes H atoms move in phase with vibrations of the carbon and metal lattices. The lattice modes are amplified by neutron scattering from the H atoms attached to lattice atoms. Uptake of hydrogen, and spillover, was greater for the Ru containing catalysts than for the Pt/C catalyst. The INS experiments have thus directly demonstrated H spillover to the carbon support of these metal catalysts.

Analysis of red and yellow ochre samples from Clearwell Caves and Catalhoyuk by vibrational spectroscopy and other techniques

Ochre samples excavated from the neolithic site at Qatalhoyuk, Turkey have been compared with "native" ochres from Clearwell Caves, UK using infrared spectroscopy backed up by Raman spectroscopy, scanning electron microscopy (with energy-dispersive X-rays (EDX) analysis), powder X-ray diffraction, diffuse reflection UV-Vis and atomic absorption spectroscopies. For the Clearwell Caves ochres, which range in colour from yellow-orange to red-brown, it is shown that the colour is related to the nature of the chromophore present and not to any differences in particle size. The darker red ochres contain predominantly haematite while the yellow ochre contains only goethite. The ochres from Qatalhoyuk contain only about one-twentieth of the levels of iron found in the Clearwell Caves ochres. The iron oxide pigment (haematite in all cases studied here) has been mixed with a soft lime plaster which also contains calcite and silicate (clay) minerals. (C) 2003 Elsevier B.V. All rights reserved.