Atomic resolution studies of oxide superlattices and ultrathin films

The ability to grow ultrathin films layer-by-layer with well-defined epitaxial relationships has allowed research groups worldwide to grow a range of artificial films and superlattices, first for semiconductors, and now with oxides. In the oxides thin film research community, there have been concerted efforts recently to develop a number of epitaxial oxide systems grown on single crystal oxide substrates that display a wide variety of novel interfacial functionality, such as enhanced ferromagnetic ordering, increased charge carrier density, increased optical absorption, etc, at interfaces. The magnitude of these novel properties is dependent upon the structure of thin films, especially interface sharpness, intermixing, defects, and strain, layering sequence in the case of superlattices and the density of interfaces relative to the film thicknesses. To understand the relationship between the interfacial thin film oxide atomic structure and its properties, atomic scale characterization is required. Transmission electron microscopy (TEM) offers the ability to study interfaces of films at high resolution. Scanning transmission electron microscopy (STEM) allows for real space imaging of materials with directly interpretable atomic number contrast. Electron energy loss spectroscopy (EELS), together with STEM, can probe the local chemical composition as well as local electronic states of transition metals and oxygen. Both techniques have been significantly improved by aberration correctors, which reduce the probe size to 1 Å, or less. Aberration correctors have thus made it possible to resolve individual atomic columns, and possibly probe the electronic structure at atomic scales. Separately, using electron probe forming lenses, structural information such as the crystal structure, strain, lattice mismatches, and superlattice ordering can be measured by nanoarea electron diffraction (NED). The combination of STEM, EELS, and NED techniques allows us to gain a fundamental understanding of the properties of oxide superlattices and ultrathin films and their relationship with the corresponding atomic and electronic structure. In this dissertation, I use the aforementioned electron microscopy techniques to investigate several oxide superlattice and ultrathin film systems. The major findings are summarized below. These results were obtained with stringent specimen preparation methods that I developed for high resolution studies, which are described in Chapter 2. The essential materials background and description of electron microscopy techniques are given in Chapter 1 and 2. In a LaMnO3-SrMnO3 superlattice, we demonstrate the interface of LaMnO3-SrMnO3 is sharper than the SrMnO3-LaMnO3 interface. Extra spectral weights in EELS are confined to the sharp interface, whereas at the rougher interface, the extra states are either not present or are not confined to the interface. Both the structural and electronic asymmetries correspond to asymmetric magnetic ordering at low temperature. In a short period LaMnO3-SrTiO3 superlattice for optical applications, we discovered a modified band structure in SrTiO3 ultrathin films relative to thick films and a SrTiO3 substrate, due to charge leakage from LaMnO3 in SrTiO3. This was measured by chemical shifts of the Ti L and O K edges using atomic scale EELS. The interfacial sharpness of LaAlO3 films grown on SrTiO3 was investigated by the STEM/EELS technique together with electron diffraction. This interface, when prepared under specific conditions, is conductive with high carrier mobility. Several suggestions for the conductive interface have been proposed, including a polar catastrophe model, where a large built-in electric field in LaAlO3 films results in electron charge transfer into the SrTiO3 substrate. Other suggested possibilities include oxygen vacancies at the interface and/or oxygen vacancies in the substrate. The abruptness of the interface as well as extent of intermixing has not been thoroughly investigated at high resolution, even though this can strongly influence the electrical transport properties. We found clear evidence for cation intermixing through the LaAlO3-SrTiO3 interface with high spatial resolution EELS and STEM, which contributes to the conduction at the interface. We also found structural defects, such as misfit dislocations, which leads to increased intermixing over coherent interfaces.

Soft computing based parameter identification in pavements and geomechanical systems

Accurate estimation of road pavement geometry and layer material properties through the use of proper nondestructive testing and sensor technologies is essential for evaluating pavement’s structural condition and determining options for maintenance and rehabilitation. For these purposes, pavement deflection basins produced by the nondestructive Falling Weight Deflectometer (FWD) test data are commonly used. The nondestructive FWD test drops weights on the pavement to simulate traffic loads and measures the created pavement deflection basins. Backcalculation of pavement geometry and layer properties using FWD deflections is a difficult inverse problem, and the solution with conventional mathematical methods is often challenging due to the ill-posed nature of the problem. In this dissertation, a hybrid algorithm was developed to seek robust and fast solutions to this inverse problem. The algorithm is based on soft computing techniques, mainly Artificial Neural Networks (ANNs) and Genetic Algorithms (GAs) as well as the use of numerical analysis techniques to properly simulate the geomechanical system. A widely used pavement layered analysis program ILLI-PAVE was employed in the analyses of flexible pavements of various pavement types; including full-depth asphalt and conventional flexible pavements, were built on either lime stabilized soils or untreated subgrade. Nonlinear properties of the subgrade soil and the base course aggregate as transportation geomaterials were also considered. A computer program, Soft Computing Based System Identifier or SOFTSYS, was developed. In SOFTSYS, ANNs were used as surrogate models to provide faster solutions of the nonlinear finite element program ILLI-PAVE. The deflections obtained from FWD tests in the field were matched with the predictions obtained from the numerical simulations to develop SOFTSYS models. The solution to the inverse problem for multi-layered pavements is computationally hard to achieve and is often not feasible due to field variability and quality of the collected data. The primary difficulty in the analysis arises from the substantial increase in the degree of non-uniqueness of the mapping from the pavement layer parameters to the FWD deflections. The insensitivity of some layer properties lowered SOFTSYS model performances. Still, SOFTSYS models were shown to work effectively with the synthetic data obtained from ILLI-PAVE finite element solutions. In general, SOFTSYS solutions very closely matched the ILLI-PAVE mechanistic pavement analysis results. For SOFTSYS validation, field collected FWD data were successfully used to predict pavement layer thicknesses and layer moduli of in-service flexible pavements. Some of the very promising SOFTSYS results indicated average absolute errors on the order of 2%, 7%, and 4% for the Hot Mix Asphalt (HMA) thickness estimation of full-depth asphalt pavements, full-depth pavements on lime stabilized soils and conventional flexible pavements, respectively. The field validations of SOFTSYS data also produced meaningful results. The thickness data obtained from Ground Penetrating Radar testing matched reasonably well with predictions from SOFTSYS models. The differences observed in the HMA and lime stabilized soil layer thicknesses observed were attributed to deflection data variability from FWD tests. The backcalculated asphalt concrete layer thickness results matched better in the case of full-depth asphalt flexible pavements built on lime stabilized soils compared to conventional flexible pavements. Overall, SOFTSYS was capable of producing reliable thickness estimates despite the variability of field constructed asphalt layer thicknesses.

Three-dimensional kinetic simulations of active suspensions: effect of chemotaxis and steric interaction

In this work, the effects of chemotaxis and steric interactions in active suspensions are analyzed by extending the kinetic model proposed by Saintillan and Shelley [1, 2]. In this model, a conservation equation for the active particle configuration is coupled to the Stokes equation for the flow arising from the force dipole exerted by the particles on the fluid. The fluid flow equations are solved spectrally and the conservation equation is solved by second-order finite differencing in space and second-order Adams-Bashforth time marching. First, the dynamics in suspensions of oxytactic run-and-tumble bacteria confined in thin liquid films surrounded by air is investigated. These bacteria modify their tumbling behavior by making temporal comparisons of the oxygen concentration, and, on average, swim towards high concentrations of oxygen. The kinetic model proposed by Saintillan and Shelley [1, 2] is modified to include run-and-tumble effects and oxygentaxis. The spatio-temporal dynamics of the oxygen and bacterial concentration are analyzed. For small film thicknesses, there is a weak migration of bacteria to the boundaries, and the oxygen concentration is high inside the film as a result of diffusion; both bacterial and oxygen concentrations quickly reach steady states. Above a critical film thickness (approximately 200 micron), a transition to chaotic dynamics is observed and is characterized by turbulent-like 3D motion, the formation of bacterial plumes, enhanced oxygen mixing and transport into the film, and hydrodynamic velocities of magnitudes up to 7 times the single bacterial swimming speed. The simulations demonstrate that the combined effects of hydrodynamic interactions and oxygentaxis create collective three-dimensional instabilities which enhances oxygen availability for the bacteria. Our simulation results are consistent with the experimental findings of Sokolov et al. [3], who also observed a similar transition with increasing film thickness. Next, the dynamics in concentrated suspensions of active self-propelled particles in a 3D periodic domain are analyzed. We modify the kinetic model of Saintillan and Shelley [1, 2] by including an additional nematic alignment torque proportional to the local concentration in the equation for the rotational velocity of the particles, causing them to align locally with their neighbors (Doi and Edwards [4]). Large-scale three- dimensional simulations show that, in the presence of such a torque both pusher and puller suspensions are unstable to random fluctuations and are characterized by highly nematic structures. Detailed measures are defined to quantify the degree and direction of alignment, and the effects of steric interactions on pattern formation will be presented. Our analysis shows that steric interactions have a destabilizing effect in active suspensions.