Analysis of nonlinear dynamics of vocal folds using high-speed video observation and biomechanical modeling

Primary voice production occurs in the larynx through vibrational movements carried out by vocal folds. However, many problems can affect this complex system resulting in voice disorders. In this context, time-frequency-shape analysis based on embedding phase space plots and nonlinear dynamics methods have been used to evaluate the vocal fold dynamics during phonation. For this purpose, the present work used high-speed video to record the vocal fold movements of three subjects and extract the glottal area time series using an image segmentation algorithm. This signal is used for an optimization method which combines genetic algorithms and a quasi-Newton method to optimize the parameters of a biomechanical model of vocal folds based on lumped elements (masses, springs and dampers). After optimization, this model is capable of simulating the dynamics of recorded vocal folds and their glottal pulse. Bifurcation diagrams and phase space analysis were used to evaluate the behavior of this deterministic system in different circumstances. The results showed that this methodology can be used to extract some physiological parameters of vocal folds and reproduce some complex behaviors of these structures contributing to the scientific and clinical evaluation of voice production. (C) 2010 Elsevier Inc. All rights reserved.

Dynamics of the 3:1 resonant planetary systems

Many of the discovered exoplanetary systems are involved inside mean-motion resonances. In this work we focus on the dynamics of the 3:1 mean-motion resonant planetary systems. Our main purpose is to understand the dynamics in the vicinity of the apsidal corotation resonance (ACR) which are stationary solutions of the resonant problem. We apply the semi-analytical method (Michtchenko et al., 2006) to construct the averaged three-body Hamiltonian of a planetary system near a 3:1 resonance. Then we obtain the families of ACR, composed of symmetric and asymmetric solutions. Using the symmetric stable solutions we observe the law of structures (Ferraz-Mello,1988), for different mass ratio of the planets. We also study the evolution of the frequencies of σ1, resonant angle, and Δω, the secular angle. The resonant domains outside the immediate vicinity of ACR are studied using dynamical maps techniques. We compared the results obtained to planetary systems near a 3:1 MMR, namely 55 Cnc b-c, HD 60532 b-c and Kepler 20 b-c.

Stability and dynamics of terrestrial planets in binary star systems: application on Alpha Centauri

Binary stars are frequent in the universe, with about 50% of the known main sequence stars being located at a multiple star system (Abt, 1979). Even though, they are universally thought as second rate sites for the location of exo-planets and the habitable zone, due to the difficulties of detection and high perturbation that could prevent planet formation and long term stability. In this work we show that planets in binary star systems can have regular orbits and remain on the habitable zone. We introduce a stability criterium based on the solution of the restricted three body problem and apply it to describe the short period planar and three-dimentional stability zones of S-type orbits around each star of the Alpha Centauri system. We develop as well a semi-analytical secular model to study the long term dynamics of fictional planets in the habitable zone of those stars and we verify that planets on the habitable zone would be in regular orbits with any eccentricity and with inclination to the binary orbital plane up until 35 degrees. We show as well that the short period oscillations on the semi-major axis is 100 times greater than the Earth's, but at all the time the planet would still be found inside the Habitable zone.

Applying many-body expansion aiming to make ab in initio molecular dynamics more efficient

In this present work we present a methodology that aims to apply the many-body expansion to decrease the computational cost of ab initio molecular dynamics, keeping acceptable accuracy on the results. We implemented this methodology in a program which we called ManBo. In the many-body expansion approach, we partitioned the total energy E of the system in contributions of one body, two bodies, three bodies, etc., until the contribution of the Nth body [1-3]: E = E1 + E2 + E3 + …EN. The E1 term is the sum of the internal energy of the molecules; the term E2 is the energy due to interaction between all pairs of molecules; E3 is the energy due to interaction between all trios of molecules; and so on. In Manbo we chose to truncate the expansion in the contribution of two or three bodies, both for the calculation of the energy and for the calculation of the atomic forces. In order to partially include the many-body interactions neglected when we truncate the expansion, we can include an electrostatic embedding in the electronic structure calculations, instead of considering the monomers, pairs and trios as isolated molecules in space. In simulations we made we chose to simulate water molecules, and use the Gaussian 09 as external program to calculate the atomic forces and energy of the system, as well as reference program for analyzing the accuracy of the results obtained with the ManBo. The results show that the use of the many-body expansion seems to be an interesting approach for reducing the still prohibitive computational cost of ab initio molecular dynamics. The errors introduced on atomic forces in applying such methodology are very small. The inclusion of an embedding electrostatic seems to be a good solution for improving the results with only a small increase in simulation time. As we increase the level of calculation, the simulation time of ManBo tends to largely decrease in relation to a conventional BOMD simulation of Gaussian, due to better scalability of the methodology presented. References [1] E. E. Dahlke and D. G. Truhlar; J. Chem. Theory Comput., 3, 46 (2007). [2] E. E. Dahlke and D. G. Truhlar; J. Chem. Theory Comput., 4, 1 (2008). [3] R. Rivelino, P. Chaudhuri and S. Canuto; J. Chem. Phys., 118, 10593 (2003).

Currents controlling sedimentation, deposits recording palaeo-hydrodynamics – lessons from the continental shelf-slope system off SE South America

The continental margin off SE South America hosts one of the world’s most energetic hydrodynamic regimes but also the second largest drainage system of the continent. Both, the ocean current system as well as the fluvial runoff are strongly controlled by the atmospheric circulation modes over the region. The distribution pattern of particular types of sediments on shelf and slope and the long-term built-up of depositional elements within the overall margin architecture are, thus, the product of both, seasonal to millennial variability as well as long-term environmental trends. This talk presents how the combination of different methodological approaches can be used to obtain a comprehensive picture of the variability of a shelf and upper-slope hydrodynamic system during Holocene times. The particular methods applied are: (a) Margin-wide stratigraphic information to elucidate the role of sea level for the oceanographic and sedimentary systems since the last glacial maximum; (b) Palaeoceanographic sediment proxies combined with palaeo-temperature indicating isotopes of bivalve shells to trace lateral shifts in the coastal oceanography (particularly of the shelf front) during the Holocene; (c) Neodymium isotopes to identify the shelf sediment transport routes resulting from the current regime; (d) Sedimentological/geochemical data to show the efficient mechanism of sand export from the shelf to the open ocean; (e) Diatom assemblages and sediment element distributions indicating palaeo-salinity and the changing marine influence to illustrate the Plata runoff history. Sea level has not only controlled the overall configuration of the shelf but also the position of the main sediment routes from the continent towards the ocean. The shelf front has shifted frequently since the last glacial times probably resulting from both, changes in the Westerly Winds intensity and in the shelf width itself. Remarkable is a southward shift of this front during the past two centuries possibly related to anthropogenic influences on the atmosphere. The oceanographic regime with its prominent hydrographic boundaries led to a clear separation of sedimentary provinces since shelf drowning. It is especially the shelf front which enhances shelf sediment export through a continuous high sand supply to the uppermost slope. Finally, the Plata River does not continuously provide sediment to the shelf but shows significant climate-related changes in discharge during the past centuries. Starting from these findings, three major fields of research should, in general, be further developed in future: (i) The immediate interaction of the hydrodynamic and sedimentary systems to close the gaps between deposit information and modern oceanographic dynamics; (ii) Material budget calculations for the marginal ocean system in terms of material fluxes, storage/retention capacities, and critical thresholds; (iii) The role of human activity on the atmospheric, oceanographic and solid material systems to unravel natural vs. anthropogenic effects and feedback mechanisms

Effect of resonant tunneling on exciton dynamics in coupled dot-well nanostructures

Excitonic dynamics in a hybrid dot-well system composed of InAs quantum dots (QDs) and an InGaAs quantum well (QW) is studied by means of femtosecond pump-probe reflection and continuous wave (cw) photoluminescence (PL) spectroscopy. The system is engineered to bring the QW ground exciton state into resonance with the third QD excited state. The resonant tunneling rate is varied by changing the effective barrier thickness between the QD and QW layers. This strongly affects the exciton dynamics in these hybrid structures as compared to isolated QW or QD systems. Optically measured decay times of the coupled system demonstrate dramatically different response to temperature change depending on the strength of the resonant tunneling or coupling strength. This reflects a competition between purely quantum mechanical and thermodynamical processes.

Numerical simulations of magma chamber dynamics at Campi Flegrei, and associated seismicity, deformation and gravity changes

Understanding the complex relationships between quantities measured by volcanic monitoring network and shallow magma processes is a crucial headway for the comprehension of volcanic processes and a more realistic evaluation of the associated hazard. This question is very relevant at Campi Flegrei, a volcanic quiescent caldera immediately north-west of Napoli (Italy). The system activity shows a high fumarole release and periodic ground slow movement (bradyseism) with high seismicity. This activity, with the high people density and the presence of military and industrial buildings, makes Campi Flegrei one of the areas with higher volcanic hazard in the world. In such a context my thesis has been focused on magma dynamics due to the refilling of shallow magma chambers, and on the geophysical signals detectable by seismic, deformative and gravimetric monitoring networks that are associated with this phenomenologies. Indeed, the refilling of magma chambers is a process frequently occurring just before a volcanic eruption; therefore, the faculty of identifying this dynamics by means of recorded signal analysis is important to evaluate the short term volcanic hazard. The space-time evolution of dynamics due to injection of new magma in the magma chamber has been studied performing numerical simulations with, and implementing additional features in, the code GALES (Longo et al., 2006), recently developed and still on the upgrade at the Istituto Nazionale di Geofisica e Vulcanologia in Pisa (Italy). GALES is a finite element code based on a physico-mathematical two dimensional, transient model able to treat fluids as multiphase homogeneous mixtures, compressible to incompressible. The fundamental equations of mass, momentum and energy balance are discretised both in time and space using the Galerkin Least-Squares and discontinuity-capturing stabilisation technique. The physical properties of the mixture are computed as a function of local conditions of magma composition, pressure and temperature.The model features enable to study a broad range of phenomenologies characterizing pre and sin-eruptive magma dynamics in a wide domain from the volcanic crater to deep magma feeding zones. The study of displacement field associated with the simulated fluid dynamics has been carried out with a numerical code developed by the Geophysical group at the University College Dublin (O’Brien and Bean, 2004b), with whom we started a very profitable collaboration. In this code, the seismic wave propagation in heterogeneous media with free surface (e.g. the Earth’s surface) is simulated using a discrete elastic lattice where particle interactions are controlled by the Hooke’s law. This method allows to consider medium heterogeneities and complex topography. The initial and boundary conditions for the simulations have been defined within a coordinate project (INGV-DPC 2004-06 V3_2 “Research on active volcanoes, precursors, scenarios, hazard and risk - Campi Flegrei”), to which this thesis contributes, and many researchers experienced on Campi Flegrei in volcanological, seismic, petrological, geochemical fields, etc. collaborate. Numerical simulations of magma and rock dynamis have been coupled as described in the thesis. The first part of the thesis consists of a parametric study aimed at understanding the eect of the presence in magma of carbon dioxide in magma in the convection dynamics. Indeed, the presence of this volatile was relevant in many Campi Flegrei eruptions, including some eruptions commonly considered as reference for a future activity of this volcano. A set of simulations considering an elliptical magma chamber, compositionally uniform, refilled from below by a magma with volatile content equal or dierent from that of the resident magma has been performed. To do this, a multicomponent non-ideal magma saturation model (Papale et al., 2006) that considers the simultaneous presence of CO2 and H2O, has been implemented in GALES. Results show that the presence of CO2 in the incoming magma increases its buoyancy force promoting convection ad mixing. The simulated dynamics produce pressure transients with frequency and amplitude in the sensitivity range of modern geophysical monitoring networks such as the one installed at Campi Flegrei . In the second part, simulations more related with the Campi Flegrei volcanic system have been performed. The simulated system has been defined on the basis of conditions consistent with the bulk of knowledge of Campi Flegrei and in particular of the Agnano-Monte Spina eruption (4100 B.P.), commonly considered as reference for a future high intensity eruption in this area. The magmatic system has been modelled as a long dyke refilling a small shallow magma chamber; magmas with trachytic and phonolitic composition and variable volatile content of H2O and CO2 have been considered. The simulations have been carried out changing the condition of magma injection, the system configuration (magma chamber geometry, dyke size) and the resident and refilling magma composition and volatile content, in order to study the influence of these factors on the simulated dynamics. Simulation results allow to follow each step of the gas-rich magma ascent in the denser magma, highlighting the details of magma convection and mixing. In particular, the presence of more CO2 in the deep magma results in more ecient and faster dynamics. Through this simulations the variation of the gravimetric field has been determined. Afterward, the space-time distribution of stress resulting from numerical simulations have been used as boundary conditions for the simulations of the displacement field imposed by the magmatic dynamics on rocks. The properties of the simulated domain (rock density, P and S wave velocities) have been based on data from literature on active and passive tomographic experiments, obtained through a collaboration with A. Zollo at the Dept. of Physics of the Federici II Univeristy in Napoli. The elasto-dynamics simulations allow to determine the variations of the space-time distribution of deformation and the seismic signal associated with the studied magmatic dynamics. In particular, results show that these dynamics induce deformations similar to those measured at Campi Flegrei and seismic signals with energies concentrated on the typical frequency bands observed in volcanic areas. The present work shows that an approach based on the solution of equations describing the physics of processes within a magmatic fluid and the surrounding rock system is able to recognise and describe the relationships between geophysical signals detectable on the surface and deep magma dynamics. Therefore, the results suggest that the combined study of geophysical data and informations from numerical simulations can allow in a near future a more ecient evaluation of the short term volcanic hazard.

Vortex dynamics in the high-temperature superconductors YBa 2 Cu 3 O 7 and Bi 2 Sr 2 CaCu 2 O 8 + delta in low- and high-dissipative transport measurements

This thesis investigates phenomena of vortex dynamics in type II superconductors depending on the dimensionality of the flux-line system and the strength of the driving force. In the low dissipative regime of Bi_2Sr_2CaCu_2O_{8+delta} (BSCCO) the influence of oxygen stoichiometry on flux-line tension was examined. An entanglement crossover of the vortex system at low magnetic fields was identified and a comprehensive B-T phase diagram of solid and fluid phases derived.In YBa_2Cu_3O_7 (YBCO) extremely long (>100 mm) high-quality measurement bridges allowed to extend the electric-field window in transport measurements by up to three orders of magnitude. Complementing analyses of the data conclusively produced dynamic exponents of the glass transition z~9 considerably higher than theoretically predicted and previously reported. In high-dissipative measurements a voltage instability appearing in the current-voltage characteristics of type II superconductors was observed for the first time in BSCCO and shown to result from a Larkin-Ovchinnikov flux-flow vortex instability under the influence of quasi-particle heating. However, in an analogous investigation of YBCO the instability was found to appear only in the temperature and magnetic-field regime of the vortex-glass state. Rapid-pulse measurements fully confirmed this correlation of vortex glass and instability in YBCO and revealed a constant rise time (~µs).

The language of action. How language translates the dynamics of our actions.

The general aim of the thesis was to investigate how and to what extent the characteristics of action organization are reflected in language, and how they influence language processing and understanding. Even though a huge amount of research has been devoted to the study of the motor effects of language, this issue is very debated in literature. Namely, the majority of the studies have focused on low-level motor effects such as effector-relatedness of action, whereas only a few studies have started to systematically investigate how specific aspects of action organization are encoded and reflected in language. After a review of previous studies on the relationship between language comprehension and action (chapter 1) and a critical discussion of some of them (chapter 2), the thesis is composed by three experimental chapters, each devoted to a specific aspect of action organization. Chapter 3 presents a study designed with the aim to disentangle the effective time course of the involvement of the motor system during language processing. Three kinematics experiments were designed in order to determine whether and, at which stage of motor planning and execution effector-related action verbs influence actions executed with either the same or a different effector. Results demonstrate that the goal of an action can be linguistically re-activated, producing a modulation of the motor response. In chapter 4, a second study investigates the interplay between the role of motor perspective (agent) and the organization of action in motor chains. More specifically, this kinematics study aims at deepening how goal can be translated in language, using as stimuli simple sentences composed by a pronoun (I, You, He/She) and a verb. Results showed that the perspective activated by the pronoun You reflects the motor pattern of the “agent” combined with the chain structure of the verb. These data confirm an early involvement of the motor system in language processing, suggesting that it is specifically modulated by the activation of the agent’s perspective. In chapter 5, the issue of perspective is specifically investigated, focusing on its role in language comprehension. In particular, this study aimed at determining how a specific perspective (induced for example by a personal pronoun) modulates motor behaviour during and after language processing. A classical compatibility effect (the Action-sentence compatibility effect) has been used to this aim. In three behavioural experiments the authors investigated how the ACE is modulated by taking first or third person perspective. Results from these experiments showed that the ACE effect occurs only when a first-person perspective is activated by the sentences used as stimuli. Overall, the data from this thesis contributed to disentangle several aspects of how action organization is translated in language, and then reactivated during language processing. This constitutes a new contribution to the field, adding lacking information on how specific aspects such as goal and perspective are linguistically described. In addition, these studies offer a new point of view to understand the functional implications of the involvement of the motor system during language comprehension, specifically from the point of view of our social interactions.

The Dynamics of a Small Model Glass Former as Viewed from Its Potential Energy Landscape

Despite intensive research during the last decades, thetheoreticalunderstanding of supercooled liquids and the glasstransition is stillfar from being complete. Besides analytical investigations,theso-called energy-landscape approach has turned out to beveryfruitful. In the literature, many numerical studies havedemonstratedthat, at sufficiently low temperatures, all thermodynamicquantities can be predicted with the help of the propertiesof localminima in the potential-energy-landscape (PEL). The main purpose of this thesis is to strive for anunderstanding ofdynamics in terms of the potential energy landscape. Incontrast to the study of static quantities, this requirestheknowledge of barriers separating the minima.Up to now, it has been the general viewpoint that thermallyactivatedprocesses ('hopping') determine the dynamics only belowTc(the critical temperature of mode-coupling theory), in thesense that relaxation rates follow from local energybarriers.As we show here, this viewpoint should be revisedsince the temperature dependence of dynamics is governed byhoppingprocesses already below 1.5Tc.At the example of a binary mixture of Lennard-Jonesparticles (BMLJ),we establish a quantitative link from the diffusioncoefficient,D(T), to the PEL topology. This is achieved in three steps:First, we show that it is essential to consider wholesuperstructuresof many PEL minima, called metabasins, rather than singleminima. Thisis a consequence of strong correlations within groups of PELminima.Second, we show that D(T) is inversely proportional to theaverageresidence time in these metabasins. Third, the temperaturedependenceof the residence times is related to the depths of themetabasins, asgiven by the surrounding energy barriers. We further discuss that the study of small (but not toosmall) systemsis essential, in that one deals with a less complex energylandscapethan in large systems. In a detailed analysis of differentsystemsizes, we show that the small BMLJ system consideredthroughout thethesis is free of major finite-size-related artifacts.

Structure and dynamics of supramolecular assemblies studied by advanced solid-state NMR spectroscopy

Ziel der vorliegenden Arbeit ist die Aufklärung von Struktur und Dynamik komplexer supramolekularer Systeme mittels Festkörper NMR Spektroskopie. Die Untersuchung von pi-pi Wechselwirkungen, welche einen entscheidenden Einfluss auf die strukturellen und dynamischen Eigenschaften supra- molekularer Systeme haben, hilft dabei, die Selbst- organisationsprozesse dieser komplexen Materialien besser zu verstehen. Mit dipolaren 1H-1H and 1H-13C Wiedereinkopplungs NMR Methoden unter schnellem MAS können sowohl 1H chemische Verschiebungen als auch dipolare 1H-1H und 1H-13C Kopplungen untersucht werden, ohne dass eine Isotopenmarkierung erforderlich ist. So erhält man detaillierte Informationen über die Struktur und die Beweglichkeit einzelner Molekül- segmente. In Verbindung mit sogenannten nucleus independent chemical shift (NICS) maps (berechnet mit ab-initio Methoden) lassen sich Abstände von Protonen relativ zu pi-Elektronensystemen bestimmen und so Strukturvorschläge ableiten. Mit Hilfe von homo- und heteronuklearen dipolaren Rotationsseitenbandenmustern könnenaußerdem Ordnungs- parameter für verschiedene Molekülsegmente bestimmt werden. Die auf diese Weise gewonnenen Informationen über die strukturellen und dynamischen Eigenschaften supramolekularer Systeme tragen dazu bei, strukturbestimmende Molekül- einheiten und Hauptordnungsphänomene zu identifizieren sowie lokale Wechselwirkungen zu quantifizieren, um so den Vorgang der Selbstorganisation besser zu verstehen.

Multi-level models and infrastructures for simulating biological system development

The hierarchical organisation of biological systems plays a crucial role in the pattern formation of gene expression resulting from the morphogenetic processes, where autonomous internal dynamics of cells, as well as cell-to-cell interactions through membranes, are responsible for the emergent peculiar structures of the individual phenotype. Being able to reproduce the systems dynamics at different levels of such a hierarchy might be very useful for studying such a complex phenomenon of self-organisation. The idea is to model the phenomenon in terms of a large and dynamic network of compartments, where the interplay between inter-compartment and intra-compartment events determines the emergent behaviour resulting in the formation of spatial patterns. According to these premises the thesis proposes a review of the different approaches already developed in modelling developmental biology problems, as well as the main models and infrastructures available in literature for modelling biological systems, analysing their capabilities in tackling multi-compartment / multi-level models. The thesis then introduces a practical framework, MS-BioNET, for modelling and simulating these scenarios exploiting the potential of multi-level dynamics. This is based on (i) a computational model featuring networks of compartments and an enhanced model of chemical reaction addressing molecule transfer, (ii) a logic-oriented language to flexibly specify complex simulation scenarios, and (iii) a simulation engine based on the many-species/many-channels optimised version of Gillespie’s direct method. The thesis finally proposes the adoption of the agent-based model as an approach capable of capture multi-level dynamics. To overcome the problem of parameter tuning in the model, the simulators are supplied with a module for parameter optimisation. The task is defined as an optimisation problem over the parameter space in which the objective function to be minimised is the distance between the output of the simulator and a target one. The problem is tackled with a metaheuristic algorithm. As an example of application of the MS-BioNET framework and of the agent-based model, a model of the first stages of Drosophila Melanogaster development is realised. The model goal is to generate the early spatial pattern of gap gene expression. The correctness of the models is shown comparing the simulation results with real data of gene expression with spatial and temporal resolution, acquired in free on-line sources.

The coupling between electron transfer and Protein/Solvent Dynamics in Photosynthetic Reaction Centers: Spectroscopic Studies in Amorphous Matrices

We investigated at the molecular level protein/solvent interactions and their relevance in protein function through the use of amorphous matrices at room temperature. As a model protein, we used the bacterial photosynthetic reaction center (RC) of Rhodobacter sphaeroides, a pigment protein complex which catalyzes the light-induced charge separation initiating the conversion of solar into chemical energy. The thermal fluctuations of the RC and its dielectric conformational relaxation following photoexcitation have been probed by analyzing the recombination kinetics of the primary charge-separated (P+QA-) state, using time resolved optical and EPR spectroscopies. We have shown that the RC dynamics coupled to this electron transfer process can be progressively inhibited at room temperature by decreasing the water content of RC films or of RC-trehalose glassy matrices. Extensive dehydration of the amorphous matrices inhibits RC relaxation and interconversion among conformational substates to an extent comparable to that attained at cryogenic temperatures in water-glycerol samples. An isopiestic method has been developed to finely tune the hydration level of the system. We have combined FTIR spectral analysis of the combination and association bands of residual water with differential light-minus-dark FTIR and high-field EPR spectroscopy to gain information on thermodynamics of water sorption, and on structure/dynamics of the residual water molecules, of protein residues and of RC cofactors. The following main conclusions were reached: (i) the RC dynamics is slaved to that of the hydration shell; (ii) in dehydrated trehalose glasses inhibition of protein dynamics is most likely mediated by residual water molecules simultaneously bound to protein residues and sugar molecules at the protein-matrix interface; (iii) the local environment of cofactors is not involved in the conformational dynamics which stabilizes the P+QA-; (iv) this conformational relaxation appears to be rather delocalized over several aminoacidic residues as well as water molecules weakly hydrogen-bonded to the RC.

Quantum effects and dynamics in hydrogen-bonded systems: a first-principles approach to spectroscopic experiments

Computer simulations have become an important tool in physics. Especially systems in the solid state have been investigated extensively with the help of modern computational methods. This thesis focuses on the simulation of hydrogen-bonded systems, using quantum chemical methods combined with molecular dynamics (MD) simulations. MD simulations are carried out for investigating the energetics and structure of a system under conditions that include physical parameters such as temperature and pressure. Ab initio quantum chemical methods have proven to be capable of predicting spectroscopic quantities. The combination of these two features still represents a methodological challenge. Furthermore, conventional MD simulations consider the nuclei as classical particles. Not only motional effects, but also the quantum nature of the nuclei are expected to influence the properties of a molecular system. This work aims at a more realistic description of properties that are accessible via NMR experiments. With the help of the path integral formalism the quantum nature of the nuclei has been incorporated and its influence on the NMR parameters explored. The effect on both the NMR chemical shift and the Nuclear Quadrupole Coupling Constants (NQCC) is presented for intra- and intermolecular hydrogen bonds. The second part of this thesis presents the computation of electric field gradients within the Gaussian and Augmented Plane Waves (GAPW) framework, that allows for all-electron calculations in periodic systems. This recent development improves the accuracy of many calculations compared to the pseudopotential approximation, which treats the core electrons as part of an effective potential. In combination with MD simulations of water, the NMR longitudinal relaxation times for 17O and 2H have been obtained. The results show a considerable agreement with the experiment. Finally, an implementation of the calculation of the stress tensor into the quantum chemical program suite CP2K is presented. This enables MD simulations under constant pressure conditions, which is demonstrated with a series of liquid water simulations, that sheds light on the influence of the exchange-correlation functional used on the density of the simulated liquid.

The Influence of Baseline Metabolic Rate and Exercise Training on Cardio-Respiratory Dynamics

Speeding the VO2 kinetics results in a reduction of the O2 deficit. Two factors might determine VO2 kinetics: oxygen delivery to muscle (Tschakovsky and Hughson 1999) and a muscle 'metabolic inertia' (Grassi et al. 1996). Therefore, in study 1 we investigated VO2 kinetics and cardiovascular system adaptations during step exercise transitions in different regions of the moderate domain. In study 2 we investigated muscle oxygenation and cardio-pulmonary adaptations during step exercise tests before, after and over a period of training. Study 1 methods: Seven subjects (26 ± 8 yr; 176 ± 5 cm; 69 ± 6 kg) performed 4 types of step transition from rest (0-50W; 0-100W) or elevate baseline (25-75W; 25-125W). GET and VO2max were assessed before testing. O2 uptake and were measured during testing. Study 2 methods: 10 subjects (25 ± 4 yr; 175 ± 9 cm; 71 ± 12 kg) performed a step transition test (0 to 100 W) before, after and during 4 weeks of endurance training (ET). VO2max and GET were assessed before and after of ET (40 minutes, 3 times a week, 60% O2max). VO2 uptake, Q and deoxyheamoglobin were measured during testing. Study 1 results: VO2 τ and the functional gain were slower in the upper regions of the moderate domain. Q increased more abruptly during rest to work condition. Q τ was faster than VO2 τ for each exercise step. Study 2 results: VO2 τ became faster after ET (25%) and particularly after 1 training session (4%). Q kinetics changed after 4 training sessions nevertheless it was always faster than VO2 τ. An attenuation in ∆[HHb] /∆VO2 was detectible. Conclusion: these investigations suggest that muscle fibres recruitment exerts a influence on the VO2 response within the moderate domain either during different forms of step transition or following ET.