Calculation of Phase Diagrams and First-Principles Study of Germanium Impacts on Phosphorus Distribution in Czochralski Silicon

Loading...
Thumbnail Image

Access rights

openAccess

URL

Journal Title

Journal ISSN

Volume Title

A1 Alkuperäisartikkeli tieteellisessä aikakauslehdessä

Date

2021-08

Major/Subject

Mcode

Degree programme

Language

en

Pages

17
4272–4288

Series

Journal of Electronic Materials, Volume 50, issue 8

Abstract

Phosphorus (P) is one of the most widely used donor dopants for fabricating a low-resistivity silicon (Si) substrate. However, its volatile nature and the relatively small equilibrium segregation coefficient in Si at the melting temperature of Si impede the efficient and effective growth of low-resistivity Czochralski (CZ) Si single crystal. The primary objective of this work is to theoretically perceive the influence of germanium co-doping on the heavily P-doped Si crystal by means of CALculation of PHase Diagrams (CALPHAD) approaches and density functional theory (DFT) calculations. Phase equilibria at the Si-rich corner of the Si-Ge-P system has been thermodynamically extrapolated based on robust thermodynamic descriptions of involved binary systems, where Si-P and Ge-P have been re-assessed in this work. Phase diagram calculation results indicate that at a given P concentration (e.g. 0.33 at.% P) Ge co-doping lowers the solidification temperature of the Si(Ge, P) alloys, as well as the relevant equilibrium segregation coefficients of P in the doped Si. DFT calculations simulated the formation of (i) monovacancy in Si as well as (ii) solutions of Si(P) and Si(Ge) with one dopant substitutionally inserted in 64- and 216-atom Si cubic supercells. Binding energies were calculated and compared for Ge-Ge, Ge-P and P-P bonds positioning at the first nearest-neighbors (1NN) to the third nearest-neighbors (3NN). P-P bonds have the largest bonding energy from 1NN to 3NN configurations. The climbing image nudged elastic band method (CL-NEB) was utilized to calculate the energy barriers of P 1NN jump in the 64-atom Si cubic supercell with/without a neighboring Ge atom. With Ge present, a higher energy barrier for P 1NN jump was obtained than that without involving Ge. This indicates that Ge can impede the P diffusion in Si matrix.

Description

Funding Information: This work is part of the POSITION-II project funded by the ECSEL Joint Undertaking under grant number Ecsel-783132-Position-II-2017-IA. www.position-2.eu . The authors would like to acknowledge the Innovation Funding Agency Business Finland for their financial support. Dr. Jari Paloheimo and Dr. Heikki Holmberg from Okmetic Oy, Finland, are acknowledged for useful discussions. One of the authors, Dr. Xiaoma Tao, acknowledges the support from National Natural Science Foundation of China [Grants Nos. 51661003 and 51961007], and the Guangxi Natural Science Foundation [Grant No. 2018GXNSFAA281254]. Publisher Copyright: © 2021, The Author(s). Copyright: Copyright 2021 Elsevier B.V., All rights reserved. | openaire: EC/H2020/783132/EU//POSITION-II

Keywords

CALPHAD, density functional theory, equilibrium segregation coefficient, formation energy, Germanium co-doping, heavily P-doped Si

Other note

Citation

Dong , H , Tao , X & Paulasto-Kröckel , M 2021 , ' Calculation of Phase Diagrams and First-Principles Study of Germanium Impacts on Phosphorus Distribution in Czochralski Silicon ' , Journal of Electronic Materials , vol. 50 , no. 8 , pp. 4272–4288 . https://doi.org/10.1007/s11664-021-08861-4