==== Front Sci Rep Sci Rep Scientific Reports 2045-2322 Nature Publishing Group UK London 78396 10.1038/s41598-020-78396-x Article Anti-oxidant copper layer by remote mode N2 plasma for low temperature copper–copper bonding Park Haesung 1 Seo Hankyeol 2 Kim Sarah Eunkyung eunkyung@seoultech.ac.kr 2 1 grid.412485.e0000 0000 9760 4919Department of Mechanical Engineering, Seoul National University of Science and Technology, 232 Gongneung-ro, Nowon-gu, Seoul, 01811 Korea 2 grid.412485.e0000 0000 9760 4919Graduate School of Nano-IT Design, Seoul National University of Science and Technology, 232 Gongneung-ro, Nowon-gu, Seoul, 01811 Korea 10 12 2020 10 12 2020 2020 10 2172017 7 2020 20 11 2020 © The Author(s) 2020Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.An anti-oxidant Cu layer was achieved by remote mode N2 plasma. Remote mode plasma treatment offers the advantages of having no defect formation, such as pinholes, by energetic ions. In this study, an activated Cu surface by Ar plasma chemically reacted with N free radicals to evenly form Cu nitride passivation over the entire Cu surface. According to chemical state analysis using XPS, Cu oxidation was effectively prevented in air, and the thickness of the Cu nitride passivation was within 3 nm. Based on statistical analysis using the DOE technique with N2 plasma variables, namely, RF power, working pressure, and plasma treatment time, we experimentally demonstrated that a lower RF power is the most effective for forming uniform Cu nitride passivation because of a lower plasma density. When the N2 plasma density reached approximately 109 cm−3 in which the remote mode was generated, high energy electrons in the plasma were significantly reduced and the amount of oxygen detected on the Cu surface was minimized. Finally, low temperature (300 °C) Cu–Cu bonding was performed with a pair of the anti-oxidant Cu layers formed by the remote mode N2 plasma. Cu atomic diffusion with new grains was observed across the bonded interface indicating significantly improved bonding quality over bare Cu–Cu bonding. Subject terms EngineeringMaterials sciencehttp://dx.doi.org/10.13039/501100003725National Research Foundation of KoreaNRF-2018R1A2B6003921http://dx.doi.org/10.13039/501100003052Ministry of Trade, Industry and Energy#20003524issue-copyright-statement© The Author(s) 2020 ==== Body Introduction In advanced packaging such as 3D packaging that vertically stacks the manufactured chips and heterogeneous packaging that places various chips in a limited area, Cu–Cu bonding as electrical connections between the chips is a key process for the high performance and reliability of the entire system1–10. In order to protect a device from thermal damage in device stacking, a low temperature (≤ 300 °C) Cu–Cu bonding process is essential. However, low temperature Cu–Cu bonding is interrupted by the native oxide film such as Cu2O which blocks Cu atomic diffusion at the Cu–Cu bonding interface11,12. Thus, Cu–Cu bonding generally requires a bonding temperature higher than 400 °C to break Cu oxide layer. Many Cu–Cu bonding studies have been reported on lowering the bonding temperature. One such approach is removing the native oxide film that has already formed on a Cu surface. This approach includes a surface activated bonding (SAB) method13,14 that breaks the Cu oxide layer and activates the Cu surface using collision from a high energy Ar beam and a wet treatment method15,16 that chemically removes the Cu oxide layer using various acids, such as sulfuric acid (H2SO4), acetic acid (CH3COOH), and citric acid (C6H8O7). Alternatively, there is a way to prevent Cu oxidation from the beginning. Temporary passivation using a self-assembled monolayer (SAM) of alkane-thiol17,18 and metal passivation using Ti, Pd, Ag, or Au19,20 allow the Cu surface to remain pristine. For high-volume manufacturing (HVM), Cu-Sn alloys have been applied to eutectic bonding thanks to its low melting point with a reasonable electrical conductivity21,22. Although Sn is used as a cap on the top of a Cu pillar bump, it is not suitable for a fine pitch (≤ 10 μm) process due to the reflow characteristics. In addition, the formation of intermetallic compounds (IMC) and the growth of Kirkendall void at the bonding interface degrade mechanical and electrical properties of the Cu/Sn/Cu interconnects. The Direct Bond Interconnect (DBI) technology23,24 is currently used in low temperature hybrid bonding applications, for example, in CMOS image sensors (CIS). That technology forms oxide-oxide bonding, and then a uniform shallow Cu recess is annealed at a low temperature (150–300 °C) to establish a Cu–Cu connection using the difference in the coefficient of thermal expansion between the oxide and Cu. However, that method requires a highly accurate chemical mechanical polishing (CMP) process control, and it is hard to apply it to ultra-fine pitch for the next generation. With this background, we propose a passivation process using a Cu nitride layer in a low temperature Cu–Cu bonding method. This technique has numerous advantages compared to the other methods introduced above. First, our passivation process is CMOS compatible because it is a dry process that uses plasma treatment. It is also inexpensive and simple because the process is performed with a conventional DC sputter used for the Cu thin film deposition by common gases such as argon and nitrogen. Moreover, Cu nitride passivation protects the Cu surface from oxidation in the atmosphere, but it is easily decomposed at a low bonding temperature because of its thermodynamic instability. This characteristic helps to form a homogeneous bonded layer in a low temperature Cu–Cu bonding process. In this study, the Cu nitride passivation was formed using a two-step Ar/N2 plasma process. Because the Cu nitride passivation can be damaged from the accompanying sputtering effect by Ar ion bombardment, the Ar + N2 mixed gas is unfavorable to the formation of uniform and thin copper nitride passivation in manufacturing. Therefore, we propose the two-step Ar/N2 plasma process: Ar plasma pre-treatment was performed first on a Cu surface to remove any contaminants and to activate Cu atoms before the N2 plasma treatment. Then, the N free radicals generated by the dissociation of N2 molecules in the N2 plasma chemically react with the activated Cu atoms at the Cu surface to form Cu nitride passivation. We report the effect of our two-step Ar/N2 plasma process on the Cu surface, the mechanism of the formation of Cu nitride passivation, and the bonding quality of low temperature Cu–Cu bonding using an anti-oxidant Cu layer. Methods Sample fabrication flow and DOE setup An 8-inch Si wafer with a 700 nm thick SiO2 layer was diced into a chip size of 1 cm2. Each piece of diced Si chip specimen was fixed by Kapton tape on an 8-inch dummy Si wafer. Using DC magnetron sputtering, a 50 nm thick Ti adhesion layer and 1 μm thick Cu thin film were deposited on the diced Si chip specimens under 5 mTorr working pressure and 2500 W DC power. After the Cu deposition, two-step Ar/N2 plasma treatment was performed in the same sputter chamber with maintained vacuum conditions. For the plasma treatment, RF power of 13.56 MHz was applied to a table chuck where the specimen was located, and it became a cathode. In this study, we focused on the effect of N2 plasma variables on the formation of Cu nitride passivation. Thus, we fixed the variables of the Ar plasma process while we evaluated the N2 plasma process under the statistically designed conditions. In the design of experiment (DOE) method, central composite design (CCD) was used for a response surface methodology (RSM) which is an efficient DOE method that estimates the regression coefficient of a quadratic polynomial model with only a few experiments25,26. A total of 20 experiments for the N2 plasma process consisted of 8 points from the 3-factor (RF power, working pressure, and treatment time) 2-level (low or high) full factorial design experiments, 6 axial points, and 6 center points. The detailed two-step Ar/N2 plasma process conditions are listed in Table 1 and a study on the Ar plasma process using the full factorial DOE method was reported in previous paper27.Table 1 Independent factors and levels (coded and actual) for CCD design (top), fixed Ar plasma treatment conditions (middle), N2 plasma treatment conditions designed by CCD (bottom). Coded factor Actual factor Coded level and actual level − α − 1 0 + 1 + α A RF power 25.68 70 135 200 244.32 B Pressure 3.48 4.5 6 7.5 8.52 C Time 31.82 100 200 300 368.18 Ar plasma treatment conditions Specimen Flow rate (sccm) RF power (W) Pressure (mTorr) Time (s) Step 1 (fixed) All 150 100 7.5 30 N2 plasma treatment conditions Specimen Flow rate (sccm) RF power (W) Pressure (mTorr) Time (s) Step 2 1 45 70 4.5 100 2 200 4.5 100 3 70 7.5 100 4 200 7.5 100 5 70 4.5 300 6 200 4.5 300 7 70 7.5 300 8 200 7.5 300 9 26 6.0 200 10 244 6.0 200 11 135 3.5 200 12 135 8.5 200 13 135 6.0 32 14 135 6.0 368 15 135 6.0 200 16 135 6.0 200 17 135 6.0 200 18 135 6.0 200 19 135 6.0 200 20 135 6.0 200 For the Cu–Cu bonding process, the above deposition procedure and two-step Ar/N2 plasma process were carried out on an 8-inch Si wafer. A pair of Cu blanket wafers were bonded at 300 °C under 700 kPa for 1 h using a thermo-compression bonding method in the SB 8e wafer bonder of SUSS Microtech. The chemical state of Cu surface and Cu–Cu bonding interface analysis methods The chemical state of the two-step Ar/N2 plasma treated Cu surface was analyzed by X-ray photoelectron spectroscopy (XPS). The peak profiles for three elements, Cu, O, and N, were measured using a Thermo Scientific K-Alpha+ XPS system with a micro-focused monochromatic Al Kα (hν = 1486.6 eV) X-ray source of 72 W at 12 kV. The bonding quality of the bonded wafer was evaluated using scanning acoustic tomography (SAT) of Hitachi FineSAT III with a 140 MHz probe. For a detailed analysis of the bonding interface, the bonded wafer was diced into 1 cm2 specimens and the diced plane was polished by Ar ion beam milling using JEOL Ltd. IB-19510CP. Then the polished clear cross section was analyzed by field emission scanning electron microscope (FE-SEM) using FEI Company Apreo S HiVac with 20,000 × magnifications. Results and discussion The chemical state analysis of the two-step Ar/N2 plasma treated Cu surface by XPS Figures 1a–c compare the XPS profiles of the two-step Ar/N2 plasma treated Cu and non-plasma treated Cu, and Table 2 shows the reference binding energy of several Cu compounds28–31. After the plasma treatment, the main peak of the Cu2p3/2 profile in Fig. 1a shifted towards a higher binding energy (from 932.48 eV to 932.58 eV) indicating that a pure metallic Cu (Cu0) combined with other elements and new substances were created, namely, Cu nitride (Cu1+) or cupric oxide (Cu2+). On the other hand, the main peak of the non-plasma treated Cu shifted towards a lower binding energy (from 932.48 eV to 932.28 eV) than pure Cu due to the formation of cuprous oxide (Cu1+), which indicated oxidation in air. In addition, a noticeable peak at 530.38 eV of the O1s profile in Fig. 1b disappeared while a new peak at 397.38 eV of the N1s profile in Fig. 1c appeared. Therefore, the chemical states of the new peak of the N1s profile and the removed peak of O1s profile are considered as a Cu nitride layer and native oxide grown on Cu, respectively. The stoichiometry of each compound was estimated to be close to Cu4N and Cu2O by atomic concentration calculation (not shown here). In order to analyze the thickness of the Cu4N, XPS depth profiling was performed for 50 s at intervals of 10 s with an etch rate of 0.1 nm/sec. The Cu2p3/2 main peaks of the two-step Ar/N2 plasma treated Cu shown in Fig. 1d and the non-plasma treated Cu shown in Fig. 1g appeared at a binding energy value for pure Cu after 30 s. The Cu4N peak of the N1s profiles shown in Fig. 1f clearly appeared until 30 s of depth profiling, which indicated that the thickness of the Cu nitride was approximately 3 nm. In contrast, there is no peak in Fig. 1i after 10 s because the non-plasma treated Cu has only chemisorbed nitrogen atoms (Nchemisorbed) on the Cu surface. The Cu2O peak of the O1s profile in non-plasma treated Cu shown in Fig. 1h was observed up to 30 s, which means that about 3 nm of native Cu oxide was formed without the plasma treatment. A peak of the O1s profile in the two-step Ar/N2 plasma treated Cu shown in Fig. 1e was observed until 10 s, however, it is considered as the peak due to chemisorbed oxygen atoms generated from a hydroxyl group (–OH) during depth profiling. Using the deconvolution technique, the binding energy, Cu–O peak area, and Cu–N peak area of the Cu2p3/2 profile and the Cu–N peak area of the N1s profile for all the specimens are plotted in Fig. 2. In the Cu2p3/2 profiles, all plasma treated Cu specimens have a Cu–N (Cu4N) peak area that is almost 10 times greater than its Cu–O (CuO) peak area. Moreover, Note that -OH peak area might be included in Cu–O peak area during the deconvolution analysis because the CuO peak is adjacent to the Cu(OH)2 peak within the narrow area of the Cu2p3/2 profile. Therefore, we considered the small area of Cu–O peak of the two-step Ar/N2 plasma treated Cu specimens to be negligible, and the surface of specimen 9 in particular to be an anti-oxidant Cu layer. In sharp contrast, the non-plasma treated Cu specimen had a significant amount of Cu–O (CuO + Cu2O) peak area. This was in agreement with the binding energy shift of the Cu2p3/2 main peak. These results prove that Cu oxide was effectively removed after the two-step Ar/N2 plasma treatment.Figure 1 XPS analysis of Cu2p3/2, O1s, and N1s: (a–c) surface profiles of two-step Ar/N2 plasma treated Cu and non-plasma treated Cu, (d–f) depth profiles of two-step Ar/N2 plasma treated Cu and (g–i) depth profiles of non-plasma treated Cu for 50 s. Table 2 Reference XPS binding energies of the Cu compounds28–31. Linked to Cu2p3/2 (eV) O1s (eV) N1s (eV) Cu 932.0–932.5 – – Cu2O 932.2–932.8 529.9–530.6 – CuO 933.6–934.4 529.4–529.8 – –OH 935.0–935.3 ~ 531.2 – Nchemisorbed – – 398.3–400.0 Cu4N 932.2–932.7 – ~ 397.4 Figure 2 XPS peak analysis by a deconvolution technique. Effect of the N2 plasma variables on the formation of Cu nitride passivation To understand the effect of the N2 plasma variables on the formation of Cu nitride passivation, the correlation between the N2 plasma variables and the Cu4N peak area (also the CuO peak area just for comparison) of the Cu2p3/2 was statistically analyzed at 95% confidence (α = 0.05) using a Minitab program. The Pareto chart of the standardized effects and main effects plot for Cu4N and CuO are shown in Fig. 3a–d. Among the three variables (A: RF power, B: working pressure, and C: treatment time), the RF power (AA) is the most dominant factor for the Cu4N formation, and the working pressure and treatment time do not have much effect. AA is a quadratic term of the variable A, meaning that the variable and response (XPS peak area) have a nonlinear correlation. From the main effects plot, the lower RF power leads to more Cu4N formation, but meanwhile, the analysis results for the CuO are the opposite because CuO production is prevented by Cu4N. Additional analysis results for the Cu4N peak area of the N1s shown in Fig. 3e and the CuO peak area of the O1s shown in Fig. 3f coincide with those for Cu2p3/2 shown in Fig. 3c,d. The effect of RF power in N2 plasma treatment is explained in detail below.Figure 3 Analysis of the response surface design at 95% confidence level using a Minitab program (a) the Pareto chart of the standardized effects for the Cu4N peak of the Cu2p3/2, and the place where the term crosses the reference line is statistically significant, (b) the Pareto chart of the standardized effects for CuO peak of the Cu2p3/2, and (c-f) the main effects plot calculated by the response for Cu4N and CuO in the Cu2p3/2, O1s, and N1s. Effect of the RF power on the N2 plasma density Depending on the applied RF power (PRF), electrons inside the plasma collide with molecules or atoms to transfer energy or dissipate into the wall. The power absorbed by electrons (Pabs) is balanced with the power consumed by generation of electrons and ions through various collisions at a steady state (Ploss). Therefore, if we ignore the power loss that occurred in the RF matching box, we can be assume that PRF = Pabs = Ploss, and the total power loss is as follows32,33. 1 Ploss=nsuBAeεT where ns, uB, A, and eεT represent the electron density at the sheath boundary, Bohm velocity (the sound speed for cold ions), interior area of the plasma chamber, and total loss energy required to generate an electron–ion pair, respectively. The ns has a relationship with the plasma density n0 as ns=hλin0 and where λi is an ion mean free path, h denotes the plasma density ratio to the radial (hR) and axial (hL) directions in the plasma chamber when the radius of the plasma chamber is R and the height is L. Equation (1) can be expressed as Ploss=PRF=hλin0uBAeεT, where hλiA=Aeff=2πR2hL+2πRLhR. Finally, the equation for the plasma density becomes32,34 2 n0=PRFuBeεTAeff The three terms of the denominator are a function of the electron temperature and are negligible for the plasma density within the designed experimental range in this research. On the other hand, the plasma density significantly increases in proportion to the RF power (26–244 W), and the calculated plasma density as a function of the RF power is plotted in Fig. 4. The blue circle marks show the major values of RF power in the CCD (26, 70, 135, 200, and 244 W) in this DOE setup. When the plasma density decreases to ~ 109 with a RF power reduction, the chemical reaction rate decreases due to the low ionization rate of the N2 plasma. Thus, a thick sheath is formed and the plasma becomes a remote mode where very few charged species reach into the substrate35. In this remote mode, many more N free radicals are generated than N2+ and N+ ions by the high ratio of electrons with low energy in the plasma. Consequently, chemical reactions for the formation of Cu nitride passivation by radicals dominate on the Cu surface rather than ion bombardment. As the plasma density increases, the plasma is converted to the direct mode and reactive ion etching (RIE) mode. These modes disturb the formation of the Cu nitride passivation because the sputtering effect caused by energetic ions accompanies it. Table 3 compares the plasma properties according to each plasma mode.Figure 4 Plasma density and mode depending on the RF power. Table 3 Plasma properties of each plasma mode. Plasma mode Remote Direct RIE Density (cm−3) ≈ 109 ~ 1010 1010–1011 Sputtering effect Almost-free Small Large Application Surface treatment (by radical) PECVD plasma etching Reactive ion etching (anisotropic) Schematic diagram Electron energy distribution in the remote mode N2 plasma The ionization, dissociation, and excitation reaction of neutral particles in the plasma depend on the kinetic energy of the electrons, and the kinetic energy of the electrons is proportional to the electron temperature. The electron temperature Te can be calculated simultaneously by Eq. (3) for the power balance equation and particle balance equation32,36. 3 KizTeuBTe=1ngdeff=1ngAeffV where Kiz is the ionization rate coefficient, gas density ng=3.3×1019p and p is the working pressure. Since the plasma chamber size is fixed, we can ignore the terms Aeff and V in Eq. (3). That is, the electron temperature is only affected by the working pressure (3.5 to 8.5 mTorr), and the electron temperature exponentially decreases as the working pressure increases. This is because when the working pressure increases at a constant plasma condition, the amounts of ions and electrons increase, and the heat energy absorbed per electron decreases. As the electron temperature along with the plasma density greatly affects the electron energy distribution function (EEDF), we can infer the chemical reactions occurring in the plasma through the EEDF. If we assume that the EEDF follows the Maxwell distribution, the equation can be written as32,37,38 4 fε=2neπ1(Te)3/2εe-ε/Te where ε is the electron energy and ne is the electron density, which is considered to be equal to the plasma density. Figure 5 shows the constructed EEDF and it demonstrates the electron energy distributions with the electron temperature and the plasma mode. The lower electron temperature, the higher distribution of electrons with lower energy, and the reactive species in the plasma significantly decrease in the remote mode with a low plasma density.Figure 5 Electron energy distribution function (EEDF) with electron temperature and the plasma mode. If the kinetic energy of electrons in the plasma is insufficient, excitation or dissociation occurs more than ionization. In our experiments, most electrons in the plasma had low electron energy (Te < 2.4 eV) because a significant amount of energy was consumed as collisional energy for ionization or dissociation of N2 molecules. Therefore, the unstable state of the N free radicals generated by low-energy electron collision reached the Cu surface, and metastable Cu nitride passivation in the form of Cu4N was produced rather than Cu3N. Table 4 lists the energy required for dissociation and ionization of Ar, N2, and Cu36,39.Table 4 Main reactions for each element in the plasma and required electron energy36,39. Element Reactions in the plasma Electron energy (eV) Ar Ar (g) + e− → Ar+ (g) + 2e− 15.8 N2 N2 (g) + e− → 2 N (g) + e− 9.8 N2 (g) + e− → N2+ (g) + 2e− (g) 15.3 N2 (g) + e− → N+ (g) + N (g) + 2e− 24.3 Cu Ar+ (g) + Cu (s) → Ar+ (g) + Cu+ (g) + e− 11.24 The mechanism of the Cu4N formation and Cu–Cu bonding To study the mechanism of the Cu4N formation, we calculated the standard enthalpy of formation. The standard enthalpy of formation of Cu4N has not been reported yet, so we derived it according to the Born-Haber cycle using the enthalpy of Cu2O and Cu3N from the literature40,41. The standard enthalpies of formation for several Cu compounds (Cu2O, Cu3N, and Cu4N) are calculated in Fig. 6a–c. The equation for the standard enthalpy of formation is expressed as 5 ΔHatomization+ΔHionizationenergy+ΔHelectronaffinity+ΔHlatticeformation=ΔHformation Figure 6 The standard enthalpy of formation according to the Born-Haber cycle (1 eV = 96.485 kJ/mol) (a) Cu2O, (b) Cu3N, and (c) Cu4N. As demonstrated in Fig. 6b, the mechanism of Cu3N formation is explained as follows. In order for solid Cu to be ionized, first the metallic bond must be broken and then it becomes a Cu atom (Cu (s) → Cu (g)), and the atomization energy of Cu is 3.51 eV. The strong triple bond of the N2 molecule is also broken to produce a N free radical (N2 (g) → 2 N (g)), which requires 9.8 eV of dissociation energy. With the first ionization energy of 7.73 eV, the Cu atom removes a valence electron and then becomes a Cu+ ion. At the same time, the neighboring N free radical that reached the Cu surface obtained three electrons, one per Cu and then became a N3− ion (N (g) → N3− (g)). Finally, three positive Cu ions and a negative N ion formed Cu3N through ion bonding with an exothermic reaction, and the standard enthalpy of formation of Cu3N reported by M. Akira is 0.77 eV41. Here, the total energy of the entire chemical reaction is conserved, so the lattice formation enthalpy is − 59.80 eV. The mechanism of Cu4N formation can be explained in the same way as shown in Fig. 6c. Cu4N has one more Cu atom than Cu3N, but the lattice formation enthalpy of Cu4N is assumed to be the same as Cu3N because the Cu atom is interstitially doped into the empty space of Cu3N, which has a cubic anti-ReO3 crystal structure42. Therefore, the standard enthalpy of formation of Cu4N is 13.23 eV, which is quite high compared to Cu2O (− 1.75 eV) as shown in Fig. 6a. Cu2O has a negative enthalpy, so it is chemically stable and hard to be decomposed. It makes Cu–Cu bonding difficult at low temperature because the Cu2O prevents the diffusion of Cu at the bonding interface. However, when the Cu4N is in a metastable state, it can be decomposed at a lower temperature than the bonding temperature, which results in pure Cu–Cu bonding. Hence, this mechanism provides spontaneous Cu–Cu bonding at low temperature. The entire bonding process and mechanism, including the formation and decomposition of Cu nitride passivation, are shown in Fig. 7.Figure 7 Low temperature Cu–Cu bonding mechanism using an anti-oxidant Cu layer by remote mode N2 plasma in the two-step Ar/N2 plasma process, including the formation and decomposition of Cu nitride passivation. The evaluation of the bonded interface of anti-oxidant Cu–Cu bonding A pair of Cu wafers of two types were bonded at 300 °C under 700 kPa for 1 h using a thermo-compression bonding method. Figure 8 shows the SAT and FE-SEM images for both the bare Cu–Cu bonding using non-plasma treated Cu and the anti-oxidant Cu–Cu bonding using remote mode N2 plasma treated Cu in the two-step Ar/N2 plasma process. In the SAT analysis, note that any unbonded voids appear as white color due to the reflection of the ultrasonic signal. As shown in Fig. 8 (left), bare Cu–Cu bonding has an extremely poor bonded interface, which is due to Cu surface oxidation; moreover, the cross-sectional image of the bonded interface in the FE-SEM image clearly shows an unbonded interface. On the other hand, anti-oxidant Cu–Cu bonding demonstrates an improved bonding quality in most of the wafer except for the central region. This unbonded area appears to be caused by trapped N2 gas when the Cu nitride passivation on the entire wafer is thermally decomposed during the bonding process. In addition, there was Cu atomic diffusion with new grains across the bonded interface. Considering that this is whole Cu blanket wafer bonding without a CMP process, it demonstrates a significant improvement and a high potential for low temperature Cu–Cu bonding.Figure 8 SAT and FE-SEM evaluation of the Cu–Cu interface bonded at 300 °C: (left) bonding with non-plasma treated Cu wafers and (right) bonding with two-step plasma treated Cu wafers. Conclusion A low temperature Cu–Cu bonding process is essential in advanced 3D device stacking technologies to protect the devices from thermal damage, and the key point of Cu–Cu bonding is preventing oxidation on a Cu surface especially for die-to-wafer stacking. In this study, we protected the Cu surface from oxidation by the formation of Cu nitride passivation using a two-step Ar/N2 plasma treatment. After the Cu surface activation by Ar plasma pre-treatment, N free radicals in the N2 plasma treatment reacted with the Cu surface to form Cu nitride passivation. According to a chemical state analysis of the two-step Ar/N2 plasma treated Cu surface by XPS and a calculation of the standard enthalpy of formation, the Cu nitride layer effectively prevented the formation of native Cu oxide and it promotes spontaneous Cu–Cu bonding. We statistically analyzed the effect of N2 plasma variables on the formation of Cu nitride passivation by the RSM using a CCD in DOE methodology. Among the N2 plasma variables, RF power was the most effective factor for the formation of Cu nitride passivation, and a lower RF power generated the remote mode N2 plasma. An anti-oxidant Cu layer was obtained in the remote mode N2 plasma with low electron temperature and low plasma density. A pair of anti-oxidant Cu wafers were bonded at 300 °C, and the bonding quality was significantly improved over bare Cu–Cu bonding. In the future, we will optimize the N2 plasma treatment conditions to provide improved Cu–Cu bonding quality at a lower temperature and perform electrical and mechanical evaluations. Publisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations. These authors contributed equally: Haesung Park and Hankyeol Seo. Acknowledgements This research was supported by the Basic Science Research Program through the National Research Foundation of Korea funded by the Ministry of Science and ICT (NRF-2018R1A2B6003921) and also was partially supported by the Ministry of Trade, Industry & Energy (#20003524) and Korea Semiconductor Research Consortium support program for the development of the future semiconductor device. Author contributions All authors discussed the results. H.P. contributed to design of plasma treatment conditions, analyzed the effects of plasma variables, and XPS, SAT data analysis. H.S. contributed specimen fabrication and XPS, FE-SEM data analysis. S.E.K devised research ideas and advised on experimental data. Competing interests The authors declare no competing interests. ==== Refs References 1. Li Y Goyal D 3D microelectronic packaging: from fundamentals to applications Springer Ser. Adv. Microelectron. 2017 57 1 15 10.1007/978-3-319-44586-1_1 2. List, R. S., Webb, C. & Kim, S. E. 3D wafer stacking technology. In Proceedings of the Advanced Metallization Conference. 29–36 (2002). 3. Ko C Chen K Low temperature bonding technology for 3D integration Microelectron. Reliab. 2012 52 302 311 10.1016/j.microrel.2011.03.038 4. Morrow, P., et al. Wafer level 3D interconnects via Cu bonding. In Proceedings of the Advanced Metallization Conference. 125–130 (2004). 5. Panigrahy AK Chen K Low temperature Cu–Cu bonding technology in three-dimensional integration: an extensive review J. Electron. Packag. 2018 140 010801 10.1115/1.4038392 6. Das S Chandrakasan AP Reif R Calibration of Rent's rule models for three-dimensional integrated circuits IEEE VLSL Syst. 2004 12 359 366 10.1109/TVLSI.2004.825833 7. Chen, F. C., Chen, M. F., Chiou, W. C. & Yu, D. H. System on integrated chips (SoICTM) for 3D heterogeneous integration. In Proceedings of the IEEE 69th Electronic Components and Technology Conference (ECTC). 594–599 (2019). 8. Hilton A Temple DS Wafer-level vacuum packaging of smart sensors MDPI Sens. 2016 16 1 33 9. Chua SL Tan CS Cu-Cu die to die surface activated bonding in atmospheric environment using Ar and Ar/N2 plasma Electrochem. Soc. 2016 75 109 116 10. Juang JY Lu CL Chen KJ Chen CC Hsu PN Chen C Tu KN Copper-to-copper direct bonding on highly (111)-oriented nanotwinned copper in no-vacuum ambient Sci. Rep. 2018 8 1 11 10.1038/s41598-018-32280-x 29311619 11. Liu J Mou Y Peng Y Sun Q Chen M Novel Cu-Ag composite nanoparticle paste for low temperature bonding Mater. Lett. 2019 248 78 81 10.1016/j.matlet.2019.03.133 12. Lim DF Wei J Leong KC Tan CS Cu passivation for enhanced low temperature (≤ 300 °C) bonding in 3D integration Microelectron. Eng. 2013 106 144 148 10.1016/j.mee.2013.01.032 13. Takagi H Kikuchi K Maeda R Chung TR Suga T Surface activated bonding of silicon wafers at room temperature Appl. Phys. Lett. 1996 68 2222 2224 10.1063/1.115865 14. He R Fujino M Yamauchi A Wang Y Suga T Combined surface activated bonding technique for low-temperature Cu/dielectric hybrid bonding ECS J. Solid State Sci. Technol. 2016 5 419 424 10.1149/2.0201607jss 15. Jangam, S., et al. Fine-Pitch (≤ 10 μm) direct Cu–Cu interconnects using in-situ formic acid vapor treatment. In Proceedings of the IEEE 69th Electronic Components and Technology Conference (ECTC). 620–627 (2019). 16. Fan J Lim DF Tan CS Effects of surface treatment on the bonding quality of wafer-level Cu-to-Cu thermo-compression bonding for 3D integration J. Micromech. Microeng. 2013 23 045025 10.1088/0960-1317/23/4/045025 17. Tan CS Lim DF Ang XF Wei J Leong KC Low temperature Cu-Cu thermo-compression bonding with temporary passivation of self-assembled monolayer and its bond strength enhancement Microelectron. Reliab. 2012 52 321 324 10.1016/j.microrel.2011.04.003 18. Peng, L., Li, H. Y., Lim, D. F., Gao, S. & Tan, C. S. Thermal reliability of fine pitch Cu-Cu bonding with self assembled monolayer (SAM) passivation for wafer-on-wafer 3D-Stacking. In Proceedings of the IEEE 61th Electronic Components and Technology Conference (ECTC). 22–26 (2011). 19. Panigrahi AK Ghosh T Vanjari SRK Gingh SG Oxidation resistive, CMOS compatible copper-based alloy ultrathin films as a superior passivation mechanism for achieving 150 °C Cu-Cu wafer on wafer thermocompression bonding IEEE Trans. Electron Dev. 2017 64 1239 1245 10.1109/TED.2017.2653188 20. Bonam S Panigrahi AK Kumar CH Vanjari SRK Singh SG Interface and reliability analysis of Au-passivated Cu–Cu fine-pitch thermocompression bonding for 3-D IC applications IEEE Trans. Compon. Packag. Manuf. Technol. 2019 9 1227 1234 10.1109/TCPMT.2019.2912891 21. Sun L Chen MH Zhang L Microstructure evolution and grain orientation of IMC in Cu-Sn TLP bonding solder joints J. Alloys Compd. 2019 786 677 687 10.1016/j.jallcom.2019.01.384 22. Wang J Wang Q Wu Z Wang D Cai J Solid-state-diffusion bonding for wafer-level fine-pitch Cu/Sn/Cu interconnect in 3-D integration IEEE Trans. Compon. Packag. Manuf. Technol. 2016 7 19 26 10.1109/TCPMT.2016.2626794 23. Gao, G., et al. Development of low temperature direct bond interconnect technology for die-to-wafer and die-to-die applications-stacking, yield improvement, reliability assessment. In Proceedings of the IEEE 2018 International Wafer Level Packaging Conference (IWLPC). 1–7 (2018). 24. Enquist, P., Fountain, G., Petteway, C., Hollingsworth, A. & Grady, H. Low cost of ownership scalable copper Direct Bond Interconnect 3D IC technology for three dimensional integrated circuit applications. In Proceedings of the IEEE 2009 International Conference on 3D System Integration. 1–6 (2009). 25. Allen TT Introduction to Engineering Statistics and Lean Six Sigma: Statistical Quality Control and Design of Experiments and SYSTEMS 2019 3 London Springer 335 374 26. Omulo G Banadda N Kabenge I Seay J Optimizing slow pyrolysis of banana peels wastes using response surface methodology Environ. Eng. Res. 2019 24 354 361 10.4491/eer.2018.269 27. Park H Kim SE Two-step plasma treatment on copper surface for low-temperature Cu thermo-compression bonding IEEE Trans. Compon. Packag. Manuf. Technol. 2020 10 332 338 10.1109/TCPMT.2019.2928323 28. Wagner, C. D. Handbook of x-ray photoelectron spectroscopy: A reference book of standard data for use in x-ray photoelectron spectroscopy. In Perkin-Elmer 82–83 (1979). 29. Kim JP Park ES Hong TE Bae JS Ha MG Jin JS Jeong ED Hong KS Electric properties and chemical bonding states of pn-junction p-CuO/n-Si by sol-gel method J. Ceram. Process. Res. 2012 13 s96 s99 30. Siriwardane RV Poston JA Characterization of copper oxides, iron oxides, and zinc copper ferrite desulfurization sorbents by X-ray photoelectron spectroscopy and scanning electron microscopy Appl. Surf. Sci. 1993 68 65 80 10.1016/0169-4332(93)90216-X 31. Venkata Subba Reddy K Uthanna S Effect of sputtering power on the physical properties of Cu3 N films formed by DC magnetron sputtering Synth. React. Inorg. Met. Org. Nano-Met Chem. 2007 37 393 395 10.1080/15533170701465697 32. Lieberman MA Lichtenberg AJ Principles of Plasma Discharges and Materials Processing 2005 2 New York Wiley 33. Balcon N Aanesland A Boswell R Pulsed RF discharges, glow and filamentary mode at atmospheric pressure in argon Plasma Sources Sci. Technol. 2007 16 217 225 10.1088/0963-0252/16/2/002 34. Berenguer C Katsonis K Plasma reactors and plasma thrusters modeling by Ar complete global models Int. J. Aerosp. Eng. 2012 2012 1 18 10.1155/2012/740869 35. Hanyaloglu BF Aydil ES Low temperature plasma deposition of silicon nitride from silane and nitrogen plasmas J. Vac. Sci. Technol. Vac. Surf. Films 1998 16 2794 2803 10.1116/1.581424 36. El-Sayed NM Electron collision rates and energy loss in argon and nitrogen glow discharge plasma Arab. J. Nucl. Sci. Appl. 2018 51 29 38 37. Khalilpour H Foroutan G The effects of electron energy distribution function on the plasma sheath structure in the presence of charged nanoparticles J. Plasma Phys. 2020 86 1 17 10.1017/S0022377820000161 38. Ramamurthi B Economou DJ Kaganovich ID Effect of electron energy distribution function on power deposition and plasma density in an inductively coupled discharge at very low pressures Plasma Sources Sci. Technol. 2003 12 302 312 10.1088/0963-0252/12/3/302 39. Hänninen T Silicon nitride based coatings grown by reactive magnetron sputtering Linköping Stud. Sci. Technol. Diss. 2018 1901 12 17 40. Wagman DD The NBS Tables of chemical thermodynamic properties: selected values for inorganic and C1 and C2 organic substances in SI units J. Phys. Chem. Ref. Data 1982 11 154 159 41. Miura A Takei T Kumada N Synthesis of Cu3 N from CuO and NaNH2 J. Asian Ceram. Soc. 2014 2 326 328 10.1016/j.jascer.2014.08.007 42. Lu N Ji A Cao Z Nearly constant electrical resistance over large temperature range in Cu3 NMx (M = Cu, Ag, Au) compounds Sci. Rep. 2013 3 1 6