Mathematical modelling of LiFePO4 cathodes

LiFePO4 is a commercially available battery material with good theoretical discharge capacity, excellent cycle life and increased safety compared with competing Li-ion chemistries. It has been the focus of considerable experimental and theoretical scrutiny in the past decade, resulting in LiFePO4 cathodes that perform well at high discharge rates. This scrutiny has raised several questions about the behaviour of LiFePO4 material during charge and discharge. In contrast to many other battery chemistries that intercalate homogeneously, LiFePO4 can phase-separate into highly and lowly lithiated phases, with intercalation proceeding by advancing an interface between these two phases. The main objective of this thesis is to construct mathematical models of LiFePO4 cathodes that can be validated against experimental discharge curves. This is in an attempt to understand some of the multi-scale dynamics of LiFePO4 cathodes that can be difficult to determine experimentally. The first section of this thesis constructs a three-scale mathematical model of LiFePO4 cathodes that uses a simple Stefan problem (which has been used previously in the literature) to describe the assumed phase-change. LiFePO4 crystals have been observed agglomerating in cathodes to form a porous collection of crystals and this morphology motivates the use of three size-scales in the model. The multi-scale model developed validates well against experimental data and this validated model is then used to examine the role of manufacturing parameters (including the agglomerate radius) on battery performance. The remainder of the thesis is concerned with investigating phase-field models as a replacement for the aforementioned Stefan problem. Phase-field models have recently been used in LiFePO4 and are a far more accurate representation of experimentally observed crystal-scale behaviour. They are based around the Cahn-Hilliard-reaction (CHR) IBVP, a fourth-order PDE with electrochemical (flux) boundary conditions that is very stiff and possesses multiple time and space scales. Numerical solutions to the CHR IBVP can be difficult to compute and hence a least-squares based Finite Volume Method (FVM) is developed for discretising both the full CHR IBVP and the more traditional Cahn-Hilliard IBVP. Phase-field models are subject to two main physicality constraints and the numerical scheme presented performs well under these constraints. This least-squares based FVM is then used to simulate the discharge of individual crystals of LiFePO4 in two dimensions. This discharge is subject to isotropic Li+ diffusion, based on experimental evidence that suggests the normally orthotropic transport of Li+ in LiFePO4 may become more isotropic in the presence of lattice defects. Numerical investigation shows that two-dimensional Li+ transport results in crystals that phase-separate, even at very high discharge rates. This is very different from results shown in the literature, where phase-separation in LiFePO4 crystals is suppressed during discharge with orthotropic Li+ transport. Finally, the three-scale cathodic model used at the beginning of the thesis is modified to simulate modern, high-rate LiFePO4 cathodes. High-rate cathodes typically do not contain (large) agglomerates and therefore a two-scale model is developed. The Stefan problem used previously is also replaced with the phase-field models examined in earlier chapters. The results from this model are then compared with experimental data and fit poorly, though a significant parameter regime could not be investigated numerically. Many-particle effects however, are evident in the simulated discharges, which match the conclusions of recent literature. These effects result in crystals that are subject to local currents very different from the discharge rate applied to the cathode, which impacts the phase-separating behaviour of the crystals and raises questions about the validity of using cathodic-scale experimental measurements in order to determine crystal-scale behaviour.

Interfacial microstructure in Y1Ba2Cu3O7-δ high Tc superconductors

The microstructures of YBa2Cu3O7-δ ceramics prepared from freeze dried powders and containing an excess of CuO have been studied by analytical electron microscopy. Special attention has been paid to the interfacial microstructure. It was found that a liquid phase formed during sintering between 890°C and 920°C and this promoted grain growth and densification. Both clean grain boundaries and boundaries containing an amorphous intergranular film, which was rich in Cu, have been observed. Both CuO and BaCuO2 were present as secondary phases.

A nanometre-scale non-periodic structural variation in high temperature superconducting ceramics and the implications for properties

Non-periodic structural variation has been found in the high Tc cuprates, YBa2Cu3O7-x and Hg0.67Pb0.33Ba2Ca2Cu 3O8+δ, by image analysis of high resolution transmission electron microscope (HRTEM) images. We use two methods for analysis of the HRTEM images. The first method is a means for measuring the bending of lattice fringes at twin planes. The second method is a low-pass filter technique which enhances information contained by diffuse-scattered electrons and reveals what appears to be an interference effect between domains of differing lattice parameter in the top and bottom of the thin foil. We believe that these methods of image analysis could be usefully applied to the many thousands of HRTEM images that have been collected by other workers in the high temperature superconductor field. This work provides direct structural evidence for phase separation in high Tc cuprates, and gives support to recent stripes models that have been proposed to explain various angle resolved photoelectron spectroscopy and nuclear magnetic resonance data. We believe that the structural variation is a response to an opening of an electronic solubility gap where holes are not uniformly distributed in the material but are confined to metallic stripes. Optimum doping may occur as a consequence of the diffuse boundaries between stripes which arise from spinodal decomposition. Theoretical ideas about the high Tc cuprates which treat the cuprates as homogeneous may need to be modified in order to take account of this type of structural variation.

Rapid synthesis of porous honeycomb Cu/Pd through a hydrogen-bubble templating method

A rapid electrochemical method based on using a clean hydrogen-bubble template to form a bimetallic porous honeycomb Cu/Pd structure has been investigated. The addition of palladium salt to a copper-plating bath under conditions of vigorous hydrogen evolution was found to influence the pore size and bulk concentration of copper and palladium in the honeycomb bimetallic structure. The surface was characterised by X-ray photoelectron spectroscopy, which revealed that the surface of honeycomb Cu/Pd was found to be rich with a Cu/Pd alloy. The inclusion of palladium in the bimetallic structure not only influenced the pore size, but also modified the dendritic nature of the internal wall structure of the parent copper material into small nanometre-sized crystallites. The chemical composition of the bimetallic structure and substantial morphology changes were found to significantly influence the surface-enhanced Raman spectroscopic response for immobilised rhodamine B and the hydrogen-evolution reaction. The ability to create free-standing films of this honeycomb material may also have many advantages in the areas of gas- and liquid-phase heterogeneous catalysis.

Nanostructured copper oxide semiconductors : a perspective on materials, synthesis methods and applications

The oxides of copper (CuxO) are fascinating materials due to their remarkable optical, electrical, thermal and magnetic properties. Nanostructuring of CuxO can further enhance the performance of this important functional material and provide it with unique properties that do not exist in its bulk form. Three distinctly different phases of CuxO, mainly CuO, Cu2O and Cu4O3, can be prepared by numerous synthesis techniques including, vapour deposition and liquid phase chemical methods. In this article, we present a review of nanostructured CuxO focusing on their material properties, methods of synthesis and an overview of various applications that have been associated with nanostructured CuxO.

Preparation and in vitro release of zein microparticles loaded with prednisolone for oral delivery

Zein has been proposed as a polymer for targeted-drug delivery via the oral route. Zein microparticles were loaded with prednisolone and evaluated as an oral delivery system. Microparticles were formulated using phase separation. Starting quantities of zein and prednisolone, along with the agitation method and temperature were found to significantly impact drug loading and loading efficiency. Vortex mixing produced the highest drug loading and loading efficiency. Drug release was measured in simulated conditions of the stomach and small intestine using the microparticles made with the method that best improved drug loading. In simulated stomach and small intestine conditions, prednisolone release reached almost 70 over 3 and 4h, respectively. While a clinically relevant dose may be delivered using c. 100mg of zein microparticles, prednisolone release from the microparticles indicates that they may not be suited as a controlled-or targeted-delivery system.

Preparation of macroporous methacrylate monolithic material with convective flow properties for bioseparation : Investigating the kinetics of pore formation and hydrodynamic performance

The preparation of macroporous methacrylate monolithic material with controlled pore structures can be carried out in an unstirred mould through careful and precise control of the polymerisation kinetics and parameters. Contemporary synthesis conditions of methacrylate monolithic polymers are based on existing polymerisation schemes without an in-depth understanding of the dynamics of pore structure and formation. This leads to poor performance in polymer usage thereby affecting final product recovery and purity, retention time, productivity and process economics. The unique porosity of methacrylate monolithic polymer which propels its usage in many industrial applications can be controlled easily during its preparation. Control of the kinetics of the overall process through changes in reaction time, temperature and overall composition such as cross-linker and initiator contents allow the fine tuning of the macroporous structure and provide an understanding of the mechanism of pore formation within the unstirred mould. The significant effect of temperature of the reaction kinetics serves as an effectual means to control and optimise the pore structure and allows the preparation of polymers with different pore size distributions from the same composition of the polymerisation mixture. Increasing the concentration of the cross-linking monomer affects the composition of the final monoliths and also decreases the average pore size as a result of pre-mature formation of highly cross-linked globules with a reduced propensity to coalesce. The choice and concentration of porogen solvent is also imperative. Different porogens and porogen mixtures present different pore structure output. Example, larger pores are obtained in a poor solvent due to early phase separation.

Atomic resolution of dielectric with low temperature coefficient

Ba(6-3x)Nd(8+2x)Ti(18)O(54) (BNTl14) is a high permittivity dielectric with low temperature coefficient (Tcf). Low coefficient of change of dielectric permittivity with temperature (Tcf) is an unusual materials property. The research is aimed at discovering how atomic structure relates to temperature coefficient. Sub-Ångström scanning transmission electron microscopy (STEM) is used to measure mixed occupancy of Nd and Ba in atomic columns. It was expected that phase separation would occur to accommodate mixing of dissimilar ions. However no evidence of phase separation was found. There is a good image match between experiment and high angle annular dark field (HAADF) simulation. Vacancies and excess Ba ions appear to be randomly arranged on the available sites which would result in distortion of TiO6 octahedra. The low Tcf may arise from TiO6 distortion.

Supramolecular ordering in oligothiophene−fullerene monolayers

Scanning tunneling microscopy (STM) of monolayers comprising oligothiophene and fullerene molecular semiconductors reveals details of their molecular-scale phase separation and ordering with potential implications for the design of organic electronic devices, in particular future bulk heterojunction solar cells. Prochiral terthienobenzenetricarboxylic acid (TTBTA) self-assembles at the solution/graphite interface into either a porous chicken wire network linked by dimeric hydrogen bonding associations of COOH groups (R22(8)) or a close-packed network linked in a novel hexameric hydrogen bonding motif (R66(24)). Analysis of high-resolution STM images shows that the chicken wire phase is racemically mixed, whereas the close-packed phase is enantiomerically pure. The cavities of the chicken wire structure can efficiently host C60 molecules, which form ordered domains with either one, two, or three fullerenes per cavity. The observed monodisperse filling and long-range co-alignment of fullerenes is described in terms of a combination of an electrostatic effect and the commensurability between the graphite and molecular network, which leads to differentiation of otherwise identical adsorption sites in the pores.

In situ investigation of explosive crystallization in a-Ge: Experimental determination of the interface response function using dynamic transmission electron microscopy

The crystallization of amorphous semiconductors is a strongly exothermic process. Once initiated the release of latent heat can be sufficient to drive a self-sustaining crystallization front through the material in a manner that has been described as explosive. Here, we perform a quantitative in situ study of explosive crystallization in amorphous germanium using dynamic transmission electron microscopy. Direct observations of the speed of the explosive crystallization front as it evolves along a laser-imprinted temperature gradient are used to experimentally determine the complete interface response function (i.e., the temperature-dependent front propagation speed) for this process, which reaches a peak of 16 m/s. Fitting to the Frenkel-Wilson kinetic law demonstrates that the diffusivity of the material locally/immediately in advance of the explosive crystallization front is inconsistent with those of a liquid phase. This result suggests a modification to the liquid-mediated mechanism commonly used to describe this process that replaces the phase change at the leading amorphous-liquid interface with a change in bonding character (from covalent to metallic) occurring in the hot amorphous material.

Glass forming ability and stability: Ternary Cu bearing Ti, Zr, Hf alloys

Four Cu bearing alloys of nominal composition Zr25Ti25Cu50, Zr34Ti16Cu50, Zr25Hf25Cu50 and Ti25Hf25Cu50 have been rapidly solidified in order to produce ribbons. All the alloys become amorphous after meltspinning. In the Zr34Ti16Cu50 alloy localized precipitation of cF24 Cu5Zr phase can be observed in the amorphous matrix. The alloys show a tendency of phase separation at the initial stages of crystallization. The difference in crystallization behavior of these alloys with Ni bearing ternary alloys can be explained by atomic size, binary heat of mixing and Mendeleev number. It has been observed that both Laves and Anti-Laves phase forming compositions are suitable for glass formation. The structures of the phases, precipitated during rapid solidification and crystallization can be viewed in terms of Bernal deltahedra and Frank-Kasper polyhedra.

Novel Mass Spectrometric Analysis Methods for Anabolic Androgenic Steroids in Sports Drug Testing

The feasibility of different modern analytical techniques for the mass spectrometric detection of anabolic androgenic steroids (AAS) in human urine was examined in order to enhance the prevalent analytics and to find reasonable strategies for effective sports drug testing. A comparative study of the sensitivity and specificity between gas chromatography (GC) combined with low (LRMS) and high resolution mass spectrometry (HRMS) in screening of AAS was carried out with four metabolites of methandienone. Measurements were done in selected ion monitoring mode with HRMS using a mass resolution of 5000. With HRMS the detection limits were considerably lower than with LRMS, enabling detection of steroids at low 0.2-0.5 ng/ml levels. However, also with HRMS, the biological background hampered the detection of some steroids. The applicability of liquid-phase microextraction (LPME) was studied with metabolites of fluoxymesterone, 4-chlorodehydromethyltestosterone, stanozolol and danazol. Factors affecting the extraction process were studied and a novel LPME method with in-fiber silylation was developed and validated for GC/MS analysis of the danazol metabolite. The method allowed precise, selective and sensitive analysis of the metabolite and enabled simultaneous filtration, extraction, enrichment and derivatization of the analyte from urine without any other steps in sample preparation. Liquid chromatographic/tandem mass spectrometric (LC/MS/MS) methods utilizing electrospray ionization (ESI), atmospheric pressure chemical ionization (APCI) and atmospheric pressure photoionization (APPI) were developed and applied for detection of oxandrolone and metabolites of stanozolol and 4-chlorodehydromethyltestosterone in urine. All methods exhibited high sensitivity and specificity. ESI showed, however, the best applicability, and a LC/ESI-MS/MS method for routine screening of nine 17-alkyl-substituted AAS was thus developed enabling fast and precise measurement of all analytes with detection limits below 2 ng/ml. The potential of chemometrics to resolve complex GC/MS data was demonstrated with samples prepared for AAS screening. Acquired full scan spectral data (m/z 40-700) were processed by the OSCAR algorithm (Optimization by Stepwise Constraints of Alternating Regression). The deconvolution process was able to dig out from a GC/MS run more than the double number of components as compared with the number of visible chromatographic peaks. Severely overlapping components, as well as components hidden in the chromatographic background could be isolated successfully. All studied techniques proved to be useful analytical tools to improve detection of AAS in urine. Superiority of different procedures is, however, compound-dependent and different techniques complement each other.

Solubility of ibuprofen in aqueous ethanol at low ethanol contents

The solubility of ibuprofen was measured in water (W) - ethanol (E) mixtures from 0 to 50% w/w ethanol at 10, 25 and 40 °C by the dissolution method using UV spectrophotometry to determine the ibuprofen concentrations. The UV calibration for ibuprofen in different water - ethanol mixtures showed Beer - Lambert linearity, however the slopes differed, which indicated the structure of the drug is influenced by the solvent system i.e. the water - ethanol ratio. The ibuprofen solubility in water (zero ethanol) is low (~ 50 ppm) but increases near exponentially with increasing ethanol content. At 40 °C, there is phase separation between 34% and 63% w/w E/(E+W). The solubility data will be used to select precipitation crystallizer conditions to directly produce free flowing ibuprofen particles (<5 m) for developing a dry powder inhaler (DPI) formulation for lung delivery.