ATOMISTIC MODELING OF DIFFUSION IN ALUMINUM
|
|
- Moris Tyler
- 5 years ago
- Views:
Transcription
1 ATOMISTIC MODELING OF DIFFUSION IN ALUMINUM S. GRABOWSKI, K. KADAU and P. ENTEL Theoretische Physik, Gerhard-Mercator-Universität Duisburg, Duisburg, Germany (Received...) Abstract We present molecular-dynamics simulations of self-diffusion in Al. In order to facilitate the description of elastic and vibrational properties as well as vacancy migration, an embedded-atom method potential was used in the simulations. This potential was specifically designed to reproduce the T = 0 K equation of state of Al obtained by ab initio total-energy calculations. We show that the temperature dependent self-diffusion coefficient obeys an Arrhenius law and that the resulting dynamical migration energy is slightly larger than the static migration energy obtained by using classical rate theory. Keywords: Self-diffusion, Vacany migration, Molecular-dynamics simulation 1. INTRODUCTION Molecular-dynamics can be a useful tool in the study of solid-state reactions provided one is able to cover a sufficiently large time range. By numerically solving the classical equations of motion, the time evolution of a system containing thousands of atoms can be simulated up to the nano-second range. This, in principle, is sufficient to calculate the temperature dependent diffusion coefficient D(T) and to investigate individual atomic jump processes in systems which have a natural high mobility of atoms. Comprehension of mechanisms underlying atomic diffusion in metals can be gained from direct observation of individual atomic jumps occurring during the simulation. Depending on the accurateness of interaction potentials between atoms specific aspects of diffusion can be studied. Of interest are, for example, deviations of D(T) from an Arrhenius law as observed for several metals. Such deviations can be investigated by analyzing the contributions of different jump mechanisms (Lorenzi and Ercolessi, 1992). Another interesting aspect is the dynamical reduction of the migration energy for self-diffusion due to low lying 1
2 phonon modes in bcc metals (Köhler and Herzig, 1988). This effect can be studied by a close inspection of the trajectories of the jumping atom and the motional behavior of surrounding atoms (Mikhin and Osetsky, 1993). Furthermore the validity and limitations of theoretical concepts like classical rate theory can be investigated without essential simplification of the N-body potential energy surface (Jacucci, 1984; Schott, Fähnle and Seeger, 1997). Here we report first results of molecular-dynamics simulations in which we employed an embedded-atom method potential specially designed for a reliable description of self-diffusion in Al. 2. COMPUTATIONAL METHODS 2.1. Construction of an optimized EAM potential The key ingredient of any molecular-dynamics simulation is the force law describing the interactions between the atoms. Embedded-atom method (EAM) potentials of the form proposed by Daw and Baskes (1984) are appropriate to describe lattice defects and lattice dynamics in metals. In the EAM case, the potential energy of a single atom is considered as the energy required to embed the atom into the local electron density provided by the other atoms plus a core-core repulsion potential which takes the form of a pair potential. The potential energy is written as E = i F(ρ i ) + i<j Φ (r ij ), (1) where the summation is over all atoms and r ij is the distance between the atoms i and j. The function F(ρ i ) represents the embedding energy of atom i depending on the background charge density ρ i = j i ρ at j (r ij ) (2) which is the sum over all atomic contributions ρ at j (r ij) from the neighboring atoms. We use double-ζ wave functions Ψ of Clementi and Roetti (1974) for the calculation of the atomic contributions ρ at (r) = N 3s Ψ 3s (r) 2 + N 3p Ψ 3p (r) 2 ρ c. (3) Here N 3s = 2.0 and N 3p = 1.0 are the numbers of s- and p-symmetry valence electrons of Al. In order to get a finite interaction range, ρ at (r) is truncated at a cutoff radius R c and the constant shift ρ c = ρ at (R c ) is applied to obtain a continuous behavior of the embedding charge density at r = R c. The second term in Eq. 1 represents a screened Coulomb potential, Φ(r) = 2 Z(r)Z(r) r (4) 2
3 with the effective charge Z(r). Here the factor 2 stems from the atomic units used in this paper. The function F(ρ) and the effective charge Z(r) are represented by cubic spline functions. In order to construct an EAM potential which allows for a reasonable description of defect and lattice vibrational properties, the parameters of these functions have been determined by a weighted least square fit to experimental data of Al. It turned out that a minimum set of experimental data is given by the lattice constant a 0, sublimation energy E s, elastic constants C 11, C 12 and C 44, vacancy formation energy Ev f, and some selected phonon frequencies of Al. These experimental data are listed in Table II. However, EAM potentials which reproduce these experimental data in the simulations, do not automatically guarantee a proper description of locally compressed atomic configurations occurring during the diffusional jumps of atoms. The static vacancy migration energy calculated from an EAM potential which has been fitted to the data in Table II, is about a factor two larger than the corresponding experimental value of (0.61 ± 0.03) ev (Erhart, 1990). The description of mono-vacancy diffusion on EAM basis can be improved if the pair potential Φ(r) is more carefully adjusted for interatomic distances which are smaller than the equilibrium lattice constant, i.e. for 0.85 < r/a 0 < 1.0, because the strong repulsive forces between the jumping atom and its nearest neighbors in the unrelaxed saddle point configuration are mainly due to Φ(r). By additionally fitting the EAM potential to the T = 0 K equation of state (see Fig. 1) calculated on ab initio basis, the resulting static vacancy migration energy is much closer to the experimental value. The potential used in our simulations was fitted to the equation of state for lattice constants in the range 0.9 < r/a 0 < 1.05, whereby the compressed lattice constant of 0.9 a 0 corresponds to an external pressure of 40 GPa. If the potential is fitted to even smaller lattice constants, i.e. higher external pressures, then the static migration energy can be reduced to 0.65 ev, which is very close to the experimental value cited above. But in this case we failed to describe the vibrational and elastic properties with comparable accuracy. In order to exclude such unphysical dynamical behavior, a subtle choice has to be made of which compressional states have to be considered in the fit. The effective pair potential corresponding to the fitted EAM potential in the case of Al is shown in Fig Molecular-dynamics In order to calculate the temperature dependence of the self-diffusion coefficient D(T), the EAM potential defined in Tab. I was used in the moleculardynamics study [for details of the method we refer to Allen and Tildesley (1987)] of fcc elementary cells (2048 lattice sites) with periodic boundary conditions. In order to observe a sufficient rate of elementary diffusion jumps, a rather high concentration of vacancies was used corresponding to 16 vacancies in the simulation box. The vacancies were initially arranged 3
4 on a bcc lattice with lattice constant 4 a (with a the lattice constant of the underlying fcc lattice). For each temperature T the lattice constant a was chosen to correspond to its zero pressure value. The classical equations of motion were integrated with the velocity Verlet algorithm with a time step of 1.5 fs and the temperature was controlled by the Nosé-Hoover thermostat (Nosé, 1984; Hoover, 1985). The box volume was kept fixed throughout each run Diffusion coefficient and migration energy The self-diffusion coefficient D can be calculated from the Einstein relation D = lim t 2 r(t) 6t and the long-time limit of the mean square displacement defined by 2 r(t) = 1 N (5) N [r i (t) r i (0) ] 2, (6) i=1 where r i (t) is the time dependent position of atom i; N is the number of atoms in the simulation box (Rapaport, 1995). In order to calculate D(T) we performed molecular-dynamics runs for temperatures between 770 K and 950 K. For each temperature the mean square displacement was sampled over a time interval of length 3.75 ns and D(T) was calculated from the slope of the linear least mean-square fit to 2 r(t). Examples are plotted in Fig DISCUSSION An Arrhenius plot of the resulting diffusion coefficient D(T) is shown in Fig. 4. From an exponential least mean-square fit to our D(T) data we obtain the dynamical migration energy E m = 0.94 ev as well as the prefactor D 0 = m 2 s 1, i.e. the diffusion coefficient can be described by the Arrhenius law ( D(T) = D 0 exp E ) m. (7) k B T From static calculations using classical rate theory (Vineyard, 1957; Flynn, 1972) we obtain the migration energy E m = 0.75 ev which is 0.19 ev smaller than our calculated dynamical migration energy. In order to estimate the error of the dynamical migration energy E m obtained from our D(T) data, we calculated the 99% confidence limit for E m within random sampling theory of regression (Spiegel, 1976). The resulting error corresponds to ev i.e., the dynamical migration energy calculated from D(T) is by at least 0.1 ev larger than the static value. 4
5 4. SUMMARY The aim of this work was to specifically design an EAM potential which allows for a reasonable description of vacancy diffusion in Al. It turned out that by an additional fit of the EAM potential to the T = 0 K equation of state, the vacancy migration energy calculated within the framework of classical rate theory could be adjusted to the experimental value. For the potential defined in Table I the static migration energy obtained is 0.75 ev. The simulated temperature dependence of the diffusion coefficient exhibits an Arrhenius law behavior and the dynamical migration energy obtained from the slope of the Arrhenius plot in the temperature range K is given by 0.94 ± ev, which is at least 0.1 ev larger than the static value. We do not have an obvious explanation for the difference between the static and the dynamic value of the vacancy migration energy. It is a well known fact that classical rate theory yields an upper limit for the vacancy jump rate and, therefore, overestimates the diffusion coefficient because return jumps are not taken into account (Hänggi, Talkner and Borkovec, 1990). According to the analytical and numerical investigations of Lorenzi et al. (1987) and Marchese et al. (1987), classical rate theory should be an excellent approximation for fcc model crystals with Lennard-Jones-type of interactions because return jumps do not lower the vacancy jump rate by more than 10%. But this error of classical rate theory is much too small to explain the obseved difference between the static and the dynamic migration energy obtained in the case of Al. Contributions to the diffusion coefficient D due to other jump processes like correlated double jumps cannot be excluded at least for the highest simulation temperatures. Therefore, a slight increase of the activation energy obtained from the Arrhenius plot can be expected because migration energies of correlated double jumps are approximately twice the migration energy of single vacancy jumps (Lorenzi and Ercolessi, 1992). Probably an intricate combination of systematical errors has occurred which could originate from vacancy-vacancy interactions due to the large concentration of vacancies used here, periodic boundary conditions and insufficient supercell size as well as insufficient statistics for single vacancy jumps at low temperatures. Nonetheless this is the first observation of elementary diffusion jumps seen in molecular-dynamics simulations for Al. References Allen, M. P. and D. J. Tildesley (1987). Computer Simulation of Liquids. Clarendon Press, Oxford. Clementi, E. and C. Roetti (1974). Atomic Data and Nuclear Data Tables, volume 14. Academic Press, New York. Daw, M. S. (1990). The embedded atom method: A review. In R. M. Nieminen, 5
6 M. J. Puska and M. J. Manninen, editors, Many-Atom Interactions in Solids, Springer Proceedings in Physics 48. Springer, Berlin. Daw, M. S. and M. I. Baskes (1984). Embedded-atom method: Derivation and application to impurities, surfaces and other defects in metals. Phys. Rev. B, 29, Erhart, P. (1990). Properties and interactions of atomic defects in metals and alloys. In H. Ullmaier, editor, Landolt-Börnstein, volume 25 of New Series Group III. Springer, Berlin. Flynn, C. P. (1972). Point Defects and Diffusion. Clarendon Press, Oxford. Gui, J. (1994). Embedded-atom method study of the effect of the order degree on the lattice parameters of Cu-based shape-memory alloys. J. Phys.: Condens. Matter, 6, Hänggi, P., P. Talkner and M. Borkovec (1990). Reaction-rate theory: Fifty years after Kramers. Rev. Mod. Phys., 62, 251. Hoover, W. G. (1985). Canonical dynamics: Equilibrium phase-space distributions. Phys. Rev. A, 31, Jacucci, G. (1984). Defect calculations beyond the harmonic model. In G. E. Murch and A. S. Nowick, editors, Diffusion in Crystalline Solids, 431. Academic Press. Köhler, U. and C. Herzig (1988). On the correlation between self-diffusion and the low-frequency LA 2/3 111 phonon mode in b.c.c. metals. Phil. Mag. A, 58, 769. Kresse, G. and J. Furthmüller (1996). Efficient iterative schemes for ab initio total-energy calculations using a plane wave basis set. Phys. Rev. B, 54, Lorenzi, G. D. and F. Ercolessi (1992). Multiple jumps and vacancy diffusion in a face-centered-cubic metal. Europhys. Lett., 20, 349. Lorenzi, G. D., G. Jacucci and C. P. Flynn (1987). Jump rate of the fcc vacancy in the short-memory-augmented-rate-theory approximation. I. Difference Monte Carlo sampling for the Vineyard rate. Phys. Rev. B, 36, Marchese, M., G. Jacucci and C. P. Flynn (1987). Jump rate of the fcc vacancy in the short-memory-augmented-rate-theory approximation. I. dynamical conversion coefficient and isotope-effect factor. Phys. Rev. B, 36, Mikhin, A. G. and Y. N. Osetsky (1993). On anomalous self-diffusion in body centered cubic metals: A computer simulation study. J. Phys.: Condens. Matter, 5, Nosé, S. (1984). A molecular dynamics method for simulations in the canonical ensemble. Mol. Phys., 52, 255. Rapaport, D. C. (1995). The Art of Molecular Dynamics Simulation. Cambridge University Press, Cambridge. 6
7 Rubini, S. and P. Ballone (1993). Quasiharnonic and molecular-dynamics study of the martensitic transformation in Ni-Al alloys. Phys. Rev. B, 48, 99. Schott, V., M. Fähnle and A. Seeger (1997). Molecular-dynamics study of selfdiffusion in Na: Validity of transition-state theory. Phys. Rev. B, 56, Spiegel, M. R. (1976). Theory and Problems of Statistics. McGraw-Hill. Vineyard, G. H. (1957). Frequency factors and isotope effects in solid state rate processes. J. Phys. Chem. Solids, 3,
8 Table I: Parameters determining the embedding energy function F(ρ) and the effective charge function Z(r). Energies are given in units of Ry and charges in units of e. The equilibrium charge density is given by ρ 0 = a.u., the equilibrium lattice constant is given by a 0 = a.u. and the cut off radius R c is 10.1 a.u.. ρ/ρ 0 F F r/a 0 Z Z
9 Table II: Experimental data used in the fit of the EAM potential and the corresponding calculated values. Phonon frequencies are given in mev. The calculated values correspond to an EAM potential which has also been fitted to the ab initio data shown in Fig.1. Experimental phonon frequencies have been taken from Rubini and Ballone (1993) and the other experimental data have been taken from Gui (1994). Expt. Calc. a 0 (a.u.) E s (ev) E f v (ev) C 11 (10 12 dyn/cm 2 ) C 12 (10 12 dyn/cm 2 ) C 44 (10 12 dyn/cm 2 ) q = 2π a 0 [100] q = 2π a 0 [110] q = 2π a 0 [111]
10 0.21 E pot (V ) (Ry) V (a.u.) Figure 1: Potential energy per atom E pot (V ) as a function of the volume V per atom. Curves for two different EAM potentials fitted to the experimental data listed in Table II are shown. The E pot (V ) curve which corresponds to an EAM potential fitted to the experimental data of Table II, is represented by a dashed line and the full line represents the E pot (V ) curve corresponding to an EAM potential with an additional fit to ab initio data, represented by full circles. Open and full circles represent E pot (V ) calculated by the Vienna ab initio Simulation Package (Kresse and Furthmüller, 1996) using ultrasoft Vanderbilt pseudopotentials. The ab initio data are shifted so that their minimum energy and the corresponding lattice constant coincide with the experimental values of E s and a 0 in Table II. 10
11 F(ρ) (mry) 300 φ (r) (mry) ρ (ρ 0 ) r (a 0 ) F eff (ρ) (mry) 195 φ eff (r) (mry) ρ (ρ 0 ) r (a 0 ) Figure 2: Embedding energy, F(ρ), and pair potential, Φ(r), for the EAM potential defined in Table I. The corresponding effective embedding energy, F eff (ρ), and effective pair potential, Φ eff (r), are also shown. According to Daw (1990) the effective potentials are defined by Φ eff (r) = Φ(r) + 2 c ρ(r) and F eff (ρ) = F(ρ) c ρ with c = F (ρ 0 ), where ρ 0 is the equilibrium embedding charge density from Table I. Both EAM potentials, defined by the pairs (F, Φ) and (F eff, Φ eff ), are equivalent. The second minimum of F eff, approximately at 2.0 ρ 0, results from a fit of the potential to the equation of state for high external pressures, i.e., smaller lattice constants. 11
12 16 < r 2 (t) > (a.u.) K 910 K 870 K 810 K t (ps) Figure 3: Mean square displacement, 2 r(t), for different temperatures. Linear least mean-square fits to these data are represented by dashed lines. The simulation time of 3.75 ns corresponds to time steps (a time step correponds to 1.5 fs). 12
13 10 11 D (T) (m 2 s 1 ) / k B T (ev 1 ) Figure 4: Arrhenius plot for the Diffusion coefficient (full circles) calculated from the mean-square displacements by using the Einstein relation. The exponential least-mean square fit to these data is given by the straight line. Simulation temperatures are between 770 and 950 K. From the result of the exponential least mean-square fit to D(T) we obtain the dynamical migration energy E m = 0.94 ev as well as the prefactor D 0 = m 2 s 1. Error bars correspond to standard deviations calculated according to σ/d = 1/ n with n the number of diffusional jumps in the simulation time. 13
ELEMENTARY DIFFUSION PROCESSES IN AL-CU-ZN ALLOYS: AN AB INITIO STUDY
ELEMENTARY DIFFUSION PROCESSES IN AL-CU-ZN ALLOYS: AN AB INITIO STUDY S. GRABOWSKI AND P. ENTEL Theoretische Tieftemperaturphysik, Gerhard-Mercator-Universität Duisburg, 47048 Duisburg, Germany E-mail:
More informationAb initio molecular dynamics simulation on temperature-dependent properties of Al Si liquid alloy
INSTITUTE OF PHYSICSPUBLISHING JOURNAL OFPHYSICS: CONDENSED MATTER J. Phys.: Condens. Matter 16 (4) 57 514 PII: S953-8984(4)7691-8 Ab initio molecular dynamics simulation on temperature-dependent properties
More informationPotentials, periodicity
Potentials, periodicity Lecture 2 1/23/18 1 Survey responses 2 Topic requests DFT (10), Molecular dynamics (7), Monte Carlo (5) Machine Learning (4), High-throughput, Databases (4) NEB, phonons, Non-equilibrium
More informationAn EAM potential for the dynamical simulation of Ni-Al alloys
J. At. Mol. Sci. doi: 10.4208/jams.022310.031210a Vol. 1, No. 3, pp. 253-261 August 2010 An EAM potential for the dynamical simulation of Ni-Al alloys Jian-Hua Zhang, Shun-Qing Wu, Yu-Hua Wen, and Zi-Zhong
More informationFrom Atoms to Materials: Predictive Theory and Simulations
From Atoms to Materials: Predictive Theory and Simulations Week 3 Lecture 4 Potentials for metals and semiconductors Ale Strachan strachan@purdue.edu School of Materials Engineering & Birck anotechnology
More informationEffective potentials for quasicrystals from ab-initio data
Effective potentials for quasicrystals from ab-initio data Peter Brommer and Franz Gähler Institut für Theoretische und Angewandte Physik Universität Stuttgart August 31, 2005 Abstract Classical effective
More informationOrdering and correlation of cluster orientations in CaCd 6
Philosophical Magazine, Vol. 00, No. 00, DD Month 200x, 1 5 Ordering and correlation of cluster orientations in CaCd 6 PETER BROMMER, FRANZ GÄHLER and MAREK MIHALKOVIČ Institut für Theoretische und Angewandte
More informationPBS: FROM SOLIDS TO CLUSTERS
PBS: FROM SOLIDS TO CLUSTERS E. HOFFMANN AND P. ENTEL Theoretische Tieftemperaturphysik Gerhard-Mercator-Universität Duisburg, Lotharstraße 1 47048 Duisburg, Germany Semiconducting nanocrystallites like
More informationClassical potentials for metals
Classical potentials for metals About 80 % of all elements are metals. The crystal structures of the elements are distributed as follows: FCC 15 HCP 26 BCC 16 All other metals 13 So if we can describe
More informationMolecular Dynamics Simulation Study of Transport Properties of Diatomic Gases
MD Simulation of Diatomic Gases Bull. Korean Chem. Soc. 14, Vol. 35, No. 1 357 http://dx.doi.org/1.51/bkcs.14.35.1.357 Molecular Dynamics Simulation Study of Transport Properties of Diatomic Gases Song
More informationThermal and Mechanical Properties of Pt-Rh Alloys
arxiv:cond-mat/9611241v1 [cond-mat.mtrl-sci] 28 Nov 1996 Thermal and Mechanical Properties of Pt-Rh Alloys G. Dereli, T. Çağın, M. Uludoğan, M. Tomak Department of Physics, Middle East Technical University,
More informationModeling the sputter deposition of thin film photovoltaics using long time scale dynamics techniques
Loughborough University Institutional Repository Modeling the sputter deposition of thin film photovoltaics using long time scale dynamics techniques This item was submitted to Loughborough University's
More informationAB INITIO MOLECULAR-DYNAMICS SIMULATIONS OF ADSORPTION OF DYE MOLECULES AT SURFACES
AB INITIO MOLECULAR-DYNAMICS SIMULATIONS OF ADSORPTION OF DYE MOLECULES AT SURFACES M. SUGIHARA, H. MEYER, AND P. ENTEL Theoretische Tieftemperaturphysik, Universität Duisburg, 47048 Duisburg, Germany
More informationVIRTUAL LATTICE TECHNIQUE AND THE INTERATOMIC POTENTIALS OF ZINC-BLEND-TYPE BINARY COMPOUNDS
Modern Physics Letters B, Vol. 16, Nos. 5 & 6 (2002) 187 194 c World Scientific Publishing Company VIRTUAL LATTICE TECHNIQUE AND THE INTERATOMIC POTENTIALS OF ZINC-BLEND-TYPE BINARY COMPOUNDS LIU YING,
More informationSupporting Information Tuning Local Electronic Structure of Single Layer MoS2 through Defect Engineering
Supporting Information Tuning Local Electronic Structure of Single Layer MoS2 through Defect Engineering Yan Chen, 1,2,,$, * Shengxi Huang, 3,6, Xiang Ji, 2 Kiran Adepalli, 2 Kedi Yin, 8 Xi Ling, 3,9 Xinwei
More informationIAP 2006: From nano to macro: Introduction to atomistic modeling techniques and application in a case study of modeling fracture of copper (1.
IAP 2006: From nano to macro: Introduction to atomistic modeling techniques and application in a case study of modeling fracture of copper (1.978 PDF) http://web.mit.edu/mbuehler/www/teaching/iap2006/intro.htm
More informationElastic constants and the effect of strain on monovacancy concentration in fcc hard-sphere crystals
PHYSICAL REVIEW B 70, 214113 (2004) Elastic constants and the effect of strain on monovacancy concentration in fcc hard-sphere crystals Sang Kyu Kwak and David A. Kofke Department of Chemical and Biological
More informationDepartment of Engineering Mechanics, SVL, Xi an Jiaotong University, Xi an
The statistical characteristics of static friction J. Wang, G. F. Wang*, and W. K. Yuan Department of Engineering Mechanics, SVL, Xi an Jiaotong University, Xi an 710049, China * E-mail: wanggf@mail.xjtu.edu.cn
More informationProperties of liquid nickel: A critical comparison of EAM and MEAM calculations
PHYSICAL REVIEW B, VOLUME 65, 0409 Properties of liquid nickel: A critical comparison of EAM and MEAM calculations F. J. Cherne and M. I. Baskes Los Alamos National Laboratory, MST-8 Structure and Property
More informationSupporting Online Material (1)
Supporting Online Material The density functional theory (DFT) calculations were carried out using the dacapo code (http://www.fysik.dtu.dk/campos), and the RPBE (1) generalized gradient correction (GGA)
More informationModeling aluminum using molecular dynamics simulation
Journal of Materials and Environmental Sciences ISSN : 2028-2508 CODEN : JMESCN Copyright 2017, University of Mohammed Premier Oujda Morocco J. Mater. Environ. Sci., 2018, Volume 9, Issue 1, Page 93-99
More informationIntroduction to molecular dynamics
1 Introduction to molecular dynamics Yves Lansac Université François Rabelais, Tours, France Visiting MSE, GIST for the summer Molecular Simulation 2 Molecular simulation is a computational experiment.
More informationComplementary approaches to high T- high p crystal structure stability and melting!
Complementary approaches to high T- high p crystal structure stability and melting! Dario ALFÈ Department of Earth Sciences & Department of Physics and Astronomy, Thomas Young Centre@UCL & London Centre
More informationThe Pennsylvania State University. The Graduate School. College of Earth and Mineral Sciences A FIRST-PRINCIPLES METHODOLOGY FOR
The Pennsylvania State University The Graduate School College of Earth and Mineral Sciences A FIRST-PRINCIPLES METHODOLOGY FOR DIFFUSION COEFFICIENTS IN METALS AND DILUTE ALLOYS A Dissertation in Materials
More informationSCIENCE CHINA Physics, Mechanics & Astronomy
SCIENCE CHINA Physics, Mechanics & Astronomy Article April 2012 Vol.55 No.4: 614 618 doi: 10.1007/s11433-012-4679-8 Stability and diffusion properties of self-interstitial atoms in tungsten: a first-principles
More information* Theoretische Physik II, Universitat Dortmund, Dortmund, Germany
JOURNAL DE PHYSIQUE IV Colloque C8, suppl6ment au Journal de Physique III, Volume 5, dkembre 1995 Structural Phase Transformation and Phonon Softening in Iron-Based Alloys H.C. Herper, E. Hoffmann, P.
More informationAccuracy and transferability of GAP models for tungsten
Accuracy and transferability of GAP models for tungsten Wojciech J. Szlachta Albert P. Bartók Gábor Csányi Engineering Laboratory University of Cambridge 5 November 214 Motivation Number of atoms 1 1 2
More informationAb initio phonon calculations in mixed systems
Ab initio phonon calculations in mixed systems Andrei Postnikov apostnik@uos.de Outline: Experiment vs. ab initio theory Ways of theory: linear response and frozen phonon approaches Applications: Be x
More informationDiffusion in multicomponent solids. Anton Van der Ven Department of Materials Science and Engineering University of Michigan Ann Arbor, MI
Diffusion in multicomponent solids nton Van der Ven Department of Materials Science and Engineering University of Michigan nn rbor, MI Coarse graining time Diffusion in a crystal Two levels of time coarse
More informationComputer simulation of symmetrical tilt grain boundaries in noble metals with MAEAM
Vol 16 No 1, January 2007 c 2007 Chin. Phys. Soc. 1009-1963/2007/16(01)/0210-07 Chinese Physics and IOP Publishing Ltd Computer simulation of symmetrical tilt grain boundaries in noble metals with MAEAM
More informationSupplementary Materials
Supplementary Materials Atomistic Origin of Brittle Failure of Boron Carbide from Large Scale Reactive Dynamics Simulations; Suggestions toward Improved Ductility Qi An and William A. Goddard III * Materials
More informationSTRUCTURAL AND MECHANICAL PROPERTIES OF AMORPHOUS SILICON: AB-INITIO AND CLASSICAL MOLECULAR DYNAMICS STUDY
STRUCTURAL AND MECHANICAL PROPERTIES OF AMORPHOUS SILICON: AB-INITIO AND CLASSICAL MOLECULAR DYNAMICS STUDY S. Hara, T. Kumagai, S. Izumi and S. Sakai Department of mechanical engineering, University of
More informationQuantum Monte Carlo simulation of spin-polarized tritium
Higher-order actions and their applications in many-body, few-body, classical problems Quantum Monte Carlo simulation of spin-polarized tritium I. Bešlić, L. Vranješ Markić, University of Split, Croatia
More informationDiffusion mechanisms in Cu grain boundaries
PHYSICAL REVIEW B VOLUME 62, NUMBER 6 Diffusion mechanisms in Cu grain boundaries 1 AUGUST 2000-II Mads R. So rensen Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545 Yuri
More informationFirst Principles Calculation of Defect and Magnetic Structures in FeCo
Materials Transactions, Vol. 47, No. 11 (26) pp. 2646 to 26 Special Issue on Advances in Computational Materials Science and Engineering IV #26 The Japan Institute of Metals First Principles Calculation
More informationIndex. A Ab-initio, 8 Accomodation, 7 Accuracy, 37 hardware, 37 software, 37 Ackland, 36, 55, 157 Andersen, 44, 47 Austenite, 3
Index A Ab-initio, 8 Accomodation, 7 Accuracy, 37 hardware, 37 software, 37 Ackland, 36, 55, 157 Andersen, 44, 47 Austenite, 3 B B2, 151 154, 156 Bain, 9, 10, 72, 73 parameters, 9 transformation, 9, 72
More informationMolecular Dynamics and Accelerated Molecular Dynamics
Molecular Dynamics and Accelerated Molecular Dynamics Arthur F. Voter Theoretical Division National Laboratory Lecture 3 Tutorial Lecture Series Institute for Pure and Applied Mathematics (IPAM) UCLA September
More informationRustam Z. Khaliullin University of Zürich
Rustam Z. Khaliullin University of Zürich Molecular dynamics (MD) MD is a computational method for simulating time evolution of a collection of interacting atoms by numerically integrating Newton s equation
More informationCharged point defects in semiconductors and the supercell approximation
INSTITUTE OF PHYSICS PUBLISHING JOURNAL OF PHYSICS: CONDENSED MATTER J. Phys.: Condens. Matter 14 (2002) 2637 2645 PII: S0953-8984(02)30267-4 Charged point defects in semiconductors and the supercell approximation
More informationTheory of Electromigration
Theory of Electromigration Electromigration is the transport of material in a conductor under the influence of an applied electric field. All conductors are susceptible to electromigration, therefore it
More informationCHAPTER 4 IMPERFECTIONS IN SOLIDS PROBLEM SOLUTIONS
4-1 CHAPTER 4 IMPERFECTIONS IN SOLIDS PROBLEM SOLUTIONS Vacancies and Self-Interstitials 4.1 In order to compute the fraction of atom sites that are vacant in copper at 1357 K, we must employ Equation
More informationSMARTMET project: Towards breaking the inverse ductility-strength relation
SMARTMET project: Towards breaking the inverse ductility-strength relation B. Grabowski, C. Tasan SMARTMET ERC advanced grant 3.8 Mio Euro for 5 years (Raabe/Neugebauer) Adaptive Structural Materials group
More informationDocument Version Publisher s PDF, also known as Version of Record (includes final page, issue and volume numbers)
Parametrization of modified embedded-atom-method potentials for Rh, Pd, Ir, and Pt based on density functional theory calculations, with applications to surface properties Beurden, van, P.; Kramer, G.J.
More informationIntroduction to Simulation - Lectures 17, 18. Molecular Dynamics. Nicolas Hadjiconstantinou
Introduction to Simulation - Lectures 17, 18 Molecular Dynamics Nicolas Hadjiconstantinou Molecular Dynamics Molecular dynamics is a technique for computing the equilibrium and non-equilibrium properties
More informationInteratomic potentials with error bars. Gábor Csányi Engineering Laboratory
Interatomic potentials with error bars Gábor Csányi Engineering Laboratory What makes a potential Ingredients Desirable properties Representation of atomic neighbourhood smoothness, faithfulness, continuity
More informationChromium Cluster on Defected Graphene
Chromium Cluster on Defected Graphene Yuhang Liu June 29, 2017 Abstract In this work, diffusion process of Cr atoms on two types of defected graphene and structure and magnetic properties of Cr cluster
More informationReactive potentials and applications
1.021, 3.021, 10.333, 22.00 Introduction to Modeling and Simulation Spring 2011 Part I Continuum and particle methods Reactive potentials and applications Lecture 8 Markus J. Buehler Laboratory for Atomistic
More informationMD Thermodynamics. Lecture 12 3/26/18. Harvard SEAS AP 275 Atomistic Modeling of Materials Boris Kozinsky
MD Thermodynamics Lecture 1 3/6/18 1 Molecular dynamics The force depends on positions only (not velocities) Total energy is conserved (micro canonical evolution) Newton s equations of motion (second order
More informationIntroduction to Molecular Dynamics
Introduction to Molecular Dynamics Dr. Kasra Momeni www.knanosys.com Overview of the MD Classical Dynamics Outline Basics and Terminology Pairwise interacting objects Interatomic potentials (short-range
More informationSupplementary Information
Supplementary Information Supplementary Figure 1: Electronic Kohn-Sham potential profile of a charged monolayer MoTe 2 calculated using PBE-DFT. Plotted is the averaged electronic Kohn- Sham potential
More informationAb initio Molecular Dynamics and Elastic Properties of TiC and TiN Nanoparticles
Mat. Res. Soc. Symp. Proc. Vol. 704 2002 Materials Research Society Ab initio Molecular Dynamics and Elastic Properties of TiC and TiN Nanoparticles A. V. Postnikov and P. Entel Theoretical Low-Temperature
More informationFirst-principles calculations and tight-binding molecular dynamics simulations of the palladium-hydrogen system
First-principles calculations and tight-binding molecular dynamics simulations of the palladium-hydrogen system A. Shabaev and D. A. Papaconstantopoulos George Mason University, Fairfax, Virginia 22030,
More informationSupplementary Information for: Controlling Cellular Uptake of Nanoparticles with ph-sensitive Polymers
Supplementary Information for: Controlling Cellular Uptake of Nanoparticles with ph-sensitive Polymers Hong-ming Ding 1 & Yu-qiang Ma 1,2, 1 National Laboratory of Solid State Microstructures and Department
More informationStructure and Curie temperature of Y 2 Fe 17 x Cr x
Vol. 46 No. 4 SCIENCE IN CHINA (Series G) August 2003 Structure and Curie temperature of Y 2 Fe 17 x Cr x HAO Shiqiang ( ) 1 & CHEN Nanxian ( ) 1,2 1. Department of Physics, Tsinghua University, Beijing
More informationInvestigations of the effects of 7 TeV proton beams on LHC collimator materials and other materials to be used in the LHC
Russian Research Center Kurchatov Institute Investigations of the effects of 7 ev proton beams on LHC collimator materials and other materials to be used in the LHC A.I.Ryazanov Aims of Investigations:
More informationAnalysis of the ultrafast dynamics of the silver trimer upon photodetachment
J. Phys. B: At. Mol. Opt. Phys. 29 (1996) L545 L549. Printed in the UK LETTER TO THE EDITOR Analysis of the ultrafast dynamics of the silver trimer upon photodetachment H O Jeschke, M E Garcia and K H
More informationMatSci 331 Homework 4 Molecular Dynamics and Monte Carlo: Stress, heat capacity, quantum nuclear effects, and simulated annealing
MatSci 331 Homework 4 Molecular Dynamics and Monte Carlo: Stress, heat capacity, quantum nuclear effects, and simulated annealing Due Thursday Feb. 21 at 5pm in Durand 110. Evan Reed In this homework,
More informationSTUDY ON CORROSION P HENOMENA OF STEELS IN PB-BI FLOW
11 th International Conference on Nuclear Engineering Tokyo, JAPAN, April -3, 3 ICONE11-36375 STUDY ON CORROSION P HENOMENA OF STEELS IN PB-BI FLOW Yingxia Qi Research Laboratory for Nuclear Reactors Tokyo
More informationOn angular momentum balance for particle systems with periodic boundary conditions
On angular momentum balance for particle systems with periodic boundary conditions Vitaly A. Kuzkin 1 Institute for Problems in Mechanical Engineering RAS Saint Petersburg State Polytechnical University
More informationAnharmonicity and Quantum Effect in Thermal Expansion of an Invar Alloy
Anharmonicity and Quantum Effect in Thermal Expansion of an Invar Alloy 1 Institute for Molecular Science (IMS), Okazaki The Graduate University for Advanced Studies (Sokendai) Toshihiko YOKOYAMA 1,, Keitaro
More informationA COMPARATIVE STUDY OF OXYGEN VACANCY MIGRATION PATHWAYS IN CRYSTALLINE POLYMORPHS OF SILICA
A COMPARATIVE STUDY OF OXYGEN VACANCY MIGRATION PATHWAYS IN CRYSTALLINE POLYMORPHS OF SILICA L. René Corrales Pacific Northwest National Laboratory Richland, WA 99352 Corresponding author: rene.corrales@pnl.gov
More informationMolecular-dynamics study of surface segregation in liquid semiconductor alloys
PHYSICAL REVIEW B VOLUME 56, NUMBER 19 15 NOVEMBER 1997-I Molecular-dynamics study of surface segregation in liquid semiconductor alloys Wenbin Yu and D. Stroud* Department of Physics, The Ohio State University,
More informationMolecular dynamics simulations of EXAFS in germanium
Cent. Eur. J. Phys. 93 2011 710-715 DOI: 10.2478/s11534-010-0074-0 Central European Journal of Physics Molecular dynamics simulations of EXAFS in germanium Research Article Janis Timoshenko Alexei Kuzmin
More informationTitle atoms in a graphite intercalation c. Citation PHYSICAL REVIEW B (2008), 78(21)
Title First-principles approach to chemic atoms in a graphite intercalation c Author(s) Toyoura, Kazuaki; Koyama, Yukinori; Oba, Fumiyasu; Tanaka, Isao Citation PHYSICAL REVIEW B (28), 78(21) Issue Date
More informationEmbedded atom approach for gold silicon system from ab initio molecular dynamics simulations using the force matching method
Bull. Mater. Sci., Vol. 9, No. 5, September 6, pp. 9 47. DOI.7/s4-6-65- c Indian Academy of Sciences. Embedded atom approach for gold silicon system from ab initio molecular dynamics simulations using
More informationAb Ini'o Molecular Dynamics (MD) Simula?ons
Ab Ini'o Molecular Dynamics (MD) Simula?ons Rick Remsing ICMS, CCDM, Temple University, Philadelphia, PA What are Molecular Dynamics (MD) Simulations? Technique to compute statistical and transport properties
More informationComputer simulation methods (2) Dr. Vania Calandrini
Computer simulation methods (2) Dr. Vania Calandrini in the previous lecture: time average versus ensemble average MC versus MD simulations equipartition theorem (=> computing T) virial theorem (=> computing
More informationSnO 2 Physical and Chemical Properties due to the Impurity Doping
, March 13-15, 2013, Hong Kong SnO 2 Physical and Chemical Properties due to the Impurity Doping Richard Rivera, Freddy Marcillo, Washington Chamba, Patricio Puchaicela, Arvids Stashans Abstract First-principles
More informationSecond-moment interatomic potential for aluminum derived from total-energy calculations and molecular dynamics application
Computational Materials Science 17 (2000) 224±229 www.elsevier.com/locate/commatsci Second-moment interatomic potential for aluminum derived from total-energy calculations and molecular dynamics application
More informationSurface stress and relaxation in metals
J. Phys.: Condens. Matter 12 (2000) 5541 5550. Printed in the UK PII: S0953-8984(00)11386-4 Surface stress and relaxation in metals P M Marcus, Xianghong Qian and Wolfgang Hübner IBM Research Center, Yorktown
More informationFirst-principles calculations of bulk and interfacial thermodynamic properties for fcc-based Al-Sc alloys
PHYSICAL REVIEW B VOLUME 57, NUMBER 18 1 MAY 1998-II First-principles calculations of bulk and interfacial thermodynamic properties for fcc-based Al-Sc alloys Mark Asta, S. M. Foiles, and A. A. Quong*
More informationA tight-binding molecular dynamics study of phonon anharmonic effects in diamond and graphite
J. Phys.: Condens. Matter 9 (1997) 7071 7080. Printed in the UK PII: S0953-8984(97)83513-8 A tight-binding molecular dynamics study of phonon anharmonic effects in diamond and graphite G Kopidakis, C Z
More informationTwo-Dimensional Spin-Polarized Hydrogen at Zero
J Low Temp Phys (2013) 171:685 692 DOI 10.1007/s10909-012-0756-7 Two-Dimensional Spin-Polarized Hydrogen at Zero Temperature L. Vranješ Markić J. Boronat Received: 13 July 2012 / Accepted: 16 September
More informationStudy of Ni/Al Interface Diffusion by Molecular Dynamics Simulation*
Engineering, 2011, 3, 227-232 doi:10.4236/eng.2011.33026 Published Online March 2011 (http://www.scirp.org/journal/eng) Study of Ni/Al Interface Diffusion by Molecular Dynamics Simulation* Chunguang Zhang,
More informationSupporting information. Realizing Two-Dimensional Magnetic Semiconductors with. Enhanced Curie Temperature by Antiaromatic Ring Based
Supporting information Realizing Two-Dimensional Magnetic Semiconductors with Enhanced Curie Temperature by Antiaromatic Ring Based Organometallic Frameworks Xingxing Li and Jinlong Yang* Department of
More informationarxiv:cond-mat/ v1 10 Jun 1994 K. M. Rabe
October 2, 2018 Phase transitions in BaTiO 3 from first principles W. Zhong and David Vanderbilt Department of Physics and Astronomy, Rutgers University, Piscataway, NJ 08855-0849 arxiv:cond-mat/9406049v1
More informationSTUDY OF INTEGRATION ALGORITHM AND TIME STEP ON MOLECULAR DYNAMIC SIMULATION. Janusz Bytnar, Anna Kucaba-Piętal
Business and Engineering Applications of Intelligent and Information Systems 11 STUDY OF INTEGRATION ALGORITHM AND TIME STEP ON MOLECULAR DYNAMIC SIMULATION Janusz Bytnar, Anna Kucaba-Piętal Abstract:
More informationStructural properties of low-density liquid alkali metals
PRAMANA c Indian Academy of Sciences Vol. 65, No. 6 journal of December 2005 physics pp. 1085 1096 Structural properties of low-density liquid alkali metals A AKANDE 1, G A ADEBAYO 1,2 and O AKINLADE 2
More informationCalculation of thermal expansion coefficient of Fe3Al with the addition of transition metal elements Abstract Introduction Methodology
Calculation of thermal expansion coefficient of Fe Al with the addition of transition metal elements Tatiana Seletskaia, Leonid Muratov, Bernard Cooper, West Virginia University, Physics Departement, Hodges
More informationMolecular Dynamics Simulation of Nanometric Machining Under Realistic Cutting Conditions Using LAMMPS
Molecular Dynamics Simulation of Nanometric Machining Under Realistic Cutting Conditions Using LAMMPS Rapeepan Promyoo Thesis Presentation Advisor: Dr. Hazim El-Mounayri Department of Mechanical Engineering
More informationSIMULATION OF NANO-SCALE CUTTING WITH MOLECULAR DYNAMICS
American Journal of Nanotechnology 5 (2): 17-26, 2014 ISSN 1949-0216 2014 Angelos P. Markopoulos et al., This open access article is distributed under a Creative Commons Attribution (CC-BY) 3.0 license
More information1 Adsorption of NO 2 on Pd(100) Juan M. Lorenzi, Sebastian Matera, and Karsten Reuter,
Supporting information: Synergistic inhibition of oxide formation in oxidation catalysis: A first-principles kinetic Monte Carlo study of NO+CO oxidation at Pd(100) Juan M. Lorenzi, Sebastian Matera, and
More informationJournal of Theoretical Physics
1 Journal of Theoretical Physics Founded and Edited by M. Apostol 53 (2000) ISSN 1453-4428 Ionization potential for metallic clusters L. C. Cune and M. Apostol Department of Theoretical Physics, Institute
More informationForce Fields in Molecular Mechanics
Force Fields in Molecular Mechanics Rajarshi Guha (9915607) and Rajesh Sardar (9915610) March 21, 2001 1 Introduction With the advent of computers chemists have realized the utility of carrying out simulations
More informationOne-shot free energy calculations for crystalline materials
Instability One-shot free energy calculations for crystalline materials Thesis for the Degree Master of Science TOMMY ANDERSSON Department of Applied Physics Chalmers University of Technology Gothenburg,
More informationDIFFUSION AND GROWTH ON STRAINED SURFACES
DIFFUSION AND GROWTH ON STRAINED SURFACES A.C. SCHINDLER AND D.E. WOLF Theoretische Physik, FB 10, Gerhard Mercator Universität Duisburg, 47048 Duisburg, Germany E-mail: A.Schindler@Uni-Duisburg.de Using
More informationMicroscopical and Microanalytical Methods (NANO3)
Microscopical and Microanalytical Methods (NANO3) 06.11.15 10:15-12:00 Introduction - SPM methods 13.11.15 10:15-12:00 STM 20.11.15 10:15-12:00 STS Erik Zupanič erik.zupanic@ijs.si stm.ijs.si 27.11.15
More informationCOMPARISON OF PROCESS OF DIFFUSION OF INTERSTITIAL OXYGEN ATOMS AND INTERSTITIAL HYDROGEN MOLECULES IN SILICON AND
COMPARISON OF PROCESS OF DIFFUSION OF INTERSTITIAL OXYGEN ATOMS AND INTERSTITIAL HYDROGEN MOLECULES IN SILICON AND GERMANIUM CRYSTALS: QUANTUMCHEMICAL SIMULATION Vasilii Gusakov Joint Institute of Solid
More informationThe effect of mutual angular misalignment in the quantized sliding of solid lubricants
Facoltà di Scienze e Tecnologie Laurea Triennale in Fisica The effect of mutual angular misalignment in the quantized sliding of solid lubricants Relatore: Prof. Nicola Manini Correlatore: Prof. Rosario
More informationStructure and energetics of Cu-Au alloys
Home Search Collections Journals About Contact us My IOPscience Structure and energetics of Cu-Au alloys This content has been downloaded from IOPscience. Please scroll down to see the full text. 2000
More informationElectronic Supplementary Information
Electronic Supplementary Material (ESI) for Physical Chemistry Chemical Physics. This journal is the Owner Societies 2014 Electronic Supplementary Information Computational investigation of structural
More informationLi = 1.6, rc Sn = 2.3, and rc S = 1.7 in bohr units. The
: Simulations of Its Structure and Electrolyte Properties N. A. W. Holzwarth Department of Physics, Wake Forest University, Winston-Salem, NC 2719, USA First principles simulations show that the ground
More informationLi 4 SnS 4 : Simulations of Its Structure and Electrolyte Properties
: Simulations of Its Structure and Electrolyte Properties N. A. W. Holzwarth Department of Physics, Wake Forest University, Winston-Salem, NC, USA Introduction Recently, there has been significant progress
More informationStructural properties of low density liquid alkali metals
PRAMANA c Indian Academy of Sciences Vol. xx, No. x journal of xxxxxx 2005 physics pp. 1 12 Structural properties of low density liquid alkali metals A AKANDE 1, G A ADEBAYO 1,2 and O AKINLADE 2 1 The
More informationUniform-acceptance force-biased Monte Carlo: A cheap way to boost MD
Dresden Talk, March 2012 Uniform-acceptance force-biased Monte Carlo: A cheap way to boost MD Barend Thijsse Department of Materials Science and Engineering Delft University of Technology, The Netherlands
More informationSolid-State Diffusion and NMR
Solid-State Diffusion and NMR P. Heitjans, S. Indris, M. Wilkening University of Hannover Germany Diffusion Fundamentals, Leipzig, 3 Sept. 005 Introduction Diffusivity in Solids as Compared to Liquids
More informationPhysics 211B : Problem Set #0
Physics 211B : Problem Set #0 These problems provide a cross section of the sort of exercises I would have assigned had I taught 211A. Please take a look at all the problems, and turn in problems 1, 4,
More informationModeling Materials. Continuum, Atomistic and Multiscale Techniques. gg CAMBRIDGE ^0 TADMOR ELLAD B. HHHHM. University of Minnesota, USA
HHHHM Modeling Materials Continuum, Atomistic and Multiscale Techniques ELLAD B. TADMOR University of Minnesota, USA RONALD E. MILLER Carleton University, Canada gg CAMBRIDGE ^0 UNIVERSITY PRESS Preface
More informationStructural Calculations phase stability, surfaces, interfaces etc
Structural Calculations phase stability, surfaces, interfaces etc Keith Refson STFC Rutherford Appleton Laboratory September 19, 2007 Phase Equilibrium 2 Energy-Volume curves..................................................................
More informationLattice Strain Due to an Atomic Vacancy
Int. J. Mol. Sci. 2009, 10, 2798-2808; doi:10.3390/ijms10062798 OPEN ACCESS Article International Journal of Molecular Sciences ISSN 1422-0067 www.mdpi.com/journal/ijms Lattice Strain Due to an Atomic
More information