Crystal Structure and Formation Energy of ε-carbide Using First Principles Calculations

Similar documents
DETERMINATION OF MECHANICAL PROPERTIES OF NANOSTRUCTURES WITH COMPLEX CRYSTAL LATTICE USING MOMENT INTERACTION AT MICROSCALE

Classical Mechanics. From Molecular to Con/nuum Physics I WS 11/12 Emiliano Ippoli/ October, 2011

C. Bulutay Topics on Semiconductor Physics. In This Lecture: Electronic Bandstructure: General Info

QUANTUM CHEMISTRY. Hückel Molecular orbital Theory Application PART I PAPER:2, PHYSICAL CHEMISTRY-I

fiziks Institute for NET/JRF, GATE, IIT JAM, M.Sc. Entrance, JEST, TIFR and GRE in Physics

Period #2 Notes: Electronic Structure of Atoms

- 5 - TEST 2. This test is on the final sections of this session's syllabus and. should be attempted by all students.

Mechanism of Roughness-induced CO 2 Microbubble Nucleation in Polypropylene Foaming

Breakdown of the Landau-Ginzburg-Wilson paradigm at quantum phase transitions

5.04 Principles of Inorganic Chemistry II

Lecture 5. Selected aspects of thermodynamics for adsorption, diffusion and desorption

PHY4605 Introduction to Quantum Mechanics II Spring 2005 Final exam SOLUTIONS April 22, 2005

Electron Theory of Permanent Magnets

Tunable and sizable band gap in silicene by surface adsorption

Electron Correlation Methods

Lecture 6. Thermodynamics and Kinetics for Adsorption, Diffusion and Desorption

Summary of equations chapters 7. To make current flow you have to push on the charges. For most materials:

THE INTERVAL LATTICE BOLTZMANN METHOD FOR TRANSIENT HEAT TRANSFER IN A SILICON THIN FILM

Massachusetts Institute of Technology Quantum Mechanics I (8.04) Spring 2005 Solutions to Problem Set 6

Energy Bands Energy Bands and Band Gap. Phys463.nb Phenomenon

A - INTRODUCTION AND OVERVIEW

Problems for HW X. C. Gwinn. November 30, 2009

MAGNETIC AND ELECTRONIC PROPERTIES OF SIMPLE, TRANSITION, AND RARE EARTH LIQUID METALS

Department of Electrical and Computer Engineering, Cornell University. ECE 4070: Physics of Semiconductors and Nanostructures.

This final is a three hour open book, open notes exam. Do all four problems.

Vector potential quantization and the photon wave-particle representation

Li ion diffusion mechanism in the crystalline electrolyte γ-li 3 PO 4

Point Lattices: Bravais Lattices

Effect of Pressure on the Structural and Electronic Properties of ZnS

Physics Graduate Prelim exam

Chapter 3 The Schrödinger Equation and a Particle in a Box

Vibrational Relaxation of HF (v=3) + CO

QUB XRD Course. The crystalline state. The Crystalline State

Physics 202, Lecture 14

EFFECTS OF ALLOYING ELEMENTS ON ELASTIC PROPERTIES OF Al BY FIRST-PRINCIPLES CALCULATIONS

The solutions of the single electron Hamiltonian were shown to be Bloch wave of the form: ( ) ( ) ikr

Analogy Between Particle in a Box and Jahn Teller Effect

ADVANCEMENT OF THE CLOSELY COUPLED PROBES POTENTIAL DROP TECHNIQUE FOR NDE OF SURFACE CRACKS

Brazilian Journal of Physics ISSN: Sociedade Brasileira de Física Brasil

Data Provided: A formula sheet and table of physical constants is attached to this paper. SOLID STATE PHYSICS

2.57/2.570 Midterm Exam No. 1 March 31, :00 am -12:30 pm

Chem 3502/4502 Physical Chemistry II (Quantum Mechanics) 3 Credits Spring Semester 2006 Christopher J. Cramer. Lecture 25, March 29, 2006

Week 10: Line Integrals

Metal Organic Frameworks: Forming Nano-Compounds and Photonic Interfaces

DISTRIBUTION OF SUB AND SUPER HARMONIC SOLUTION OF MATHIEU EQUATION WITHIN STABLE ZONES

Nuclear Reactions. Shape, interaction, and excitation structures of nuclei scattering expt. cf. Experiment by Rutherford (α scatt.

Quantum Theory Project. Tutorial 45 th Sanibel Symposium Plenary Sessions VII XVII (Tuesday Friday) So Hirata Samuel B. Trickey March 7, 2005

Structural Effect of Thioureas on the Detection of Chemical Warfare Agent Simulants

Light and Optics Propagation of light Electromagnetic waves (light) in vacuum and matter Reflection and refraction of light Huygens principle

a * a (2,1) 1,1 0,1 1,1 2,1 hkl 1,0 1,0 2,0 O 2,1 0,1 1,1 0,2 1,2 2,2

An Approximate DFT Method: The Density-Functional Tight-Binding (DFTB) Method

Physical Chemistry (II) CHEM Lecture 25. Lecturer: Hanning Chen, Ph.D. 05/01/2018

SUPPLEMENTARY INFORMATION

Topological Quantum Compiling

Some parameters of varicaps with gradient base area based on Shottky barrier

The Predom module. Predom calculates and plots isothermal 1-, 2- and 3-metal predominance area diagrams. Predom accesses only compound databases.

Appendix A Light Absorption, Dispersion and Polarization

Quantum Mechanics Qualifying Exam - August 2016 Notes and Instructions

Lecture: P1_Wk2_L5 Contact Mechanics. Ron Reifenberger Birck Nanotechnology Center Purdue University 2012

Name Solutions to Test 3 November 8, 2017

MAC-solutions of the nonexistent solutions of mathematical physics

Influence of Carbon Vacancies on CO Chemisorption on TiC(001): A Theoretical Study

Quantum Analogs Chapter 4 Student Manual

A Brief Note on Quasi Static Thermal Stresses In A Thin Rectangular Plate With Internal Heat Generation

Jack Simons, Henry Eyring Scientist and Professor Chemistry Department University of Utah

The heat budget of the atmosphere and the greenhouse effect

Ch. 24 Molecular Reaction Dynamics 1. Collision Theory 2. Diffusion-Controlled Reaction

Topological Insulators in 2D and 3D

Practice Problems Solution

Data Provided: A formula sheet and table of physical constants is attached to this paper. SOLID STATE PHYSICS

LECTURE NOTES ON PATH INTEGRALS

Data Provided: A formula sheet and table of physical constants is attached to this paper. SOLID STATE PHYSICS

Nuclear Reactions. Shape, interaction, and excitation structures of nuclei scattering expt. cf. Experiment by Rutherford (a scatt.

Chem 130 Second Exam

Physics 24 Exam 1 February 18, 2014

Crystalline Structures The Basics

Crystals. Fig From Principles of Electronic Materials and Devices, Third Edition, S.O. Kasap ( McGraw-Hill, 2005)

Supporting Information

Variational and other methods

Bypassing no-go theorems for consistent interactions in gauge theories

Physics 201 Lab 3: Measurement of Earth s local gravitational field I Data Acquisition and Preliminary Analysis Dr. Timothy C. Black Summer I, 2018

Overview. Before beginning this module, you should be able to: After completing this module, you should be able to:

13.4 Work done by Constant Forces

Numerical Study Of Coated Electrical Contacts

25 Which of the following summarises the change in wave characteristics on going from infra-red to ultraviolet in the electromagnetic spectrum?

General Model of Diffusion of Interstitial Oxygen in Silicon and Germanium Crystals

SOLUTION OF QUADRATIC NONLINEAR PROBLEMS WITH MULTIPLE SCALES LINDSTEDT-POINCARE METHOD. Mehmet Pakdemirli and Gözde Sarı

University of Washington Department of Chemistry Chemistry 453 Winter Quarter 2010 Homework Assignment 4; Due at 5p.m. on 2/01/10

Strategy: Use the Gibbs phase rule (Equation 5.3). How many components are present?

INTRODUCTION. The three general approaches to the solution of kinetics problems are:

Photonic irradiation induces nano-structures an EPR study

On the Linear Stability of Compound Capillary Jets

F k B T = 1 2 K. θ q e iq x, (2) q 2 θq. d 2 1

4 The dynamical FRW universe

Physics 1402: Lecture 7 Today s Agenda

Candidates must show on each answer book the type of calculator used.

Phys 6321 Final Exam - Solutions May 3, 2013

Chapter 9 Many Electron Atoms

University of Washington Department of Chemistry Chemistry 453 Winter Quarter 2009

PC-SAFT: Theory and Application

Transcription:

2010 Crystl Structure nd Formtion Energy of ε-crbide Using First Principles Clcultions Je Hoon Jng, In Gee Kim, Dong Woo Suh nd H. K. D. H. hdeshi Computtionl Metllurgy Lbortory Grdute Institute of Ferrous Technology Pohng University of Science nd Technology

Introduction Mrtensite (α ) ε-crbide η-crbide χ-crbide Cementite (θ) Fe 2.4 C Fe 2 C Fe 2.5 C Fe 3 C Silicon promotes the formtion of ε-crbide below 520 K. S. S. Nyk et.l, Mterils Science nd Engineering. 498, pp.442-456(2008) 900, 100s 250, 30 s 1.0wt% Si : No ε-crbide 1.7wt% Si : ε-crbide M s = 302(1.0 wt%si), 293 (1.7 wt%si) 200, 20s ε-crbide forms without redistribution of Si. S. J. rnrd, G. D. W. Smith, Proceedings of the solid-solid phse trnsformtion., pp.881(1981) No initil prtitioning of Si between ε, θ nd mrtensite S. S. bu, K. Hono, T. Skuri, Metl Mter. Trns. 25 (1994) p. 499

FLPW method E. Wimmer, H. Krkuer, M. Weinert, nd. J. Freemn, Phys. Rev. 28, 864 (1981) nd references therein. M. Weinert, E. Wimmer, nd. J. Freemn, Phys. Rev. 26, 4571 (1982). Muffin-Tin (MT) Sphere Region Interstitil Region $ v G ( r) = c # G( K, r, Wve Function Expnsion! k v k + G " Kmx k, ) ) G( k, r) i( k+ G) ( r "% e = $ v v [ ul r ul "# ' r Y lm lm ( ) + & lm ( )]& lm, r! Interstitil (*,), r! MT sphere v

Clcultion Prmeters E. Wimmer, H. Krkuer, M. Weinert, nd. J. Freemn, Phys. Rev. 28, 864 (1981) nd references therein. M. Weinert, E. Wimmer, nd. J. Freemn, Phys. Rev. 26, 4571 (1982). Perdew, J. P., urke, K., Ernzerhof, M., Phys. Rev. Let. 77,3865 (1996) ll-electron Full-potentil LPW method Generlized Grdient pproximtion for Exchnge-Correltion Potentil Plne-Wve Cutoff : 21 Ry Str-Function Cutoff : 340 Ry k-points : 88 Fe 2.4 C, 365 (Fe 11 M)C 5 Muffin-tin Sphere : Fe, Si, l, Mn (2.04.u.), C (1.30.u.) Mixing Method : royden

Epsilon Crbide C 2 C 3 h C 2 C 3 Formul Unit Structure Spce group c c/ Fe 2.4 C(Fe 2 C~Fe 3 C) hexgonl P6 3 22 or P6 3 /mmc 4.767 Å 4.354 Å 0.913 C 3 C 2 S. Nkgur, J. Phys. Soc. Jpn, 14 (1959) 186.

Fe 3 C nd Fe 2 C (ε-crbide( -crbide) 4.661 Å ( 2.3%) 4.785 Å (+0.3%) c 4.294 Å ( 1.4%) c 4.321 Å ( 0.8%) c/ 0.9213 c/ 0.903 x 0.3190 x 0.3302 # E = E( Fe C2) " 6! E(Fe) " 2! E(C) 8 6 = 5.09kJ/mol # E = E Fe C ) " 6! E(Fe) " 3! E(C) 9 ( 6 3 = 7.00kJ/mol

Fe 2.4 C (ε-crbide( -crbide) Fe 2.4 c 4.740 Å ( 0.6%) 8.631 Å ( 0.9%) # E = E( Fe C5 ) " 12! E(Fe) " 5! E(C) 17 12 = 6.24kJ/mol

Results 0.0-0.5-1.0-1.5-2.0-2.5 Lttice Chnge(%) 0.5 Fe3C Fe2.4C Fe2C c System (Å) c (Å) c h / h Mesured, ε 4.767 8.708 1.582 Fe 3 C 4.661( 2.3%) 8.588 ( 1.4%) 1.596 Fe 2.4 C 4.740( 0.6%) 8.631( 0.9%) 1.577 Fe 2 C 4.785(+0.3%) 8.642( 0.8%) 1.564

Si,, l nd Mn Substitution Si l Mn 4.7303Å ( 0.2%) 4.742Å (+0.0%) 4.738 Å ( 0.1%) c 8.5901 Å ( 0.5%) c 8.685 Å (+0.6%) c 8.664 Å (+0.4%)!E = 9.08kJ/mol!E = 4.98kJ/mol!E = 4.40kJ/mol

Results Lttice Chnge(%) 0.6 0.4 0.2 0.0-0.2-0.4 pure Si l Mn c System (Å) c (Å) c h / h Mesured, ε 4.767 8.708 1.582 Clculted, ε 4.740 8.631 1.577 Si substituted 4.730( 0.2%) 8.590( 0.5%) 1.573 l substituted 4.742(+0.0%) 8.685(+0.6%) 1.586 Mn substituted 4.738( 0.1%) 8.664(+0.4%) 1.584

Results Formtion Energy (kj/mol) 10 8 6 4 2 0 epsilon cementite pure Si l Mn E (kj/mol) pure-crbide Si substituted l substituted Mn substituted ε-crbide 6.24 9.08(+2.84) 4.98( 1.26) 4.40( 1.84) cementite 5.38 7.70(+2.32) 4.53( 0.85) 5.07( 0.31)

Orienttion Reltionship H. K. D. H. hdeshi, inite (1992) (101) α (1011) ε (211) α (1010) ε [011] α [0001] ε (111) α (1210) ε C 2 C 3 C 2 C 3 C 3 C 2 C [1210] 1 ε h 4. 7% Mrtensite with 1 wt% C = 2.85 Å c = 2.98 Å [111] α (0001) ε 9.2% (011) α

Summry First Principles Clcultion cn be pplied for hypotheticl crystl structure. Si ddition increses the formtion energy of θ nd ε- crbide. The formtion energy clcultion : ε-crbide θ The role of silicon in trnsition of crbide : Reducing the misfit Mngnese ddition : stble ε-crbide

Thnk You!!

Equilibrium Clcultion µ pr with silicon µ equi without silicon Crbon θ, ε-crbide

ppendix - DFT Hohenberg-Kohn Theorem : The ground stte property is functionl of electron density. E [ n] = " V ( r)! n( r) dr + < # T + Vee # >! E! n [ n] ext = V [ n]! EGS, E[ ngs] EGS E = Kohn-Shm Eqution : Introducing the non-intercting fictitious prticle. E 1 = " ext C 2 " r [ n] V ( r)! n( r) dr + T[ n] + V ( r)! n( r) d + Exc[ n] ext ( r)! T +! n [ n] + V C ( r)! Exc +! n [ n] = µ! E! n [ n]! T[ n] = + v eff r! n ( ) = µ ( 1 2 % & * ) + v #! = " eff i i! i ' 2 $ ( r) ( r) ( r) v! ( ) V ( ) VC( ) xc eff r = ext r + r + n N ( r ) = " ( r)! i= 1 i 2 E! n [ n]

First-Principles Clcultion Quntum Mechnics tomic Position Schrödinger Eqution Electron Density Physicl Principles Structure Moleculr Crystl Density Defect Surfce Interfce Thermodynmic Internl Energy Entropy Free Energy Formtion Energy Phse Trnsform Phse digrms Mechnicl - Compressibility Elstic constnts ulk modulus Vibrtions Etc. Electron Ioniztion Energy bnd nd gps Mgnetic

Formtion Energy Formtion energy (kj/mol) 300 250 200 150 100 50 0-50 G xs =ΩX(1 X) (Fe,Si)3C (Fe,l)3C (Fe,Mn)3C 0 10 20 30 40 50 60 70 80 Concentrtion (t %)

Ternry Phse Digrm t 773K Equilibrium Mn Si l 0.2 0.2 0.2 θ+γ θ+γ θ+γ Fe 0.8 Fe 0.8 Fe 0.8 Pr-equilibrium 0.2 Mn θ+γ 0.2 Si γ θ+γ 0.2 l θ+γ Fe 0.8 Fe 0.8 Fe 0.8

Equilibrium Phse Digrm t 723K γ-θ-ε γ-θ-ε γ γ2+ ε γ2+ ε γ+ ε γ+ ε γ γ+ θ+ ε γ+ θ θ+ ε γ+γ2+ ε γ+ θ γ+ θ+ ε

Gibbs Free Energy of Crbide H(ε) H(θ) G(ε) G(θ) G(θ) G(ε) H(θ) H(ε) Fe-C System Mn-C System

Gibbs Free Energy of l-crbide H(ε) H(θ) G(θ) G(ε) H(θ) H(ε) G(ε) G(θ) Si-C System l-c System

Density Functionl Theory full-reltivistic sclr-reltivistic non-reltivisitic 2 [# 1?? + V ( r) + V ( r) + V ( r)]ø ( r) = Ø ( r ) 2 non-periodic periodic symmetrized rel-spce ll electron, full-potentil ll electron, sphericl potentil pseudopotentil (vlence electron s) jellium pproximtion ext C xc! " i! locl density pprox. (LD) generlized grident pprox. (GG) non-spinpolrized spinpolrized, vector-spin density LD+U, OEP hybrid functionls current functionls rel-spce grid plne-wves (PW) non-liner methods - PW: ugmented method - KKR-GF linerized methods - LPW - SW - LMTO liner combintion of t omic orbitls (LCO) - tight-binding - Guss-O - Slter type-o - numericl O lugel, S. ihlmyer, G., Computtionl nnoscience: Do it yourself, 31:85-129 (2006)