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.

Topics in equilibrium and nonequilibrium thermodynamics: computing crystalline free energies and engineering Maxwell’s demon.

This dissertation covers two separate topics in statistical physics. The first part of the dissertation focuses on computational methods of obtaining the free energies (or partition functions) of crystalline solids. We describe a method to compute the Helmholtz free energy of a crystalline solid by direct evaluation of the partition function. In the many-dimensional conformation space of all possible arrangements of N particles inside a periodic box, the energy landscape consists of localized islands corresponding to different solid phases. Calculating the partition function for a specific phase involves integrating over the corresponding island. Introducing a natural order parameter that quantifies the net displacement of particles from lattices sites, we write the partition function in terms of a one-dimensional integral along the order parameter, and evaluate this integral using umbrella sampling. We validate the method by computing free energies of both face-centered cubic (FCC) and hexagonal close-packed (HCP) hard sphere crystals with a precision of $10^{-5}k_BT$ per particle. In developing the numerical method, we find several scaling properties of crystalline solids in the thermodynamic limit. Using these scaling properties, we derive an explicit asymptotic formula for the free energy per particle in the thermodynamic limit. In addition, we describe several changes of coordinates that can be used to separate internal degrees of freedom from external, translational degrees of freedom. The second part of the dissertation focuses on engineering idealized physical devices that work as Maxwell's demon. We describe two autonomous mechanical devices that extract energy from a single heat bath and convert it into work, while writing information onto memory registers. Additionally, both devices can operate as Landauer's eraser, namely they can erase information from a memory register, while energy is dissipated into the heat bath. The phase diagrams and the efficiencies of the two models are solved and analyzed. These two models provide concrete physical illustrations of the thermodynamic consequences of information processing.

Coupling of magnetic, strain and electric polarization fields in the structure of multiferroic material

The comprehensive study on the coupling of magnetism, electrical polarization and the crystalline lattice with the off-stoichiometric effects in self-doped multiferroic hexagonal h-LuMnxO3±δ (0.92≤x≤1.12) ceramic oxides was carried out for the PhD work. There is a complex coupling of the three ferroic degrees. The cancelation of the magnetic moments of ions in the antiferromagnetic order, electric polarization with specific vortex/antivortex topology and lattice properties have pushed researchers to find out ways to disclose the underlying physics and chemistry of magneto-electric and magneto-elastic couplings of h-RMnO3 multiferroic materials. In this research work, self-doping of Lu-sites or Mn-sites of h-LuMnxO3±δ ceramics prepared via solid state route was done to pave a way for deeper understanding of the antiferromagnetic transition, the weak ferromagnetism often reported in the same crystalline lattices and the ferroelectric properties coupled to the imposed lattice changes. Accordingly to the aim of the PhD thesis, the objectives set for the sintering study in the first chapter on experimental results were two. First, study of sintering off-stoichiometric samples within conditions reported in the bibliography and also extracted from the phase diagrams of the LuMnxO3±δ, with a multiple firings ending with a last high temperature step at 1300ºC for 24 hours. Second, explore longer annealing times of up to 240 hours at the fixed temperature of 1300 ºC in a search for improving the properties of the solid solution under study. All series of LuMnxO3±δ ceramics for each annealing time were characterized to tentatively build a framework enabling comparison of measured properties with results of others available in literature. XRD and Rietveld refinement of data give the evolution the lattice parameters as a function to x. Shrinkage of the lattice parameters with increasing x values was observed, the stability limit of the solid solution being determined by analysis of lattice parameters. The evolution of grain size and presence of secondary phases have been investigated by means of TEM, SEM, EDS and EBSD techniques. The dependencies of grain growth and regression of secondary phases on composition x and time were further characterized. Magnetic susceptibility of samples and magnetic irreversibility were extensively examined in the present work. The dependency of magnetic susceptibility, Neel ordering transition and important magnetic parameters are determined and compared to observation in other multiferroics in the following chapter of the thesis. As a tool of high sensitivity to detect minor traces of the secondary phase hausmannite, magnetic measurements are suggested for cross-checking of phase diagrams. Difficulty of previous studies on interpreting the magnetic anomaly below 43 K in h-RMnO3 oxides was discussed and assigned to the Mn3O4 phase, with supported of the electron microscopy. Magneto-electric coupling where AFM ordering is coupled to dielectric polarization is investigated as a function of x and of sintering condition via frequency and temperature dependent complex dielectric constant measurements in the final chapter of the thesis. Within the limits of solid solubility, the crystalline lattice of off-stoichiometric ceramics was shown to preserve the magneto-electric coupling at TN. It represents the first research work on magneto-electric coupling modified by vacancy doping to author’s knowledge. Studied lattices would reveal distortions at the atomic scale imposed by local changes of x dependent on sintering conditions which were widely inspected by using TEM/STEM methods, complemented with EDS and EELS spectroscopy all together to provide comprehensive information on cross coupling of distortions, inhomogeneity and electronic structure assembled and discussed in a specific chapter. Internal interfaces inside crystalline grains were examined. Qualitative explanations of the measured magnetic and ferroelectric properties were established in relation to observed nanoscale features of h-LuMnxO3±δ ceramics. Ferroelectric domains and topological defects are displayed both in TEM and AFM/PFM images, the later technique being used to look at size, distribution and switching of ferroelectric domains influenced by vacancy doping at the micron scale bridging to complementary TEM studies on the atomic structure of ferroelectric domains. In support to experimental study, DFT simulations using Wien2K code have been carried out in order to interpret the results of EELS spectra of O K-edge and to obtain information on the cation hybridization to oxygen ions. The L3,2 edges of Mn is used to access the oxidation state of the Mn ions inside crystalline grains. In addition, rehybridization driven ferroelectricity is also evaluated by comparing the partial density of states of the orbitals of all ions of the samples, also the polarization was calculated and correlated to the off-stoichiometric effect.

Extraction and purification of theobromine using ionic liquids

Theobromine is an alkaloid present in cocoa and it is used in the treatment of atherosclerosis, hypertension, angina, among others. Due to its importance, the aim of this work consists on the development of an efficient and sustainable technology for the extraction of theobromine from cocoa beans. For the development of a purification technique for theobromine extracted from cocoa, aqueous biphasic systems (ABS) composed of ionic liquids (ILs) were initially studied to infer on the most promising systems. Cholinium-based ILs, based on a non-toxic and biocompatible cation, were used combined with two polymers (PPG 400 and PEG 400) and an inorganic salt (K3PO4). The respective phase diagrams at 298 K and atmospheric pressure were determined, as well as their extraction efficiencies for theobromine. The results obtained indicate that K3PO4 has a greater ability to induce the formation of ABS compared to PEG 400 and PPG 400. ABS consisting of K3PO4 also have a high potential for the extraction of theobromine, with extraction efficiencies ranging between 96.4 and 99.9 %. Based on the most promising ILs for the purification step, they were further used in aqueous solution to extract theobromine from cocoa beans, with extraction yields ranging between 4.5% and 6.5 wt%. Finally, ABS were applied to the aqueous solutions containing theobromine from the cocoa extract, with extraction efficiencies ranging between 96.7 and 99.0%.

Conditions for observing IR chemiluminescence and the role of condensation in metal oxidation reactions

Infrared chemiluminescence (IRCL) studies of cw metal oxidation reactions wherein metal atoms entrained in a carrier gas were mixed with an oxidizer by means of a nozzle system are described. One goal of the work was to determine the vibrational distribution of the product molecule produced by the chemical reaction. In order to observe IRCL it was important to operate the system at the appropriate P-T point in the phase diagram of both the metal and metal salt, otherwise rapid condensation quenched any IRCL that was present. If the nucleation rate was greater 1010 3 than ~ cm-sec-I, then only "black body" radiation could be seen from the reaction. Most of the studies were on the Li/I2 system which is unique in that the phase diagrams of Li and LiI in the P-T ranges of interest are almost identical. This property permitted a relatively easy control with respect to condensation and the measurement of IRCL in the 10-28 um range for the excited LiI molecule.

Uso da glicerina oriunda da síntese do biodiesel em sistemas microemulsionados visando aplicação de herbicidas

The program PROBIODIESEL from the Ministry of Science and Technology has substantially increased glycerine, obtained as a sub-product of biodiesel production process, making it necessary to seek alternatives for the use of this co-product. On the other hand, herbicides although play a role of fundamental importance in the agricultural production system in force, have been under growing concern among the various segments of society because of their potential environmental risk. In this work, we used glycerin in microemulsion systems for application of herbicides, to improve efficiency and lower environmental pollution caused by the loss of those products to the environment. To obtain the systems of microemulsinados were used Unitol L90 NP and Renex 40 as surfactants, butanol as co-surfactant, coconut oil as oil phase and aqueous phase as we used solutions of glycerin + water. Through the determination of phase diagrams, the microemulsion region was found in the system E (L90 Unitol, coconut oil and glycerin + water 1:1). Three points were chosen to the aqueous phase rich in characterization and application in the solubilization of glyphosate and atrazine. Three experiments were performed in Horta, Department of Plant Sciences, Plant Science Sector, UFERSA, Mossoró-RN. The first experiment was conducted in randomized complete blocks with 20 treatments and four replications. The treatments consisted of five doses of the herbicide glyphosate (0.0, 0.45, 0.9, 1.35 and 1.8 L ha-1) diluted with four sauces: C1, C2, C3 (microemulsions) and C4 (water). The phytotoxicity of Brachiaria brizantha was measured at 7, 14, 28 and 60 DAA (days after application). At 60 DAA, we evaluated the biomass of plants. The second experiment was developed in randomized complete blocks with 20 treatments and four repetitions. The treatments consisted of five doses of the herbicide atrazine (0.0, 0.4, 0.8, 1.6 and 2.4 L ha-1) diluted with four sauces: C1, C2, C3 (microemulsions) and C4 (water). The phytotoxicity on Zea mays and Talinum paniculatum was evaluated at 2, 7, 20 DAA. The experiment III was developed in randomized complete blocks with 16 treatments and three repetitions. The treatments consisted of 16 combinations among the constituents of the microemulsion: Unitol L90 surfactant (0.0, 1.66, 5.0, 15 %) and glycerin (0.0, 4.44, 13.33 and 40.0 %). The phytotoxicity on Zea mays was evaluated at 1, 7 and 14 DAA. At 14 DAA, we evaluated the biomass of plants. The control plants using the microemulsions was lower than in the water due to the poisoning caused by the initial microemulsions in the leaves of the plants, a fact that hinders the absorption and translocation of the herbicide. There was no toxicity in Zea mays plants caused by the herbicide, however, were highly intoxicated by microemulsions. T. paniculatum was better controlled in spraying with the microemulsions, regardless of the dose of the herbicide. The glycerine did not cause plant damage. Higher poisoning the plants are caused by tensoactive Unitol L90 and higher rates occur with the use of higher concentrations of surfactant and glycerin, or microemulsion. The microemulsions used hampered the action of glyphosate in controlling B. brizantha and caused severe poisoning in corn, and these poisonings attributed mainly to the action of surfactant

Effects of solid poly (ethylene glycols) addition to the solutions of aniline or N,N-dimethylaniline with water:experimental measurements and modelling

In this work, the liquid-liquid and solid-liquid phase behaviour of ten aqueous pseudo-binary and three binary systems containing polyethylene glycol (PEG) 2050, polyethylene glycol 35000, aniline, N,N-dimethylaniline and water, in the temperature range 298.15-350.15 K and at ambient pressure of 0.1 MPa, was studied. The obtained temperature-composition phase diagrams showed that the only functional co-solvent was PEG2050 for aniline in water, while PEG35000 even showed a clear anti-solvent effect in the N,N-dimethylaniline aqueous system. The experimental solid-liquid equilibria (SLE) data have been correlated by the non-random two-liquid (NRTL) model, and the correlation results are in accordance with the experimental results.

STUDY ON FUSIBILITY DIAGRAMS CONTAINING LICL AND PHASE-DIAGRAM OF LICL-KCL-LIF TERNARY BY DIRECT CHLORINATION

In order to avoid the hygroscopicity of LiCl specimem, the method of directly chlorinating Li_2CO_3 with NH_4Cl was successfully introducing into the thermal analysis of the system containing LiCl. The three fusibility diagram of LiCl-KCl, LiClNaCl, LiCl-LiF were determined using the method. The results are in agreement wish the values reported in the literature, and phase diagram of LiCl-KCl-LiF ternary system was constructed based on these results. Temperature of the ternary eutectic, composed of 57.3mol％...