Development of the ReaxFF Reactive Force Field for Cu/Si Systems with Application to Copper Cluster Formation during Cu Diffusion Inside Silicon
Abstract
Transition-metal impurities such as nickel, copper, and iron in solid-state materials like silicon have a significant impact on the electrical performance of integrated circuits and solar cells. To study the impact of copper impurities inside bulk silicon on the electrical properties of the material, one needs to understand the configurational space of copper atoms incorporated inside the silicon lattice. In this work, we developed a ReaxFF reactive force field and used it to perform molecular dynamics simulations on models with up to 762 atoms to study the various configurations of individual and crystalline clusters of copper atoms inside bulk silicon by examining copper's diffusional behavior in silicon. The ReaxFF Cu/Si parameter set was developed by training against density functional theory (DFT) data, including the energy barrier for an individual Cu atom traveling inside a silicon lattice. We found that the diffusion of copper atoms is dependent on temperature. Moreover, we show that individual copper atoms start to form clusters inside bulk silicon at temperatures above 500 K. Our simulation results provide a comprehensive understanding of the effects of temperature and copper concentration on the formation of copper clusters inside a silicon lattice. Finally, the stress–strain relationship of Cu/Si compounds under uniaxial tensile loading has been obtained. Our results indicate a decrease in the elastic modulus with increasing Cu-impurity concentration. We observe spontaneous microcracking of the Si during the stress–strain tests as a consequence of the formation of a small Cu cluster adjacent to the Si surface.
Additional Information
© 2021 American Chemical Society. Received 11 May 2021. Revised1 3 August 2021. Published online 30 August 2021. This work was supported by NSF DMR grant #DMR-1842952. KR, MT and ACTvD acknowledge the support by a grant from the U.S. Army Research Laboratory through the Collaborative Research Alliance (CRA) for Multi-Scale Multidisciplinary Modeling of Electronic Materials (MSME) under Cooperative Agreement No. W911NF- 12-2-0023. Computations for this research were performed on the Pennsylvania State University's Institute for Cyber ScienceAdvanced Cyber Infrastructure (ICS-ACI). W.A.G. and J.E.M. were supported by NSF (CBET-1805022). The authors declare no competing financial interest.Attached Files
Supplemental Material - jp1c04178_si_001.pdf
Files
Name | Size | Download all |
---|---|---|
md5:0a9186293501bb9d86b46103b8a421e0
|
63.7 kB | Preview Download |
Additional details
- Eprint ID
- 110639
- DOI
- 10.1021/acs.jpcc.1c04178
- Resolver ID
- CaltechAUTHORS:20210831-184319829
- NSF
- DMR-1842952
- Army Research Office (ARO)
- W911NF-12-2-0023
- NSF
- CBET-1805022
- Created
-
2021-08-31Created from EPrint's datestamp field
- Updated
-
2021-10-19Created from EPrint's last_modified field
- Other Numbering System Name
- WAG
- Other Numbering System Identifier
- 1486