E1: Diffusion I
Monday PM, November 26, 2007
Constitution B (Sheraton)
9:00 AM - **E1.1
Diffusion in Intermetallics Using On-the-fly Kinetic Monte Carlo.
Murray Daw 1 , Erdi Bleda 1 , Xing Gao 1 Show Abstract
1 Dept of Physics, Clemson, Clemson, South Carolina, United States
We have calculated the vacancy-assisted tracer diffusivities in these binary intermetallic alloys using on-the-fly kinetic monte carlo: Cu-Au, Au-Ag, and Cu-Ni. This set encompases behavior from weak segregation to moderately strong ordering. The energetics are based on the Embedded Atom Method, and the saddlepoints are found on-the-fly (that is, beginning from scratch with each valley.) Using this technique we have calculated the tracer diffusivities over a wide range of temperature and composition. We evaluate the strengths and weaknesses of the approach. Where possible we compare to experimental data.
9:30 AM - E1.2
Kinetic Monte Carlo Simulations of Defect Diffusion in hcp alpha-Zirconium.
Cristina Arevalo 1 3 , Jose Perlado 1 , Maria Caturla 2 Show Abstract
1 , Fusion Nuclear Institute, madrid Spain, 3 , Escuela Técnica Superior de Ingenieros, Sevilla Spain, 2 , Universidad de Alicante, Alicante Spain
The study of point defect clustering in hexagonal-close-packed (hcp) metals is dominated by a consideration of the geometry of the hcp lattice and lattice parameters ratio (c/a). Because of this crystallographic anisotropy, defect anisotropic diffusion is expected (jump distances and jump rates depend on jump directions). This study has focused on hcp alpha-Zirconium (c/a 1.594, lower than ideal 1.633 and similar to Titanium, c/a 1.587). We have created a new and original model for the understanding of the microscopic evolution (defect diffusion) in hcp metals, using a kinetic Monte Carlo (kMC) simulation technique. This technique allows us to understand the evolution of damage accumulation, due to either neutron or electron irradiation, for long times (hours-months). Several multiscale modelling simulations steps have been used in order to understand the microscopical fission reactor cladding behaviour. We have focused on zirconium alloys claddings (Zircaloy-4 and Zr-2.5%Nb). The first step we have simulated has been the neutron spectra. We have obtained spectra for current pressure water reactors (PWR) and high burn-up advanced reactors. We have also obtained neutron spectra in several burnt steps and we have represented the isotopic variation in the cladding (looking for Helium and Hydrogen apparition in the system). Taking these results as input data we have reproduced the Primary Knock-on Atom (PKA) spectra. From those data a systematic analysis of primary damage has been obtained using binary collisions code SRIM for high energy recoils, in order to get distribution of cascades and subcascades for these recoil energies. With these data we have studied the evolution of the microstructure during irradiation under environment conditions of 600K, dose rate 1e-6 dpa/s and final dose of 0.5 dpa. We have selected the molecular dynamics (MD) simulations cascades that fit best to PKA spectra and we have used them as input data on defect energetics and cascade damage. We have considered isotropic motion for vacancies and we have studied how the accumulation of damage is affected considering from one dimension to three dimension movement for interstitials.
9:45 AM - E1.3
Kinetic Monte-Carlo Simulation of the Effect of Strain Field on Hydrogen Diffusion in hcp Zr.
Kenichi Nakashima 1 , Naoki Soneda 1 , Akiyoshi Nomoto 1 , Misako Iwasawa 1 Show Abstract
1 Materials Science Research Laboratory, Central Research Institute of Electric Power Industry, Komae, Tokyo, Japan
Hydride formation and associated embrittlement is very important for the integrity of the nuclear fuel cladding and the safety of the nuclear power plants. However, the mechanism of the hydride formation ahead of a crack tip of the cladding has not been well understood especially at the atomic level. In the fuel cladding, relatively strong strain field is produced by not only the hoop stress due to internal pressure but also the hydrides or the stress concentration at the crack tip. Such a strain field should affect the diffusion of the point defects produced by irradiation and solute atoms such as hydrogen, and thus, the macroscopic phenomena such as hydride formation might be affected by the results of cumulative effect of migration energy changes.In this paper, diffusion of a hydrogen atom in a strain field in hcp zirconium crystal is studied by means of kinetic Monte Carlo (kMC) computer simulation method. In order to consider the strain effect, we first calculate migration energies of a hydrogen atom under a wide range of given uniform strain field using molecular dynamics. The computation box of the kMC simulations is subdivided into small (1nm cube) cells, and the strain field of the box is modeled by assigning an appropriate constant strain field to each of the small cells. Then, effect of strain field to hydrogen atom diffusion is considered by choosing appropriate migration energy for the hydrogen jump depending on the strain field of the cell in which the hydrogen is located.In the present simulations, we place edge dislocations as a source of strong strain field in the computation box where hydrogen atoms are also uniformly distributed in the box as an initial condition. Then the crystal is heated up to a given temperature (600K) to let the hydrogen atoms diffuse in the strain field in the crystal. No interaction between the hydrogen atoms is considered in the present simulations. General tendency of the hydrogen diffusion is that hydrogen atoms agglomerate in the tension strain region, and the population of the hydrogen in the compression region decreases. Steady state is reached after 1ms. This result clearly demonstrates that the strain fields ahead of crack tips or hydrides provide locations where hydrogen atoms prefer to agglomerate.
10:00 AM - **E1.4
The Lattice Monte Carlo Method for Addressing Mass and Thermal Transport Problems in Materials.
Irina Belova 1 , Graeme Murch 1 Show Abstract
1 School of Engineering, The University of Newcastle, Callaghan, New South Wales, Australia
10:30 AM - E1.5
In-diffusion and Out-diffusion of Oxygen in Ag-MgO Composites: Analysis with Finite Element and Monte Carlo Methods.
Irina Belova 1 , Andreas Oechsner 2 1 , Graeme Murch 1 Show Abstract
1 School of Engineering, The University of Newcastle, Callaghan, New South Wales, Australia, 2 Department of Applied Mechanics, Technical University of Malaysia, Skudai, Johor , Malaysia
The in-diffusion of oxygen during the formation of Ag-MgO composites results in oxygen segregation to the metal-ceramic interface and subsequent weakening of the interface bonding. Removal of oxygen by out-diffusion can be achieved at high temperatures. These diffusion processes are simulated using finite element and Lattice Monte Carlo methods.
E2: Diffusion II
Monday PM, November 26, 2007
Constitution B (Sheraton)
11:15 AM - E2.1
Quantitative Modeling of Self-Interstitial Diffusion in Silicon.
N. Modine 1 Show Abstract
1 1132, Sandia National Labs, Albuquerque, New Mexico, United States
Predictive modeling of the early-time transient annealing of radiation damage in electronic devices requires a detailed, quantitative understanding of the behavior of the fundamental defects in the device material. The isolated self-interstitial in silicon is extremely difficult to observe experimentally, and therefore accurate theoretical results should be very valuable. We apply electronic structure calculations based on the Kohn-Sham Density Functional Theory (DFT) in concert with Kinetic Monte-Carlo (KMC) techniques to study diffusion of the silicon self-interstitial as a function of the majority and minority carrier populations. We perform carefully converged DFT calculations to identify the structures that are locally stable (stable or metastable) for each charge state and the reaction pathways and transition states for transformation between these configurations. The resulting transition barriers are incorporated into a KMC model to determine the diffusion rate as a function of temperature and carrier concentrations. For non-equilibrium carrier concentrations, fluctuations in the charge state of the defect lead to non-Arrhenius behavior via a modified Bourgoin-Corbett mechanism. The rather complex behavior of the calculated diffusion coefficient will be explained in terms of fundamental physical mechanisms. The effects of transient trapping of the silicon self-interstitial in a complex with the oxygen impurity will also be discussed.Sandia is a multiprogram laboratory operated by Sandia Corporation, a Lockheed Martin Company, for the United States Department of Energy’s National Nuclear Security Administration under contract DE- AC04-94AL85000.
11:30 AM - E2.2
Diffusion of H and C through Fe Alloys.
Donald Johnson 1 , Emily Carter 2 3 Show Abstract
1 Chemistry, Princeton University, Princeton, New Jersey, United States, 2 Mechanical and Aerospace Engineering, Princeton University, Princeton, New Jersey, United States, 3 Program in Applied and Computaitonal Mathematics, Princeton University, Princeton, New Jersey, United States
Many industrial applications, such as automotive engines or high-pressure vessels for oil recovery, require strong materials to withstand high temperature or high pressure environments. Under these harsh conditions, corrosion of steel by various reactive agents is a great concern, even under reducing atmospheres. For example, attack by H2S leads to incorporation of hydrogen into steel and subsequent embrittlement. Similarly, at high temperature, attack by CO and subsequent incorporation of C into steel leads to metal dusting and carburization. In this work, we investigate whether formation of an Fe alloy thin film on steel can curb corrosion in reducing atmospheres by limiting C and H diffusion into and through steel. In particular, we employ density functional theory to calculate rate constants and energy barriers for the diffusion of H and C through bulk FeAl and bulk Fe3Si alloys. We also investigate the reaction pathways of H/C diffusing from the alloy surface to subsurface layers, in order to determine the overall kinetics of hydrogen and carbon incorporation into these Fe alloys.
11:45 AM - **E2.3
Determination of Multicomponent Interdiffusion Coefficients and its Applications.
Yongho Sohn 1 Show Abstract
1 AMPAC-MMAE, University of Central Florida, Orlando, Florida, United States
Solid-state diffusion is a subject of great interest for its intellectual merit and practical applications in materials and coatings for advanced energy generation systems such as gas turbines and nuclear reactors. This talk will highlight the importance of multicomponent-multiphase interdiffusion with specific examples from materials and coatings used in gas turbine engines and metallic nuclear fuels. Results and analysis from both laboratory experiments and field applications are presented to emphasize the cross-fertilization of science and applications. Examples will include, for gas turbine applications, ternary interdiffusion in Ni3Al (L12) with Ir or Ta ternary alloying additions, ternary and quaternary interdiffusion in austenitic NiCr and FeNiCr alloys with Al, Si, Ge or Pd alloying additions for improved oxidation resistance, and interdiffusion between oxidation resistant coatings and superalloy substrates. Diffusion controlled degradation of metallic nuclear fuels due to interdiffusion with cladding and thermotransport will be presented with appropriate analysis for the determination of critical thermo-kinetic parameters that in turn are used for fuel performance prediction. This work was financially supported by CAREER award from National Science Foundation (NSF-CAREER DMR-0238356), U.S. Department of Energy (DE-FC26-02NT41431) subcontract through Clemson University (No. 01-01-SR103), and U.S. Department of Energy (DE-AC07-05ID14517) subcontract through Idaho National Laboratory (0051953) with technical assistance from Oak Ridge National Laboratory and National Institute of Materials Science (NIMS) of Japan.
12:15 PM - E2.4
Oxygen Diffusion in δ-Bi2O3: A Molecular Dynamics and Density Functional Theory Approach.
Dilpuneet Aidhy 1 , Susan Sinnott 1 , Juan Nino 1 , Eric Wachsman 1 , Simon Phillpot 1 Show Abstract
1 Department of Materials Science and Engineering, University of Florida, Gainesville, Florida, United States
Among all the fluorite-based electrolytes proposed for use in solid oxide fuel cells, δ-Bi2O3 has the highest oxygen ion conductivity. To achieve this fluorite crystal structure at relatively lower temperatures, δ-Bi2O3 is doped with different lanthanides . However, doping decreases the oxygen conductivity; a direct proportional dependence of oxygen conductivity on atomic polarizability and radius are observed to be the two contributing factors [1, 2]. Here, the effects of ionic polarizability and ionic radius on oxygen conductivity are separated using molecular dynamics (MD) simulations. Using MD simulations, artificial cations are generated to separate the two factors’ effects by keeping one constant and varying the other. In the pure δ-Bi2O3, we find that high bismuth polarizability is important to achieve high oxygen conductivity; the ionic radius is not important. We have further observed that neither the oxygen radius nor its polarizability plays a role in its high diffusion. In doped δ-Bi2O3, elastic stresses are developed which effect the oxygen diffusion. Using both MD and density functional theory, we also present a detailed structural and oxygen diffusion mechanism comparison between the pure and doped δ-Bi2O3. This work is supported by NASA under the grant NAG3–2930 and by DOE through the High Temperature Electrochemistry Center (HiTeC) at the University of Florida, Contract No. DE-AC05-76RL01830. 1.Jiang N., Wachsman E. D., J. Am. Cer. Soc., 82 (11) (1999), 3057.2.Wachsman E. D., Ionics, 7 (2001), 1.
12:30 PM - E2.5
Ab initio Study of Surface and Surface Oxygen Diffusion Properties of LaMnO3.
Yueh-Lin Lee 1 , Dane Morgan 1 Show Abstract
1 Materials Science Program, University of Wisconsin-Madison, Madison, Wisconsin, United States
(La,Sr)MnO3 (LSM) is the cathode catalyst generally used in commercial solid oxide fuel cells (SOFCs). However, SOFCs are significantly limited by the oxygen reduction reaction (ORR) kinetics on this material, particularly at lower temperatures. To understand the relationship between surface structure and ORR in LSM, we adopt an atomistic modeling approach based on ab initio density functional theory (DFT) simulations. Initial work has focused on understanding the energetics of surfaces, near surface vacancies, and oxygen transport in LaMnO3 (LMO). The energies for LMO surfaces obtained by simply truncating the bulk orthorhombic phase structure are calculated to be 64 meV/Å2 (100), 76 meV/Å2 (110), and 75 meV/Å2 (111). Vacancy formation energies associated with the most stable (100) surface are found to deviate dramatically from their bulk value near the surface. This will lead to orders of magnitude differences in vacancy concentrations near the surfaces, potentially strongly impacting the ORR mechanisms. The electronic correlation effects in LMO were explored with GGA+U and the relative vacancy formation energies were found to be largely insensitive to the choice of U. The oxygen vacancy and binding energetics have been incorporated into a thermokinetic model to predict the oxygen coverage on the surface, near surface vacancy concentrations, and oxygen diffusivity on the surface and in near surface layers.
12:45 PM - E2.6
Cation Migration Processes in Magnesium Aluminate Spinel.
Samuel Murphy 1 , Jonathan Ball 1 , Grimes Robin 1 , Blas Uberuaga 2 , Kurt Sickafus 2 , Roger Smith 3 Show Abstract
1 Materials, Imperial College London, London United Kingdom, 2 , Los Alamos National Laboratory, Los Alamos, New Mexico, United States, 3 Department of Mathematical Science, Loughborough University, Loughborough United Kingdom
Magnsium aluminate spinel has demonstrated a significant resistance to degradation of its mechanical properties under neutron irradiation. As a consequence, spinel has been considered for numerous applications within the nuclear industry. Spinel’s high tolerance to some types of radiation is currently believed to be a result of a high interstitial – vacancy recombination rate coupled with its ability to tolerate a significant amount of disorder between its cation sublattices. Here we report on atomistic simulations, based on pair potentials, to investigate the behaviour of residual point defects within the spinel matrix. Isolated interstitial defects of all species in spinel were found to form complex split interstitial structures. For example, the oxygen interstitial rather than being a simple Oi defect, forms an Oi-VMg-Oi structure aligned along the 〈1-10〉direction. Both the Al3+ and Mg2+ interstitial ions split across the tetrahedral magnesium site, aligned along〈110〉.Migration is an important process to consider when investigating the interstitial – vacancy recombination rate as the interstitials and vacancies may have to migrate over significant distances within the lattice in order to recombine. The present studies suggest an interesting new concept; Al3+ prefers to migrate via a mechanism rather than via a vacancy mechanism on its own sublattice or via an interstitial mechanism. Aluminium interstitial migration is limited as the Al3+ interstitial ion collapses to form an AlMg antisite specie and a split Mgi-VMg-Mgi. The predicted activstion energies for Mg2+ to be transported via a vacancy mechanism, on its native sublattice, or via an interstitial mechanism are found to be essentially the same.Previous experimental and theoretical investigations of radiation damage in spinel have indicated that the number of antisite defects are increased after interaction with radiation. As mentioned previously spinel can tolerate a significant amount of disorder on the cation sublattices, however, the barriers to formation as well as annihilation of antisite defects via a vacancy mechanism are found to be very high. Effectively decoupling the two sublattices and ensuring these processes are not promoted as a result of thermal activation.
E3: Microstructural Evolution
Monday PM, November 26, 2007
Constitution B (Sheraton)
2:30 PM - **E3.1
Nucleation and Growth of Thin Films and Minerals.
John Harding 1 Show Abstract
1 Engineering Materials, University of Sheffield, Sheffield United Kingdom
The modelling of nucleation and growth remains a basic problem in materials. It is common to a range of fields: from biomineralization and nanotechnology to corrosion and scale formation in boilers. Crystal growth is a multiscale problem in both length-scales and time-scales. It is therefore necessary to use a range of simulation methods from full-scale ab initio methods to coarse-grained or mean field techniques to capture the important physical phenomena. Frequently, the substrate can control the phase, shape and properties of the material that is grown on it. It is possible to design substrates that induce unusual phases and shapes because these are the most stable option at an early stage of the growth and thereafter become “frozen in” because the activation barrier to a thermodynamically stable state is prohibitively high. Unusual diffusion mechanisms, particularly cooperative mechanisms  are often present at surfaces and interfaces and can determine the evolution of the system. We will use a number of examples from the growth of biominerals , thin films (in particular wurtzite structures)  and growth of metals on oxides  to illustrate how a range of simulation methods can be used to understand these phenomena M.Y. Lavrentiev, N.L. Allan, J.H. Harding, D.J. Harris and J.A. Purton, Comp. Mater. Sci. 36 (2006) 54-59 J.H. Harding and D.M. Duffy; J. Mater. Chem. 16 (2006) 1105-1112 C.L. Freeman, F. Claeyssens, N.L. Allan and J.H. Harding, Phys. Rev. Lett. 96 (2006) 066102. J.A. Venables, L. Giordano and J.H. Harding; J. Phys. Cond.Mat.18 (2006) S411-S427
3:00 PM - E3.2
Heterogeneous Nucleation in Alloys using MCCask, a Parallel Monte Carlo Code in the Off-lattice Transmutation Ensemble.
Babak Sadigh 1 , Enrique Martinez 1 , Alfredo Caro 1 , Luis Zepeda-Ruiz 1 , Magdalena Caro 1 , Edmundo Lopasso 2 Show Abstract
1 , Lawrence Livermore National Laboratory, Livermore, California, United States, 2 , Centro Atomico Bariloche, 8400 Bariloche Argentina
We present a massively parallel off-lattice Monte Carlo code in the transmutation ensemble. Besides the usual constrains for the Semi-Grand Canonical ensemble a new one as been introduced to allow simulations within the miscibility gap in a segregating binary solution. The second order moment of the concentration, i. e. the dispersion around the average concentration, has been added so that the concentration can be fixed even within the miscibility gap. We derive the transition matrix for this algorithm which drives the system towards equilibrium while maintaining ergodicity and macroscopic eversibility. By implementing a subdivision of each processor's domain we solve the problem of the parallelizing algorithm in a way that minimizes communications. The code follows a route from an initial state towards a state closer to thermodynamic equilibrium, as defined by the minimum of free energy for the model considered and the state variable specifications. Heterogeneous precipitation in the FeCu system is used as example of its applications. This code opens the possibility to study dirty interfaces by placing the impurities at locations dictated by free energy criteria with no approximations.
3:15 PM - E3.3
Kinetics of Precipitation in Fe-Cu Alloys: a Multiscale Modelling.
Frederic Soisson 1 , Chu Chun Fu 1 Show Abstract
1 , CEA Saclay, Gif-sur-Yvette France
The kinetics of copper precipitation in α-iron, both during thermal ageing and under irradiation, strongly depend on the details of point defects diffusion properties. They are modelled here by Monte Carlo simulations based on a diffusion model by vacancy and interstitial jumps, which takes into the dependence of point defect concentrations and migration barriers on the local atomic environment. These parameters are fitted on ab initio data, calculated within the Density Functional Theory.For thermal ageing, the simulated precipitation kinetics are in good agreement with experimental ones. The Fe-Cu system is characterized by a low mutual solubility, which results in the formation of almost pure copper precipitates; and by a large difference between the vacancy formation in bcc iron and metastable bcc copper, which leads to a strong trapping of vacancies by the precipitates. As a result, precipitates containing up to several tens of copper atoms can be much more mobile than individual copper atoms. This result is analysed with a simple model of cluster diffusion, which suggests that a similar behaviour could be observed in other alloys.Under irradiation, the evolution of the system can be modified by the coupling between point defects and copper fluxes: we discuss the effects of these phenomena on the precipitation kinetics.
3:45 PM - E3.5
Molecular Theory of Ostwald Ripening.
Victor Burlakov 1 Show Abstract
1 Materials, University of Oxford, Oxford United Kingdom
Microstructure coarsening, or Ostwald ripening (OR), is a phenomenon often observed at the late stages of many first order phase transformations in two-phase mixtures, binary alloys, clusters on surfaces, oil-water emulsions, and during epitaxial growth and synthesis of nanoparticles. During OR small clusters of atoms/molecules dissolve transferring their mass to bigger clusters. In contrast to the classical theories of OR [1-2], which are based on Gibbs-Thomson effect, we analyze the mass exchange between spherical molecular clusters considering emitted and absorbed molecular fluxes and taking into account that some of the emitted molecules return to the cluster of origin. The return probability for molecules is found to be determined by the regime of molecular transport, and is higher for bigger clusters suggesting that the return process itself may account for a new, purely kinetic ripening mechanism in the absence of the Gibbs-Thomson effect. We show that in the limit characteristic to Wagner theory  and for molecular mean free path small compared to cluster radii the ripening process is entirely due to the kinetic mechanism. The developed molecular theory reproduces major features of the classical theories and identifies parameter regions of their applicability. It also identifies an intermediate regime  with respect to the Lifshitz-Slyozov  and Wagner  theories, when the time dependence of average cluster size and the cluster size distribution are entirely different from those predicted by the classical theories. Implications of the obtained results are discussed.References M. Lifshitz, V. V. Slezov., J. Phys. Chem. Solids, 19, 35 (1961). C. Z. Wagner, Electrochem. 65, 581 (1961). V. M. Burlakov, Phys. Rev. Lett. 97, 155703 (2006).
E4: Evolution of Microstructure and Properties
Monday PM, November 26, 2007
Constitution B (Sheraton)
4:30 PM - E4.1
Quantitative Phase Field Model for Diffusion-controlled Microstructural Evolutions in the Solid State.
Arnoldo Badillo 1 , Robert Averback 1 , Pascal Bellon 1 Show Abstract
1 Materials Science and Engineering, UIUC, Urbana, Illinois, United States
Current phase field models for diffusion-controlled evolutions in the solid state are based on phenomenological kinetic equations, such the Cahn-Hilliard diffusion equation for conserved order parameters, and as a result, simulations relying on these equations lack absolute time and space scales. We propose here a new derivation of such equations, starting from an atomistic description and using a coarse graining in space and time to define the continuous variables for the phase field model. For a binary alloy where atoms interact through pairwise interactions, this approach makes it possible to derive the coarse-grained free energy of an inhomogeneous system, as well as the kinetic equation describing the evolution of its composition field under isothermal condition. It is shown that physical parameters of the coarse grained model, including the gradient energy coefficient and the atomic mobility, display dependence with the size of the cells used for coarse graining. These dependences, however, produce the desirable result that, within a given mean-field approximation, measurable quantities such as concentration profiles at equilibrium and kinetics of phase separation are independent of the coarse graining cell size. This is at variance with past results where coarse graining was applied in the space but not in the time domain. Simulations based on this new coarse graining possess physical length and time scales, which are directly inherited from the underlying microscopic model.Coarse-grained static and kinetic equations are here derived in the point and pair mean-field approximation, and their extension to higher order approximations is outlined. Examples of the application of the model are given for spinodal decomposition in three dimensions with either a direct atom-exchange mechanism or with a vacancy assisted mechanism for atomic diffusion. We also demonstrate that this model can be applied to the evolution of defect clusters in irradiated solids by including a source term for the production of point defects and point-defect clusters, a recombination term for vacancy-interstitial annihilation, and a sink term for the elimination of point defects on clusters or on extended defects such as dislocation. The presence of an absolute space scale in our simulations makes it possible to study recoil spectrum effects, that is, to compare irradiation conditions with the same average rate of production of defects but with different size distribution of defects clusters.
4:45 PM - E4.2
Combined Experimental and Computational Investigation of Microstructural Plasticity.
Luke Brewer 1 , Corbett Battaile 3 , W. Counts 2 , Thomas Buchheit 4 Show Abstract
1 Materials Characterization, Sandia National Laboratories, Albuquerque, New Mexico, United States, 3 Computational Materials Science and Engineering, Sandia National Laboratories, Albuquerque, New Mexico, United States, 2 Microstructure Physics and Metal Forming, Max-Planck-Institut für Eisenforschung GmbH , Dusseldorf Germany, 4 Microsystem Materials, Sandia National Laboratories, Albuquerque, New Mexico, United States
A great deal is known about the evolution of plastic deformation fields at the continuum level; however, there is a lack of combined theoretical and experimental investigation of plasticity at the microstructural length scale. This study examines the evolution of plastic deformation in a given microstructural field of view as a function of macroscale strain. Well-annealed tensile bars of high purity Ni were strained in steps of 0, 1%, 5%, and 10% strain. At each strain increment, the distribution of plastic deformation was mapped using electron backscatter diffraction (EBSD) orientation and misorientation maps. The misorientation maps from this data were compared to misorientation maps calculated using both a local and a novel non-local crystal plasticity-based finite element method (FEM). The degree of agreement between experiment and simulation differed upon the length scale of the comparison. The distribution of intragranular crystal rotation was qualitatively better captured by the local FEM model, while the non-local model had closer agreement with experiment in the total amount of intragranular crystal rotation. Sandia is a multiprogram laboratory operated by Sandia Corporation, a Lockheed Martin Company, for the United Stated Department of Energy (DOE) under contract DE-AC0494AL85000.
5:00 PM - E4.3
Modeling the Effect of Particle Size Distribution on Platinum Surface Area Loss in PEMFC Cathodes.
Edward Holby 1 , Dane Morgan 1 , Wenchao Sheng 2 , Yang Shao-Horn 2 Show Abstract
1 Materials Science Program, University of Wisconsin-Madison, Madison, Wisconsin, United States, 2 , Massachusetts Institute of Technology, Cambridge, Massachusetts, United States
The long-term durability of platinum catalysts in proton exchange membrane fuel cell cathodes is an important issue in achieving the Department of Energy’s 5000 hour PEMFC lifetime goal. Fuel cell efficiency is decreased as electrochemically active surface area (ECASA) of platinum is lost. It is important to better understand the mechanisms of platinum surface area loss in the cathode and if they can be mitigated in order to meet lifetime requirements. We have used a kinetic model to explore the influence of Particle Size Distribution (PSD) on the loss of ECASA under both constant potential holds and cycled potentials for carbon supported Pt catalysts (Pt/C). This model extends recent work from Darling and Meyers1, and includes the effects of platinum dissolution into the ionomer, platinum reprecipitation both on and off the carbon support and surface oxidation and passivation. All of these strongly couple to particle size, due to the Gibbs-Thomson effect. It is found that the mean diameter and functional form of the initial PSD are important considerations in the long-term loss of cathodic platinum ECASA in PEMFCs. Also, we confirm that changes in the PSD upon heat treatment help to explain the improved durability associated with heat-treating Pt/C catalysts.1. Darling, R.M. and J.P. Meyers, Kinetic Model of Platinum Dissolution in PEMFCs. Journal of the Electrochemical Society, 2003 150(11): p. A1523-A1527.
5:15 PM - E4.4
Relaxation of Biaxial Strain in Ultra-Thin Films of Face-centered Cubic Metals through Ductile Void Growth, Grain Nucleation, and Structural Phase Transitions.
Kedarnath Kolluri 1 , M. Rauf Gungor 1 , Dimitrios Maroudas 1 Show Abstract
1 Chemical Engineering, University of massachusetts, Amherst, Massachusetts, United States
Nanometer-scale-thick metallic films are used increasingly in modern technologies, from microelectronics to various areas of nanofabrication. Processing of such materials generates voids and structural defects and may induce structural transformations. The structural evolution of these materials during the relaxation of electro- and thermo-mechanical stresses underlies many mechanical reliability problems. Improvement of the reliability, functionality, and performance of nano-scale devices requires a fundamental understanding of the atomistic mechanisms of strain relaxation in thin films and the associated structural changes in order to establish links between structural evolution and mechanical behavior.In this presentation, we report a computational analysis of atomistic mechanisms of relaxation of biaxially applied tensile strain over a broad range of strain levels in free standing ultra-thin copper films with the film surface oriented normal to the  crystallographic direction; the films studied contained voids extending throughout the thickness of the film with the void axes oriented along . Our study is based on isothermal-isostrain large-scale molecular-dynamics simulations, using an embedded-atom-method parameterization to describe the interatomic interactions. Our analysis reveals four regimes of structural response as the applied biaxial strain level, ε, increases. In the first regime (ε < 2%), the films respond elastically. In the second regime (2% < ε < 8%), strain relaxation is governed by ductile void growth that is mediated by dislocation emission from the surface of the void and formation of a plastic zone around the void. Dislocation glide and dislocation-dislocation interactions lead to formation of multi-layered microtwin boundaries and Lomer-Cottrell locks in the thin-film material. In addition, dislocation arrays are nucleated at the surface of the void and lead eventually to lattice rotations and nanometer-sized grain formation in the thin film. In the third regime (8% ≤ ε < 12%), dislocations and dislocation arrays are nucleated at the surface of the thin film in addition to those that are nucleated at the void surface. The fourth regime (ε ≥ 12%) is characterized by the nucleation of a high density of dislocation arrays, leading to a phase transformation of the thin film from a face-centered cubic (fcc) lattice structure to a hexagonal close-packed (hcp) lattice structure. Throughout the strain range examined, we characterize the surface roughness of the strain-relaxed films and calculate their elastic properties. Furthermore, we examine the validity of static lattice stability criteria in capturing the phase transformation onsets in the metallic thin films and construct the underlying free-energy landscapes in order to determine accurately these transformation onsets.
5:30 PM - E4.5
Simulation of the Columnar-to-Equiaxed Transition in Alloy Solidification -- the Effect of Alloy Solidification Range, Nucleation Undercooling and Density of Nuclei on the Transition.
Huijuan Dai 1 , Hongbiao Dong 1 2 , Helen Atkinson 1 , Peter Lee 3 Show Abstract
1 Department of Engineering, University of Leicester, Leicester United Kingdom, 2 Precision Casting Facility, Rolls-Royce plc., Derby United Kingdom, 3 Department of Materials, Imperial College London, London United Kingdom
5:45 PM - E4.6
Abnormal Grain Growth of Fe-3%Si Steel Approached by Solid-State Wetting Mechanism.
Kyung-Jun Ko 1 , Pil-Ryung Cha 2 , Jong-Tae Park 3 , Jae-Kwan Kim 3 , Nong-Moon Hwang 1 Show Abstract
1 Material Science and Engineering, Seoul National University, Seoul Korea (the Republic of), 2 Advanced Materials Engineering, Kook-Min University, Seoul Korea (the Republic of), 3 POSCO Technical Research Laboratories, POSCO, Pohang Korea (the Republic of)
Abnormal grain growth (AGG) which is also called secondary recrystallization takes place in many metallic systems especially after recrystallization of deformed polycrystals. A famous example is the evolution of the strong Goss texture after secondary recrystallization in Fe-3%Si steel. Extensive efforts have been made to identify the growth advantage of Goss grains due to its scientific and technological interests. However, the Goss secondary recrystallization or AGG still remains as a puzzling phenomenon since its first report by Goss in 1935.In this study, we suggest a new approach to the growth advantage of AGG by solid-state wetting, where a grain wets or penetrates the grain boundary or the triple junction of neighboring grains just as the liquid phase wets along the grain boundary or the triple junction. If the energy sum of the two grain boundaries is lower than the energy of the other grain boundary in contact at the triple junction, the high energy grain boundary will be replaced by the two low energy grain boundaries through the wetting process. Once the solid-state wetting occurs, the triple junction, which tends to be a rate-determining step in grain growth, migrates much faster than the grain boundaries in contact and therefore, the triple junction constraint in grain growth disappears. According to the solid-state wetting mechanism, the sub-grain boundaries which have very low energy can play a critical role in inducing AGG, because they provide grains with the exclusive growth advantage by increasing the wetting probabilityAGG was approached based on the solid-state wetting mechanism by both experiments and the 3-dimensional grain growth simulation using Phase-Field model. The effect of solid-state wetting on AGG by phase field model (PFM) simulation is compared the microstructural feature of AGG such as the formation of numerous island and peninsular grains observed experimentally in Fe-3%Si alloy. The misorientation angles between island or peninsular grains and the abnormally-growing Goss grains were mostly low angles less than 15o with some high angle boundary of coincidence lattice relationship, which indicates that they have the low grain boundary energy. Very high frequency of island and peninsular grains formed at or near the growth front of abnormally growing grains could be best explained by solid-state wetting. 3-sided or 4-sided grains with negative curvatures, which are formed near the growth front of abnormally growing grains, are the two-dimensional section vertical to the triple junction wetting in three-dimensional polycrystalline structure. These 3-sided or 4-sided grains were observed near the abnormally-growing Goss grains and were shown to be connected to the abnormally-growing grain by sequential polishing and observation.
E5/T5: Joint Session: Modeling Defects in Nuclear Materials I
Tuesday AM, November 27, 2007
Constitution B (Sheraton)
9:30 AM - **E5.1/T5.1
Defects and Impurities in Nuclear Materials from First Principles.
Chu Chun Fu 1 , F. Willaime 1 Show Abstract
1 DEN/DMN, SRMP-CEA, Sacaly, Gif sur Yvette France
Ferritic steels play an central role in metallurgical and nuclear technology, in particularas structural materials for fission and future fusion nuclear reactors.Although their mechanical properties have been extensively investigated, little is knownabout the structural, electronic and magnetic properties at the origin of their macroscopicbehavior. First principles calculations within the Density Functional Theory (DFT) provide suchinformation at atomic scale, which is not directly accessible through experiments.However, the application of first principles studies in this field is rather new compared withother Solid State Physics and Material Science disciplines.The objective of this talk is to report key contributions of recent DFT studies which allowto reconcile theory and experiments:The determination of the 3D migration of self-interstitial atoms (SIAs) -- elementary defectscreated by irradiation -- induced an overall revision of the widely accepted picture of damageaccumulation under irradiation predicted by empirical potentials . The coupled ab initio andmesoscopic kinetic Monte Carlo simulation provided strong evidence to clarify long-standingcontroversial interpretations of electrical resistivity recovery experiments concerning themobility of vacancies, SIAs, and their clusters . The behavior of Carbon, one of theessential alloying element in steels, and of impurities such as Phosphorus which are responsibleof reactor pressure vessel (RPV) steels embrittlement will also be discussed in detail . C. C. Fu, F. Willaime and P. Ordejon, Phys. Rev. Lett. 92, 175503 (2004) C. C. Fu et al. Nature Mater. 4, 68 (2005) E. Meslin et al. Phys. Rev. B 75, 094303 (2007)
10:00 AM - **E5.2/T5.2
Computer Simulation of Defect Properties in Irradiated Metals.
David Bacon 1 , Alexander Barashev 1 , Andrew Calder 1 , Yuri Osetsky 2 Show Abstract
1 Engineering, University of Liverpool, Liverpool United Kingdom, 2 , Oak Ridge National Laboratory, Oak Ridge, Tennessee, United States
10:30 AM - E5.3/T5.3
Ab initio-based Radiation-induced Segregation Modeling in Fe-Ni-Cr Alloys.
Julie Tucker 1 , Todd Allen 1 , Dane Morgan 2 Show Abstract
1 Nuclear Engineering & Engineering Physics, University of Wisconsin - Madison, Madison, Wisconsin, United States, 2 Materials Science and Engineering, Univerisity of Wisconsin - Madison, Madison, Wisconsin, United States
High concentrations of point defects, such as those created in radiation environments, can cause severe material degradation as they migrate and cluster. Radiation induced segregation (RIS), the process by which the local composition of an alloy is altered near point defect sinks, is a phenomenon that has concerned the nuclear industry for decades. While substantial progress has been made in the area of RIS prediction by empirical fitting, many questions remain about the diffusion mechanisms of point defects and how they are affected by local environment changes in multi-component alloys. This research uses ab initio methods to determine diffusion coefficients associated with both vacancy and interstitial migration in Ni rich fcc Fe-Ni-Cr alloys. We find that the alloy kinetics differs significantly from that predicted by simple extrapolation of empirical fits to high-temperature data. The calculated diffusion coefficients are used to parameterize a rate theory type RIS model.
10:45 AM - E5.4/T5.4
Microstructural Evolution Under High Flux Irradiation of Dilute Fe-CuNiMnSi Alloys Studied by Atomic Kinetic Monte Carlo Model – Effect of the Self Interstitials.
Christophe Domain 1 2 , Edwige Vincent 1 2 , Raoul Ngayam Happy 2 1 , Charlotte Becquart 2 Show Abstract
1 MMC, EDF R&D, Moret sur Loing France, 2 Laboratoire de Métallurgie Physique et Génie des Matériaux, UMR 8517, Villeneuve d'Ascq France
E6/T6: Joint Session: Modeling Defects in Nuclear Materials II
Tuesday PM, November 27, 2007
Constitution B (Sheraton)
11:30 AM - **E6.1/T6.1
Multi-Scale Simulations of Ion-Solid Interactions in SiC and GaN.
Fei Gao 1 , William Weber 1 , Haiyan Xiao 2 3 , Zhouwen Rong 1 , Yanwen Zhang 1 , Ram Devanathan 1 , Lumin Wang 2 , Xiaotao Zu 3 Show Abstract
1 MS K8-93, Pacific Northwest National Laboratory, Richland, Washington, United States, 2 Department of Nuclear Engineering and Radiological Sciences, University of Michigan, Ann Arbor, Michigan, United States, 3 Department of Applied Physics, University of Electronic Science and Technology of China, Chengdu, Sichuan, China
Recent progress in multi-scale simulations of fundamental ion-solid interactions in SiC and GaN is discussed. Large-scale ab initio simulation methods (up to a few thousand atoms) have been developed for the study of ion-solid interactions in materials, and these methods are employed to investigate the properties of defects and defect clusters. Atomic structures, formation energies and binding energies of these small clusters are determined, and their relative stabilities are described. Furthermore, ab initio molecular dynamics methods have been used to calculate the threshold displacement energy surface and to simulate the primary damage states for the PKA (primary knock-on atom) energies up to 1 keV in SiC and GaN. These simulations provide significant insights into electronic effects on ion-solid interaction processes. Molecular dynamics (MD) have been used to simulate the high-energy ion-solid interactions with energies up to 50 keV, while kinetic Monte Carlo methods have been employed to investigate the recovery of defects during annealing, with input parameters determined by ab initito calculations and empirical potential MD simulations. A large number of 10 keV displacement cascades are randomly generated in a model crystal to simulate multiple ion-solid interactions and damage accumulation, as well as the mechanisms controlling the crystalline-to-amorphous transition. The amorphous-to-crystalline (a-c) transition has also been studied using MD methods, with simulation times of up to a few hundred ns, in 4H-SiC at temperatures between 1000 and 2000 K. Based on a model developed in the previous annealing simulations of 3C-SiC, the activation energy spectra for recrystallization along the three directions have been determined. In general, the activation spectra show that there are a number of activation energy peaks associated with different recrystallization processes. These activation energy values for full recrystallization are in the range of from 1.2 to 1.7 eV for the amorphous layers with the a-c interfaces along [-12-10] and [-1010] directions, and 1.1 to 2.3 eV for the amorphous layer with the a-c interfaces along  direction.
12:00 PM - E6.2/T6.2
Defect Structure and Stability in Uranium and Zirconium Nitrides.
Robin Grimes 1 , Eugene Kotomin 2
1 materials, Imperial College London, London United Kingdom, 2