Towards Closed-loop Deep Brain Stimulation: Behavior Recognition from Human STN

Deep brain stimulation (DBS) provides significant therapeutic benefit for movement disorders such as Parkinson’s disease (PD). Current DBS devices lack real-time feedback (thus are open loop) and stimulation parameters are adjusted during scheduled visits with a clinician. A closed-loop DBS system may reduce power consumption and side effects by adjusting stimulation parameters based on patient’s behavior. Thus behavior detection is a major step in designing such systems. Various physiological signals can be used to recognize the behaviors. Subthalamic Nucleus (STN) Local field Potential (LFP) is a great candidate signal for the neural feedback, because it can be recorded from the stimulation lead and does not require additional sensors. This thesis proposes novel detection and classification techniques for behavior recognition based on deep brain LFP. Behavior detection from such signals is the vital step in developing the next generation of closed-loop DBS devices. LFP recordings from 13 subjects are utilized in this study to design and evaluate our method. Recordings were performed during the surgery and the subjects were asked to perform various behavioral tasks. Various techniques are used understand how the behaviors modulate the STN. One method studies the time-frequency patterns in the STN LFP during the tasks. Another method measures the temporal inter-hemispheric connectivity of the STN as well as the connectivity between STN and Pre-frontal Cortex (PFC). Experimental results demonstrate that different behaviors create different m odulation patterns in STN and it’s connectivity. We use these patterns as features to classify behaviors. A method for single trial recognition of the patient’s current task is proposed. This method uses wavelet coefficients as features and support vector machine (SVM) as the classifier for recognition of a selection of behaviors: speech, motor, and random. The proposed method is 82.4% accurate for the binary classification and 73.2% for classifying three tasks. As the next step, a practical behavior detection method which asynchronously detects behaviors is proposed. This method does not use any priori knowledge of behavior onsets and is capable of asynchronously detect the finger movements of PD patients. Our study indicates that there is a motor-modulated inter-hemispheric connectivity between LFP signals recorded bilaterally from STN. We utilize a non-linear regression method to measure this inter-hemispheric connectivity and to detect the finger movements. Our experimental results using STN LFP recorded from eight patients with PD demonstrate this is a promising approach for behavior detection and developing novel closed-loop DBS systems.

Mathematical approach for predicting non-negative tristimulus values using the CAT02 chromatic adaptation transform

It has been reported that for certain colour samples, the chromatic adaptation transform CAT02 imbedded in the CIECAM02 colour appearance model predicts corresponding colours with negative tristimulus values (TSVs), which can cause problems in certain applications. To overcome this problem, a mathematical approach is proposed for modifying CAT02. This approach combines a non-negativity constraint for the TSVs of corresponding colours with the minimization of the colour differences between those values for the corresponding colours obtained by visual observations and the TSVs of the corresponding colours predicted by the model, which is a constrained non-linear optimization problem. By solving the non-linear optimization problem, a new matrix is found. The performance of the CAT02 transform with various matrices including the original CAT02 matrix, and the new matrix are tested using visual datasets and the optimum colours. Test results show that the CAT02 with the new matrix predicted corresponding colours without negative TSVs for all optimum colours and the colour matching functions of the two CIE standard observers under the test illuminants considered. However, the accuracy with the new matrix for predicting the visual data is approximately 1 CIELAB colour difference unit worse compared with the original CAT02. This indicates that accuracy has to be sacrificed to achieve the non-negativity constraint for the TSVs of the corresponding colours.

Kinetic study of thermal 1,3-dipolar cycloaddition of azomethine ylides using differential scanning calorimetry as monitoring window

Kinetics of 1,3-dipolar cycloaddition involving azomethine ylides, generated from thermal [1,2]-prototropy of the corresponding imino ester, employing differential scanning calorimetry (DSC), is surveyed. Glycine and phenylalanine derived imino esters have different behavior. The first one prefers reacting with itself at 75 ºC, rather than with the dipolarophile. However, the α-substituted imino ester gives the cycloadduct at higher temperatures. The thermal dynamic analysis by 1H NMR of the neat reaction mixture of the glycine derivative reveals the presence of signals corresponding to the dipole in very small proportion. The non-isothermal and isothermal DSC curves of the cycloaddition of phenylalaninate and diisobutyl fumarate are obtained from freshly prepared samples. The application of known kinetic models and mathematical multiple non-linear regressions (NLR) allow to determine and to compare Ea, lnA, reaction orders, and reaction enthalpy. Finally a rate equation for each different temperature can be established for this particular thermal cycloaddition.

Repeatability, reproducibility, and accuracy of a novel pushbroom hyperspectral system

This study analyzes the repeatability, reproducibility and accuracy of a new hyperspectral system based on a pushbroom sensor as a means of measuring spectral features and color of materials and objects. The hyperspectral system consisted of a CCD camera, a spectrograph and an objective lens. An additional linear moving system allowed the mechanical scanning of the complete scene. A uniform overhead luminaire with daylight configuration was used to irradiate the scene using d:45 geometry. We followed the guidelines of the ASTM E2214-08 Standard Practice for Specifying and Verifying the Performance of Color-Measuring Instruments that define the standards and latest multidimensional procedures. The results obtained are analyzed in-depth and compared to those recently reported by other authors for spectrophotometers and multispectral systems. It can be concluded that hyperspectral systems are reliable and can be used in the industry to perform spectral and color readings with a high spatial resolution.

The effect of extreme cold temperatures on the risk of death in the two major Portuguese cities

It is well known that meteorological conditions influence the comfort and human health. Southern European countries, including Portugal, show the highest mortality rates during winter, but the effects of extreme cold temperatures in Portugal have never been estimated. The objective of this study was the estimation of the effect of extreme cold temperatures on the risk of death in Lisbon and Oporto, aiming the production of scientific evidence for the development of a real-time health warning system. Poisson regression models combined with distributed lag non-linear models were applied to assess the exposure-response relation and lag patterns of the association between minimum temperature and all-causes mortality and between minimum temperature and circulatory and respiratory system diseases mortality from 1992 to 2012, stratified by age, for the period from November to March. The analysis was adjusted for over dispersion and population size, for the confounding effect of influenza epidemics and controlled for long-term trend, seasonality and day of the week. Results showed that the effect of cold temperatures in mortality was not immediate, presenting a 1–2-day delay, reaching maximumincreased risk of death after 6–7 days and lasting up to 20–28 days. The overall effect was generally higher and more persistent in Lisbon than in Oporto, particularly for circulatory and respiratory mortality and for the elderly. Exposure to cold temperatures is an important public health problem for a relevant part of the Portuguese population, in particular in Lisbon.

Coral reefs in a century of rapid environmental change

Coral reefs are the most diverse marine ecosystem and embrace possibly millions of plant, animal and protist species. Mutualistic symbioses are a fundamental feature of coral reefs that have been used to explain their structure, biodiversity and existence. Complex inter-relationships between hosts, habitats and symbionts belie closely coupled nutrient and community dynamics that create the circumstances for something from nothing (or the oasis in a nutrient desert). The flip side of these dynamics is a close dependency between species, which results in a series of non-linear relationships as conditions change. These responses are being highlighted as anthropogenic influences increase across the world's tropical and subtropical coastlines. Caribbean as well as Indo-Pacific coral populations are now in a serious decline in many parts of the world. This has resulted in a significant reorganization of how coral reef ecosystems function. Among the spectra of changes brought about by humans is rapid climate change. Mass coral bleaching - the loss of the dinoflagellate symbionts from reef-building corals - and mortality has affected the world's coral reefs with increasing frequency and intensity since the late 1970s. Mass bleaching events, which often cover thousands of square kilometres of coral reefs, are triggered by small increases (+1-3degreesC) in water temperature. These increases in sea temperature are often seen during warm phase weather conditions (e.g. ENSO) and are increasing in size and magnitude. The loss of living coral cover (e.g. 16% globally in 1998, an exceptionally warm year) is resulting in an as yet unspecified reduction in the abundance of a myriad of other species. Projections from general circulation models (GCM) used to project changes in global temperature indicate that conditions even under the mildest greenhouse gas emission scenarios may exceed the thermal tolerances of most reef-building coral communities. Research must now explore key issues such as the extent to which the thermal tolerances of corals and their symbionts are dynamic if bleaching and disease are linked; how the loss of high densities of reef-building coral will affect other dependent species; and, how the loss of coral populations will affect the millions of people globally who depend on coral reefs for their daily survival.

Adaptive phase measurements in linear optical quantum computation

Photon counting induces an effective non-linear optical phase shift in certain states derived by linear optics from single photons. Although this non-linearity is non-deterministic, it is sufficient in principle to allow scalable linear optics quantum computation (LOQC). The most obvious way to encode a qubit optically is as a superposition of the vacuum and a single photon in one mode-so-called 'single-rail' logic. Until now this approach was thought to be prohibitively expensive (in resources) compared to 'dual-rail' logic where a qubit is stored by a photon across two modes. Here we attack this problem with real-time feedback control, which can realize a quantum-limited phase measurement on a single mode, as has been recently demonstrated experimentally. We show that with this added measurement resource, the resource requirements for single-rail LOQC are not substantially different from those of dual-rail LOQC. In particular, with adaptive phase measurements an arbitrary qubit state a alpha/0 > + beta/1 > can be prepared deterministically.

Optimization of co-current spray drying process of sugar-rich foods. Part I - Moisture and glass transition temperature profile during drying

A steady state mathematical model for co-current spray drying was developed for sugar-rich foods with the application of the glass transition temperature concept. Maltodextrin-sucrose solution was used as a sugar-rich food model. The model included mass, heat and momentum balances for a single droplet drying as well as temperature and humidity profile of the drying medium. A log-normal volume distribution of the droplets was generated at the exit of the rotary atomizer. This generation created a certain number of bins to form a system of non-linear first-order differential equations as a function of the axial distance of the drying chamber. The model was used to calculate the changes of droplet diameter, density, temperature, moisture content and velocity in association with the change of air properties along the axial distance. The difference between the outlet air temperature and the glass transition temperature of the final products (AT) was considered as an indicator of stickiness of the particles in spray drying process. The calculated and experimental AT values were close, indicating successful validation of the model. (c) 2004 Elsevier Ltd. All rights reserved.

Recursive filtering of images with symmetric extension

Recursive filters are widely used in image analysis due to their efficiency and simple implementation. However these filters have an initialisation problem which either produces unusable results near the image boundaries or requires costly approximate solutions such as extending the boundary manually. In this paper, we describe a method for the recursive filtering of symmetrically extended images for filters with symmetric denominator. We begin with an analysis of symmetric extensions and their effect on non-recursive filtering operators. Based on the non-recursive case, we derive a formulation of recursive filtering on symmetric domains as a linear but spatially varying implicit operator. We then give an efficient method for decomposing and solving the linear implicit system, along with a proof that this decomposition always exists. This decomposition needs to be performed only once for each dimension of the image. This yields a filtering which is both stable and consistent with the ideal infinite extension. The filter is efficient, requiring less computation than the standard recursive filtering. We give experimental evidence to verify these claims. (c) 2005 Elsevier B.V. All rights reserved.

Analytical solutions of linear finite- and small-strain one-dimensional consolidation

Analytical solutions are presented for linear finite-strain one-dimensional consolidation of initially unconsolidated soil layers with surcharge loading for both one- and two-way drainage. These solutions complement earlier solutions for initially unconsolidated soil layers without surcharge and initially normally consolidated soil layers with surcharge. Small-strain solutions for the consolidation of initially unconsolidated soil layers with surcharge loading are also presented, and the relationship between the earlier solutions for initially unconsolidated soil without surcharge and the corresponding small-strain solutions, which was not addressed in the earlier work, is clarified. The new solutions for initially unconsolidated soil with surcharge loading can be applied to the analysis of low stress consolidation tests and to the partial validation of numerical solutions of non-linear finite-strain consolidation. They also clarify a formerly perplexing aspect of finite-strain solution charts first noted in numerical solutions. Copyright (C) 2004 John Wiley Sons, Ltd.

Suppressing chaos via Lyapunov-Krasovskii's method

An algorithm for suppressing the chaotic oscillations in non-linear dynamical systems with singular Jacobian matrices is developed using a linear feedback control law based upon the Lyapunov-Krasovskii (LK) method. It appears that the LK method can serve effectively as a generalised method for the suppression of chaotic oscillations for a wide range of systems. Based on this method, the resulting conditions for undisturbed motions to be locally or globally stable are sufficient and conservative. The generalized Lorenz system and disturbed gyrostat equations are exemplified for the validation of the proposed feedback control rule. (c) 2005 Elsevier Ltd. All rights reserved.

Prediction of chaotic instabilities in a dragline bucket swing

The occurrence of chaotic instabilities is investigated in the swing motion of a dragline bucket during operation cycles. A dragline is a large, powerful, rotating multibody system utilised in the mining industry for removal of overburden. A simplified representative model of the dragline is developed in the form of a fundamental non-linear rotating multibody system with energy dissipation. An analytical predictive criterion for the onset of chaotic instability is then obtained in terms of critical system parameters using Melnikov's method. The model is shown to exhibit chaotic instability due to a harmonic slew torque for a range of amplitudes and frequencies. These chaotic instabilities could introduce irregularities into the motion of the dragline system, rendering the system difficult to control by the operator and/or would have undesirable effects on dragline productivity and fatigue lifetime. The sufficient analytical criterion for the onset of chaotic instability is shown to be a useful predictor of the phenomenon under steady and unsteady slewing conditions via comparisons with numerical results. (c) 2005 Elsevier Ltd. All rights reserved.

Coupled thermal, structural and vibrational analysis of a hypersonic engine for flight test

This paper describes a relatively simple and quick method for implementing aerodynamic heating models into a finite element code for non-linear transient thermal-structural and thermal-structural-vibrational analyses of a Mach 10 generic HyShot scramjet engine. The thermal-structural-vibrational response of the engine was studied for the descent trajectory from 60 to 26 km. Aerodynamic heating fluxes, as a function of spatial position and time for varying trajectory points, were implemented in the transient heat analysis. Additionally, the combined effect of varying dynamic pressure and thermal loads with altitude was considered. This aero-thermal-structural analysis capability was used to assess the temperature distribution, engine geometry distortion and yielding of the structural material due to aerodynamic heating during the descent trajectory, and for optimising the wall thickness, nose radius of leading edge, etc. of the engine intake. A structural vibration analysis was also performed following the aero-thermal-structural analysis to determine the changes in natural frequencies of the structural vibration modes that occur at the various temperatures associated with the descent trajectory. This analysis provides a unique and relatively simple design strategy for predicting and mitigating the thermal-structural-vibrational response of hypersonic engines. (C) 2006 Elsevier SAS. All rights reserved.

The cross-entropy method for continuous multi-extremal optimization

In recent years, the cross-entropy method has been successfully applied to a wide range of discrete optimization tasks. In this paper we consider the cross-entropy method in the context of continuous optimization. We demonstrate the effectiveness of the cross-entropy method for solving difficult continuous multi-extremal optimization problems, including those with non-linear constraints.

Emergent temporal behaviour and collaborative work

Although collaboration manifestly takes place in time, the role of time in shaping the behaviour of collaborations, and collaborative systems, is not well understood. Time is more than clock-time or the subjective experience of time; its effects on systems include differential rates of change of system elements, temporally non-linear behaviour and phenomena such as entrainment and synchronization. As a system driver, it generates emergent effects shaping systems and their behaviour. In the paper we present a systems view of time, and consider the implications of such a view through the case of collaborative development of a new university timetabling system. Teasing out the key temporal phenomena using the notion of temporal trajectories helps us understand the emergent temporal behaviour and suggests a means for improving outcomes.