Novel Features for Brain-Computer Interfaces

While conventional approaches of BCI feature extraction are based on the power spectrum, we have tried using nonlinear features for classifying BCI data. In this paper, we report our test results and findings, which indicate that the proposed method is a potentially useful addition to current feature extraction techniques

Van der Waals forces in density functional theory: perturbational long-range electron interaction corrections

Long-range exchange and correlation effects, responsible for the failure of currently used approximate density functionals in describing van der Waals forces, are taken into account explicitly after a separation of the electron-electron interaction in the Hamiltonian into short- and long-range components. We propose a "range-separated hybrid" functional based on a local density approximation for the short-range exchange-correlation energy, combined with a long-range exact exchange energy. Long-range correlation effects are added by a second-order perturbational treatment. The resulting scheme is general and is particularly well-adapted to describe van der Waals complexes, like rare gas dimers.Comment: 8 pages, 1 figure, submitted to Phys. Rev.

Electron affinities of the first- and second- row atoms: benchmark ab initio and density functional calculations

A benchmark ab initio and density functional (DFT) study has been carried out on the electron affinities of the first- and second-row atoms. The ab initio study involves basis sets of $spdfgh$ and $spdfghi$ quality, extrapolations to the 1-particle basis set limit, and a combination of the CCSD(T), CCSDT, and full CI electron correlation methods. Scalar relativistic and spin-orbit coupling effects were taken into account. On average, the best ab initio results agree to better than 0.001 eV with the most recent experimental results. Correcting for imperfections in the CCSD(T) method improves the mean absolute error by an order of magnitude, while for accurate results on the second-row atoms inclusion of relativistic corrections is essential. The latter are significantly overestimated at the SCF level; for accurate spin-orbit splitting constants of second-row atoms inclusion of (2s,2p) correlation is essential. In the DFT calculations it is found that results for the 1st-row atoms are very sensitive to the exchange functional, while those for second-row atoms are rather more sensitive to the correlation functional. While the LYP correlation functional works best for first-row atoms, its PW91 counterpart appears to be preferable for second-row atoms. Among ``pure DFT'' (nonhybrid) functionals, G96PW91 (Gill 1996 exchange combined with Perdew-Wang 1991 correlation) puts in the best overall performance. The best results overall are obtained with the 1-parameter hybrid modified Perdew-Wang (mPW1) exchange functionals of Adamo and Barone [J. Chem. Phys. {\bf 108}, 664 (1998)], with mPW1LYP yielding the best results for first-row, and mPW1PW91 for second-row atoms. Indications exist that a hybrid of the type $a$ mPW1LYP + $(1-a)$ mPW1PW91 yields better results than either of the constituent functionals.Comment: Phys. Rev. A, in press (revised version, review of issues concerning DFT and electron affinities added

Predictive response-relevant clustering of expression data provides insights into disease processes

This article describes and illustrates a novel method of microarray data analysis that couples model-based clustering and binary classification to form clusters of ;response-relevant' genes; that is, genes that are informative when discriminating between the different values of the response. Predictions are subsequently made using an appropriate statistical summary of each gene cluster, which we call the ;meta-covariate' representation of the cluster, in a probit regression model. We first illustrate this method by analysing a leukaemia expression dataset, before focusing closely on the meta-covariate analysis of a renal gene expression dataset in a rat model of salt-sensitive hypertension. We explore the biological insights provided by our analysis of these data. In particular, we identify a highly influential cluster of 13 genes-including three transcription factors (Arntl, Bhlhe41 and Npas2)-that is implicated as being protective against hypertension in response to increased dietary sodium. Functional and canonical pathway analysis of this cluster using Ingenuity Pathway Analysis implicated transcriptional activation and circadian rhythm signalling, respectively. Although we illustrate our method using only expression data, the method is applicable to any high-dimensional datasets

Raman Scattered He II λ\lambda 6545 Line in the Symbiotic Star V1016 Cygni

We present a spectrum of the symbiotic star V1016 Cyg observed with the 3.6 m Canada-France-Hawaii Telescope, in order to illustrate a method to measure the covering factor of the neutral scattering region around the giant component with respect to the hot emission region around the white dwarf component. In the spectrum, we find broad wings around H$\alpha$ and a broad emission feature around 6545${\rm \AA}$ that is blended with the [N II]$\lambda$ 6548 line. These two features are proposed to be formed by Raman scattering by atomic hydrogen, where the incident radiation is proposed to be UV continuum radiation around Ly$\beta$ in the former case and He II $\lambda$ 1025 emission line arising from $n=6\to n=2$ transitions for the latter feature. We remove the H$\alpha$ wings by a template Raman scattering wing profile and subtract the [N II] $\lambda$ 6548 line using the 3 times stronger [N II] $\lambda$ 6583 feature in order to isolate the He II Raman scattered 6545 \AA line. We obtain the flux ratio $F_{6545}/F_{6560}=0.24$ of the He II $\lambda$ 6560 emission line and the 6545 \AA feature for V1016 Cyg. Under the assumption that the He II emission from this object is isotropic, this ratio is converted to the ratio $\Phi_{6545}/\Phi_{1025}=0.17$ of the number of the incident photons and that of the scattered photons. This implies that the scattering region with H I column density $N_{HI}\ge 10^{20}{\rm cm^{-2}}$ covers 17 per cent of the emission region. By combining the presumed binary period $\sim 100$ yrs of this system we infer that a significant fraction of the slow stellar wind from the Mira component is ionized and that the scattering region around the Mira extends a few tens of AU, which is closely associated with the mass loss process of the Mira component.Comment: 12 pages, 6 figures, accepted for publication in Ap

Photochemistry of Furyl- and Thienyldiazomethanes: Spectroscopic Characterization of Triplet 3-Thienylcarbene

Photolysis (λ \u3e 543 nm) of 3-thienyldiazomethane (1), matrix isolated in Ar or N2 at 10 K, yields triplet 3-thienylcarbene (13) and α-thial-methylenecyclopropene (9). Carbene 13 was characterized by IR, UV/vis, and EPR spectroscopy. The conformational isomers of 3-thienylcarbene (s-E and s-Z) exhibit an unusually large difference in zero-field splitting parameters in the triplet EPR spectrum (|D/hc| = 0.508 cm–1, |E/hc| = 0.0554 cm–1; |D/hc| = 0.579 cm–1, |E/hc| = 0.0315 cm–1). Natural Bond Orbital (NBO) calculations reveal substantially differing spin densities in the 3-thienyl ring at the positions adjacent to the carbene center, which is one factor contributing to the large difference in D values. NBO calculations also reveal a stabilizing interaction between the sp orbital of the carbene carbon in the s-Z rotamer of 13 and the antibonding σ orbital between sulfur and the neighboring carbon—an interaction that is not observed in the s-E rotamer of 13. In contrast to the EPR spectra, the electronic absorption spectra of the rotamers of triplet 3-thienylcarbene (13) are indistinguishable under our experimental conditions. The carbene exhibits a weak electronic absorption in the visible spectrum (λmax = 467 nm) that is characteristic of triplet arylcarbenes. Although studies of 2-thienyldiazomethane (2), 3-furyldiazomethane (3), or 2-furyldiazomethane (4) provided further insight into the photochemical interconversions among C5H4S or C5H4O isomers, these studies did not lead to the spectroscopic detection of the corresponding triplet carbenes (2-thienylcarbene (11), 3-furylcarbene (23), or 2-furylcarbene (22), respectively)

CCSD(T) Study of CD3-O-CD3 and CH3-O-CD3 Far-Infrared Spectra

From a vibrationally corrected 3D potential energy surface determined with highly correlated ab initio calculations (CCSD(T)), the lowest vibrational energies of two dimethyl-ether isotopologues, 12CH3–16O–12CD3 (DME-d3) and 12CD3–16O–12CD3 (DME-d6), are computed variationally. The levels that can be populated at very low temperatures correspond to the COC-bending and the two methyl torsional modes. Molecular symmetry groups are used for the classification of levels and torsional splittings. DME-d6 belongs to the G36 group, as the most abundant isotopologue 12CH3–16O–12CH3 (DME-h6), while DME-d3 is a G18 species. Previous assignments of experimental Raman and far-infrared spectra are discussed from an effective Hamiltonian obtained after refining the ab initio parameters. Because a good agreement between calculated and experimental transition frequencies is reached, new assignments are proposed for various combination bands corresponding to the two deuterated isotopologues and for the 020 → 030 transition of DME-d6. Vibrationally corrected potential energy barriers, structural parameters, and anharmonic spectroscopic parameters are provided. For the 3N – 9 neglected vibrational modes, harmonic and anharmonic fundamental frequencies are obtained using second-order perturbation theory by means of CCSD and MP2 force fields. Fermi resonances between the COC-bending and the torsional modes modify DME-d3 intensities and the band positions of the torsional overtones

Ribosomal oxygenases are structurally conserved from prokaryotes to humans

2-Oxoglutarate (2OG)-dependent oxygenases have important roles in the regulation of gene expression via demethylation of N-methylated chromatin components1,2 and in the hydroxylation of transcription factors3 and splicing factor proteins4. Recently, 2OG-dependent oxygenases that catalyse hydroxylation of transfer RNA5,6,7 and ribosomal proteins8 have been shown to be important in translation relating to cellular growth, TH17-cell differentiation and translational accuracy9,10,11,12. The finding that ribosomal oxygenases (ROXs) occur in organisms ranging from prokaryotes to humans8 raises questions as to their structural and evolutionary relationships. In Escherichia coli, YcfD catalyses arginine hydroxylation in the ribosomal protein L16; in humans, MYC-induced nuclear antigen (MINA53; also known as MINA) and nucleolar protein 66 (NO66) catalyse histidine hydroxylation in the ribosomal proteins RPL27A and RPL8, respectively. The functional assignments of ROXs open therapeutic possibilities via either ROX inhibition or targeting of differentially modified ribosomes. Despite differences in the residue and protein selectivities of prokaryotic and eukaryotic ROXs, comparison of the crystal structures of E. coli YcfD and Rhodothermus marinus YcfD with those of human MINA53 and NO66 reveals highly conserved folds and novel dimerization modes defining a new structural subfamily of 2OG-dependent oxygenases. ROX structures with and without their substrates support their functional assignments as hydroxylases but not demethylases, and reveal how the subfamily has evolved to catalyse the hydroxylation of different residue side chains of ribosomal proteins. Comparison of ROX crystal structures with those of other JmjC-domain-containing hydroxylases, including the hypoxia-inducible factor asparaginyl hydroxylase FIH and histone Nε-methyl lysine demethylases, identifies branch points in 2OG-dependent oxygenase evolution and distinguishes between JmjC-containing hydroxylases and demethylases catalysing modifications of translational and transcriptional machinery. The structures reveal that new protein hydroxylation activities can evolve by changing the coordination position from which the iron-bound substrate-oxidizing species reacts. This coordination flexibility has probably contributed to the evolution of the wide range of reactions catalysed by oxygenases

The pharmacogenetics and pharmacogenomics knowledge base: accentuating the knowledge

PharmGKB is a knowledge base that captures the relationships between drugs, diseases/phenotypes and genes involved in pharmacokinetics (PK) and pharmacodynamics (PD). This information includes literature annotations, primary data sets, PK and PD pathways, and expert-generated summaries of PK/PD relationships between drugs, diseases/phenotypes and genes. PharmGKB's website is designed to effectively disseminate knowledge to meet the needs of our users. PharmGKB currently has literature annotations documenting the relationship of over 500 drugs, 450 diseases and 600 variant genes. In order to meet the needs of whole genome studies, PharmGKB has added new functionalities, including browsing the variant display by chromosome and cytogenetic locations, allowing the user to view variants not located within a gene. We have developed new infrastructure for handling whole genome data, including increased methods for quality control and tools for comparison across other data sources, such as dbSNP, JSNP and HapMap data. PharmGKB has also added functionality to accept, store, display and query high throughput SNP array data. These changes allow us to capture more structured information on phenotypes for better cataloging and comparison of data. PharmGKB is available at www.pharmgkb.or

A multi-threshold sampling method for TOF PET signal processing

As an approach to realizing all-digital data acquisition for positron emission tomography (PET), we have previously proposed and studied a multithreshold sampling method to generate samples of a PET event waveform with respect to a few user-defined amplitudes. In this sampling scheme, one can extract both the energy and timing information for an event. In this paper, we report our prototype implementation of this sampling method and the performance results obtained with this prototype. The prototype consists of two multi-threshold discriminator boards and a time-to-digital converter (TDC) board. Each of the multi-threshold discriminator boards takes one input and provides up to 8 threshold levels, which can be defined by users, for sampling the input signal. The TDC board employs the CERN HPTDC chip that determines the digitized times of the leading and falling edges of the discriminator output pulses. We connect our prototype electronics to the outputs of two Hamamatsu R9800 photomultiplier tubes (PMTs) that are individually coupled to a 6.25 x 6.25 x 25mm{sup 3} LSO crystal. By analyzing waveform samples generated by using four thresholds, we obtain a coincidence timing resolution of about 340 ps and an {approx}18% energy resolution at 511 keV. We are also able to estimate the decay-time constant from the resulting samples and obtain a mean value of 44 ns with an {approx}9 ns FWHM. In comparison, using digitized waveforms obtained at a 20 GSps sampling rate for the same LSO/PMT modules we obtain {approx}300 ps coincidence timing resolution, {approx}14% energy resolution at 511 keV, and {approx}5 ns FWHM for the estimated decay-time constant. Details of the results on the timing and energy resolutions by using the multi-threshold method indicate that it is a promising approach for implementing digital PET data acquisition