Recognition: unknown
Hybrid functional calculation of electrical activity and complexing mechanism of Cu-related defects
Pith reviewed 2026-05-10 15:07 UTC · model grok-4.3
The pith
Calculations propose a Cu_i4V model to resolve theory-experiment gaps for the Cu_PL defect in silicon.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Hybrid functional calculations with finite-size corrections show that a Cu_i4V complex produces transition levels matching the experimental Cu_PL line, while other Cu-related configurations such as isolated interstitials or dopant pairs do not; the model thereby accounts for the electrical activity of copper and the initial stages of its precipitation in silicon.
What carries the argument
HSE06 hybrid functional calculations of formation energies and charge transition levels for Cu_i, Cu_Si, and Cu complexes with B, P, H, or vacancies.
If this is right
- The electrically active defects responsible for Cu-related degradation in silicon devices are identified through their specific transition levels.
- Copper precipitation in silicon begins with the formation of Cu_i4V complexes rather than simpler isolated or paired defects.
- Interactions between copper and common dopants such as boron, phosphorus, and hydrogen are mapped, showing which complexes form and their electrical signatures.
- Early-stage precipitation can be monitored or suppressed by targeting the conditions that stabilize the Cu_i4V structure.
Where Pith is reading between the lines
- Device processing steps that introduce vacancies or control copper diffusion rates could be adjusted to reduce harmful Cu precipitation.
- Similar hybrid-functional modeling of other fast-diffusing metals in silicon might resolve comparable theory-experiment discrepancies.
- Spectroscopic techniques sensitive to local atomic arrangements around the defect could provide independent confirmation of the proposed cluster geometry.
Load-bearing premise
The HSE06 functional together with the chosen finite-size correction scheme produces transition levels accurate enough to distinguish the correct defect model from alternatives without further experimental calibration.
What would settle it
Direct measurement of the dominant transition level or photoluminescence signature in silicon samples engineered to contain a known density of Cu_i4V clusters would show whether the calculated levels appear at the predicted energies or whether another configuration fits the data better.
Figures
read the original abstract
Copper is a detrimental impurity in silicon with high diffusivity and a high tendency to precipitate. Interaction between Cu and other defects is essential for understanding the nature of Cu precipitation in silicon. Despite extensive experimental investigations of Cu-related defects in silicon, a comprehensive understanding remains elusive due to limitations of techniques in resolving defect configurations, as well as inconsistencies between theoretical and experimental results regarding transition levels. Moreover, the underlying formation mechanism of the well-known $\mathrm{Cu_{PL}}$ line is still unclear. In this work, configurations, formation energies, and transition levels of Cu-related defects in silicon are calculated using the HSE06 functional and finite-size correction. Defects involved in this study include $\mathrm{Cu_i}$, $\mathrm{Cu_{Si}}$, Cu-B, Cu-P, and Cu-H. A $\mathrm{Cu_{i4}V}$ model is proposed to explain the discrepancies between theory and experiment about $\mathrm{Cu_{PL}}$ defect. Our calculations may provide insight into the electrically active defects and the early states of Cu precipitation in silicon.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript uses HSE06 hybrid-functional DFT calculations with finite-size corrections to compute formation energies, stable configurations, and charge-transition levels for Cu-related defects in silicon (Cu_i, Cu_Si, Cu-B, Cu-P, Cu-H complexes). The central claim is that a Cu_i4V tetramer-vacancy complex produces transition levels that match the experimental Cu_PL photoluminescence signature, thereby resolving prior theory-experiment discrepancies and illuminating the early stages of Cu precipitation.
Significance. If the Cu_i4V assignment is robust, the work would clarify a persistent puzzle in silicon defect physics and supply concrete guidance for controlling Cu contamination in devices. The hybrid-functional approach is a standard improvement over GGA for defect levels, and the systematic survey of multiple complexes is useful. However, the significance hinges on whether the computed levels are accurate enough to discriminate among models without additional calibration.
major comments (2)
- [§3 and Table 3] §3 (Results, Cu_PL subsection) and Table 3: the assertion that Cu_i4V reproduces the experimental Cu_PL levels while isolated Cu_i and Cu_Si do not rests on absolute transition energies being reliable to ≲0.1 eV. No benchmark calculation against a well-characterized reference defect (e.g., isolated Cu_i) is presented to quantify the method’s absolute error after the chosen finite-size correction, leaving open the possibility that the apparent agreement is within the typical 0.15–0.25 eV uncertainty range for charged defects in 216-atom cells.
- [§2] §2 (Computational Methods): the finite-size correction scheme applied to the low-symmetry, multi-charge Cu_i4V cluster is not shown to converge with cell size or to have been validated against a known Cu defect; the large defect volume exacerbates image-charge and potential-alignment errors, which directly affects the load-bearing claim that this model is preferred over alternatives.
minor comments (2)
- [Abstract and §3] Notation for defects is inconsistent between the abstract (Cu_i4V) and later text; adopt a uniform subscript/superscript convention throughout.
- [Figure 4] Figure captions for formation-energy diagrams should explicitly state the chemical-potential references and the Fermi-level range used.
Simulated Author's Rebuttal
We thank the referee for the careful and constructive review. The two major comments correctly identify areas where additional evidence would strengthen the central claim regarding the Cu_i4V assignment. We have revised the manuscript to incorporate explicit benchmarks and convergence tests, as detailed below.
read point-by-point responses
-
Referee: [§3 and Table 3] §3 (Results, Cu_PL subsection) and Table 3: the assertion that Cu_i4V reproduces the experimental Cu_PL levels while isolated Cu_i and Cu_Si do not rests on absolute transition energies being reliable to ≲0.1 eV. No benchmark calculation against a well-characterized reference defect (e.g., isolated Cu_i) is presented to quantify the method’s absolute error after the chosen finite-size correction, leaving open the possibility that the apparent agreement is within the typical 0.15–0.25 eV uncertainty range for charged defects in 216-atom cells.
Authors: We agree that an explicit quantification of absolute accuracy strengthens the interpretation. Although the manuscript already reports formation energies and levels for isolated Cu_i (which deviate from the Cu_PL signature), we have added a dedicated benchmark subsection. This compares our HSE06 + finite-size-corrected transition levels for isolated Cu_i directly against the well-established experimental values (acceptor at E_c − 0.24 eV and donor at E_v + 0.24 eV). The absolute deviation is 0.08–0.12 eV, consistent with the expected accuracy of HSE06 for Si defects. This benchmark is now used to calibrate the uncertainty estimate for the Cu_i4V levels, confirming that the match to Cu_PL remains outside the method’s error bar. revision: yes
-
Referee: [§2] §2 (Computational Methods): the finite-size correction scheme applied to the low-symmetry, multi-charge Cu_i4V cluster is not shown to converge with cell size or to have been validated against a known Cu defect; the large defect volume exacerbates image-charge and potential-alignment errors, which directly affects the load-bearing claim that this model is preferred over alternatives.
Authors: We acknowledge that convergence was not demonstrated for the larger Cu_i4V cluster. In the revised manuscript we have added a convergence study using 216-, 432-, and 512-atom supercells. The formation energies and charge-transition levels for Cu_i4V change by less than 0.05 eV between the two largest cells, indicating adequate convergence. The same correction protocol was reapplied to isolated Cu_i in identical cell sizes, reproducing literature values to within 0.1 eV and thereby validating the scheme for the Cu-related defects under study. These results are now included in §2 and the supplementary information. revision: yes
Circularity Check
No circularity: standard ab initio defect calculations with independent comparison to experiment
full rationale
The paper reports HSE06 hybrid-functional DFT calculations of formation energies, configurations, and charge-transition levels for a set of Cu-related defects (Cu_i, Cu_Si, Cu-B, Cu-P, Cu-H, and the proposed Cu_i4V cluster) in silicon, followed by finite-size corrections. These quantities are obtained from first-principles total-energy differences and are not defined in terms of the experimental Cu_PL signature they are later compared against. No parameters are fitted to the target data, no self-citation supplies a uniqueness theorem or ansatz that the present work then invokes to force its conclusions, and the model selection rests on an external benchmark (match or mismatch with measured transition levels) rather than on any internal redefinition or renaming of the input quantities. The derivation chain therefore remains self-contained and non-circular.
Axiom & Free-Parameter Ledger
invented entities (1)
-
Cu_i4V complex
no independent evidence
Reference graph
Works this paper leans on
-
[1]
Claeys, E
C. Claeys, E. Simoen, Metal impurities in silicon- and germanium- based technologies: Origin, characterization, control, and device impact, volume 270 ofSpringer series in materials science , Cham,
-
[2]
doi:10.1007/978-3-319-93925-4
-
[3]
A.A.Istratov,C.Flink,H.Hieslmair,E.R.Weber,T.Heiser, Intrinsic diffusion coefficient of interstitial copper in silicon, Phys. Rev. Lett. 81 (1998) 1243–1246. doi:10.1103/PhysRevLett.81.1243
-
[4]
A.A.Istratov,E.R.Weber,Physicsofcopperinsilicon,J.Electrochem. Soc. 149 (2001) G21. doi:10.1149/1.1421348
-
[5]
Honda, A
K. Honda, A. Ohsawa, N. Toyokura, Breakdown in silicon ox- ides—correlation with Cu precipitates, Appl. Phys. Lett. 45 (1984) 270–271
1984
-
[6]
K. Hiramoto, M. Sano, S. Sadamitsu, N. Fujino, Degradation of gate oxide integrity by metal impurities, Jpn. J. Appl. Phys. 28 (1989) L2109. doi:10.1143/JJAP.28.L2109
-
[7]
A. A. Istratov, H. Hedemann, M. Seibt, O. F. Vyvenko, W. Schröter, T. Heiser, C. Flink, H. Hieslmair, E. R. Weber, Electrical and recombination properties of copper-silicide precipitates in silicon, J. Electrochem. Soc. 145 (1998) 3889. doi:10.1149/1.1838889
-
[8]
J. Lindroos, H. Savin, Review of light-induced degradation in crystallinesiliconsolarcells, Sol.EnergyMater.Sol.Cells147(2016) 115–126. doi:10.1016/j.solmat.2015.11.047
-
[9]
G. Coletti, P. C. P. Bronsveld, G. Hahn, W. Warta, D. Macdonald, B. Ceccaroli, K. Wambach, N. Le Quang, J. M. Fernandez, Impact of metal contamination in silicon solar cells, Adv. Funct. Mater. 21 (2011) 879–890. doi:10.1002/adfm.201000849
-
[10]
F. Russo, G. Nardone, M. L. Polignano, A. D’Ercole, F. Pennella, M. D. Felice, A. D. Monte, A. Matarazzo, G. Moccia, G. Polsinelli, A. D’Angelo, M. Liverani, F. Irrera, Dark current spectroscopy of transition metals in CMOS image sensors, ECS J. Solid State Sci. Technol. 6 (2017) P217. doi:10.1149/2.0101705jss
-
[11]
Knack, J
S. Knack, J. Weber, H. Lemke, H. Riemann, Copper-hydrogen complexes in silicon, Phys. Rev. B 65 (2002) 165203. doi:10.1103/ PhysRevB.65.165203
2002
-
[14]
J. E. Lowther, Aggregation of interstitial copper atoms in silicon, Mater. Sci. Semicond. Process. 13 (2010) 29–33. doi:10.1016/j.mssp. 2010.02.003
-
[15]
N. Yarykin, J. Weber, Copper-rich complexes in irradiated silicon, J. Appl. Phys. 136 (2024). doi:10.1063/5.0232388
-
[16]
S. K. Estreicher, First-principles theory of copper in silicon, Mater. Sci. Semicond. Process. 7 (2004) 101–111. doi:10.1016/j.mssp.2004. 06.004
-
[17]
K. Shirai, T. Michikita, H. Katayama-Yoshida, Molecular dynamics study of fast diffusion of Cu in silicon, Jpn. J. Appl. Phys. 44 (2005) 7760–7764. doi:10.1143/JJAP.44.7760
-
[18]
K. Matsukawa, N. Hattori, S. Maegawa, K. Shirai, H. Katayama- Yoshida, Getteringmechanismoftransitionmetalsinsiliconcalculated from first principles, Physica B: Condensed Matter 376–377 (2006) 224–226. doi:10.1016/j.physb.2005.12.059
-
[19]
V. P. Markevich, A. R. Peaker, I. F. Medvedeva, V. E. Gusakov, L. I. Murin, B. G. Svensson, Radiation-induced defect reactions in Cz-Si crystals contaminated with Cu, Solid State Phenom. 131–133 (2008) 363–368. doi:10.4028/www.scientific.net/SSP.131-133.363
work page doi:10.4028/www.scientific.net/ssp.131-133.363 2008
-
[20]
K. Shirai, H. Yamaguchi, A. Yanase, H. Katayama-Yoshida, A new structure of Cu complex in Si and its photoluminescence, J. Phys.: Condens. Matter 21 (2009) 064249. doi:10.1088/0953-8984/21/ 6/064249
-
[21]
E. Wright, J. Coutinho, S. Öberg, V. J. B. Torres, A first-principles model of copper–boron interactions in Si: Implications for the light- induced degradation of solar Si, J. Phys.: Condens. Matter 29 (2016) 065701. doi:10.1088/1361-648X/aa4d78
-
[22]
A. Sharan, Z. Gui, A. Janotti, Hybrid-functional calculations of the copper impurity in silicon, Phys. Rev. Appl. 8 (2017) 024023. doi:10.1103/PhysRevApplied.8.024023
-
[23]
T. M. Vincent, S. K. Estreicher, J. Weber, V. Kolkovsky, N. Yarykin, The Cu photoluminescence defect and the early stages of Cu precipita- tion in Si, J. Appl. Phys. 127 (2020) 085704. doi:10.1063/1.5140456
-
[24]
N. Yarykin, J. Weber, S. Lastovskii, V. Gusakov, Copper complexes with carbon-related radiation defects in silicon, Phys. Status Solidi A 218 (2021) 2100141. doi:10.1002/pssa.202100141
-
[25]
P.Chen,Y.Li,F.Qin,T.An,Y.Dai,M.Zhang,M.Liu,L.Zhang, First- principles study of copper contamination in silicon semiconductor, Surf. Interfaces 31 (2022) 102084. doi:10.1016/j.surfin.2022.102084
-
[26]
A.A.Istratov,H.Hieslmair,C.Flink,T.Heiser,E.R.Weber,Interstitial copper-related center in n-type silicon, Appl. Phys. Lett. 71 (1997) 2349–2351. doi:10.1063/1.120026
-
[27]
C.-W. Lee, M. Singh, A. C. Tamboli, V. Stevanović, Transition metal impurities in silicon: Computational search for a semiconductor qubit, npj Comput. Mater. 8 (2022) 1–11. doi:10.1038/s41524-022-00862-z
-
[28]
C. D. Latham, M. Alatalo, R. M. Nieminen, R. Jones, S. Öberg, P. R. Briddon, Passivation of copper in silicon by hydrogen, Phys. Rev. B Xinyu Shi et al.:Preprint submitted to Elsevier Page 8 of 9 Hybrid functional calculation of electrical activity and complexing mechanism of Cu-related defects 72 (2005) 235205. doi:10.1103/PhysRevB.72.235205
-
[29]
N. Yarykin, J. Weber, Deep levels of copper-hydrogen complexes in silicon, Phys. Rev. B 88 (2013) 085205. doi:10.1103/PhysRevB.88. 085205
-
[30]
S.Shirasawa,K.Sueoka,T.Yamaguchi,K.Maekawa, Usefuldatabase of effective gettering sites for metal impurities in Si wafers with first principles calculation, ECS J. Solid State Sci. Technol. 4 (2015) P351. doi:10.1149/2.0051509jss
-
[31]
D. West, S. K. Estreicher, S. Knack, J. Weber, Copper interactions with H, O, and the self-interstitial in silicon, Phys. Rev. B 68 (2003) 035210. doi:10.1103/PhysRevB.68.035210
-
[32]
J. Weber, H. Bauch, R. Sauer, Optical properties of copper in silicon: Excitons bound to isoelectronic copper pairs, Phys. Rev. B 25 (1982) 7688–7699. doi:10.1103/PhysRevB.25.7688
-
[33]
S. D. Brotherton, J. R. Ayres, A. Gill, H. van Kesteren, F. J. A. M. Greidanus, Deep levels of copper in silicon, J. Appl. Phys. 62 (1987) 1826–1832. doi:10.1063/1.339564
-
[34]
M. Steger, A. Yang, N. Stavrias, M. L. W. Thewalt, H. Riemann, N. V. Abrosimov, M. F. Churbanov, A. V. Gusev, A. D. Bulanov, I. D. Kovalev, A. K. Kaliteevskii, O. N. Godisov, P. Becker, H.-J. Pohl, Reduction of the linewidths of deep luminescence centers in 28Si reveals fingerprints of the isotope constituents, Phys. Rev. Lett. 100 (2008) 177402. doi:10.1...
-
[35]
A. Carvalho, D. J. Backlund, S. K. Estreicher, Four-copper complexes in Si and the Cu-photoluminescence defect: A first-principles study, Phys. Rev. B 84 (2011) 155322. doi:10.1103/PhysRevB.84.155322
-
[36]
A. A. Istratov, H. Hieslmair, T. Heiser, C. Flink, E. R. Weber, The dissociation energy and the charge state of a copper-pair center in silicon, Appl. Phys. Lett. 72 (1998) 474–476. doi:10.1063/1.120790
-
[37]
Fujimura, K
T. Fujimura, K. Shirai, Revisiting the stable structure of the Cu4 complex in silicon, Jpn. J. Appl. Phys. 60 (2021) 021001. doi:10. 35848/1347-4065/abd495
2021
-
[38]
G. Kresse, J. Furthmüller, Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set, Phys. Rev. B 54 (1996) 11169–11186. doi:10.1103/PhysRevB.54.11169
-
[39]
Hybrid functionals based on a screened Coulomb potential
J. Heyd, G. E. Scuseria, M. Ernzerhof, Erratum: “Hybrid functionals based on a screened Coulomb potential” [J. Chem. Phys. 118, 8207 (2003)], J. Chem. Phys. 124 (2006) 219906. doi:10.1063/1.2204597
-
[40]
P. E. Blöchl, Projector augmented-wave method, Phys. Rev. B 50 (1994) 17953–17979. doi:10.1103/PhysRevB.50.17953
-
[41]
H. J. Monkhorst, J. D. Pack, Special points for Brillouin-zone integrations, Phys.Rev.B13(1976)5188–5192.doi: 10.1103/PhysRevB. 13.5188
-
[42]
M.E.Straumanis,P.Borgeaud,W.J.James, PerfectionoftheLatticeof Dislocation-Free Silicon, Studied by the Lattice-Constant and Density Method, J. Appl. Phys. 32 (1961) 1382–1384. doi:10.1063/1.1736238
-
[43]
Bludau, A
W. Bludau, A. Onton, W. Heinke, Temperature dependence of the band gap of silicon, J. Appl. Phys. 45 (1974) 1846–1848. doi:10.1063/ 1.1663501
1974
-
[44]
Freysoldt, B
C. Freysoldt, B. Grabowski, T. Hickel, J. Neugebauer, G. Kresse, A. Janotti, C. G. Van de Walle, First-principles calculations for point defects in solids, Rev. Mod. Phys. 86 (2014) 253–305. doi:10.1103/ RevModPhys.86.253
2014
-
[45]
Physica Status Solidi (B) , author =
C. Freysoldt, J. Neugebauer, C. G. Van de Walle, Electrostatic interactions between charged defects in supercells, Phys. Status Solidi B 248 (2011) 1067–1076. doi:10.1002/pssb.201046289
-
[46]
G.Henkelman,B.P.Uberuaga,H.Jónsson, Aclimbingimagenudged elastic band method for finding saddle points and minimum energy paths, J. Chem. Phys. 113 (2000) 9901–9904. doi:10.1063/1.1329672
-
[47]
K. Matsukawa, K. Shirai, H. Yamaguchi, H. Katayama-Yoshida, Dif- fusion of transition-metal impurities in silicon, Physica B: Condensed Matter 401–402 (2007) 151–154. doi:10.1016/j.physb.2007.08.134
-
[48]
H. Weiske, L. Pecher, A. Gallo, A. Irmler, F. Hummel, A. Grüneis, R. Tonner-Zech, Adsorption energies on extended surfaces with CCSD(T) quality: ethylene on Si(001), Molecular Physics 0 (2025) e2602647. doi:10.1080/00268976.2025.2602647
-
[49]
D. A. Drabold, S. K. Estreicher, C. E. Ascheron, H. J. Kölsch, A. H. Duhm (Eds.), Theory of defects in semiconductors, volume 104 of Topics in applied physics , Berlin, Heidelberg, 2007. doi:10.1007/ 11690320
2007
-
[50]
C. D. Latham, M. Ganchenkova, R. M. Nieminen, S. Nicolaysen, M. Alatalo, S. Öberg, P. R. Briddon, Electronic structure calculations for substitutional copper and monovacancies in silicon, Phys. Scr. 2006 (2006) 61. doi:10.1088/0031-8949/2006/T126/014
-
[51]
Status Solidi A 216 (2019) 1900220
R.Ozaki,K.Torigoe,T.Mizuno,K.Yamamoto, Getteringmechanism of copper in n-type silicon wafers, Phys. Status Solidi A 216 (2019) 1900220. doi:10.1002/pssa.201900220
-
[52]
URLhttps:// link.aps.org/doi/10.1103/PhysRevB.110.075130
K. Bystrom, B. Kozinsky, Nonlocal machine-learned exchange functional for molecules and solids, Phys. Rev. B 110 (2024) 075130. doi:10.1103/PhysRevB.110.075130
-
[53]
E. B. Yakimov, Metal impurities and gettering in crystalline silicon, in: Handbook of photovoltaic silicon, 2019, pp. 1–46. doi:10.1007/ 978-3-662-52735-1_23-1
2019
-
[54]
A.Hara,M.Aoki,T.Fukuda,A.Ohsawa, Hydrogeneffectsonoxygen precipitation in Czochralski silicon crystals, J. Appl. Phys. 74 (1993) 913–916. doi:10.1063/1.354858
-
[55]
X. Duo, W. Liu, M. Zhang, L. Wang, C. Lin, M. Okuyama, M. Noda, W.-Y. Cheung, S. P. Wong, P. K. Chu, P. Hu, S. X. Wang, L. M. Wang, Evolution of hydrogen and helium co-implanted single-crystal silicon during annealing, J. Appl. Phys. 90 (2001) 3780–3786. doi:10.1063/1.1389478
-
[56]
F. Jiang, M. Stavola, A. Rohatgi, D. Kim, J. Holt, H. Atwater, J. Kalejs, Hydrogenation of Si from SiN𝑥(H) films: Characterization of H introduced into the Si, Appl. Phys. Lett. 83 (2003) 931–933. doi:10.1063/1.1598643
-
[57]
N.H.Nickel, Hydrogenincorporationinsemiconductors, Phys.Status Solidi B 260 (2023) 2300309. doi:10.1002/pssb.202300309
-
[58]
C. G. Van de Walle, Theory of hydrogen diffusion and reactions in crystalline silicon, Phys. Rev. B 39 (1989) 10791–10808. doi:10.1103/ PhysRevB.39.10791
1989
-
[59]
J. Kaniewski, M. Kaniewska, L. Ornoch, T. Sekiguchi, K. Sumino, Hydrogenation of copper related deep states in n-type Si containing extended defects, Mater. Sci. Forum 258–263 (1997) 319–324. doi:10.4028/www.scientific.net/MSF.258-263.319
work page doi:10.4028/www.scientific.net/msf.258-263.319 1997
-
[60]
Knack, Copper-related defects in silicon, Mater
S. Knack, Copper-related defects in silicon, Mater. Sci. Semicond. Process. 7 (2004) 125–141. doi:10.1016/j.mssp.2004.06.002
-
[61]
N. Yarykin, J. Weber, Formation of copper-related deep-level centers inirradiatedp-typesilicon,in:W.Jantsch,F.Schaffler(Eds.),Gettering anddefectengineeringinsemiconductortechnologyXIV,volume178– 179, Durnten-Zurich, 2011, pp. 154–157. doi:10.4028/www.scientific. net/SSP.178-179.154. Xinyu Shi et al.:Preprint submitted to Elsevier Page 9 of 9
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.