Shapes and Dynamics of Blood Cells in Poiseuille and Shear Flows

The dynamics, shape, deformation, and orientation of red blood cells in microcirculation affect the rheology, flow resistance and transport properties of whole blood. This leads to important correlations of cellular and continuum scales. Furthermore, the dynamics of RBCs subject to different flow conditions and vessel geometries is relevant for both fundamental research and biomedical applications (e.g drug delivery). In this thesis, the behaviour of RBCs is investigated for different flow conditions via computer simulations. We use a combination of two mesoscopic particle-based simulation techniques, dissipative particle dynamics and smoothed dissipative particle dynamics. We focus on the microcapillary scale of several μm. At this scale, blood cannot be considered at the continuum but has to be studied at the cellular level. The connection between cellular motion and overall blood rheology will be investigated. Red blood cells are modelled as viscoelastic objects interacting hydrodynamically with a viscous fluid environment. The properties of the membrane, such as resistance against bending or shearing, are set to correspond to experimental values. Furthermore, thermal fluctuations are considered via random forces. Analyses corresponding to light scattering measurements are performed in order to compare to experiments and suggest for which situations this method is suitable. Static light scattering by red blood cells characterises their shape and allows comparison to objects such as spheres or cylinders, whose scattering signals have analytical solutions, in contrast to those of red blood cells. Dynamic light scattering by red blood cells is studied concerning its suitability to detect and analyse motion, deformation and membrane fluctuations. Dynamic light scattering analysis is performed for both diffusing and flowing cells. We find that scattering signals depend on various cell properties, thus allowing to distinguish different cells. The scattering of diffusing cells allows to draw conclusions on their bending rigidity via the effective diffusion coefficient. The scattering of flowing cells allows to draw conclusions on the shear rate via the scattering amplitude correlation. In flow, a RBC shows different shapes and dynamic states, depending on conditions such as confinement, physiological/pathological state and cell age. Here, two essential flow conditions are studied: simple shear flow and tube flow. Simple shear flow as a basic flow condition is part of any more complex flow. The velocity profile is linear and shear stress is homogeneous. In simple shear flow, we find a sequence of different cell shapes by increasing the shear rate. With increasing shear rate, we find rolling cells with cup shapes, trilobe shapes and quadrulobe shapes. This agrees with recent experiments. Furthermore, the impact of the initial orientation on the dynamics is studied. To study crowding and collective effects, systems with higher haematocrit are set up. Tube flow is an idealised model for the flow through cylindric microvessels. Without cell, a parabolic flow profile prevails. A single red blood cell is placed into the tube and subject to a Poiseuille profile. In tube flow, we find different cell shapes and dynamics depending on confinement, shear rate and cell properties. For strong confinements and high shear rates, we find parachute-like shapes. Although not perfectly symmetric, they are adjusted to the flow profile and maintain a stationary shape and orientation. For weak confinements and low shear rates, we find tumbling slippers that rotate and moderately change their shape. For weak confinements and high shear rates, we find tank-treading slippers that oscillate in a limited range of inclination angles and strongly change their shape. For the lowest shear rates, we find cells performing a snaking motion. Due to cell properties and resultant deformations, all shapes differ from hitherto descriptions, such as steady tank-treading or symmetric parachutes. We introduce phase diagrams to identify flow regimes for the different shapes and dynamics. Changing cell properties, the regime borders in the phase diagrams change. In both flow types, both the viscosity contrast and the choice of stress-free shape are important. For in vitro experiments, the solvent viscosity has often been higher than the cytosol viscosity, leading to a different pattern of dynamics, such as steady tank-treading. The stress-free state of a RBC, which is the state at zero shear stress, is still controversial, and computer simulations enable direct comparisons of possible candidates in equivalent flow conditions.

(Quantum-) critical dynamics of magnetoelectric materials

In this thesis the critical dynamics of several magnetoelectric compounds at their phase transition were examined. Mostly measurements of the dielectric properties in the frequency range of below 1 Hz up to 5 GHz were employed to evaluate the critical exponents for both magnetic field and temperature-dependent measurements. Most of the materials that are part of this work show anomalous behavior, especially at very low temperatures where quantum fluctuations are of the order of or even dominate those induced thermally. This anomalous behavior manifests in different forms. In Dy2Ti2O7 we demonstrate the existence of electric dipoles on magnetic monopoles. Here the dynamics at the critical endpoint located at 0.36K and in a magnetic field of 1T parallel to the [111] direction are of special interest. At this critical endpoint the expected critical slowing down of the dynamics could not only not be observed but instead the opposite, critical speeding-up by several orders of magnitude, could be demonstrated. Furthermore, we show that the phase diagram of Dy2Ti2O7 in this field direction can be reproduced solely from the dynamical properties, for example the resonance frequency of the observed relaxation that is connected to the monopole movement. Away from this point of the phase diagram the dynamics are slowing-down with reduction of temperature as one would expect. Additional measurements on Y2Ti2O7, a structurally identical but non-magnetic material, show only slowing down with reduction of temperature and no additional features. A possible explanation for the observed critical speeding-up is a coherent movement of magnetic monopoles close to the critical field that increases the resonance frequency by reducing the damping of the process. LiCuVO4 on the other hand behaves normally at its phase transition as long as the temperature is higher than 0.4 K. In this temperature regime the dynamics show critical slowing-down analogous to classical ferroelectric materials. This analogy extends also towards higher frequencies where the permittivity displays a ‘dispersion’ minimum that is temperature-dependent but of the order of 2 GHz. Below 0.4K the observed behavior changes drastically. Here we found no longer relaxational behavior but instead an excitation with very low energy. This low energy excitation was predicted by theory and is caused by nearly gapless soliton excitations within the 1D Cu2+ chains of LiCuVO4. Finally, in TbMnO3 the dynamics of the phase transition into the multiferroic phase was observed at roughly 27 K, a much higher temperature compared to the other materials. Here the expected critical slowing-down was observed, even though in low-frequency measurements this transition into the ferroelectric phase is overshadowed by the so-called c-axis relaxation. Therefore, only frequencies above 1MHz could be used to determine the critical exponents for both temperatureand magnetic-field-dependent measurements. This was done for both the peak frequency as well as the relaxation strength. In TbMnO3 an electromagnetic soft-mode with small optical weight causes the observed fluctuations, similar to the case of multiferroic MnWO4.

Phase Transitions in Engineered Ultracold Quantum Systems

The study of quantum degenerate gases has many applications in topics such as condensed matter dynamics, precision measurements and quantum phase transitions. We built an apparatus to create 87Rb Bose-Einstein condensates (BECs) and generated, via optical and magnetic interactions, novel quantum systems in which we studied the contained phase transitions. For our first experiment we quenched multi-spin component BECs from a miscible to dynamically unstable immiscible state. The transition rapidly drives any spin fluctuations with a coherent growth process driving the formation of numerous spin polarized domains. At much longer times these domains coarsen as the system approaches equilibrium. For our second experiment we explored the magnetic phases present in a spin-1 spin-orbit coupled BEC and the contained quantum phase transitions. We observed ferromagnetic and unpolarized phases which are stabilized by the spin-orbit coupling’s explicit locking between spin and motion. These two phases are separated by a critical curve containing both first-order and second-order transitions joined at a critical point. The narrow first-order transition gives rise to long-lived metastable states. For our third experiment we prepared independent BECs in a double-well potential, with an artificial magnetic field between the BECs. We transitioned to a single BEC by lowering the barrier while expanding the region of artificial field to cover the resulting single BEC. We compared the vortex distribution nucleated via conventional dynamics to those produced by our procedure, showing our dynamical process populates vortices much more rapidly and in larger number than conventional nucleation.

Increased CSF Aβ during the very early phase of cerebral Aβ deposition in mouse models

Abnormalities in brains of Alzheimer's disease (AD) patients are thought to start long before the first clinical symptoms emerge. The identification of affected individuals at this 'preclinical AD' stage relies on biomarkers such as decreased levels of the amyloid-β peptide (Aβ) in the cerebrospinal fluid (CSF) and positive amyloid positron emission tomography scans. However, there is little information on the longitudinal dynamics of CSF biomarkers, especially in the earliest disease stages when therapeutic interventions are likely most effective. To this end, we have studied CSF Aβ changes in three Aβ precursor protein transgenic mouse models, focusing our analysis on the initial Aβ deposition, which differs significantly among the models studied. Remarkably, while we confirmed the CSF Aβ decrease during the extended course of brain Aβ deposition, a 20-30% increase in CSF Aβ40 and Aβ42 was found around the time of the first Aβ plaque appearance in all models. The biphasic nature of this observed biomarker changes stresses the need for longitudinal biomarker studies in the clinical setting and the search for new 'preclinical AD' biomarkers at even earlier disease stages, by using both mice and human samples. Ultimately, our findings may open new perspectives in identifying subjects at risk for AD significantly earlier, and in improving the stratification of patients for preventive treatment strategies.

Nonlinear Optics and Carrier Dynamics in Nanostructured and Two-Dimensional Materials

Understanding and measuring the interaction of light with sub-wavelength structures and atomically thin materials is of critical importance for the development of next generation photonic devices.  One approach to achieve the desired optical properties in a material is to manipulate its mesoscopic structure or its composition in order to affect the properties of the light-matter interaction.  There has been tremendous recent interest in so called two-dimensional materials, consisting of only a single to a few layers of atoms arranged in a planar sheet.  These materials have demonstrated great promise as a platform for studying unique phenomena arising from the low-dimensionality of the material and for developing new types of devices based on these effects.  A thorough investigation of the optical and electronic properties of these new materials is essential to realizing their potential.  In this work we present studies that explore the nonlinear optical properties and carrier dynamics in nanoporous silicon waveguides, two-dimensional graphite (graphene), and atomically thin black phosphorus. We first present an investigation of the nonlinear response of nanoporous silicon optical waveguides using a novel pump-probe method. A two-frequency heterodyne technique is developed in order to measure the pump-induced transient change in phase and intensity in a single measurement. The experimental data reveal a characteristic material response time and temporally resolved intensity and phase behavior matching a physical model dominated by free-carrier effects that are significantly stronger and faster than those observed in traditional silicon-based waveguides.  These results shed light on the large optical nonlinearity observed in nanoporous silicon and demonstrate a new measurement technique for heterodyne pump-probe spectroscopy. Next we explore the optical properties of low-doped graphene in the terahertz spectral regime, where both intraband and interband effects play a significant role. Probing the graphene at intermediate photon energies enables the investigation of the nonlinear optical properties in the graphene as its electron system is heated by the intense pump pulse. By simultaneously measuring the reflected and transmitted terahertz light, a precise determination of the pump-induced change in absorption can be made. We observe that as the intensity of the terahertz radiation is increased, the optical properties of the graphene change from interband, semiconductor-like absorption, to a more metallic behavior with increased intraband processes. This transition reveals itself in our measurements as an increase in the terahertz transmission through the graphene at low fluence, followed by a decrease in transmission and the onset of a large, photo-induced reflection as fluence is increased.  A hybrid optical-thermodynamic model successfully describes our observations and predicts this transition will persist across mid- and far-infrared frequencies.  This study further demonstrates the important role that reflection plays since the absorption saturation intensity (an important figure of merit for graphene-based saturable absorbers) can be underestimated if only the transmitted light is considered. These findings are expected to contribute to the development of new optoelectronic devices designed to operate in the mid- and far-infrared frequency range.  Lastly we discuss recent work with black phosphorus, a two-dimensional material that has recently attracted interest due to its high mobility and direct, configurable band gap (300 meV to 2eV), depending on the number of atomic layers comprising the sample. In this work we examine the pump-induced change in optical transmission of mechanically exfoliated black phosphorus flakes using a two-color optical pump-probe measurement. The time-resolved data reveal a fast pump-induced transparency accompanied by a slower absorption that we attribute to Pauli blocking and free-carrier absorption, respectively. Polarization studies show that these effects are also highly anisotropic - underscoring the importance of crystal orientation in the design of optical devices based on this material. We conclude our discussion of black phosphorus with a study that employs this material as the active element in a photoconductive detector capable of gigahertz class detection at room temperature for mid-infrared frequencies.

Preferential location of prilocaine and etidocaine in phospholipid bilayers: A molecular dynamics study

In this work, we report a 20-ns constant pressure molecular dynamics simulation of the uncharged form of two amino-amide local anesthetics (LA). etidocaine and prilocaine, present at 1:3 LA:lipid, molar ratio inside the membrane, in the hydrated liquid crystal bilayer phase of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphatidylcholine (POPC). Both LAs induced lateral expansion and a concomitant contraction in the bilayer thickness. A decrease in the acyl chain segment order parameter, -S(CD), compared to neat bilayers, was also observed. Besides, both LA molecules got preferentially located in the hydrophobic acyl chains region, with a maximum probability at similar to 12 and similar to 10 angstrom from the center of the bilayer for prilocaine and etidocaine, respectively. (C) 2009 Elsevier B.V. All rights reserved.

Phase transitions in gene networks evolved under different selection rules

Mathematical models of gene regulation are a powerful tool for understanding the complex features of genetic control. While various modeling efforts have been successful at explaining gene expression dynamics, much less is known about how evolution shapes the structure of these networks. An important feature of gene regulatory networks is their stability in response to environmental perturbations. Regulatory systems are thought to have evolved to exist near the transition between stability and instability, in order to have the required stability to environmental fluctuations while also being able to achieve a wide variety of functions (corresponding to different dynamical patterns). We study a simplified model of gene network evolution in which links are added via different selection rules. These growth models are inspired by recent work on `explosive' percolation which shows that when network links are added through competitive rather than random processes, the connectivity phase transition can be significantly delayed, and when it is reached, it appears to be first order (discontinuous, e.g., going from no failure at all to large expected failure) instead of second order (continuous, e.g., going from no failure at all to very small expected failure). We find that by modifying the traditional framework for networks grown via competitive link addition to capture how gene networks evolve to avoid damage propagation, we also see significant delays in the transition that depend on the selection rules, but the transitions always appear continuous rather than `explosive'.

Whole-cell-analysis of live cardiomyocytes using wide-field interferometric phase microscopy.

We apply wide-field interferometric microscopy techniques to acquire quantitative phase profiles of ventricular cardiomyocytes in vitro during their rapid contraction with high temporal and spatial resolution. The whole-cell phase profiles are analyzed to yield valuable quantitative parameters characterizing the cell dynamics, without the need to decouple thickness from refractive index differences. Our experimental results verify that these new parameters can be used with wide field interferometric microscopy to discriminate the modulation of cardiomyocyte contraction dynamics due to temperature variation. To demonstrate the necessity of the proposed numerical analysis for cardiomyocytes, we present confocal dual-fluorescence-channel microscopy results which show that the rapid motion of the cell organelles during contraction preclude assuming a homogenous refractive index over the entire cell contents, or using multiple-exposure or scanning microscopy.

Transient scaling and resurgence of chimera states in networks of Boolean phase oscillators.

We study networks of nonlocally coupled electronic oscillators that can be described approximately by a Kuramoto-like model. The experimental networks show long complex transients from random initial conditions on the route to network synchronization. The transients display complex behaviors, including resurgence of chimera states, which are network dynamics where order and disorder coexists. The spatial domain of the chimera state moves around the network and alternates with desynchronized dynamics. The fast time scale of our oscillators (on the order of 100ns) allows us to study the scaling of the transient time of large networks of more than a hundred nodes, which has not yet been confirmed previously in an experiment and could potentially be important in many natural networks. We find that the average transient time increases exponentially with the network size and can be modeled as a Poisson process in experiment and simulation. This exponential scaling is a result of a synchronization rate that follows a power law of the phase-space volume.

Chemical Reaction Dynamics of the Liquid/Vapor Interface Studied by Mass Spectrometry

This thesis presents investigations of chemical reactions occurring at the liquid/vapor interface studied using novel sampling methodologies coupled with detection by mass spectrometry. Chapters 2 and 3 utilize the recently developed technique of field-induced droplet ionization mass spectrometry (FIDI-MS), in which the application of a strong electric field to a pendant microliter droplet results in the ejection of highly charged progeny droplets from the liquid surface. In Chapter 2, this method is employed to study the base-catalyzed dissociation of a surfactant molecule at the liquid/vapor interface upon uptake of ammonia from the gas phase. This process is observed to occur without significant modulation of the bulk solution pH, suggesting a transient increase in surface pH following the uptake of gaseous ammonia. Chapter 3 presents real-time studies of the oxidation of the model tropospheric organic compound glycolaldehyde by photodissociation of iron (III) oxalate complexes. The oxidation products of glycolaldehyde formed in this process are identified, and experiments in a deoxygenated environment identify the role of oxygen in the oxidation pathway and in the regeneration of iron (III) following photo-initiated reduction. Chapter 4 explores alternative methods for the study of heterogeneous reaction processes by mass spectrometric sampling from liquid surfaces. Bursting bubble ionization (BBI) and interfacial sampling with an acoustic transducer (ISAT) generate nanoliter droplets from a liquid surface that can be sampled via the atmospheric pressure interface of a mass spectrometer. Experiments on the oxidation of oleic acid by ozone using ISAT are also presented. Chapters 5 and 6 detail mechanistic studies and applications of free-radical-initiated peptide sequencing (FRIPS), a technique employing gas-phase free radical chemistry to the sequencing of peptides and proteins by mass spectrometry. Chapter 5 presents experimental and theoretical studies on the anomalous mechanism of dissociation observed in the presence of serine and threonine residues in peptides. Chapter 6 demonstrates the combination of FRIPS with ion mobility-mass spectrometry (IM-MS) for the separation of isomeric peptides.

Environmental Dynamics of Dissolved Black Carbon in Aquatic Ecosystems

Black carbon (BC), the incomplete combustion product from biomass and fossil fuel burning, is ubiquitously found in soils, sediments, ice, water and atmosphere. Because of its polyaromatic molecular characteristic, BC is believed to contribute significantly to the global carbon budget as a slow-cycling, refractory carbon pool. However, the mass balance between global BC generation and accumulation does not match, suggesting a removal mechanism of BC to the active carbon pool, most probable in a dissolved form. The presence of BC in waters as part of the dissolved organic matter (DOM) pool was recently confirmed via ultrahigh resolution mass spectrometry, and dissolved black carbon (DBC), a degradation product of charcoal, was found in marine and coastal environments. However, information on the loadings of DBC in freshwater environments and its global riverine flux from terrestrial systems to the oceans remained unclear. The main objectives of this study were to quantify DBC in diverse aquatic ecosystems and to determine its environmental dynamics. Surface water samples were collected from aquatic environments with a spatially significant global distribution, and DBC concentrations were determined by a chemical oxidation method coupled with HPLC detection. While it was clear that biomass burning was the main sources of BC, the translocation mechanism of BC to the dissolved phase was not well understood. Data from the regional studies and the developed global model revealed a strong positive correlation between DBC and dissolved organic carbon (DOC) dynamics, indicating a co-generation and co-translocation between soil OC and BC. In addition, a DOC-assistant DBC translocation mechanism was identified. Taking advantage of the DOC-DBC correlation model, a global riverine DBC flux to oceans on the order of 26.5 Mt C yr-1 (1 Mt = 1012 g) was determined, accounting for 10.6% of the global DOC flux. The results not only indicated that DOC was an important environmental intermediate for BC transfer and storage, but also provided an estimate of a major missing link in the global BC budget. The ever increasing DBC export caused by global warming will change the marine DOM quality and may have important consequences for carbon cycling in marine ecosystem.

Altered oscillatory dynamics of CA1 parvalbumin-basket cells during theta-gamma rhythmopathies of temporal lobe epilepsy

Recent reports in human demonstrate a role of theta– gamma coupling in memory for spatial episodes and a lack of coupling in people experiencing temporal lobe epilepsy, but the mechanisms are unknown. Using multisite silicon probe recordings of epileptic rats engaged in episodic-like object recognition tasks, we sought to evaluate the role of theta– gamma coupling in the absence of epileptiform activities. Our data reveal a specific association between theta– gamma (30 – 60 Hz) coupling at the proximal stratum radiatum of CA1 and spatial memory deficits. We targeted the microcircuit mechanisms with a novel approach to identify putative interneuronal types in tetrode recordings (parvalbumin basket cells in particular) and validated classification criteria in the epileptic context with neurochemical identification of intracellularly recorded cells. In epileptic rats, putative parvalbumin basket cells fired poorly modulated at the falling theta phase, consistent with weaker inputs from Schaffer collaterals and attenuated gamma oscillations, as evaluated by theta-phase decomposition of current–source density signals. We propose that theta– gamma interneuronal rhythmopathies of the temporal lobe are intimately related to episodic memory dysfunction in this condition.

Life history & environmental effects on telomere dynamics in Atlantic salmon

While much of the study of molecular biology inevitably focuses on the parts of the genome that contain active genes, there are also non-coding regions that nonetheless play an essential role in maintaining genome integrity. One such region are telomeres, which cap the ends of all eukaryotic chromosomes and play an important role in chromosome protection. Telomere loss occurs at each cell division as a result of the ‘end replication problem’ and a relatively short telomere length is indicative of poor biological state. Thus far, the majority of studies on the dynamics and role of telomeres have been biased towards certain taxa. Research to date has mostly focussed on humans, other mammals and birds. There has been far less research on the telomere dynamics of ectotherms. It is important that we do so, especially since ectothermic vertebrates do not seem to down-regulate telomerase expression in the same way as endotherms, suggesting that their telomere dynamics may be less predictable in the later life stages. The main objective of this thesis was therefore to investigate how life history and environmental effects may influence telomere dynamics in Atlantic salmon Salmo salar. I carried out carefully designed experiments, both in the laboratory and in the wild, using a longitudinal approach where possible, in order to address a number of specific questions that are connected to this central theme. In chapter 2, I demonstrate that there can be significant links between parental life history and offspring telomere dynamics. Maternal life history traits, in particular egg size, were most strongly related to offspring telomere length at the embryonic stages. Paternal life history traits, such as early life growth rate, had a greater association with offspring telomere dynamics in the later stages of development. In chapter 3, using a wild Atlantic salmon population, I found that most individuals experienced a reduction in telomere length during the migratory phase of their life cycle; however the relative rate of telomere loss was dependent on sex, with males experiencing a relatively greater loss. Unexpectedly, I also found that juvenile salmon that had the shortest telomeres at the time of outward migration, had the greatest probability of surviving through to the return migration. In chapter 4, again using a wild system involving experimental manipulations of juvenile Atlantic salmon in Scottish streams, I found that telomere length in juvenile fish was influenced by parental traits and by direct environmental effects. Faster-growing fish had shorter telomeres and there was a greater cost (in terms of reduced telomere length) if the growth occurred in a harsher environment. I also found a positive association between offspring telomere length and the growth history of their fathers (but not mothers), represented by the number of years that fathers had spent at sea. Chapter 5 explored the hypotheses that oxidative DNA damage, catalase (CAT) antioxidant activity and cell proliferation rate are underlying mechanisms linking incubation temperature and telomere dynamics in salmon embryos. No evidence was found for any such effects, but telomere lengths in salmon embryos were found to be significantly affected by the temperature of the water in which they were living. There is also evidence that telomere length significantly increases during embryonic development. In summary, this thesis has shown that a complex mix of environmental and parental effects appear to influence telomere dynamics in Atlantic salmon, with parental effects especially evident during early life stages. It also demonstrated that telomeres lengthen through the embryo stages of development before reducing once the fry begin feeding, indicating that the patterns of telomere loss commonly found in endotherms may differ in ectotherms. Reasons for this variation in telomere dynamics are presented in the final Discussion chapter of the thesis.

About stabilization of non-minimum phase systems by output feedback

This thesis work has been motivated by an internal benchmark dealing with the output regulation problem of a nonlinear non-minimum phase system in the case of full-state feedback. The system under consideration structurally suffers from finite escape time, and this condition makes the output regulation problem very hard even for very simple steady-state evolution or exosystem dynamics, such as a simple integrator. This situation leads to studying the approaches developed for controlling Non-minimum phase systems and how they affect feedback performances. Despite a lot of frequency domain results, only a few works have been proposed for describing the performance limitations in a state space system representation. In particular, in our opinion, the most relevant research thread exploits the so-called Inner-Outer Decomposition. Such decomposition allows splitting the Non-minimum phase system under consideration into a cascade of two subsystems: a minimum phase system (the outer) that contains all poles of the original system and an all-pass Non-minimum phase system (the inner) that contains all the unavoidable pathologies of the unstable zero dynamics. Such a cascade decomposition was inspiring to start working on functional observers for linear and nonlinear systems. In particular, the idea of a functional observer is to exploit only the measured signals from the system to asymptotically reconstruct a certain function of the system states, without necessarily reconstructing the whole state vector. The feature of asymptotically reconstructing a certain state functional plays an important role in the design of a feedback controller able to stabilize the Non-minimum phase system.