Analysis of temporal structure in spike trains of visual cortical area MT

The temporal structure of neuronal spike trains in the visual cortex can provide detailed information about the stimulus and about the neuronal implementation of visual processing. Spike trains recorded from the macaque motion area MT in previous studies (Newsome et al., 1989a; Britten et al., 1992; Zohary et al., 1994) are analyzed here in the context of the dynamic random dot stimulus which was used to evoke them. If the stimulus is incoherent, the spike trains can be highly modulated and precisely locked in time to the stimulus. In contrast, the coherent motion stimulus creates little or no temporal modulation and allows us to study patterns in the spike train that may be intrinsic to the cortical circuitry in area MT. Long gaps in the spike train evoked by the preferred direction motion stimulus are found, and they appear to be symmetrical to bursts in the response to the anti-preferred direction of motion. A novel cross-correlation technique is used to establish that the gaps are correlated between pairs of neurons. Temporal modulation is also found in psychophysical experiments using a modified stimulus. A model is made that can account for the temporal modulation in terms of the computational theory of biological image motion processing. A frequency domain analysis of the stimulus reveals that it contains a repeated power spectrum that may account for psychophysical and electrophysiological observations.

Some neurons tend to fire bursts of action potentials while others avoid burst firing. Using numerical and analytical models of spike trains as Poisson processes with the addition of refractory periods and bursting, we are able to account for peaks in the power spectrum near 40 Hz without assuming the existence of an underlying oscillatory signal. A preliminary examination of the local field potential reveals that stimulus-locked oscillation appears briefly at the beginning of the trial.

Nuclear Hartree-Fock calculations on parallel computers

Hartree-Fock (HF) calculations have had remarkable success in describing large nuclei at high spin, temperature and deformation. To allow full range of possible deformations, the Skyrme HF equations can be discretized on a three-dimensional mesh. However, such calculations are currently limited by the computational resources provided by traditional supercomputers. To take advantage of recent developments in massively parallel computing technology, we have implemented the LLNL Skyrme-force static and rotational HF codes on Intel's DELTA and GAMMA systems at Caltech.

We decomposed the HF code by assigning a portion of the mesh to each node, with nearest neighbor meshes assigned to nodes connected by communication· channels. This kind of decomposition is well-suited for the DELTA and the GAMMA architecture because the only non-local operations are wave function orthogonalization and the boundary conditions of the Poisson equation for the Coulomb field.

Our first application of the HF code on parallel computers has been the study of identical superdeformed (SD) rotational bands in the Hg region. In the last ten years, many SD rotational bands have been found experimentally. One very surprising feature found in these SD rotational bands is that many pairs of bands in nuclei that differ by one or two mass units have nearly identical deexcitation gamma-ray energies. Our calculations of the five rotational bands in ^(192)Hg and ^(194)Pb show that the filling of specific orbitals can lead to bands with deexcitation gamma-ray energies differing by at most 2 keV in nuclei differing by two mass units and over a range of angular momenta comparable to that observed experimentally. Our calculations of SD rotational bands in the Dy region also show that twinning can be achieved by filling or emptying some specific orbitals.

The interpretation of future precise experiments on atomic parity nonconservation (PNC) in terms of parameters of the Standard Model could be hampered by uncertainties in the atomic and nuclear structure. As a further application of the massively parallel HF calculations, we calculated the proton and neutron densities of the Cesium isotopes from A = 125 to A = 139. Based on our good agreement with experimental charge radii, binding energies, and ground state spins, we conclude that the uncertainties in the ratios of weak charges are less than 10^(-3), comfortably smaller than the anticipated experimental error.

Essays on cooperation and reciprocity

This dissertation comprises three essays that use theory-based experiments to gain understanding of how cooperation and efficiency is affected by certain variables and institutions in different types of strategic interactions prevalent in our society.

Chapter 2 analyzes indefinite horizon two-person dynamic favor exchange games with private information in the laboratory. Using a novel experimental design to implement a dynamic game with a stochastic jump signal process, this study provides insights into a relation where cooperation is without immediate reciprocity. The primary finding is that favor provision under these conditions is considerably less than under the most efficient equilibrium. Also, individuals do not engage in exact score-keeping of net favors, rather, the time since the last favor was provided affects decisions to stop or restart providing favors.

Evidence from experiments in Cournot duopolies is presented in Chapter 3 where players indulge in a form of pre-play communication, termed as revision phase, before playing the one-shot game. During this revision phase individuals announce their tentative quantities, which are publicly observed, and revisions are costless. The payoffs are determined only by the quantities selected at the end under real time revision, whereas in a Poisson revision game, opportunities to revise arrive according to a synchronous Poisson process and the tentative quantity corresponding to the last revision opportunity is implemented. Contrasting results emerge. While real time revision of quantities results in choices that are more competitive than the static Cournot-Nash, significantly lower quantities are implemented in the Poisson revision games. This shows that partial cooperation can be sustained even when individuals interact only once.

Chapter 4 investigates the effect of varying the message space in a public good game with pre-play communication where player endowments are private information. We find that neither binary communication nor a larger finite numerical message space results in any efficiency gain relative to the situation without any form of communication. Payoffs and public good provision are higher only when participants are provided with a discussion period through unrestricted text chat.

等离子体效应对类氦氖Kα线系电偶极辐射的影响

采用离子球模型，通过自洽求解Boltzmann方程和Poisson方程，得到类氦氖离子Kα线系的两条电偶极辐射光谱能量随等离子体环境的漂移．结果显示，Kα线系电偶极谱线随等离子体电子密度增大发生红移，红移量与等离子体电子密度有近似的正比关系；随着等离子体电子温度的降低，光谱红移对等离子体电子密度的敏感性增大。另外，所研究的两条谱线间的能量间隔随等离子体电子密度的增大而减小，减小量随等离子体电子密度的变化也呈现出近似的线性规律。值得注意的是，类氦氖Kα线系中两条电偶极谱线分别为互组合线与共振谱线，而其能量差

Euskal testu ahantzi bat frantses komedia batean

[EN] In this paper a short basque text included in a french comedy from XVIIth century is published and analyzed. That text appears in the second scene ofthe Raymond Poisson's comedy, "Le Poète Basque". It seems difficult to specify which variety of Basque is used by Poisson; anyway, I try to demonstrate that text shows characteristics of the dialects of Labourd and Low Navarre.

类氢氖在高温高密度等离子体中的光谱漂移

采用离子球模型，通过自治求解Poisson方程和Dirac方程，得到氖的类氢离子低能级组态的能级能量随等离子体电子温度和电子密度的变化关系，进一步研究了等离子体电子温度和电子密度对光谱漂移的影响。结果表明：光谱漂移随着等离子体电子密度的增大而增大，随着电子温度的升高而减小；谱线精细结构分裂随着电子密度的增大而减小．随着电子温度的升高而增大。等离子体对束缚电子的屏蔽是决定光谱漂移的主要原因。这些变化规律不仅对等离子体光谱模拟结果产生影响，而且使实验上观测光谱的相对或绝对漂移成为可能．从而为高密度等离子体诊断

Seismic reflection experiments imaging the physical nature of crustal structures in southern California

Seismic reflection methods have been extensively used to probe the Earth's crust and suggest the nature of its formative processes. The analysis of multi-offset seismic reflection data extends the technique from a reconnaissance method to a powerful scientific tool that can be applied to test specific hypotheses. The treatment of reflections at multiple offsets becomes tractable if the assumptions of high-frequency rays are valid for the problem being considered. Their validity can be tested by applying the methods of analysis to full wave synthetics.

Three studies illustrate the application of these principles to investigations of the nature of the crust in southern California. A survey shot by the COCORP consortium in 1977 across the San Andreas fault near Parkfield revealed events in the record sections whose arrival time decreased with offset. The reflectors generating these events are imaged using a multi-offset three-dimensional Kirchhoff migration. Migrations of full wave acoustic synthetics having the same limitations in geometric coverage as the field survey demonstrate the utility of this back projection process for imaging. The migrated depth sections show the locations of the major physical boundaries of the San Andreas fault zone. The zone is bounded on the southwest by a near-vertical fault juxtaposing a Tertiary sedimentary section against uplifted crystalline rocks of the fault zone block. On the northeast, the fault zone is bounded by a fault dipping into the San Andreas, which includes slices of serpentinized ultramafics, intersecting it at 3 km depth. These interpretations can be made despite complications introduced by lateral heterogeneities.

In 1985 the Calcrust consortium designed a survey in the eastern Mojave desert to image structures in both the shallow and the deep crust. Preliminary field experiments showed that the major geophysical acquisition problem to be solved was the poor penetration of seismic energy through a low-velocity surface layer. Its effects could be mitigated through special acquisition and processing techniques. Data obtained from industry showed that quality data could be obtained from areas having a deeper, older sedimentary cover, causing a re-definition of the geologic objectives. Long offset stationary arrays were designed to provide reversed, wider angle coverage of the deep crust over parts of the survey. The preliminary field tests and constant monitoring of data quality and parameter adjustment allowed 108 km of excellent crustal data to be obtained.

This dataset, along with two others from the central and western Mojave, was used to constrain rock properties and the physical condition of the crust. The multi-offset analysis proceeded in two steps. First, an increase in reflection peak frequency with offset is indicative of a thinly layered reflector. The thickness and velocity contrast of the layering can be calculated from the spectral dispersion, to discriminate between structures resulting from broad scale or local effects. Second, the amplitude effects at different offsets of P-P scattering from weak elastic heterogeneities indicate whether the signs of the changes in density, rigidity, and Lame's parameter at the reflector agree or are opposed. The effects of reflection generation and propagation in a heterogeneous, anisotropic crust were contained by the design of the experiment and the simplicity of the observed amplitude and frequency trends. Multi-offset spectra and amplitude trend stacks of the three Mojave Desert datasets suggest that the most reflective structures in the middle crust are strong Poisson's ratio (σ) contrasts. Porous zones or the juxtaposition of units of mutually distant origin are indicated. Heterogeneities in σ increase towards the top of a basal crustal zone at ~22 km depth. The transition to the basal zone and to the mantle include increases in σ. The Moho itself includes ~400 m layering having a velocity higher than that of the uppermost mantle. The Moho maintains the same configuration across the Mojave despite 5 km of crustal thinning near the Colorado River. This indicates that Miocene extension there either thinned just the basal zone, or that the basal zone developed regionally after the extensional event.

Geometric integration applied to moving mesh methods and degenerate Lagrangians

Moving mesh methods (also called r-adaptive methods) are space-adaptive strategies used for the numerical simulation of time-dependent partial differential equations. These methods keep the total number of mesh points fixed during the simulation, but redistribute them over time to follow the areas where a higher mesh point density is required. There are a very limited number of moving mesh methods designed for solving field-theoretic partial differential equations, and the numerical analysis of the resulting schemes is challenging. In this thesis we present two ways to construct r-adaptive variational and multisymplectic integrators for (1+1)-dimensional Lagrangian field theories. The first method uses a variational discretization of the physical equations and the mesh equations are then coupled in a way typical of the existing r-adaptive schemes. The second method treats the mesh points as pseudo-particles and incorporates their dynamics directly into the variational principle. A user-specified adaptation strategy is then enforced through Lagrange multipliers as a constraint on the dynamics of both the physical field and the mesh points. We discuss the advantages and limitations of our methods. The proposed methods are readily applicable to (weakly) non-degenerate field theories---numerical results for the Sine-Gordon equation are presented.

In an attempt to extend our approach to degenerate field theories, in the last part of this thesis we construct higher-order variational integrators for a class of degenerate systems described by Lagrangians that are linear in velocities. We analyze the geometry underlying such systems and develop the appropriate theory for variational integration. Our main observation is that the evolution takes place on the primary constraint and the 'Hamiltonian' equations of motion can be formulated as an index 1 differential-algebraic system. We then proceed to construct variational Runge-Kutta methods and analyze their properties. The general properties of Runge-Kutta methods depend on the 'velocity' part of the Lagrangian. If the 'velocity' part is also linear in the position coordinate, then we show that non-partitioned variational Runge-Kutta methods are equivalent to integration of the corresponding first-order Euler-Lagrange equations, which have the form of a Poisson system with a constant structure matrix, and the classical properties of the Runge-Kutta method are retained. If the 'velocity' part is nonlinear in the position coordinate, we observe a reduction of the order of convergence, which is typical of numerical integration of DAEs. We also apply our methods to several models and present the results of our numerical experiments.

Surface-enhanced resonance Raman scattering spectroscopy of single R6G molecules

Surface-enhanced resonance Raman scattering (SERRS) of Rhodamine 6G (R6G) adsorbed on colloidal silver clusters has been studied. Based on the great enhancement of the Raman signal and the quench of the fluorescence, the SERRS spectra of R6G were recorded for the samples of dye colloidal solution with different concentrations. Spectral inhomogeneity behaviours from single molecules in the dried sample films were observed with complementary evidences, such as spectral polarization, spectral diffusion, intensity fluctuation of vibrational lines and even "breathing" of the molecules. Sequential spectra observed from a liquid sample with an average of 0.3 dye molecules in the probed volume exhibited the expected Poisson distribution for actually measuring 0, 1 or 2 molecules. Difference between the SERRS spectra of R6G excited by linearly and circularly polarized light were experimentally measured.

Mechanistic studies of reactions at the single-molecule level using microfluidics with applications in molecular diagnostics

Motivated by needs in molecular diagnostics and advances in microfabrication, researchers started to seek help from microfluidic technology, as it provides approaches to achieve high throughput, high sensitivity, and high resolution. One strategy applied in microfluidics to fulfill such requirements is to convert continuous analog signal into digitalized signal. One most commonly used example for this conversion is digital PCR, where by counting the number of reacted compartments (triggered by the presence of the target entity) out of the total number of compartments, one could use Poisson statistics to calculate the amount of input target.

However, there are still problems to be solved and assumptions to be validated before the technology is widely employed. In this dissertation, the digital quantification strategy has been examined from two angles: efficiency and robustness. The former is a critical factor for ensuring the accuracy of absolute quantification methods, and the latter is the premise for such technology to be practically implemented in diagnosis beyond the laboratory. The two angles are further framed into a “fate” and “rate” determination scheme, where the influence of different parameters is attributed to fate determination step or rate determination step. In this discussion, microfluidic platforms have been used to understand reaction mechanism at single molecule level. Although the discussion raises more challenges for digital assay development, it brings the problem to the attention of the scientific community for the first time.

This dissertation also contributes towards developing POC test in limited resource settings. On one hand, it adds ease of access to the tests by incorporating massively producible, low cost plastic material and by integrating new features that allow instant result acquisition and result feedback. On the other hand, it explores new isothermal chemistry and new strategies to address important global health concerns such as cyctatin C quantification, HIV/HCV detection and treatment monitoring as well as HCV genotyping.

Edge function method applied to thin plates resting on Winkler foundation

The Edge Function method formerly developed by Quinlan⁽²⁵⁾ is applied to solve the problem of thin elastic plates resting on spring supported foundations subjected to lateral loads the method can be applied to plates of any convex polygonal shapes, however, since most plates are rectangular in shape, this specific class is investigated in this thesis. The method discussed can also be applied easily to other kinds of foundation models (e.g. springs connected to each other by a membrane) as long as the resulting differential equation is linear. In chapter VII, solution of a specific problem is compared with a known solution from literature. In chapter VIII, further comparisons are given. The problems of concentrated load on an edge and later on a corner of a plate as long as they are far away from other boundaries are also given in the chapter and generalized to other loading intensities and/or plates springs constants for Poisson's ratio equal to 0.2

Formal methods for control synthesis in partially observed environments : application to autonomous robotic manipulation

Modern robots are increasingly expected to function in uncertain and dynamically challenging environments, often in proximity with humans. In addition, wide scale adoption of robots requires on-the-fly adaptability of software for diverse application. These requirements strongly suggest the need to adopt formal representations of high level goals and safety specifications, especially as temporal logic formulas. This approach allows for the use of formal verification techniques for controller synthesis that can give guarantees for safety and performance. Robots operating in unstructured environments also face limited sensing capability. Correctly inferring a robot's progress toward high level goal can be challenging.

This thesis develops new algorithms for synthesizing discrete controllers in partially known environments under specifications represented as linear temporal logic (LTL) formulas. It is inspired by recent developments in finite abstraction techniques for hybrid systems and motion planning problems. The robot and its environment is assumed to have a finite abstraction as a Partially Observable Markov Decision Process (POMDP), which is a powerful model class capable of representing a wide variety of problems. However, synthesizing controllers that satisfy LTL goals over POMDPs is a challenging problem which has received only limited attention.

This thesis proposes tractable, approximate algorithms for the control synthesis problem using Finite State Controllers (FSCs). The use of FSCs to control finite POMDPs allows for the closed system to be analyzed as finite global Markov chain. The thesis explicitly shows how transient and steady state behavior of the global Markov chains can be related to two different criteria with respect to satisfaction of LTL formulas. First, the maximization of the probability of LTL satisfaction is related to an optimization problem over a parametrization of the FSC. Analytic computation of gradients are derived which allows the use of first order optimization techniques.

The second criterion encourages rapid and frequent visits to a restricted set of states over infinite executions. It is formulated as a constrained optimization problem with a discounted long term reward objective by the novel utilization of a fundamental equation for Markov chains - the Poisson equation. A new constrained policy iteration technique is proposed to solve the resulting dynamic program, which also provides a way to escape local maxima.

The algorithms proposed in the thesis are applied to the task planning and execution challenges faced during the DARPA Autonomous Robotic Manipulation - Software challenge.

Índice de massa corporal, percepção corporal e morbidade psiquiátrica menor entre funcionários de uma universidade pública do Rio de Janeiro

A presente investigação teve como principal objetivo explorar a associação entre o índice de massa corporal, percepção corporal e morbidade psiquiátrica menor (MPM) em uma população de funcionários técnico-administrativos de uma Universidade Pública do Rio de Janeiro. Foi realizado um estudo seccional (tipo censo) entre 3892 funcionários da Universidade com idade entre 20 e 81 anos. A morbidade psiquiátrica menor foi definida a partir de um instrumento utilizado para rastreamento de distúrbios psiquiátricos menor na população geral (GHQ-12). Encontrou-se nesta população uma prevalência de morbidade psiquiátrica menor de 22,3% para homens e 34,5% para mulheres. Elegeu-se para análise dos dados o método de regressão logística multivariada, considerando no modelo o IMC e a percepção corporal como variáveis explicativas e morbidade psiquiátrica menor como variável desfecho. Não foi encontrada associação entre o IMC e morbidade psiquiátrica menor, mesmo após o ajuste por sexo, idade, atividade física, renda familiar e morbidade referida. A percepção corporal, entretanto, se mostrou associada a presença de morbidade psiquiátrica menor. Sentir-se fora do peso considerado como ideal esteve associado a presença de sofrimento psíquico em homens e mulheres.

Ground-state energy shifts of H-like Ti under dense and hot plasma conditions

In this study, by adopting the ion sphere model, the self-consistent. field method is used with the Poisson-Boltzmann equation and the Dirac equation to calculate the ground-state energies of H-like Ti at a plasma electron density from 10(22) cm(-3) to 10(24) cm(-3) and the electron temperature from 100 eV to 3600 eV. The ground-state energy shifts of H-like Ti show different trends with the electron density and the electron temperature. It is shown that the energy shifts increase with the increase in the electron density and decrease with the increase in the electron temperature. The energy shifts are sensitive to the electron density, but only sensitive to the low electron temperature. In addition, an accurately fitting formula is obtained to fast estimate the ground-state energies of H-like Ti. Such fitted formula can also be used to estimate the critical electron density of pressure ionization for the ground state of H-like Ti.

Análise das alterações decorrentes do envelhecimento na mecânica respiratória pela técnica de oscilações forçadas

O crescimento do percentual de idosos na população ocorre mundialmente tornando necessário conhecer o impacto do processo de envelhecimento, neste contexto, do sistema respiratório. O desconhecimento do impacto do envelhecimento associado a diferentes graus de exposição a poluentes e a presença de comorbidade(s) dificulta a diagnose das pneumopatias acarretando aos idosos piora da qualidade de vida. São vantagens da Técnica de Oscilações Forçadas (FOT): alto potencial de aplicação em idosos, fácil realização, análise detalhada da mecânica respiratória, desempenho de papel complementar, bem como de alternativa na impossibilidade de realização dos exames tradicionais. Foi realizado um estudo experimental comparativo que objetivou investigar o impacto do envelhecimento no sistema respiratório pela FOT e pela espirometria entre grupos de diferentes faixas etárias, sendo a idade a variável independente e as variáveis dependentes, os parâmetros oscilométricos resistência em regime contínuo (R0) e das vias aéreas centrais (Rm), inclinação da resistência (S), frequência de ressonância (fr), reatância média (Xm), complacência dinâmica (Cdin,sr) e os parâmetros espirométricos (VEF1, CVF, VEF1/CVF e FEF/CVF). Foram realizados entrevista, exame clínico, radiografia torácica, avaliação da mecânica respiratória pela FOT e da função pulmonar pela espirometria. 255 indivíduos com idades entre 20 e 86 anos foram entrevistados. Destes, 175 foram excluídos, restando os 80 voluntários analisados, que foram divididos em 6 grupos de acordo com a faixa etária (A: 20 a 29 anos; B: 30 a 39 anos; C: 40 a 49 anos; D: 50 a 59 anos; E: 60 a 69 anos; F: 70 anos ou mais). Foram utilizados os testes de Shapiro-Wilkins, na avaliação da normalidade dos dados biométricos em cada grupo, Oneway ANOVA, na comparação entre os grupos, e Tukey HSD na comparação entre as classes subjacentes. A análise da associação entre duas variáveis foi realizada inicialmente pela regressão univariada entre os parâmetros oscilométricos, a idade e a altura. A regressão múltipla entre os parâmetros oscilométricos, idade e altura foi realizada em conjunto. Foi realizada a análise de confundimento ou modificação de efeito sobre o parâmetro altura na relação entre a idade e os parâmetros oscilométricos. A correção pelo fator altura foi realizada quando sua análise apresentava fator de confundimento. Quanto aos parâmetros resistivos, não foram observadas alterações em R0 e Rm com o envelhecimento enquanto que o declínio observado em S é discreto e não-significativo. Em relação aos parâmetros reativos, verificouse que Cdin,sr e Xm diminuem enquanto que fr aumenta com o processo de envelhecimento. Todas estas alterações são significativas. Todavia, a diminuição da Cdin,sr não apresenta relação com a idade e sim com a altura, que constituiu modificação do efeito. Nos demais parâmetros oscilométricos, a altura constituiu fator de confundimento. Quanto à espirometria, observou-se declínio significativo do VEF1, do VEF1/CVF e da CVF. O índice FEF/CVF apresentou declínio nãosignificativo. Concluindo, a resistência do sistema respiratório e a complacência dinâmica não se modificam enquanto a homogeneidade do sistema respiratório diminui com o processo de envelhecimento.