Fluctuations in and out of equilibrium: Thermalization, quantum measurements and Coulomb disorder

This purely theoretical thesis covers aspects of two contemporary research fields: the non-equilibrium dynamics in quantum systems and the electronic properties of three-dimensional topological insulators. In the first part we investigate the non-equilibrium dynamics in closed quantum systems. Thanks to recent technologies, especially from the field of ultracold quantum gases, it is possible to realize such systems in the laboratory. The focus is on the influence of hydrodynamic slow modes on the thermalization process. Generic systems in equilibrium, either classical or quantum, in equilibrium are described by thermodynamics. This is characterized by an ensemble of maximal entropy, but constrained by macroscopically conserved quantities. We will show that these conservation laws slow down thermalization and the final equilibrium state can be approached only algebraically in time. When the conservation laws are violated thermalization takes place exponential in time. In a different study we calculate probability distributions of projective quantum measurements. Newly developed quantum microscopes provide the opportunity to realize new measurement protocols which go far beyond the conventional measurements of correlation functions. The second part of this thesis is dedicated to a new class of materials known as three-dimensional topological insulators. Also here new experimental techniques have made it possible to fabricate these materials to a high enough quality that their topological nature is revealed. However, their transport properties are not fully understood yet. Motivated by unusual experimental results in the optical conductivity we have investigated the formation and thermal destruction of spatially localized electron- and hole-doped regions. These are caused by charged impurities which are introduced into the material in order to make the bulk insulating. Our theoretical results are in agreement with the experiment and can explain the results semi-quantitatively. Furthermore, we study emergent lengthscales in the bulk as well as close to the conducting surface.

DIRECT NUMERICAL SIMULATION OF INCOMPRESSIBLE MULTIPHASE FLOW WITH PHASE CHANGE

Phase change problems arise in many practical applications such as air-conditioning and refrigeration, thermal energy storage systems and thermal management of electronic devices. The physical phenomenon in such applications are complex and are often difficult to be studied in detail with the help of only experimental techniques. The efforts to improve computational techniques for analyzing two-phase flow problems with phase change are therefore gaining momentum. The development of numerical methods for multiphase flow has been motivated generally by the need to account more accurately for (a) large topological changes such as phase breakup and merging, (b) sharp representation of the interface and its discontinuous properties and (c) accurate and mass conserving motion of the interface. In addition to these considerations, numerical simulation of multiphase flow with phase change introduces additional challenges related to discontinuities in the velocity and the temperature fields. Moreover, the velocity field is no longer divergence free. For phase change problems, the focus of developmental efforts has thus been on numerically attaining a proper conservation of energy across the interface in addition to the accurate treatment of fluxes of mass and momentum conservation as well as the associated interface advection. Among the initial efforts related to the simulation of bubble growth in film boiling applications the work in \cite{Welch1995} was based on the interface tracking method using a moving unstructured mesh. That study considered moderate interfacial deformations. A similar problem was subsequently studied using moving, boundary fitted grids \cite{Son1997}, again for regimes of relatively small topological changes. A hybrid interface tracking method with a moving interface grid overlapping a static Eulerian grid was developed \cite{Juric1998} for the computation of a range of phase change problems including, three-dimensional film boiling \cite{esmaeeli2004computations}, multimode two-dimensional pool boiling \cite{Esmaeeli2004} and film boiling on horizontal cylinders \cite{Esmaeeli2004a}. The handling of interface merging and pinch off however remains a challenge with methods that explicitly track the interface. As large topological changes are crucial for phase change problems, attention has turned in recent years to front capturing methods utilizing implicit interfaces that are more effective in treating complex interface deformations. The VOF (Volume of Fluid) method was adopted in \cite{Welch2000} to simulate the one-dimensional Stefan problem and the two-dimensional film boiling problem. The approach employed a specific model for mass transfer across the interface involving a mass source term within cells containing the interface. This VOF based approach was further coupled with the level set method in \cite{Son1998}, employing a smeared-out Heaviside function to avoid the numerical instability related to the source term. The coupled level set, volume of fluid method and the diffused interface approach was used for film boiling with water and R134a at the near critical pressure condition \cite{Tomar2005}. The effect of superheat and saturation pressure on the frequency of bubble formation were analyzed with this approach. The work in \cite{Gibou2007} used the ghost fluid and the level set methods for phase change simulations. A similar approach was adopted in \cite{Son2008} to study various boiling problems including three-dimensional film boiling on a horizontal cylinder, nucleate boiling in microcavity \cite{lee2010numerical} and flow boiling in a finned microchannel \cite{lee2012direct}. The work in \cite{tanguy2007level} also used the ghost fluid method and proposed an improved algorithm based on enforcing continuity and divergence-free condition for the extended velocity field. The work in \cite{sato2013sharp} employed a multiphase model based on volume fraction with interface sharpening scheme and derived a phase change model based on local interface area and mass flux. Among the front capturing methods, sharp interface methods have been found to be particularly effective both for implementing sharp jumps and for resolving the interfacial velocity field. However, sharp velocity jumps render the solution susceptible to erroneous oscillations in pressure and also lead to spurious interface velocities. To implement phase change, the work in \cite{Hardt2008} employed point mass source terms derived from a physical basis for the evaporating mass flux. To avoid numerical instability, the authors smeared the mass source by solving a pseudo time-step diffusion equation. This measure however led to mass conservation issues due to non-symmetric integration over the distributed mass source region. The problem of spurious pressure oscillations related to point mass sources was also investigated by \cite{Schlottke2008}. Although their method is based on the VOF, the large pressure peaks associated with sharp mass source was observed to be similar to that for the interface tracking method. Such spurious fluctuation in pressure are essentially undesirable because the effect is globally transmitted in incompressible flow. Hence, the pressure field formation due to phase change need to be implemented with greater accuracy than is reported in current literature. The accuracy of interface advection in the presence of interfacial mass flux (mass flux conservation) has been discussed in \cite{tanguy2007level,tanguy2014benchmarks}. The authors found that the method of extending one phase velocity to entire domain suggested by Nguyen et al. in \cite{nguyen2001boundary} suffers from a lack of mass flux conservation when the density difference is high. To improve the solution, the authors impose a divergence-free condition for the extended velocity field by solving a constant coefficient Poisson equation. The approach has shown good results with enclosed bubble or droplet but is not general for more complex flow and requires additional solution of the linear system of equations. In current thesis, an improved approach that addresses both the numerical oscillation of pressure and the spurious interface velocity field is presented by featuring (i) continuous velocity and density fields within a thin interfacial region and (ii) temporal velocity correction steps to avoid unphysical pressure source term. Also I propose a general (iii) mass flux projection correction for improved mass flux conservation. The pressure and the temperature gradient jump condition are treated sharply. A series of one-dimensional and two-dimensional problems are solved to verify the performance of the new algorithm. Two-dimensional and cylindrical film boiling problems are also demonstrated and show good qualitative agreement with the experimental observations and heat transfer correlations. Finally, a study on Taylor bubble flow with heat transfer and phase change in a small vertical tube in axisymmetric coordinates is carried out using the new multiphase, phase change method.

Single particle motion of polymeric, colloidal, and biological materials

Small particles and their dynamics are of widespread interest due both to their unique properties and their ubiquity. Here, we investigate several classes of small particles: colloids, polymers, and liposomes. All these particles, due to their size on the order of microns, exhibit significant similarity in that they are large enough to be visualized in microscopes, but small enough to be significantly influenced by thermal (or Brownian) motion. Further, similar optical microscopy and experimental techniques are commonly employed to investigate all these particles. In this work, we develop single particle tracking techniques, which allow thorough characterization of individual particle dynamics, observing many behaviors which would be overlooked by methods which time or ensemble average. The various particle systems are also similar in that frequently, the signal-to-noise ratio represented a significant concern. In many cases, development of image analysis and particle tracking methods optimized to low signal-to-noise was critical to performing experimental observations. The simplest particles studied, in terms of their interaction potentials, were chemically homogeneous (though optically anisotropic) hard-sphere colloids. Using these spheres, we explored the comparatively underdeveloped conjunction of translation and rotation and particle hydrodynamics. Developing off this, the dynamics of clusters of spherical colloids were investigated, exploring how shape anisotropy influences the translation and rotation respectively. Transitioning away from uniform hard-sphere potentials, the interactions of amphiphilic colloidal particles were explored, observing the effects of hydrophilic and hydrophobic interactions upon pattern assembly and inter-particle dynamics. Interaction potentials were altered in a different fashion by working with suspensions of liposomes, which, while homogeneous, introduce the possibility of deformation. Even further degrees of freedom were introduced by observing the interaction of particles and then polymers within polymer suspensions or along lipid tubules. Throughout, while examination of the trajectories revealed that while by some measures, the averaged behaviors accorded with expectation, often closer examination made possible by single particle tracking revealed novel and unexpected phenomena.

Computational modeling of laser-induced delamination testing of thin films

Thin film adhesion often determines microelectronic device reliability and it is therefore essential to have experimental techniques that accurately and efficiently characterize it. Laser-induced delamination is a novel technique that uses laser-generated stress waves to load thin films at high strain rates and extract the fracture toughness of the film/substrate interface. The effectiveness of the technique in measuring the interface properties of metallic films has been documented in previous studies. The objective of the current effort is to model the effect of residual stresses on the dynamic delamination of thin films. Residual stresses can be high enough to affect the crack advance and the mode mixity of the delimitation event, and must therefore be adequately modeled to make accurate and repeatable predictions of fracture toughness. The equivalent axial force and bending moment generated by the residual stresses are included in a dynamic, nonlinear finite element model of the delaminating film, and the impact of residual stresses on the final extent of the interfacial crack, the relative contribution of shear failure, and the deformed shape of the delaminated film is studied in detail. Another objective of the study is to develop techniques to address issues related to the testing of polymeric films. These type of films adhere well to silicon and the resulting crack advance is often much smaller than for metallic films, making the extraction of the interface fracture toughness more difficult. The use of an inertial layer which enhances the amount of kinetic energy trapped in the film and thus the crack advance is examined. It is determined that the inertial layer does improve the crack advance, although in a relatively limited fashion. The high interface toughness of polymer films often causes the film to fail cohesively when the crack front leaves the weakly bonded region and enters the strong interface. The use of a tapered pre-crack region that provides a more gradual transition to the strong interface is examined. The tapered triangular pre-crack geometry is found to be effective in reducing the stresses induced thereby making it an attractive option. We conclude by studying the impact of modifying the pre-crack geometry to enable the testing of multiple polymer films.

Mechanism of pH-dependent activation of the sodium-proton antiporter NhaA

Data files to accompany the article in Nature Communications, in press.

Multiple sulfur isotope fractionations in inorganic aqueous systems

New constraints on isotope fractionation factors in inorganic aqueous sulfur systems based on theoretical and experimental techniques relevant to studies of the sulfur cycle in modern environments and the geologic rock record are presented in this dissertation. These include theoretical estimations of equilibrium isotope fractionation factors utilizing quantum mechanical software and a water cluster model approach for aqueous sulfur compounds that span the entire range of oxidation state for sulfur. These theoretical calculations generally reproduce the available experimental determinations from the literature and provide new constraints where no others are available. These theoretical calculations illustrate in detail the relationship between sulfur bonding environment and the mass dependence associated with equilibrium isotope exchange reactions involving all four isotopes of sulfur. I additionally highlight the effect of isomers of protonated compounds (compounds with the same chemical formula but different structure, where protons are bound to either sulfur or oxygen atoms) on isotope partitioning in the sulfite (S4+) and sulfoxylate (S2+) systems, both of which are key intermediates in oxidation-reduction processes in the sulfur cycle. I demonstrate that isomers containing the highest degree of coordination around sulfur (where protonation occurs on the sulfur atom) have a strong influence on isotopic fractionation factors, and argue that isomerization phenomenon should be considered in models of the sulfur cycle. Additionally, experimental results of the reaction rates and isotope fractionations associated with the chemical oxidation of aqueous sulfide are presented. Sulfide oxidation is a major process in the global sulfur cycle due largely to the sulfide-producing activity of anaerobic microorganisms in organic-rich marine sediments. These experiments reveal relationships between isotope fractionations and reaction rate as a function of both temperature and trace metal (ferrous iron) catalysis that I interpret in the context of the complex mechanism of sulfide oxidation. I also demonstrate that sulfide oxidation is a process associated with a mass dependence that can be described as not conforming to the mass dependence typically associated with equilibrium isotope exchange. This observation has implications for the inclusion of oxidative processes in environmental- and global-scale models of the sulfur cycle based on the mass balance of all four isotopes of sulfur. The contents of this dissertation provide key reference information on isotopic fractionation factors in aqueous sulfur systems that will have far-reaching applicability to studies of the sulfur cycle in a wide variety of natural settings.

Path independent limiting criteria in sheet metal forming

The Forming Limit Diagram (FLD) is a conventional failure diagram to estimate necking limits of sheet metal for proportional loading conditions. Previous studies reveal that the FLD is not suitable for predicting the influence of nonlinear strain paths. This paper presents methodical comparison among all common available strain path independent strain/stress based limiting criteria. All the strain path independent strain based limiting criteria (Traditional Failure Diagram (TFD), Extended Forming Limit Diagram (XFLD), Extended Stress Ratio Based Forming Limit Diagram (ESRFLD), and Polar Effective Plastic StrainDiagram (PEPSD)) and stress based limiting criteria (Traditional Stress based Failure Diagram (TFSD), Stress Based Forming Limit Diagram (FLSD), Stress Ratio and Stress Based Forming Limit Diagram (SRFLSD), Extended Stress Based Forming Limit Diagram (XFLSD), and Polar Effective Stress Diagram (PESSD)) are approximately path-independent for smaller amount of pre-straining and path dependent for large pre-straining conditions. From advance image correlation technique precisely determination of local strains near necked area is possible today. However direct determination of local stresses near necked area is not possible. Therefore, local stresses and equivalent stress are determined by employing both yield criterion and strain-hardening law. Similarly equivalent strain is calculated by the use of yield criterion. Therefore, the choice of yield criterion has an impact on the results for TFD, XFLD, ESRFLD and PEPSD. However, selections of both yield criterion and strain-hardening law have substantial influence on the results for TFSD, FLSD, SRFLSD, XFLSD and PESSD. The inherent calculation error can be minimized by generation of experimental data for each material and then selection of representable yield criterion and strain-hardening law. Improvement of experimental techniques and generation of rigorous material data bank in various strain paths may help researchers to diagnose and resolve the issue. TFD, TFSD and XFLSD have inherent variables to take care the effect of through thickness stress, however rigorous experimental verification is needed before the field application.

Mechanical and Electrical Properties of Single-walled Carbon Nanotubes Synthesized by Chemical Vapor Deposition

Despite the tremendous application potentials of carbon nanotubes (CNTs) proposed by researchers in the last two decades, efficient experimental techniques and methods are still in need for controllable production of CNTs in large scale, and for conclusive characterizations of their properties in order to apply CNTs in high accuracy engineering. In this dissertation, horizontally well-aligned high quality single-walled carbon nanotubes (SWCNTs) have been successfully synthesized on St-cut quartz substrate by chemical vapor deposition (CVD). Effective radial moduli (Eradial) of these straight SWCNTs have been measured by using well-calibrated tapping mode and contact mode atomic force microscopy (AFM). It was found that the measured Eradial decreased from 57 to 9 GPa as the diameter of the SWCNTs increased from 0.92 to 1.91 nm. The experimental results were consistent with the recently reported theoretical simulation data. The method used in this mechanical property test can be easily applied to measure the mechanical properties of other low-dimension nanostructures, such as nanowires and nanodots. The characterized sample is also an ideal platform for electrochemical tests. The electrochemical activities of redox probes Fe(CN)63-/4-, Ru(NH3)63+, Ru(bpy)32+ and protein cytochrome c have been studied on these pristine thin films by using aligned SWCNTs as working electrodes. A simple and high performance electrochemical sensor was fabricated. Flow sensing capability of the device has been tested for detecting neurotransmitter dopamine at physiological conditions with the presence of Bovine serum albumin. Good sensitivity, fast response, high stability and anti-fouling capability were observed. Therefore, the fabricated sensor showed great potential for sensing applications in complicated solution.

Experimental investigation of stratified combustion: Application of laser diagnostic techniques

In order to better understand the stratified combustion, the propagation of flame through stratified mixture field in laminar and turbulent flow conditions has been studied by using combined PIV/PLIF techniques. A great emphasis was placed on developing methods to improve the accuracy of local measurements of flame propagation. In particular, a new PIV approach has been developed to measure the local fresh gas velocity near preheat zone of flame front. To improve the resolution of measurement, the shape of interrogation window has been continuously modified based on the local flame topology and gas expansion effect. Statistical analysis of conditioned local measurements by the local equivalence ratio of flames allows the characterization of the properties of flame propagation subjected to the mixture stratification in laminar and turbulent flows, especially the highlight of the memory effect.

Experimental investigation of turbine stator well rim seal, reingestion and interstage seal flows using gas concentration techniques and displacement measurements

Gas turbine engine performance requires effective and reliable internal cooling over the duty cycle of the engine. Life predictions for rotating components subject to the main gas path temperatures are vital. This demands increased precision in the specification of the internal air system flows which provide turbine stator well cooling and sealing. This in turn requires detailed knowledge of the flow rates through rim seals and interstage labyrinth seals. Knowledge of seal movement and clearances at operating temperatures is of great importance when prescribing these flows. A test facility has been developed at the University of Sussex, incorporating a two stage turbine rated at 400 kW with an individual stage pressure ratio of 1.7:1. The mechanical design of the test facility allows internal cooling geometry to be rapidly re-configured, while cooling flow rates of between 0.71 CW, ENT and 1.46 C W, ENT, may be set to allow ingress or egress dominated cavity flows. The main annulus and cavity conditions correspond to in cavity rotational Reynolds numbers of 1.71×106< Reφ <1.93×106. Displacement sensors have been used to establish hot running seal clearances over a range of stator well flow conditions, allowing realistic flow rates to be calculated. Additionally, gas seeding techniques have been developed, where stator well and main annulus flow interactions are evaluated by measuring changes in gas concentration. Experiments have been performed which allow rim seal and re-ingestion flows to be quantified. It will be shown that this work develops the measurement of stator well cooling flows and provides data suitable for the validation of improved thermo-mechanical and CFD codes, beneficial to the engine design process. Copyright © 2011 by Rolls-Royce plc.

Experimental investigation of negative pressure intrusion techniques of acetabular cementation in total hip arthroplasty

The main mode of failure of the acetabular component in total hip arthroplasty is aseptic loosening. Successive generations of cementation techniques have evolved to alleviate this problem.

Switched Diversity Techniques for Indoor Off-body Communications Channels: An Experimental Analysis and Modeling

This paper investigates the potential improvement in signal reliability for indoor off-body communications channels operating at 5.8 GHz using switched diversity techniques. In particular we investigate the performance of switch-and-stay combining (SSC), switch-and-examine combining (SEC) and switch-and-examine combining with post-examining selection (SECps) schemes which utilize multiple spatially separated antennas at the base station. During the measurements a test subject, wearing an antenna on his chest, performed a number of walking movements towards and then away from a uniform linear array. It was found that all of the considered diversity schemes provided a worthwhile signal improvement. However, the performance of the diversity systems varied according to the switching threshold that was adopted. To model the fading envelope observed at the output of each of the combiners, we have applied diversity specific equations developed under the assumption of Nakagami-$m$ fading. As a measure of the goodness-of-fit, the Kullback-Leibler divergence between the empirical and theoretical probability density functions (PDFs) was calculated and found to be close to 0. To assist with the interpretation of the goodness-of-fit achieved in this study, the standard deviation, $\sigma$, of a zero-mean, $\sigma^2$ variance Gaussian PDF used to approximate a zero-mean, unit variance Gaussian PDF is also presented. These were generally quite close to 1 indicating that the theoretical models provided an adequate fit to the measured data.