Solving Schrödinger equation for two dimensional potentials using supersymmetry

The formalism of supersymmetric Quantum Mechanics can be extended to arbitrary dimensions. We introduce this formalism and explore its utility to solve the Schrödinger equation for a bidimensinal potential. This potential can be applied in several systems in physical and chemistry context , for instance, it can be used to study benzene molecule.

Statistical Model to DNA Melting

We explore a DNA statistical model to obtain information about the behavior of the thermodynamics quantities. Special attention is given to the thermal denaturation of this macromolecule.

Methods for construction and analysis of computational models in systems biology : applications to the modelling of the heat shock response and the self-assembly of intermediate filaments

Systems biology is a new, emerging and rapidly developing, multidisciplinary research field that aims to study biochemical and biological systems from a holistic perspective, with the goal of providing a comprehensive, system- level understanding of cellular behaviour. In this way, it addresses one of the greatest challenges faced by contemporary biology, which is to compre- hend the function of complex biological systems. Systems biology combines various methods that originate from scientific disciplines such as molecu- lar biology, chemistry, engineering sciences, mathematics, computer science and systems theory. Systems biology, unlike “traditional” biology, focuses on high-level concepts such as: network, component, robustness, efficiency, control, regulation, hierarchical design, synchronization, concurrency, and many others. The very terminology of systems biology is “foreign” to “tra- ditional” biology, marks its drastic shift in the research paradigm and it indicates close linkage of systems biology to computer science. One of the basic tools utilized in systems biology is the mathematical modelling of life processes tightly linked to experimental practice. The stud- ies contained in this thesis revolve around a number of challenges commonly encountered in the computational modelling in systems biology. The re- search comprises of the development and application of a broad range of methods originating in the fields of computer science and mathematics for construction and analysis of computational models in systems biology. In particular, the performed research is setup in the context of two biolog- ical phenomena chosen as modelling case studies: 1) the eukaryotic heat shock response and 2) the in vitro self-assembly of intermediate filaments, one of the main constituents of the cytoskeleton. The range of presented approaches spans from heuristic, through numerical and statistical to ana- lytical methods applied in the effort to formally describe and analyse the two biological processes. We notice however, that although applied to cer- tain case studies, the presented methods are not limited to them and can be utilized in the analysis of other biological mechanisms as well as com- plex systems in general. The full range of developed and applied modelling techniques as well as model analysis methodologies constitutes a rich mod- elling framework. Moreover, the presentation of the developed methods, their application to the two case studies and the discussions concerning their potentials and limitations point to the difficulties and challenges one encounters in computational modelling of biological systems. The problems of model identifiability, model comparison, model refinement, model inte- gration and extension, choice of the proper modelling framework and level of abstraction, or the choice of the proper scope of the model run through this thesis.

Potentials of biogas production from young bulls manure fed with different diets

Physical and chemical characteristics of manure are modified by different animal production systems. In cattle feeding system for young bulls there is an inversion of the proportion between forage and concentrate. In other words, the animals receive a smaller amount of forage compared to the traditional system. These changes in the manure characteristics involve changes in the treatment systems. The aim of this study was to determine the potential production of biogas of batch digesters fed with manure from young bulls that received two diets containing different proportions between forage and concentrate, with or without inoculums and submitted to three levels of temperature (25, 35 and 40(0)C). The evaluated parameters were total solids (TS) and volatile solids (VS) reduction and biogas potentials production. The digesters fed with manure from animals that received the diet 2 (80%C + 20% R) showed the largest reductions of TS and VS. About the potentials of biogas production there was interaction between the factors diet and inoculums, but no effects of temperatures. The treatment content manure from animals fed with diet 2 without inoculums presented the greatest potential of biogas production per kg of TS added (0.2123 m³).

Multivariate statistical analysis to support the minimum streamflow regionalization

ABSTRACT This study aimed to develop a methodology based on multivariate statistical analysis of principal components and cluster analysis, in order to identify the most representative variables in studies of minimum streamflow regionalization, and to optimize the identification of the hydrologically homogeneous regions for the Doce river basin. Ten variables were used, referring to the river basin climatic and morphometric characteristics. These variables were individualized for each of the 61 gauging stations. Three dependent variables that are indicative of minimum streamflow (Q7,10, Q90 and Q95). And seven independent variables that concern to climatic and morphometric characteristics of the basin (total annual rainfall – Pa; total semiannual rainfall of the dry and of the rainy season – Pss and Psc; watershed drainage area – Ad; length of the main river – Lp; total length of the rivers – Lt; and average watershed slope – SL). The results of the principal component analysis pointed out that the variable SL was the least representative for the study, and so it was discarded. The most representative independent variables were Ad and Psc. The best divisions of hydrologically homogeneous regions for the three studied flow characteristics were obtained using the Mahalanobis similarity matrix and the complete linkage clustering method. The cluster analysis enabled the identification of four hydrologically homogeneous regions in the Doce river basin.

Cognitive Auditory Evoked Potentials in Investigation of Hearing Discrimination

Preattentive perception of occasional deviating stimuli in the stream of standard stimuli can be recorded with cognitive event-related potential (ERP) mismatch negativity (MMN). The earlier detection of stimuli at the auditory cortex can be examined with N1 and P2 ERPs. The MMN recording does not require co-operation, it correlates with perceptual threshold, and even complex sounds can be used as stimuli. The aim of this study was to examine different aspects that should be considered when measuring discrimination of hearing with ERPs. The MMN was found to be stimulusintensity- dependent. As the intensity of sine wave stimuli was increased from 40 to 80 dB HL, MMN mean amplitudes increased. The effect of stimulus frequency on the MMN was studied so that the pitch difference would be equal in each stimulus block according to the psychophysiological mel scale or the difference limen of frequency (DLF). However, the blocks differed from each other. The contralateral white noise masking (50 dB EML) was found to attenuate the MMN amplitude when the right ear was stimulated. The N1 amplitude was attenuated and, in contrast, P2 amplitude was not affected by contralateral white noise masking. The perception and production of vowels by four postlingually deafened patients with a cochlear implant were studied. The MMN response could be elicited in the patient with the best vowel perception abilities. The results of the studies show that concerning the MMN recordings, the stimulus parameters and recording procedure design have a great influence on the results.

Neural mechanisms underlying conscious and unconscious vision. Evidence from event-related potentials and transcranial magnetic stimulation

Vision affords us with the ability to consciously see, and use this information in our behavior. While research has produced a detailed account of the function of the visual system, the neural processes that underlie conscious vision are still debated. One of the aims of the present thesis was to examine the time-course of the neuroelectrical processes that correlate with conscious vision. The second aim was to study the neural basis of unconscious vision, that is, situations where a stimulus that is not consciously perceived nevertheless influences behavior. According to current prevalent models of conscious vision, the activation of visual cortical areas is not, as such, sufficient for consciousness to emerge, although it might be sufficient for unconscious vision. Conscious vision is assumed to require reciprocal communication between cortical areas, but views differ substantially on the extent of this recurrent communication. Visual consciousness has been proposed to emerge from recurrent neural interactions within the visual system, while other models claim that more widespread cortical activation is needed for consciousness. Studies I-III compared models of conscious vision by studying event-related potentials (ERP). ERPs represent the brain’s average electrical response to stimulation. The results support the model that associates conscious vision with activity localized in the ventral visual cortex. The timing of this activity corresponds to an intermediate stage in visual processing. Earlier stages of visual processing may influence what becomes conscious, although these processes do not directly enable visual consciousness. Late processing stages, when more widespread cortical areas are activated, reflect the access to and manipulation of contents of consciousness. Studies IV and V concentrated on unconscious vision. By using transcranial magnetic stimulation (TMS) we show that when early visual cortical processing is disturbed so that subjects fail to consciously perceive visual stimuli, they may nevertheless guess (above chance-level) the location where the visual stimuli were presented. However, the results also suggest that in a similar situation, early visual cortex is necessary for both conscious and unconscious perception of chromatic information (i.e. color). Chromatic information that remains unconscious may influence behavioral responses when activity in visual cortex is not disturbed by TMS. Our results support the view that early stimulus-driven (feedforward) activation may be sufficient for unconscious processing. In conclusion, the results of this thesis support the view that conscious vision is enabled by a series of processing stages. The processes that most closely correlate with conscious vision take place in the ventral visual cortex ~200 ms after stimulus presentation, although preceding time-periods and contributions from other cortical areas such as the parietal cortex are also indispensable. Unconscious vision relies on intact early visual activation, although the location of visual stimulus may be unconsciously resolved even when activity in the early visual cortex is interfered with.

Statistical analysis of survey-based event history data with application to modeling of unemployment duration

Longitudinal surveys are increasingly used to collect event history data on person-specific processes such as transitions between labour market states. Surveybased event history data pose a number of challenges for statistical analysis. These challenges include survey errors due to sampling, non-response, attrition and measurement. This study deals with non-response, attrition and measurement errors in event history data and the bias caused by them in event history analysis. The study also discusses some choices faced by a researcher using longitudinal survey data for event history analysis and demonstrates their effects. These choices include, whether a design-based or a model-based approach is taken, which subset of data to use and, if a design-based approach is taken, which weights to use. The study takes advantage of the possibility to use combined longitudinal survey register data. The Finnish subset of European Community Household Panel (FI ECHP) survey for waves 1–5 were linked at person-level with longitudinal register data. Unemployment spells were used as study variables of interest. Lastly, a simulation study was conducted in order to assess the statistical properties of the Inverse Probability of Censoring Weighting (IPCW) method in a survey data context. The study shows how combined longitudinal survey register data can be used to analyse and compare the non-response and attrition processes, test the missingness mechanism type and estimate the size of bias due to non-response and attrition. In our empirical analysis, initial non-response turned out to be a more important source of bias than attrition. Reported unemployment spells were subject to seam effects, omissions, and, to a lesser extent, overreporting. The use of proxy interviews tended to cause spell omissions. An often-ignored phenomenon classification error in reported spell outcomes, was also found in the data. Neither the Missing At Random (MAR) assumption about non-response and attrition mechanisms, nor the classical assumptions about measurement errors, turned out to be valid. Both measurement errors in spell durations and spell outcomes were found to cause bias in estimates from event history models. Low measurement accuracy affected the estimates of baseline hazard most. The design-based estimates based on data from respondents to all waves of interest and weighted by the last wave weights displayed the largest bias. Using all the available data, including the spells by attriters until the time of attrition, helped to reduce attrition bias. Lastly, the simulation study showed that the IPCW correction to design weights reduces bias due to dependent censoring in design-based Kaplan-Meier and Cox proportional hazard model estimators. The study discusses implications of the results for survey organisations collecting event history data, researchers using surveys for event history analysis, and researchers who develop methods to correct for non-sampling biases in event history data.

Ca2+ dependence of gluconeogenesis stimulation by glucagon at different cytosolic NAD+-NADH redox potentials

The influence of Ca2+ on hepatic gluconeogenesis was measured in the isolated perfused rat liver at different cytosolic NAD+-NADH potentials. Lactate and pyruvate were the gluconeogenic substrates and the cytosolic NAD+-NADH potentials were changed by varying the lactate to pyruvate ratios from 0.01 to 100. The following results were obtained: a) gluconeogenesis from lactate plus pyruvate was not affected by Ca2+-free perfusion (no Ca2+ in the perfusion fluid combined with previous depletion of the intracellular pools); gluconeogenesis was also poorly dependent on the lactate to pyruvate ratios in the range of 0.1 to 100; only for a ratio equal to 0.01 was a significantly smaller gluconeogenic activity observed in comparison to the other ratios. b) In the presence of Ca2+, the increase in oxygen uptake caused by the infusion of lactate plus pyruvate at a ratio equal to 10 was the most pronounced one; in Ca2+-free perfusion the increase in oxygen uptake caused by lactate plus pyruvate infusion tended to be higher for all lactate to pyruvate ratios; the most pronounced difference was observed for a lactate/pyruvate ratio equal to 1. c) In the presence of Ca2+ the effects of glucagon on gluconeogenesis showed a positive correlation with the lactate to pyruvate ratios; for a ratio equal to 0.01 no stimulation occurred, but in the 0.1 to 100 range stimulation increased progressively, producing a clear parabolic dependence between the effects of glucagon and the lactate to pyruvate ratio. d) In the absence of Ca2+ the relationship between the changes caused by glucagon in gluconeogenesis and the lactate to pyruvate ratio was substantially changed; the dependence curve was no longer parabolic but sigmoidal in shape with a plateau beginning at a lactate/pyruvate ratio equal to 1; there was inhibition at the lactate to pyruvate ratios of 0.01 and 0.1 and a constant stimulation starting with a ratio equal to 1; for the lactate to pyruvate ratios of 10 and 100, stimulation caused by glucagon was much smaller than that found when Ca2+ was present. e) The effects of glucagon on oxygen uptake in the presence of Ca2+ showed a parabolic relationship with the lactate to pyruvate ratios which was closely similar to that found in the case of gluconeogenesis; the only difference was that inhibition rather than stimulation of oxygen uptake was observed for a lactate to pyruvate ratio equal to 0.01; progressive stimulation was observed in the 0.1 to 100 range. f) In the absence of Ca2+ the effects of glucagon on oxygen uptake were different; the dependence curve was sigmoidal at the onset, with a well-defined maximum at a lactate to pyruvate ratio equal to 1; this maximum was followed by a steady decline at higher ratios; at the ratios of 0.01 and 0.1 inhibition took place; oxygen uptake stimulation caused by glucagon was generally lower in the absence of Ca2+ except when the lactate to pyruvate ratio was equal to 1. The results of the present study demonstrate that stimulation of gluconeogenesis by glucagon depends on Ca2+. However, Ca2+ is only effective in helping gluconeogenesis stimulation by glucagon at highly negative redox potentials of the cytosolic NAD+-NADH system. The triple interdependence of glucagon-Ca2+-NAD+-NADH redox potential reveals highly complex interrelations that can only be partially understood at the present stage of knowledge

Clinical assessment of the effect of digital filtering on the detection of ventricular late potentials

Ventricular late potentials are low-amplitude signals originating from damaged myocardium and detected on the body surface by ECG filtering and averaging. Digital filters present in commercial equipment may interfere with the ability of arrhythmia stratification. We compared 40-Hz BiSpec (BI) and classical 40- to 250-Hz band-pass Butterworth bidirectional (BD) filters in terms of impact on time domain variables and diagnostic properties. In a transverse retrospective age-adjusted case-control study, 221 subjects with sinus rhythm without bundle branch block were divided into three groups after signal-averaged ECG acquisition: GI (N = 40), clinically normal controls, GII (N = 158), subjects with coronary heart disease without sustained monomorphic ventricular tachycardia (SMVT), and GIII (N = 23), subjects with heart disease and documented SMVT. Conventional variables analyzed from vector magnitude data after averaging to 0.3 µV final noise were obtained by application of each filter to the averaged signal, and evaluated in pairs by numerical comparison and by diagnostic agreement assessment, using conventional and optimized thresholds of normality. Significant differences were found between BI and BD variables in all groups, with diagnostic results showing significant disagreement between both filters [kappa value of 0.61 (P<0.05) for GII and 0.31 for GIII (P = NS)]. Sensitivity for SMVT was lower with BI than with BD (65.2 vs 91.3%, respectively, P<0.05). Filters provided significantly different numerical and diagnostic results and the BI filter showed only limited clinical application to risk stratification of ventricular arrhythmia.

Topographic abnormality of slow cortical potentials in schizophrenia

A recent study from our laboratory has provided evidence for the generation of slow potentials occurring in anticipation to task-performance feedback stimuli, in multiple association cortical areas, consistently including two prefrontal areas. In the present study, we intended to determine whether these slow potentials would indicate some abnormality (topographic) in schizophrenic patients, and thus serve as an indication of abnormal association cortex activity. We recorded slow potentials while subjects performed a paired-associates memory task. A 123-channel EEG montage and common average reference were used for 20 unmedicated schizophrenic (mean duration of illness: 11.3 ± 9.2 years; mean number of previous hospitalizations: 1.2 ± 1.9) and 22 healthy control subjects during a visual paired-associates matching task. For the topographic analysis, we used a simple index of individual topographic deviation from normality, corrected for absolute potential intensities. Slow potentials were observed in all subjects. Control subjects showed a simple spatial pattern of voltage extrema (left central positive and right prefrontal negative), whereas schizophrenic patients presented a more complex, fragmented pattern. Topographic deviation was significantly different between groups (P < 0.001). The increased topographic complexity in schizophrenics could be visualized in grand averages computed across subjects. Increased topographic complexity could also be seen when grand averages were computed for subgroups of patients assembled either according to task-performance (high versus low) or by their scores on psychopathological scales. There was no significant correlation between topographic deviation and psychopathology scores. We conclude that the slow potential topographic abnormalities of schizophrenia indicate an abnormality in the configuration of large-scale electrical activity in association cortices.

Statistical segmentation methods and color variance analysis of retinal images

In this research, the effectiveness of Naive Bayes and Gaussian Mixture Models classifiers on segmenting exudates in retinal images is studied and the results are evaluated with metrics commonly used in medical imaging. Also, a color variation analysis of retinal images is carried out to find how effectively can retinal images be segmented using only the color information of the pixels.

Comparison of the reliability of multifocal visual evoked cortical potentials generated by pattern reversal and pattern pulse stimulation

This study compared the effectiveness of the multifocal visual evoked cortical potentials (mfVEP) elicited by pattern pulse stimulation with that of pattern reversal in producing reliable responses (signal-to-noise ratio >1.359). Participants were 14 healthy subjects. Visual stimulation was obtained using a 60-sector dartboard display consisting of 6 concentric rings presented in either pulse or reversal mode. Each sector, consisting of 16 checks at 99% Michelson contrast and 80 cd/m² mean luminance, was controlled by a binary m-sequence in the time domain. The signal-to-noise ratio was generally larger in the pattern reversal than in the pattern pulse mode. The number of reliable responses was similar in the central sectors for the two stimulation modes. At the periphery, pattern reversal showed a larger number of reliable responses. Pattern pulse stimuli performed similarly to pattern reversal stimuli to generate reliable waveforms in R1 and R2. The advantage of using both protocols to study mfVEP responses is their complementarity: in some patients, reliable waveforms in specific sectors may be obtained with only one of the two methods. The joint analysis of pattern reversal and pattern pulse stimuli increased the rate of reliability for central sectors by 7.14% in R1, 5.35% in R2, 4.76% in R3, 3.57% in R4, 2.97% in R5, and 1.78% in R6. From R1 to R4 the reliability to generate mfVEPs was above 70% when using both protocols. Thus, for a very high reliability and thorough examination of visual performance, it is recommended to use both stimulation protocols.

STATISTICAL ANALYSYS OF THE SCFE OF A BRAZILAN MINERAL COAL

The influence of some process variables on the productivity of the fractions (liquid yield times fraction percent) obtained from SCFE of a Brazilian mineral coal using isopropanol and ethanol as primary solvents is analyzed using statistical techniques. A full factorial 23 experimental design was adopted to investigate the effects of process variables (temperature, pressure and cosolvent concentration) on the extraction products. The extracts were analyzed by the Preparative Liquid Chromatography-8 fractions method (PLC-8), a reliable, non destructive solvent fractionation method, especially developed for coal-derived liquids. Empirical statistical modeling was carried out in order to reproduce the experimental data. Correlations obtained were always greater than 0.98. Four specific process criteria were used to allow process optimization. Results obtained show that it is not possible to maximize both extract productivity and purity (through the minimization of heavy fraction content) simultaneously by manipulating the mentioned process variables.