This dissertation is dedicated to develop models and methods to bridge atomistic and continuum scales of deformation processes in asphalt-aggregate interfacial composite materials systems. The deformation and failure behaviors, e.g. nanoscale strength, deformation, stiffness, and adhesion/cohesion at asphalt-aggregate interfaces are all evaluated by means of atomistic simulations. The atomistic modeling approach is employed to simulate mechanical properties, which is connected by their common dependence on the nanoscale bonding and their sensitive dependences on mechanics and moisture sensitivity. Specifically, CVFF-aug forcefield is employed in the atomistic calculations to study the fundamental failure processes that appear at the interface as a result of a mechanical deformation. There are five primary aspects to this dissertation. First, the multiscale features of asphalt concrete materials are characterized by using nanoscale characterization & fabrication devices, e.g. High Resolution Optical Microscope (HROM), Environmental Scanning Electron Microscope (ESEM), Transmission Electron Microscope (TEM), Focused Ion Beam (FIB), and Atomistic Force Microscope (AFM). Second, based on the multiscale devices characterization of the interfaces, a 2-layer atomistic bitumen-rock interface structure is constructed. Interface structure evolution under uniaxial tension is performed with various deformation rates. Comparison is made between both theoretical and experimental characterizations of interface configuration. Molecular dynamics (MD) simulations are used to investigate potential relationships between interface structure and morphology. Influences of deformation rate and temperature factors are discussed in terms of interface region stress-strain relation and loading time duration. Third, molecular dynamics simulations are also performed to provide a characterization of atomic scale mechanical behaviors for a 3-layer confined shear structure which leads to interfacial shear failure. In addition, atomistic static simulation approach is employed to calculate a couple of mineral crystals' elastic constants. Furthermore, molecular dynamics simulations are also used to predict the static, thermodynamic, and mechanical properties of three asphalt molecular models. Fourth, the high performance parallel computing technology is extensively employed throughout this dissertation. In addition to use the large-scale MD program, LAMMPS, the author developed a high performance parallel distributive computing program, MPI_multistress, to implement the multiscale understanding/predicting of materials mechanical behaviors. Finally, this research also focuses on the evaluation of the susceptibility of aggregates and asphalts to moisture damage through understanding the nano-mechanisms that influence adhesive bond between aggregates and asphalt, as well as the cohesive strength and moisture susceptibility of the specific asphalt-aggregate interfaces. Surface energy theory and pull-out simulation are used to compute the adhesive bond strength between the aggregates and asphalt, as well as the cohesive bond strength within the binder. In general, this dissertation has focused on the development of nanoscale modeling methods to assess asphalt-aggregate interfacial atomistic deformation and failure behaviors, as well as moisture effects on asphalt mixture strength. Simulation results provide valuable insights into mechanistic details of nanoscale interactions, particularly under conditions of various deformation rates and different temperatures. The results obtained show that a reasonable agreement between the theoretical and pavement industry observations is satisfactory. We conclude that the theoretical calculations presented here are useful in asphalt concrete industry for predicting the mechanical properties of asphalt-aggregate interfaces, which are difficult to obtain experimentally because of their small size. / Ph. D.

The metallic binary-alloy fuel Uranium-Zirconium is important for use in the new generation of advanced fast reactors. Uranium-Zirconium goes through a phase transition at higher temperatures to a (gamma) Body Centered Cubic (BCC) phase. The BCC high temperature phase is particularly important since it corresponds to the temperature range in which the fast reactors will operate. A semi-empirical Modified Embedded Atom Method (MEAM) potential is presented for Uranium-Zirconium. This is the first interatomic potential created for the U-Zr system. The bulk physical properties of the Uranium-Zirconium binary alloy were reproduced using Molecular Dynamics (MD) and Monte Carlo (MC) simulations with the MEAM potential. The simulation of bulk metallic alloy separation and ordering phenomena on the atomic scale using iterative MD and MC simulations with interatomic potentials has never been done before. These simulations will help the fundamental understanding of complex phenomena in the metallic fuels. This is a large step in making a computationally acceptable fuel performance code, able to replicate and predict fuel behavior.

No description available.

Irradiation induces a high concentration of defects in the structural materials of nuclear reactors, which are typically of body-centered cubic Iron (BCC Fe) and its alloys. The primary effect of irradiation is hardening which is caused by the blocking of dislocations with defects and defect clusters like point defects, self-interstitial loops, and voids. The dislocation-defect interactions are atomistic in nature due to the very small length and time scales involved, i.e., of the order of nanometers and picoseconds. To predict the effect of dislocation-defect interactions on the macroscopic mechanical and plastic behavior of the material, it is critically important to develop robust coupling schemes by which accurate atomic level physics of the rate-limiting kinetic processes can be informed into a coarse-grained model such as crystal plasticity. In this thesis we will develop an atomistically informed constitutive model. Relevant atomistic processes are identified from molecular dynamics simulations. The respective unit process studies are conducted using atomistic reaction pathway sampling methods like Nudged Elastic Band method. Stress-dependent activation energies and activation volumes are computed for various rate-liming unit processes like thermally activated dislocation motion via kinkpair nucleation, dislocation pinning due to self interstitial atom, etc. Constitutive laws are developed based on transition state theory, that informs the atomistically determined activation parameters into a coarse-grained crystal plasticity model. The macroscopic deformation behavior predicted by the crystal plasticity model is validated with experimental results and the characteristic features explained in the light of atomistic knowledge of the constituting kinetics. We also investigate on unique irradiation induced defects such as stacking fault tetrahedra, that are formed under non-irradiated condition. This thesis also includes our work on materials with internal interfaces that can resist irradiation induced damage. Overall, the research presented in this thesis involves the implementation and development of novel computational paradigm that encompasses computational approaches of various length and time scales towards robust predictions of the mechanical behavior of irradiated materials.

<p>Ni-base superalloys are among the highest temperature capable alloys and are used pervasively throughout the transportation, energy, and nuclear industries. However, their microstructures have been largely limited to containing the γ´ (cubic) and γ´´ (tetragonal) phases to enable high strength at elevated temperatures, and this fixation has restricted alloy development opportunities. In the past three decades, a new set of alloys, strengthened by the γ´´´ (orthorhombic) phase, was developed by Haynes International. The alloys exhibit comparable strength to existing Ni-based superalloys and show a 25% decrease in the thermal expansion coefficient, designed for tighter clearances (thus improving engine efficiency) and help to reduce thermally induced fatigue from engine cycling. </p> <p>The newest iteration of such alloys, HAYNES^® 244^®, has a nominal composition of Ni-22.5Mo-8Cr-6W (wt.%), and each alloying element is used to help precipitate the γ´´´-Ni₂(Cr, Mo, W) phase. The deformation mechanisms of this material are currently unknown. Previous studies investigating the predecessor alloy, HAYNES^® 242^® alloy, showed deformation twinning to be the dominant deformation mechanism during mechanical testing, but the physical phenomena responsible for this mode of deformation were not clearly elucidated. As a result, the primary motivation of this project is to understand the deformation behavior of the 244 alloy from the atomistic level and upwards. </p> <p>This work details efforts to elucidate these deformation mechanisms using an integrated computational and experimental approach. First-principles calculations were performed to determine the entire generalized stacking fault energy (GSFE) surface and slip pathways of the γ´´´ phase for dislocation slip. The various planar defects that could form from dislocation slip were predicted to provide significant barriers for dislocation motion due to their very high planar defect energies (~1000 mJ/m²), likely precluding shearing of the precipitates. We incorporated these results into phase field dislocation dynamics (PFDD) to simulate dislocation-precipitate interactions of finite size. These results showed that the planar defect energies of the γ´´´ phase largely govern the deformation behavior and critical resolved shear stress for precipitate shearing, regardless of precipitate shape, size, or orientation. Extensive mechanical testing conducted from room temperature up to 760 ºC over strain rates ranging from 10^-9 s^-1 to 10^-4 s^-1 combined with transmission electron microscopy validated the predicted deformation structures of creep and tensile samples. Shearing of individual precipitates by intrinsic and extrinsic stacking faults, as well as extensive deforming twinning, was observed. The integrated GSFE and PFDD simulations showed that the precipitates would resist dislocation shearing and favor twinning as the preferred deformation mechanism at all temperatures and strain rates investigated. These results provide pathways for microstructural and composition modification to further increase the strength of γ´´´ strengthened alloys in the future.</p> <p><br></p>

Amorphous and bulk amorphous metallic alloys are an intriguing class of structural materials and possess a range of interesting properties, including near theoretical strength, high hardness, extremely low damping characteristics, excellent wear properties, high corrosion resistance, low shrinkage during cooling and almost perfect as-cast surfaces with good potential for forming and shaping. In this study, new Ti-based bulk amorphous alloys are tried to be modeled and synthesized. For that purpose, electronic theory of alloys in the pseudo potential approximation was used as a tool for understanding the theory lying beneath the bulk glass forming ability (BGFA). The results from this approach were evaluated both separately and together with the other theories supposed by our colleagues. Glass forming parameters of ordering energy, &amp / #916 / HM, viscosity, mismatch entropy, Rc was calculated for various Ti-based binary and ternary and the change in these parameters in both cases was evaluated. The results of the theoretical calculations of glass forming parameters has shown good relation with the literature data that the predicted alloying elements, i.e. Mo, Hf, Zr, B, Fe, Sn, and Be, to increase GFA for Ti2Ni binary system were generally used in the production of Ti-based bulk amorphous alloys. In the second part of this thesis, new Ti-based compositions with high GFA were tried to be synthesized with light of these results and encouraging conclusions were drawn. The production of these alloys were made with centrifugal casting method which is relatively a new technique for producing such alloys and the characterization of these alloys were made with metallographic, X-ray and thermal means.

In recent years, the use of first-principles based atomistic modeling technique has become extremely popular to gain better insights on the various locally modulated electronic properties of nano materials and structures. Atomistic modeling offers the benefit of predicting crystal structures, visualizing orbital distribution and electron density, as well as understanding material properties which are hard to access experimentally. The single layer MoS2 has emerged as a suitable choice for the next generation nano devices, owing to its distinctive electrical, optical and mechanical properties like, better electrostatics, increased photo luminescence, higher mechanical flexibility, etc. The realization of decananometer scale digital switches with the single layer MoS2 as the channel may provide many significant advantages such as, high On/Off current ratio, excellent electrostatic control of the gate, low leakage, etc. However, there are quite a few critical issues such as, forming low resistance source/drain contacts, achieving higher effective mobility, ensuring large scale controlled growth, etc. which need to be addressed for successful implementation of the atomically thin transistors in integrated circuits. Recent experimental demonstration showing the coexistence of metallic and semiconducting phases in the same monolayer MoS2, has attracted much attention for its use in ultra-low contact resistance-MoS2 transistors. Howbeit, the electronic structures of the metallic-to-semiconducting phase boundaries, which appear to dictate the carrier injection in such transistors, are not yet well understood. In this work, we first develop the geometrically optimized atomistic models of the 2H-1T′ hetero-phase structures with two distinct phase boundaries, β and γ. We then apply density functional theory to calculate the electronic structures for those optimized geometries. Furthermore, we employ non equilibrium Green’s function formalism to evaluate the transmission spectra and the local density of states in order to assess the Schottky barrier nature of the phase boundaries. Nonetheless, the symmetry of the source-channel and drain-channel junction, is a unique property of a metal-oxide semiconductor field effect transistor (MOSFET), which needs to be preserved while realizing sub-10 nm channel length devices using advanced technology. Employing experimental-findings-driven atomistic modeling technique, we demonstrate that such symmetry might not be preserved in an atomically thin phase-engineered MoS2- based MOSFET. It originates from the two distinct atomic patterns at phase boundaries (β and β*) when the semiconducting phase (channel) is sandwiched between the two metallic phases (source and drain). Next, using first principles based quantum transport calculations we demonstrate that due to the clusterization of “Mo” atoms in 1T′ MoS2, the transmission along the zigzag direction is significantly higher than that in the armchair direction. Moreover, to achieve excellent impedance matching with various metal contacts (such as, “Au”, “Pd”, etc.), we further develop the atomistic models of metal-1T′ MoS2 edge contact geometries and compute their resistance values. Other than the charge carrier transport, analysing the heat transport across the channel is also crucial in designing the ultra-thin next generation transistors. Hence, in this thesis work, we have investigated the electro-thermal transport properties of single layer MoS2, in quasi ballistic regime. Besides the perfect monolayer in its pristine form, we have also considered various line defects which have been experimentally observed in mechanically exfoliated MoS2 samples. Furthermore, a comprehensive study on the phonon thermal conductivity of a suspended monolayer MoS2, has been incorporated in this thesis. The studies presented in this thesis could be useful for understanding the carrier transport in atomically thin devices and designing the ultra-thin next generation transistors.

No description available.

Le dioxyde d’uranium UO2 est le combustible standard dans les réacteurs nucléaires à eau pressurisée (REP). Durant le fonctionnement du réacteur les pastilles combustibles subissent des contraintes thermiques et mécaniques. Pour cette raison il est très important de bien connaître les propriétés de ce système à la fois dans les conditions de fonctionnement normales et accidentelles (300 à 2000K). Lors des réactions de fission de l’uranium, des gaz rares comme le xénon sont produits à l’intérieur du combustible. En raison de leur faible solubilité, ces gaz vont former des bulles intra- et inter- granulaires dans l’UO2. La présence de ces bulles dans le combustible a un impact sur les propriétés macroscopiques de ce dernier. A l'échelle nanométrique, les bulles intragranulaires prennent la forme d’un octaèdre facetté, essentiellement suivant les directions (111) et (100). Devant la complexité de l’étude de la stabilité de cet octaèdre, nous avons décomposé le problème afin de pouvoir l’étudier de façon plus systématique et de découpler les différents effets. Dans un premier temps, nous avons déterminé la stabilité des surfaces planes (111) et (100) de l’UO2 et les modifications de microstructure engendrées par leur relaxation. Dans un deuxième temps, nous avons caractérisé les isothermes d’adsorption du xénon sur ces surfaces relaxées, en les comparant à ceux de l’incorporation dans une boîte vide pour identifier les effets de surface. Une attention particulière a été portée sur la microstructure du xénon dans ces systèmes. Finalement, nous avons effectué une analyse des propriétés mécaniques (profils de pression et de contrainte au voisinage des surfaces). / Uranium dioxide UO2 is the standard fuel in nuclear pressurized water reactors (PWR). During the operation of the reactor the fuel pellets undergo thermal and mechanical stresses. For this reason it is very important to understand these thermomechanical properties of this system both in normal operation conditions and accidental situations (300 to 2000K). During fission reactions of uranium, rare gases such as xenon are produced within the fuel. Due to their low solubility, these gases will either be released or form intra- and inter-granular bubbles inside the UO2. The presence of these bubbles in the fuel has an impact on the thermomechanical properties of the latter. We focus in this thesis on the study of intragranularbubbles and their impact on the thermomechanical properties of UO2 , through modeling at the atomic scale. At this scale, intragranular bubbles take the shape of an octahedron, presenting mainly (111) and (100) facets. Given the complexity of the study of the stability of this octahedron, we have simplified the problem in order to study it in a more systematic way and to decouple the various effects. First, the stability of (100) and (111) extended surfaces of UO2 and microscructural modifications generated by their relaxation were studied. In a second step, we dermined adsorption isotherms of xenon on these relaxed surfaces, and compared them to the incorporation ones inside an empty box in order to isolate surface effects. A specific attention has been given to the microstructure of xenon in these systems. Finally, an analysis of the mechanical properties (pressure and stress profiles near by the surface).

Ce travail théorique vise à étudier, via les méthodes Premiers Principes, les propriétés des complexes de métaux de transitions, left[Mleft(dmbpyright)_{3}right]^{n+}nCi^{-} pour un usage en batterie. Pour cette étude ab-initio, les composés mono et bi-nucléaires ont été retenus. La pertinance de notre modélisation a été validée sur les composés mononucléaires. Nous nous sommes interessé au complexes de Fe, Ru et Cu pour lesquels une validation expérimentale était possible. Notre étude a principalement consisté à faire varier les degrés de liberté que nous possédons pour optimiser le voltage et la cinétique de chargement des batteries. Pour cela, nous avons fait varier le TM = Fe, Ru, et Cu, la nature des contre-ions Ci^{-}=PF_{6}^{-}, TFSI^{-} et ClO_{4}^{-} en interaction avec le polymère lors du processus de charge, ainsi que la longeur de la chaîne alkyl qui sépare les deux monomers dans le cas des composés binucléaires. Le composé à base de Fe avec une chaîne -left(CH_{2}right)_{n=6}- a été retenu comme le meilleur candidat pour une application batterie. Le composé à base Ru montre un comportement proche de celui du Fe, quant-au complexe de Cu, il présente des changements de géométrie locale sous chargement trop importants, le rendant peu apte à conduire à une cinétique efficace. Cette étude nous a permis de déterminer que l'approximation PBE était le meilleur choix possible pour modéliser nos complexes dans les conditions de fonctionnement en batterie (dans le champ créé par les contre-ions) et que l'approximation PBE0, généralement utilisée dans la littérature, ne pouvait rendre compte de la physico-chimie de nos composés dans de telles conditions.De surcroît, nous avons dévelopé pour le complexe de Fe, un potentiel atomistique de type “Champ de forces” de manière à pouvoir aborder les aspects dynamiques impliquant de plus grandes tailles de boîte de simulation. Ici, nous modélisons une structure 3D, totalement réticulée à partir de nos monomères à base de Fe. Nous nous sommes servi de la base de donnés DFT que nous avions généré (énergies, géométries, état de spin et fréquences vibrationnelles calculées) pour ajuster les paramètres entrant dans l'écriture du modèle. La construction de la géométrie initiale du polymère 3D a nécessité l'écriture d'un code de calcul visant à produire un arrangement complétement réticulé et à assigner les charges effectives issues des calculs DFT. Ce modèle nous a permis de déterminer les coefficients de diffusion des contre-ions pour les états totalement chargé et non-chargé. Un calcul plus ambitieux vise à déterminer les chemins de diffusion des contre-ions lors d'un processus de chargement en considérant un seul centre de degré d'oxydation 3+ au centre du polymère 3D, pour lequel les centres actifs possèdent un degré d'oxidation 2+. Les contre-ions assurent la neutralité globale.Keyword: Polymer, Electrochemistry, Li-ion Battery, DFT, Force Field development, 3D structure, Atomistic modeling / Abstract Standard redox potentials for mono and bi-nuclear transition metal (TM) complexes left[Mleft(dmbpyright)_{3}right]^{n+}nCi^{-}, have been investigated using First Principles Calculation. Three metal centers are investigated: Fe, Ru, and Cu. Our modeling is validated on mono-nuclear compounds. This approach consists in determining the best small polymer (bi-nuclear) made out of these monomers for a battery application. For that, we varied the three available degrees of freedom i.e., the nature of the central TM atom (Fe, Ru, and Cu), counter-ions Ci=PF_{6}^{-}, TFSI^{-} and ClO_{4}^{-} in interaction with the polymer, and the alkyl chain -left(CH_{2}right)_{n}- of length n that connects both mono-nuclear in the bi-nuclear compound. The Iron compound with -left(CH_{2}right)_{n=6}- is found to be the best candidate. The left[Culeft(dmbpyright)_{2}right]^{n+}nCi^{-} complex shows too much structure deformation upon loading, making it less reliable for cathode material. Moreover, we studied two XC functional, PBE and PBE0 and found, for three complexes PBE approximation retains the ligand field picture whereas PBE0 functional induces an exaggerated and unexpected band dispersion by dissolving the ligand field picture expected for the octahedral environment of the TM in the studied complexes. These findings validate that hybrid functional for which it was designed to localize and cancel self-interaction error does not work for all system. More particularly, the PBE0 approximation fails to model the three complexes (Fe, Ru, and Cu) in functional conditions (in the field made by the counter-ions).Abstract Further, we have developed an atomistic potential relying on the Force Field scheme for the Iron complex in order to study the dynamical properties of this compound at larger simulation scale (3D reticulated polymerization made of our Fe complex monomers). We made an intensive use of our DFT data (energies, geometries, spin-state configurations and calculated vibrational properties) to develop the required parameters entering the model. Moreover, computational techniques (written python language) were developed specifically to create a 3D structure of transition metal complexes satisfying the condition to be fully reticulated. Bounding conditions had to be designed and a procedure aiming at fixing reliable and physical effective charges on each atom of the simulation cell (compatible with DFT results) were developed. Our first simulations have been attached to calculate the diffusion coefficients of the counter-ions in both the fully loaded and unloaded states. A more ambitious and realistic calculation aims at investigating the paths of the counter-ions when one single center starts to be loaded in an unloaded environment.Abstract Keyword: Polymer, Electrochemistry, Li-ion Battery, DFT, Force Field development, 3D structure, Atomistic modeling