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.

Role of charge in destabilizing AlH4 and BH4 complex anions for hydrogen storage applications : Ab initio density functional calculations

NaAlH4 and LiBH4 are potential candidate materials for mobile hydrogen storage applications, yet they have the drawback of being highly stable and desorbing hydrogen only at elevated temperatures. In this letter, ab initio density functional theory calculations reveal how the stabilities of the AlH4 and BH4 complex anions will be affected by reducing net anionic charge. Tetrahedral AlH4 and BH4 complexes are found to be distorted with the decrease of negative charge. One H-H distance becomes smaller and the charge density will overlap between them at a small anion charge. The activation energies to release of H2 from AlH4 and BH4 complexes are thus greatly decreased. We demonstrate that point defects such as neutral Na vacancies or substitution of a Na atom with Ti on the NaAlH4(001) surface can potentially cause strong distortion of neighboring AlH4 complexes and even induce spontaneous dehydrogenation. Our results help to rationalize the conjecture that the suppression of charge transfer to AlH4 and BH4 anion as a consequence of surface defects should be very effective for improving the recycling performance of H2 in NaAlH4 and LiBH4. The understanding gained here will aid in the rational design and development of hydrogen storage materials based on these two systems.

Nuclear instruments and methods in physics research section B : beam interactions with materials and atoms

The impact-induced deposition of Al13 clusters with icosahedral structure on Ni(0 0 1) surface was studied by molecular dynamics (MD) simulation using Finnis–Sinclair potentials. The incident kinetic energy (Ein) ranged from 0.01 to 30 eV per atom. The structural and dynamical properties of Al clusters on Ni surfaces were found to be strongly dependent on the impact energy. At much lower energy, the Al cluster deposited on the surface as a bulk molecule. However, the original icosahedral structure was transformed to the fcc-like one due to the interaction and the structure mismatch between the Al cluster and Ni surface. With increasing the impinging energy, the cluster was deformed severely when it contacted the substrate, and then broken up due to dense collision cascade. The cluster atoms spread on the surface at last. When the impact energy was higher than 11 eV, the defects, such as Al substitutions and Ni ejections, were observed. The simulation indicated that there exists an optimum energy range, which is suitable for Al epitaxial growth in layer by layer. In addition, at higher impinging energy, the atomic exchange between Al and Ni atoms will be favourable to surface alloying.

Energy window effect in chemisorption of C36 on diamond (001)-(/2×1) surface

In this paper, the collision of a C36, with D6h symmetry, on diamond (001)-(/2×1) surface was investigated using molecular dynamics (MD) simulation based on the semi-empirical Brenner potential. The incident kinetic energy of the C36 ranges from 20 to 150 eV per cluster. The collision dynamics was investigated as a function of impact energy Ein. The C36 cluster was first impacted towards the center of two dimers with a fixed orientation. It was found that when Ein was lower than 30 eV, C36 bounces off the surface without breaking up. Increasing Ein to 30-45 eV, bonds were formed between C36 and surface dimer atoms, and the adsorbed C36 retained its original free-cluster structure. Around 50-60 eV, the C36 rebounded from the surface with cage defects. Above 70 eV, fragmentation both in the cluster and on the surface was observed. Our simulation supported the experimental findings that during low-energy cluster beam deposition small fullerenes could keep their original structure after adsorption (i.e. the memory effect), if Ein is within a certain range. Furthermore, we found that the energy threshold for chemisorption is sensitive to the orientation of the incident C36 and its impact position on the asymmetric surface.

Deformation and failure of graphene sheet and graphene-polymer interface

With a monolayer honeycomb-lattice of sp2-hybridized carbon atoms, graphene has demonstrated exceptional electrical, mechanical and thermal properties. One of its promising applications is to create graphene-polymer nanocomposites with tailored mechanical and physical properties. In general, the mechanical properties of graphene nanofiller as well as graphene-polymer interface govern the overall mechanical performance of graphene-polymer nanocomposites. However, the strengthening and toughening mechanisms in these novel nanocomposites have not been well understood. In this work, the deformation and failure of graphene sheet and graphene-polymer interface were investigated using molecular dynamics (MD) simulations. The effect of structural defects on the mechanical properties of graphene and graphene-polymer interface was investigated as well. The results showed that structural defects in graphene (e.g. Stone-Wales defect and multi-vacancy defect) can significantly deteriorate the fracture strength of graphene but may still make full utilization of corresponding strength of graphene and keep the interfacial strength and the overall mechanical performance of graphene-polymer nanocomposites.

Molecular dynamics simulation of fracture strength and morphology of defective graphene

Different types of defects can be introduced into graphene during material synthesis, and significantly influence the properties of graphene. In this work, we investigated the effects of structural defects, edge functionalisation and reconstruction on the fracture strength and morphology of graphene by molecular dynamics simulations. The minimum energy path analysis was conducted to investigate the formation of Stone-Wales defects. We also employed out-of-plane perturbation and energy minimization principle to study the possible morphology of graphene nanoribbons with edge-termination. Our numerical results show that the fracture strength of graphene is dependent on defects and environmental temperature. However, pre-existing defects may be healed, resulting in strength recovery. Edge functionalization can induce compressive stress and ripples in the edge areas of graphene nanoribbons. On the other hand, edge reconstruction contributed to the tensile stress and curved shape in the graphene nanoribbons.

Serine protease inhibition and mitochondrial dysfunction associated with cisplatin resistance in human tumor cell lines: targets for therapy

Indicators of mitochondrial function were studied in two different cell culture models of cis-diamminedichloroplatinum-II (CDDP) resistance: the intrinsically resistant human ovarian cancer cell line CI-80-13S, and resistant clones (HeLa-S1a and HeLa-S1b) generated by stable expression of the serine protease inhibitor—plasminogen activator inhibitor type-2 (PAI-2), in the human cervical cancer cell line HeLa. In both models, CDDP resistance was associated with sensitivity to killing by adriamycin, etoposide, auranofin, bis[1,2-bis(diphenylphosphino)ethane]gold(I) chloride {[Au(DPPE)2]Cl}, CdCl2 and the mitochondrial inhibitors rhodamine-123 (Rhl23), dequalinium chloride (DeCH), tetraphenylphosphonium (TPP), and ethidium bromide (EtBr) and with lower constitutive levels of ATP. Unlike the HeLa clones, CI-80-13S cells were additionally sensitive to chloramphenicol, 1-methyl-4-phenylpyridinium ion (MPP+), rotenone, thenoyltrifluoroacetone (TTFA), and antimycin A, and showed poor reduction of 1-[4,5-dimethylthiazol-2-yl]-2,5-diphenyltetrazolium bromide (MTT), suggesting a deficiency in NADH dehydrogenase and/or succinate dehydrogenase activities. Total platinum uptake and DNA-bound platinum were slightly lower in CI-80-13S than in sensitive cells. The HeLa-S1a and HeLa-S1b clones, on the other hand, showed poor reduction of triphenyltetrazolium chloride (TTC), indicative of low cytochrome c oxidase activity. Total platinum uptake by HeLa-S1a was similar to HeLa, but DNA-bound platinum was much lower than for the parent cell line. The mitochondria of CI-80-13S and HeLa-S1a showed altered morphology and were fewer in number than those of JAM and HeLa. In both models, CDDP resistance was associated with less platinum accumulation and with mitochondrial and membrane defects, brought about one case with expression of a protease inhibitor which is implicated in tumor progression. Such markers may identify tumors suitable for treatment with gold phosphine complexes or other mitochondrial inhibitors.

In situ observation of structural changes in polycrystalline silver catalysts by environmental scanning electron microscopy

Morphology changes induced in polycrystalline silver catalysts as a result of heating in either oxygen, water or oxygen-methanol atmospheres have been investigated by environmental scanning electron microscopy (ESEM), FT-Raman spectroscopy and temperature programmed desorption (TPD). The silver catalyst of interest consisted of two distinct particle types, one of which contained a significant concentration of sub-surface hydroxy species (in addition to surface adsorbed atomic oxygen). Heating the sample to 663 K resulted in the production of 'pin-holes' in the silver structure as a consequence of near-surface explosions caused by sub-surface hydroxy recombination. Furthermore, 'pin-holes' were predominantly found in the vicinity of surface defects, such as platelets and edge structures. Reaction between methanol and oxygen also resulted in the formation of 'pin-holes' in the silver surface, which were inherently associated with the catalytic process. A reaction mechanism is suggested that involves the interaction of methanol with sub-surface oxygen species to form sub-surface hydroxy groups. The sub-surface hydroxy species subsequently erupt through the silver surface to again produce 'pin-holes'.

A combined environmental scanning electron microscopy and Raman microscopy study of methanol oxidation on silver(I) oxide

The techniques of environmental scanning electron microscopy (ESEM) and Raman microscopy have been used to respectively elucidate the morphological changes and nature of the adsorbed species on silver(I) oxide powder, during methanol oxidation conditions. Heating Ag2O in either water vapour or oxygen resulted firstly in the decomposition of silver(I) oxide to polycrystalline silver at 578 K followed by sintering of the particles at higher temperature. Raman spectroscopy revealed the presence of subsurface oxygen and hydroxyl species in addition to surface hydroxyl groups after interaction with water vapour. Similar species were identified following exposure to oxygen in an ambient atmosphere. This behaviour indicated that the polycrystalline silver formed from Ag2O decomposition was substantially more reactive than silver produced by electrochemical methods. The interaction of water at elevated temperatures subsequent to heating silver(I) oxide in oxygen resulted in a significantly enhanced concentration of subsurface hydroxyl species. The reaction of methanol with Ag2O at high temperatures was interesting in that an inhibition in silver grain growth was noted. Substantial structural modification of the silver(I) oxide material was induced by catalytic etching in a methanol/air mixture. In particular, "pin-hole" formation was observed to occur at temperatures in excess of 773 K, and it was also recorded that these "pin- holes" coalesced to form large-scale defects under typical industrial reaction conditions. Raman spectroscopy revealed that the working surface consisted mainly of subsurface oxygen and surface Ag=O species. The relative lack of sub-surface hydroxyl species suggested that it was the desorption of such moieties which was the cause of the "pin-hole" formation.

In situ imaging of catalytic etching on silver during methanol oxidation conditions by environmental scanning electron microscopy

Polycrystalline silver is used to catalytically oxidise methanol to formaldehyde. This paper reports the results of extensive investigations involving the use of environmental scanning electron microscopy (ESEM) to monitor structural changes in silver during simulated industrial reaction conditions. The interaction of oxygen, nitrogen, and water, either singly or in combination, with a silver catalyst at temperatures up to 973 K resulted in the appearance of a reconstructed silver surface. More spectacular was the effect an oxygen/methanol mixture had on the silver morphology. At a temperature of ca. 713 K pinholes were created in the vicinity of defects as a consequence of subsurface explosions. These holes gradually increased in size and large platelet features were created. Elevation of the catalyst temperature to 843 K facilitated the wholescale oxygen induced restructuring of the entire silver surface. Methanol reacted with subsurface oxygen to produce subsurface hydroxyl species which ultimately formed water in the subsurface layers of silver. The resultant hydrostatic pressure forced the silver surface to adopt a "hill and valley" conformation in order to minimise the surface free energy. Upon approaching typical industrial operating conditions widespread explosions occurred on the catalyst and it was also apparent that the silver surface was extremely mobile under the applied conditions. The interaction of methanol alone with silver resulted in the initial formation of pinholes primarily in the vicinity of defects, due to reaction with oxygen species incorporated in the catalyst during electrochemical synthesis. However, dramatic reduction in the hole concentration with time occurred as all the available oxygen became consumed. A remarkable correlation between formaldehyde production and hole concentration was found.

In situ Raman studies of the selective oxidation of methanol to formaldehyde and ethene to ethylene oxide on a polycrystalline silver catalyst

The combined techniques of in situ Raman microscopy and scanning electron microscopy (SEM) have been used to study the selective oxidation of methanol to formaldehyde and the ethene epoxidation reaction over polycrystalline silver catalysts. The nature of the oxygen species formed on silver was found to depend critically upon the exact morphology of the catalyst studied. Bands at 640, 780 and 960 cm-1 were identified only on silver catalysts containing a significant proportion of defects. These peaks were assigned to subsurface oxygen species situated in the vicinity of surface dislocations, AgIII=O sites formed on silver atoms modified by the presence of subsurface oxygen and O2 - species stabilized on subsurface oxygen-modified silver sites, respectively. The selective oxidation of methanol to formaldehyde was determined to occur at defect sites, where reaction of methanol with subsurface oxygen initially produced subsurface OH species (451 cm-1) and adsorbed methoxy species. Two distinct forms of adsorbed ethene were identified on oxidised silver sites. One of these was created on silver sites modified by the interaction of subsurface oxygen species, and the other on silver crystal planes containing a surface coverage of atomic oxygen species. The selective oxidation of ethene to ethylene oxide was achieved by the reaction between ethene adsorbed on modified silver sites and electrophilic AgIII=O species, whereas the combustion reaction was perceived to take place by the reaction of adsorbed ethene with nucleophilic surface atomic oxygen species. Defects were determined to play a critical role in the epoxidation reaction, as these sites allowed the rapid diffusion of oxygen into subsurface positions, and consequently facilitated the formation of the catalytically active AgIII=O sites.

Monitoring of complex systems of interacting dynamic systems

Increases in functionality, power and intelligence of modern engineered systems led to complex systems with a large number of interconnected dynamic subsystems. In such machines, faults in one subsystem can cascade and affect the behavior of numerous other subsystems. This complicates the traditional fault monitoring procedures because of the need to train models of the faults that the monitoring system needs to detect and recognize. Unavoidable design defects, quality variations and different usage patterns make it infeasible to foresee all possible faults, resulting in limited diagnostic coverage that can only deal with previously anticipated and modeled failures. This leads to missed detections and costly blind swapping of acceptable components because of one’s inability to accurately isolate the source of previously unseen anomalies. To circumvent these difficulties, a new paradigm for diagnostic systems is proposed and discussed in this paper. Its feasibility is demonstrated through application examples in automotive engine diagnostics.

Investigation between the S377G3 GATA-4 polymorphism and migraine

Migraine is a common and painful neurological disorder, with genetic and environmental components. Several conditions have been shown to be comorbid with migraine, notably a cardiac malformation affecting the interatrial septum and leading to patent foramen ovale (PFO). Mutations in the development regulatory gene GATA-4, located on human chromosome 8p23.1-p22, have been found to be responsible for some cases of congenital heart defects including PFO. To determine whether the GATA-4 gene is involved in migraine, the present study performed an association analysis of a common GATA-4 variant that results in a change of amino acid (S377G), in a large case/control population (275 unrelated Caucasian migraineurs versus 275 control individuals). The results showed that there was no significant association for this polymorphism between migraine and controls (χ² = 0.84, P = 0.66). Thus it appears that the GATA-4 (S377G) mutation does not play a significant role in common migraine susceptibility.

The search for migraine genes : an overview of current knowledge

Migraine is a complex familial condition that imparts a significant burden on society. There is evidence for a role of genetic factors in migraine, and elucidating the genetic basis of this disabling condition remains the focus of much research. In this review we discuss results of genetic studies to date, from the discovery of the role of neural ion channel gene mutations in familial hemiplegic migraine (FHM) to linkage analyses and candidate gene studies in the more common forms of migraine. The success of FHM regarding discovery of genetic defects associated with the disorder remains elusive in common migraine, and causative genes have not yet been identified. Thus we suggest additional approaches for analysing the genetic basis of this disorder. The continuing search for migraine genes may aid in a greater understanding of the mechanisms that underlie the disorder and potentially lead to significant diagnostic and therapeutic applications.

Interactions between undifferentiated and osteogenically differentiated mesenchymal stromal cells during osteogenesis

The body of work presented in this dissertation has demonstrated that the interactions between donor cells and host cells are critical for bone repair and regeneration. The donor cells secrete VEGF which activates the downstream PI3K/Akt signaling pathway, ultimately leading to host cell recruitment and robust bone regeneration. The findings from this dissertation may provide a scientific rationale for the development of novel therapeutic strategies in the treatment and management of bone defects.