Latent variable models for neural data analysis

The brain is perhaps the most complex system to have ever been subjected to rigorous scientific investigation. The scale is staggering: over 10^11 neurons, each making an average of 10^3 synapses, with computation occurring on scales ranging from a single dendritic spine, to an entire cortical area. Slowly, we are beginning to acquire experimental tools that can gather the massive amounts of data needed to characterize this system. However, to understand and interpret these data will also require substantial strides in inferential and statistical techniques. This dissertation attempts to meet this need, extending and applying the modern tools of latent variable modeling to problems in neural data analysis.

It is divided into two parts. The first begins with an exposition of the general techniques of latent variable modeling. A new, extremely general, optimization algorithm is proposed - called Relaxation Expectation Maximization (REM) - that may be used to learn the optimal parameter values of arbitrary latent variable models. This algorithm appears to alleviate the common problem of convergence to local, sub-optimal, likelihood maxima. REM leads to a natural framework for model size selection; in combination with standard model selection techniques the quality of fits may be further improved, while the appropriate model size is automatically and efficiently determined. Next, a new latent variable model, the mixture of sparse hidden Markov models, is introduced, and approximate inference and learning algorithms are derived for it. This model is applied in the second part of the thesis.

The second part brings the technology of part I to bear on two important problems in experimental neuroscience. The first is known as spike sorting; this is the problem of separating the spikes from different neurons embedded within an extracellular recording. The dissertation offers the first thorough statistical analysis of this problem, which then yields the first powerful probabilistic solution. The second problem addressed is that of characterizing the distribution of spike trains recorded from the same neuron under identical experimental conditions. A latent variable model is proposed. Inference and learning in this model leads to new principled algorithms for smoothing and clustering of spike data.

Response of multiple simultaneously recorded macaque area LIP neurons in a memory saccade task

Cells in the lateral intraparietal cortex (LIP) of rhesus macaques respond vigorously and in spatially-tuned fashion to briefly memorized visual stimuli. Responses to stimulus presentation, memory maintenance, and task completion are seen, in varying combination from neuron to neuron. To help elucidate this functional segmentation a new system for simultaneous recording from multiple neighboring neurons was developed. The two parts of this dissertation discuss the technical achievements and scientific discoveries, respectively.

Technology. Simultanous recordings from multiple neighboring neurons were made with four-wire bundle electrodes, or tetrodes, which were adapted to the awake behaving primate preparation. Signals from these electrodes were partitionable into a background process with a 1/f-like spectrum and foreground spiking activity spanning 300-6000 Hz. Continuous voltage recordings were sorted into spike trains using a state-of-the-art clustering algorithm, producing a mean of 3 cells per site. The algorithm classified 96% of spikes correctly when tetrode recordings were confirmed with simultaneous intracellular signals. Recording locations were verified with a new technique that creates electrolytic lesions visible in magnetic resonance imaging, eliminating the need for histological processing. In anticipation of future multi-tetrode work, the chronic chamber microdrive, a device for long-term tetrode delivery, was developed.

Science. Simultaneously recorded neighboring LIP neurons were found to have similar preferred targets in the memory saccade paradigm, but dissimilar peristimulus time histograms, PSTH). A majority of neighboring cell pairs had a difference in preferred directions of under 45° while the trial time of maximal response showed a broader distribution, suggesting homogeneity of tuning with het erogeneity of function. A continuum of response characteristics was present, rather than a set of specific response types; however, a mapping experiment suggests this may be because a given cell's PSTH changes shape as well as amplitude through the response field. Spike train autocovariance was tuned over target and changed through trial epoch, suggesting different mechanisms during memory versus background periods. Mean frequency-domain spike-to-spike coherence was concentrated below 50 Hz with a significant maximum of 0.08; mean time-domain coherence had a narrow peak in the range ±10 ms with a significant maximum of 0.03. Time-domain coherence was found to be untuned for short lags (10 ms), but significantly tuned at larger lags (50 ms).

Attosecond pulse extreme-ultraviolet photoionization in a two-color laser field

Attosecond-pulse extreme-ultraviolet (XUV) photoionization in a two-color laser field is investigated. Attosecond pulse trains with different numbers of pulses are examined, and their strong dependence on photoelectronic spectra is found. Single-color driving-laser-field-assisted attosecond XUV photoionization cannot determine the number of attosecond pulses from the photoelectronic energy spectrum that are detected orthogonally to the beam direction and the electric field vector of the linearly polarized laser field. A two-color-field-assisted XUV photoionization scheme is proposed for directly determining the number of attosecond pulses from a spectrum detected orthogonally. (C) 2005 Optical Society of America.

脉冲啁啾对于阿秒脉冲的影响

采用单电子近似和软核势模型，通过数值求解一维含时薛定谔方程，理论研究了当脉冲分别带有正、负啁啾的情况下所产生的阿秒脉冲，分析了不同脉冲啁啾特性对阿秒脉冲的强度和宽度的影响，研究结果表明，无论是正啁啾还是负啁啾，随着啁啾量的增加，都将使激光脉冲由产生单个阿秒脉冲趋向于产生阿秒脉冲链，正啁啾和负啁啾对于阿秒脉冲宽度的影响是不同的，负啁啾对于阿秒脉冲宽度影响很小，适当的负啁啾有利于缩小阿秒脉冲的宽度；而正啁啾脉冲产生的阿秒脉冲较无啁啾时展宽，且随着啁啾量的增加，其阿秒脉冲宽度迅速增大。

影响单电子非线性汤姆孙散射因素的研究

应用电子汤姆孙散射的经典理论，通过理论分析和计算机模拟，研究了超短超强激光脉冲作用下电子产生的辐射脉冲的性质．计算表明，在这种情况下，电子的辐射通常以阿秒脉冲列的形式出现．讨沦了不同激光场参数(包括激光强度、脉宽、初相位和偏振态)、不同电子初始状态(初始速度和位置)对辐射脉冲的时间和空间特性的影响．通常在相对论光强条件下，激光强度越大，电子辐射越强，脉宽越窄，中心频率越大，并且方向性越好；电子在线偏振激光中产生的辐射效率，比在同样强度下圆偏振激光中产生的效率更高；无论入射光是线偏振光，还是圆偏振光，辐射场

Alternative trajectories for decarbonizing urban transportation

4 p.

Trends in Income and Price Elasticities of Transport Demand (1850-2010)

40 p.

Long Run Demand for Energy Services: the Role of Economic and Technological Development

36 p.

Decarbonising urban transportation

18 p.

Modelización acústica del interior de un Tren de Alta Velocidad

[ES]El presente Trabajo de Fin de Grado, titulado ‘Modelización Acústica del Interior de un Tren de Alta Velocidad’, tiene como objetivo el análisis acústico a bajas frecuencias del habitáculo de un coche de un tren de alta velocidad. La temática ha sido elegida debido a la creciente presencia de trenes de alta velocidad en nuestro entorno, en los cuales el Control Acústico o de Ruidos es un aspecto muy importante para el confort del medio de transporte. Dentro de los ruidos por los que se ve afectado el tren, los de baja frecuencia son los más difíciles de tratar y controlar, y es por ello que van a constituir la línea de trabajo de este proyecto. La estructura de lo que se va a tratar en este proyecto es la siguiente: En primer lugar, se analizará el contexto en que se encuadra este proyecto y las razones que han llevado a su realización. Posteriormente, se explicarán los fundamentos teóricos que hay detrás de los análisis que se van a realizar. Más adelante se pasará al análisis del caso práctico que se ha elegido para ilustrar el proyecto: el habitáculo de un coche del tren Serie 120 de Renfe, construido por CAF [1, 2, 3]. Se realizarán tanto análisis modales como de respuesta forzada. Dicho ejemplo servirá de base para asentar las conclusiones y proponer la aplicación que se les puede dar a éstas, así como las líneas de investigación para las cuales este proyecto puede ser un punto de partida.

飞秒激光诱导波纹状微突起结构

利用飞秒激光振荡器产生的脉冲对镀有铬层的玻璃和石英基片进行微加工,发现两种样品表面均有波纹状的微突起结构产生。这些微突起结构离开样品表面的高度为10～300 nm不等,并且随着激光功率的增大而增加,在一定功率下达到饱和状态。它们的形貌、尺寸和高度取决于入射飞秒激光的能流以及飞秒脉冲的参数。通过化学方法证明了这些微突起结构是由玻璃和石英的主要成分SiO2组成的,并非样品表面的铬元素。此外,通过选取适当的飞秒激光功率和样品加工速度,制作了两种不同周期和线宽的光栅结构,显示出飞秒激光振荡器良好的加工性能。

Estimação de velocidade de composições usando processamento de vídeo

Esta dissertação tem como objetivo estudar um método para estimar a velocidade de veículos ferroviários usando processamento de vídeo. O sistema proposto consiste na utilização de câmeras ao longo das vias férreas que permitam não só o monitoramento do tráfego ferroviário, mas cujo vídeo captado possa ser utilizado para a obtenção de estimativas para a velocidade instantânea dos trens que por ela trafegam. Tal sistema seria utilizado independentemente dos sistemas de controle já utilizados pela operadora do sistema ferroviário, permitindo que os controladores possam ter uma segunda análise no caso de falha da primeira, assim como um modelo que permita avaliar a velocidade instantânea do veículo ferroviário ao longo do percurso. Os algoritmos de rastreamento empregados para esse fim abordaram diferentes métodos. Confrontaram-se os resultados obtidos com os algoritmos propostos com os dados empíricos de forma a determinar aquele com melhor resposta dada as características do sistema. O algoritmo que apresentou os melhores resultados emprega um único bloco de referência para todos os quadros comparados. A métrica de similaridade responsável por determinar quais blocos são mais ou menos similares dentro do universo de busca estipulado é a soma de diferenças absolutas (SAD, Sum of Absolute Differences). O tempo de processamento requerido por cada um dos métodos de rastreamento estudados também foi considerado nas avaliações de resultados apresentadas. Uma comparação realizada entre as velocidades coletadas e aquelas informadas pelo sistema de controle mostraram que os resultados obtidos com o sistema atual, utilizando a sinalização apenas por circuito de via apresenta resultados pouco confiáveis com erros bastante significativos. Os resultados obtidos com o sistema proposto apresentaram erros menores quando comparados àqueles obtidos pelo sistema vigente, apresentando-se assim como uma solução viável e de baixo custo quando comparada às técnicas atualmente empregadas para a medida de velocidade de trens.

Determinação de compostos monoaromáticos voláteis provenientes de emissões de veículos leves, em estacionamentos subterrâneos.

Este trabalho visa determinar a contribuição das emissões evaporativas provenientes dos veículos leves de passageiro, para a degradação da qualidade do ar atmosférico. O objetivo principal é avaliar as concentrações compostos monoaromáticos voláteis Benzeno, Tolueno, Etilbenzeno e Xilenos (BTEX) em ambientes confinados, sendo este realizado em um local que caracterize a realidade da frota veicular da Região metropolitana do Rio de Janeiro. As amostras foram coletadas em um estacionamento subterrâneo de um Shopping Center da zona norte do Rio de Janeiro, através do sistema de amostragem ativa, utilizando cartucho de carvão ativo como adsorvente. As amostras foram extraídas com solvente orgânico e analisadas posteriormente por Cromatografia gasosa acoplada à espectrometria de massas (CGEM). As médias dos resultados obtidos foram 52,7 g.m-3 para o benzeno, 203,6 g.m-3 para o tolueno, 44,6 g.m-3 para o etilbenzeno, 115,7 g.m-3 para os xilenos, sendo o tolueno o composto encontrado em maior abundância. Esses resultados foram comparados com resultados encontrados na literatura de emissões veiculares em ambientes confinados como garagens e túneis. Foi investigada a correlação com as emissões do veículo em movimento, obtidas através de estudos previamente realizados em um túnel de grande circulação e as emissões obtidas no estacionamento subterrâneo. Através desses dados ficou demonstrada diferença das fontes de emissão.

Welding of Al alloy plates using a novel high-average-power pulsed CO2 laser

A novel high-average-power pulsed CO2 laser with a unique electrode structure is presented. The operation of a 5-kW transverse-flow CO2 laser with the preionized pulse-train switched technique results in pulsation of the laser power, and the average laser power is about 5 kW. The characteristic of this technique is switching the preionized pulses into pulse trains so as to use the small preionized power (hundreds of watts) to control the large main-discharge power (tens of kilowatts). By this means, the cost and the complexity of the power supply are greatly reduced. The welding of LF2, LF21, LD2, and LY12 aluminum alloy plates has been successfully achieved using this laser. (c) 2005 Society of Photo-Optical Instrumentation Engineers.

Welding of Al alloy plates using a novel high-average-power pulsed CO2 laser

A novel high-average-power pulsed CO2 laser with a unique electrode structure is presented. The operation of a 5-kW transverse-flow CO2 laser with the preionized pulse-train switched technique results in pulsation of the laser power, and the average laser power is about 5 kW. The characteristic of this technique is switching the preionized pulses into pulse trains so as to use the small preionized power (hundreds of watts) to control the large main-discharge power (tens of kilowatts). By this means, the cost and the complexity of the power supply are greatly reduced. The welding of LF2, LF21, LD2, and LY12 aluminum alloy plates has been successfully achieved using this laser. (c) 2005 Society of Photo-Optical Instrumentation Engineers.