pith. machine review for the scientific record. sign in

arxiv: 2604.11675 · v1 · submitted 2026-04-13 · ❄️ cond-mat.mtrl-sci

Recognition: unknown

Hybrid functional calculation of electrical activity and complexing mechanism of Cu-related defects

Authors on Pith no claims yet

Pith reviewed 2026-05-10 15:07 UTC · model grok-4.3

classification ❄️ cond-mat.mtrl-sci
keywords copper defectssiliconHSE06Cu_PLtransition levelsdefect complexesprecipitationelectrical activity
0
0 comments X

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.

The paper uses hybrid functional calculations to determine the structures, formation energies, and charge transition levels of copper-related defects and complexes in silicon. It identifies mismatches between earlier theoretical results and experimental data for the Cu_PL photoluminescence center, then introduces a specific copper interstitial-vacancy cluster as the structure that aligns with observed levels. This identification clarifies which defects are electrically active and how copper begins to aggregate during precipitation. The work matters because copper is a fast-diffusing contaminant that harms silicon device performance, so pinpointing its defect forms aids control of contamination effects.

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

These are editorial extensions of the paper, not claims the author makes directly.

  • 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

Figures reproduced from arXiv: 2604.11675 by An-An Sun, Hao Hu, Meng Chen, Shang-Peng Gao, Siqing Shen, Xinyu Shi, Yongli Liang, Zirui He.

Figure 1
Figure 1. Figure 1: Energy barriers for Cu+ i and Cu0 i diffusion in silicon. Atomic configurations of Cui at initial state (T-site), saddle point (H-site) and final state (T-site) are shown in the insets. 0. 0 0. 2 0. 4 0. 6 0. 8 1 . 0 0. 5 1 . 0 1 . 5 2. 0 2. 5 3. 0 3. 5 4. 0 4. 5 5. 0  [PITH_FULL_IMAGE:figures/full_fig_p005_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Calculated formation energies of isolated defects, including vacancy, Cui , and CuSi. However, previous DFT results show varied results on the exact symmetries of CuSi with different charge states [21, 27, 49]. Our calculations find that structural relaxations with different initial symmetries yield slightly different geometries, but the corresponding formation energy differences are consistently smaller t… view at source ↗
Figure 3
Figure 3. Figure 3: Configurations of (a) Cui -BSi and (b) CuSi-PSi. The red atom is copper, the green atom is boron, and the blue atom is phosphorus. passivates various kinds of deep-level defects, and creates its own electrically and optically active centers [30]. Isolated Hi has only a single donor state and a single acceptor state in the band gap, with a (+1∕ − 1) transition level at 𝐸v + 0.81 eV. H− i and H 0 i prefer a … view at source ↗
Figure 4
Figure 4. Figure 4: Configurations of Cu-H complex. The vertical and horizontal rightward direction corresponds to [001] and [110] direction separatedly. (a) CuSi-Hi with 𝐶2𝑣 symmetry. (b) CuSi-Hi with 𝐶𝑠 symmetry. (c) CuSi-Hi2 with 𝐶2 symmetry. (d) CuSi-Hi2 with 𝐶𝑠 symmetry. (e) Calculated formation energies of CuSi and CuSi-Hin (n=1, 2, 3) defects. either 𝐶𝑠 or 𝐶1 at neutral and +1 charge state. Our calculated results are i… view at source ↗
Figure 5
Figure 5. Figure 5: The calculated binding energy (𝐸b ) of Cu-H complex. The energy needed to absorb or release electrons from the environment (i.e., Si) is included in 𝐸b . 3.4. CuPL defect As is mentioned above, two possible defect complexes, namely CuSi-Cui3 and Cui4V, were proposed to explain the unresolved PL line at 1.014 eV in Cu-contaminated silicon samples [31, 32]. Here we re-examine both structures with the same sc… view at source ↗
Figure 6
Figure 6. Figure 6: (a) Energy barriers between CuSi-Cui3 and Cui4V. Atomic configurations at initial state (CuSi-Cui3), saddle point and final state (Cui4V) are shown in the insets. (b) Calculated formation energies of defects associated with CuPL. (c) Calculated formation energies of CuSi-Cuin(n=0, 1, 2, 3) defects. Xinyu Shi et al.: Preprint submitted to Elsevier Page 6 of 9 [PITH_FULL_IMAGE:figures/full_fig_p008_6.png] view at source ↗
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.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

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)
  1. [§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] §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)
  1. [Abstract and §3] Notation for defects is inconsistent between the abstract (Cu_i4V) and later text; adopt a uniform subscript/superscript convention throughout.
  2. [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

2 responses · 0 unresolved

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
  1. 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

  2. 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

0 steps flagged

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

0 free parameters · 0 axioms · 1 invented entities

Only the abstract is available, so the ledger reflects standard elements of hybrid-DFT defect calculations plus the newly proposed complex; no explicit free parameters or ad-hoc axioms are stated.

invented entities (1)
  • Cu_i4V complex no independent evidence
    purpose: To account for the observed Cu_PL photoluminescence line and resolve discrepancies in transition levels
    Introduced in the abstract as a proposed model based on the HSE06 results; no independent experimental confirmation is mentioned.

pith-pipeline@v0.9.0 · 5499 in / 1090 out tokens · 78796 ms · 2026-05-10T15:07:50.777129+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

59 extracted references · 50 canonical work pages

  1. [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. [2]

    doi:10.1007/978-3-319-93925-4

  3. [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. [4]

    A.A.Istratov,E.R.Weber,Physicsofcopperinsilicon,J.Electrochem. Soc. 149 (2001) G21. doi:10.1149/1.1421348

  5. [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

  6. [6]

    Hiramoto, M

    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. [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. [8]

    Lindroos, H

    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. [9]

    Coletti, P

    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. [10]

    Russo, G

    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. [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

  12. [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

  13. [15]

    Yarykin, J

    N. Yarykin, J. Weber, Copper-rich complexes in irradiated silicon, J. Appl. Phys. 136 (2024). doi:10.1063/5.0232388

  14. [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

  15. [17]

    Shirai, T

    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

  16. [18]

    Matsukawa, N

    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

  17. [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

  18. [20]

    Shirai, H

    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

  19. [21]

    Wright, J

    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

  20. [22]

    Sharan, Z

    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

  21. [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

  22. [24]

    Yarykin, J

    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

  23. [25]

    Interfaces 31 (2022) 102084

    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

  24. [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

  25. [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

  26. [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

  27. [29]

    Yarykin, J

    N. Yarykin, J. Weber, Deep levels of copper-hydrogen complexes in silicon, Phys. Rev. B 88 (2013) 085205. doi:10.1103/PhysRevB.88. 085205

  28. [30]

    Solid State Sci

    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

  29. [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

  30. [32]

    Weber, H

    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

  31. [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

  32. [34]

    Steger, A

    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...

  33. [35]

    Carvalho, D

    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

  34. [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

  35. [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

  36. [38]

    Kresse and J

    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

  37. [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

  38. [40]

    P. E. Blöchl, Projector augmented-wave method, Phys. Rev. B 50 (1994) 17953–17979. doi:10.1103/PhysRevB.50.17953

  39. [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

  40. [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

  41. [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

  42. [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

  43. [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

  44. [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

  45. [47]

    Matsukawa, K

    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

  46. [48]

    Weiske, L

    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

  47. [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

  48. [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

  49. [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

  50. [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

  51. [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

  52. [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

  53. [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

  54. [56]

    Jiang, M

    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

  55. [57]

    doi:10.1002/pssb.202300309

    N.H.Nickel, Hydrogenincorporationinsemiconductors, Phys.Status Solidi B 260 (2023) 2300309. doi:10.1002/pssb.202300309

  56. [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

  57. [59]

    Kaniewski, M

    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

  58. [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

  59. [61]

    Yarykin, J

    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