abstract: A hybrid molecular dynamics (MD) simulation framework is developed to emulate mechanochemical reaction of mechanophores in epoxy-based nanocomposites. Two different force fields, a classical force field and a bond order based force field are hybridized to mimic the experimental processes from specimen preparation to mechanical loading test. Ultra-violet photodimerization for mechanophore synthesis and epoxy curing for thermoset polymer generation are successfully simulated by developing a numerical covalent bond generation method using the classical force field within the framework. Mechanical loading tests to activate mechanophores are also virtually conducted by deforming the volume of a simulation unit cell. The unit cell deformation leads to covalent bond elongation and subsequent bond breakage, which is captured using the bond order based force field. The outcome of the virtual loading test is used for local work analysis, which enables a quantitative study of mechanophore activation. Through the local work analysis, the onset and evolution of mechanophore activation indicating damage initiation and propagation are estimated; ultimately, the mechanophore sensitivity to external stress is evaluated. The virtual loading tests also provide accurate estimations of mechanical properties such as elastic, shear, bulk modulus, yield strain/strength, and Poisson’s ratio of the system. Experimental studies are performed in conjunction with the simulation work to validate the hybrid MD simulation framework. Less than 2% error in estimations of glass transition temperature (Tg) is observed with experimentally measured Tgs by use of differential scanning calorimetry. Virtual loading tests successfully reproduce the stress-strain curve capturing the effect of mechanophore inclusion on mechanical properties of epoxy polymer; comparable changes in Young’s modulus and yield strength are observed in experiments and simulations. Early damage signal detection, which is identified in experiments by observing increased intensity before the yield strain, is captured in simulations by showing that the critical strain representing the onset of the mechanophore activation occurs before the estimated yield strain. It is anticipated that the experimentally validated hybrid MD framework presented in this dissertation will provide a low-cost alternative to additional experiments that are required for optimizing material design parameters to improve damage sensing capability and mechanical properties. In addition to the study of mechanochemical reaction analysis, an atomistic model of interphase in carbon fiber reinforced composites is developed. Physical entanglement between semi-crystalline carbon fiber surface and polymer matrix is captured by introducing voids in multiple graphene layers, which allow polymer matrix to intertwine with graphene layers. The hybrid MD framework is used with some modifications to estimate interphase properties that include the effect of the physical entanglement. The results are compared with existing carbon fiber surface models that assume that carbon fiber has a crystalline structure and hence are unable to capture the physical entanglement. Results indicate that the current model shows larger stress gradients across the material interphase. These large stress gradients increase the viscoplasticity and damage effects at the interphase. The results are important for improved prediction of the nonlinear response and damage evolution in composite materials. / Dissertation/Thesis / Doctoral Dissertation Mechanical Engineering 2017

An alternate family of “high” voltage (where the equilibrium voltage lies between 3.6 V and 4.2 V) polyanion cathode materials is reported in this thesis with the objective of improving specific energy density (Wh/kg) and developing a better understanding of polyanion electrochemistry. The electrochemical properties, synthesis and the structure of novel fluorosulfate materials crystallizing in the tavorite and the triplite type mineral structures are described. These materials display highest discharge voltages reported for any Fe2+/Fe3+ redox couple. LiFeSO4F was prepared in both the tavorite and the triplite polymorphs using inexpensive and scalable methods. Complete structural characterization was performed using X-ray and neutron based diffraction methods. A rapid synthesis of fluorosulfates can be achieved by using microwave heating. The local rapid heating created by the microwaves generates nanocrystalline LiFeSO4F tavorite with defects that induce significant microstrain. To date, this is unique to the microwave synthesis method. Phase transformation to the more stable triplite framework, facilitated by the lattice defects which include hydroxyl groups, is therefore easily triggered. The formation of nanocrystalline tavorite leads to nanocrystalline triplite, which greatly favors its electrochemical performance because of the inherently disordered nature of the triplite structure. Direct synthesis of the electrochemically active triplite type compound can be carried out either by extending the duration of the solvothermal reactions or by the partial substitution of Fe by Mn to produce LiFe1-xMnxSO4F. This study, overall, has led to a better understanding of the transformation of tavorite to the triplite phase. To examine Li and the Na ion conduction and their correlation with the electrochemical performance of 3-D, 2-D and 1-D ion conductors, atomistic scale simulations have been used to investigate tavorite type LiFeSO4F, NaFeSO4F, olivine type NaMPO4 (M= Fe, Mn, Fe0.5Mn0.5) and layered Na2FePO4F. These calculations predict high mobility of the Li-ion in the tavorite type LiFeSO4F but sluggish Na-ion transport in iso-structural NaFeSO4F. High mobility of the Na-ion is predicted for phosphate layered and olivine structures. Finally, the synthesis and structural details of NaMSO4F (M=Fe, Mn) and NH4MSO4F (M=Fe, Mn) are presented in the last chapter to show the structural diversity present in the fluorosulfate family.

An alternate family of “high” voltage (where the equilibrium voltage lies between 3.6 V and 4.2 V) polyanion cathode materials is reported in this thesis with the objective of improving specific energy density (Wh/kg) and developing a better understanding of polyanion electrochemistry. The electrochemical properties, synthesis and the structure of novel fluorosulfate materials crystallizing in the tavorite and the triplite type mineral structures are described. These materials display highest discharge voltages reported for any Fe2+/Fe3+ redox couple. LiFeSO4F was prepared in both the tavorite and the triplite polymorphs using inexpensive and scalable methods. Complete structural characterization was performed using X-ray and neutron based diffraction methods. A rapid synthesis of fluorosulfates can be achieved by using microwave heating. The local rapid heating created by the microwaves generates nanocrystalline LiFeSO4F tavorite with defects that induce significant microstrain. To date, this is unique to the microwave synthesis method. Phase transformation to the more stable triplite framework, facilitated by the lattice defects which include hydroxyl groups, is therefore easily triggered. The formation of nanocrystalline tavorite leads to nanocrystalline triplite, which greatly favors its electrochemical performance because of the inherently disordered nature of the triplite structure. Direct synthesis of the electrochemically active triplite type compound can be carried out either by extending the duration of the solvothermal reactions or by the partial substitution of Fe by Mn to produce LiFe1-xMnxSO4F. This study, overall, has led to a better understanding of the transformation of tavorite to the triplite phase. To examine Li and the Na ion conduction and their correlation with the electrochemical performance of 3-D, 2-D and 1-D ion conductors, atomistic scale simulations have been used to investigate tavorite type LiFeSO4F, NaFeSO4F, olivine type NaMPO4 (M= Fe, Mn, Fe0.5Mn0.5) and layered Na2FePO4F. These calculations predict high mobility of the Li-ion in the tavorite type LiFeSO4F but sluggish Na-ion transport in iso-structural NaFeSO4F. High mobility of the Na-ion is predicted for phosphate layered and olivine structures. Finally, the synthesis and structural details of NaMSO4F (M=Fe, Mn) and NH4MSO4F (M=Fe, Mn) are presented in the last chapter to show the structural diversity present in the fluorosulfate family.

Selon la théorie de Cottrell et Bilby, les dislocations à travers leur champ de contrainte interagissent avec les atomes de soluté qui s’agrègent au cœur et autour des dislocations (atmosphère de Cottrell). Ces atmosphères « bloquent » les dislocations et fragilisent le matériau. Dans cette thèse, les techniques de simulations à l’échelle atomique telles que la Dynamique Moléculaire, les simulations Monte Carlo Cinétique, Monte Carlo Métropolis ont été développées qui permettent de prendre en compte les interactions entre plusieurs centaines d’atomes de carbone et la dislocation, pour étudier la cinétique de formation ainsi que la structure d’une atmosphère de Cottrell. Par ailleurs, la technique de simulation est appliquée à deux autres problématiques: premièrement, il est connu que les atomes de C dans la ferrite se mettent en ordre (mise en ordre de Zener). La stabilité de cette phase est étudiée en fonction de la température et la concentration de C. Deuxièmement, la ségrégation des atomes de soluté dans les nano-cristaux de Ni ainsi que la stabilité des nano-cristaux avec les atomes de soluté dans les joints de grain à haute température est étudiée. / The objective of the thesis was to understand the microscopic features at the origin of ageing in metals. The originality of this contribution was the com- bination of three complementary computational techniques : (1) Metropolis Monte Carlo (MMC), (2) Atomic Kinetic Monte Carlo (AKMC), and, (3) Molecular Dynamics (MD). It consisted of four main sections : Firstly the ordering occurring in bulk alpha-iron via MMC and MD was studied. Various carbon contents and temperatures were investigated in order to obtain a “phase diagram”. Secondly, the generation of systems containing a dislocation interacting with many carbon atoms, namely a Cottrell Atmosphere, with MMC technique was described. The equilibrium structure of the atmosphere and the stress field around the atmospheres proves that the stress field around the dislocation was affected but not cancelled out by the atmosphere. Thirdly, the kinetics of the carbon migration and Cottrell atmosphere evolution were investigated via AKMC. The activation energies for carbon atom migration were calculated from the local stress field and the arrangement of the neigh- bouring carbon atoms. Lastly, an application of the combined use of MMC and MD to describe grain boundary segregation of solute atoms in fcc nickel was presented. The grain growth was inhibited due to the solute atoms in the grain boundary.