Adsorption of enol hexafluoroacetylacetone (hfacH) on nickel oxide (NiO) fcc (100) and metallic Ni fcc (100) surfaces and the stability of the adsorbate was examined using first-principles quantum mechanical simulations. It was shown that an hfacH molecule can be unstable and dissociate on an Ni metal surface. On an NiO surface; however, an hfacH molecule can be deprotonated and form a hexafluoroacetylacetonate anion (hfac−) bonded stably with positively charged Ni atoms of the surface. The results are consistent with observations of the interaction of hfacH with NiO and Ni surfaces in earlier experiments. The results also explain the mechanisms of the adsorption steps in the thermal atomic layer etching of Ni based on the cyclic processes of surface oxidation and formation of volatile organo-nickel complexes.
I. INTRODUCTION
Atomic layer etching (ALE) via formation of volatile organometallic complexes is expected to establish low-damage and atomically controlled etching processes.1–9 Hexafluoroacetylacetone (hfacH) may be used as an etchant gas for etching various metal oxides including magnetic materials.1,5,6,10 When hfacH gas molecules are exposed to a metal oxide surface at room temperature, they are expected to form organometallic complexes such as Ni(hfac)2, which is volatile at higher surface temperatures. In this process, hydrogen atoms removed from the hfacH also form volatile water molecules and remove oxygen from the surface, whereby the etching of a metal oxide proceeds thermally.1,2,4,7
As an application of such a process with hfacH molecules, a thermal ALE process of metal can be constructed by combining the two steps of oxidation and of etching of the oxidized layer by the formation of volatile organometallic complexes. In this ALE process, it is important that the etchant organic gas molecules do not form volatile organometallic complexes on the metallic surface so that the metal surface is not etched by the etchant gas. In this way, the etching of the oxidized layer stops when the metal surface is exposed to the etchant gas, which is termed a self-limiting step. In the ALE process, the cycle comprising the oxidation step and the self-limiting oxide etching step is repeated a predefined number of times so that the metal surface can be etched to a desired depth in a controlled manner.1,7,11
The self-limited nature of the oxide etching step arises from the fact that hfacH is unstable and dissociates on an Ni surface, which has been confirmed experimentally.1,7 However, the mechanisms of this stability, explaining the reasons why an hfacH molecule decomposes on an Ni surface while it forms Ni(hfac)2 on an NiO surface, have not been well understood prior to this study. As the first step to understand these phenomena, therefore, we examine herein the stability of hfacH molecules adsorbed on Ni and NiO surfaces using first-principles quantum mechanical (QM) simulations.
The goal of this study is to present qualitative insight into the mechanisms of the adsorption of hfacH on Ni and NiO surfaces. More specifically, our aim is to answer why an hfacH molecule tends to be adsorbed stably on an NiO surface, whereas it tends to decompose on an Ni surface. As mentioned above, the adsorption processes play a crucial role in establishing the self-limiting etching step of NiO by hfacH exposure at an elevated temperature. Although we use first-principles QM simulations to achieve this goal, a precise evaluation of the reaction energies for adsorption (i.e., adsorption energies) on a real Ni or NiO surface is beyond the scope of the present study. This is owing to the highly complex nature of the QM simulation for such systems. Therefore, as is shown in the subsequent sections, we employ simple models for the hfacH adsorption processes and analyze the model systems using QM simulation.
II. MODELING
The simulation of this study was performed with gaussian 09 Revision D.0112 software. Our QM simulation is based on the density functional theory (DFT) with the B3LYP functional13–15 and the 6-311G (d) AO basis set.16 The QM simulation was used to evaluate the total energy of a system consisting of multiple atoms (e.g., a single molecule or a system of a molecule interacting with a surface) in its ground state, when the positions of all nuclei were fixed in space. The QM simulation was also used to obtain the most stable atomic configuration (geometric configuration) of a system at zero temperature. The latter is called geometry optimization. Partial geometry optimization can be performed with the positions of some selected atoms being fixed.
The accuracy of system energy evaluated by such a simulation depends on the assumptions used in the numerical techniques. In the DFT framework, the most suitable functionals need to be selected for quantitative evaluation of the system energy. However, it cannot be known a priori which functionals are the most suitable for the system one wants to calculate. It is generally believed that B3LYP offers reasonable computational results for organic systems as well as systems containing transition metals such as Ni.13–15 This is why we used B3LYP in our simulation presented in this study. However, there have been many studies comparing various functionals for DFT calculations,14,15,17–20 which typically show that the selection of best functionals depends on the system and the physical quantities (such as reaction energies, bond lengths, band structures, vibrational frequency shifts, etc.) to be evaluated. Discussion on the most appropriate selection of DFT functionals is beyond the scope of the present work and deferred to a future study.
Similarly, we used the 6-311G(d) basis set as it is widely used in DFT calculations in general.16 It should be noted, however, that the 6-311G(d) basis set does not include the diffuse functions, the use of which could improve the representation of highly polarized atoms. Furthermore, we did not include the dispersion interaction (i.e., van der Waals interaction) in our simulation for the sake of simplicity. The dispersion interactions have much lower attractive interaction energies than those of typical covalent bond interactions but are of much longer range. In this study, we assume that the effects of dispersion interactions on adsorption energies are relatively minor and do not alter our conclusion on the adsorption mechanisms of hfacH on Ni and NiO surfaces. However, the inclusion of dispersion interactions is necessary for the quantitative evaluation of adsorption energies and barriers.
Owing to the high cost of QM simulations in general, it is practically impossible to employ, as a computational model system, a realistic system comprising an hfacH molecule incident upon a surface material that consists of a large number of atoms. Therefore, to represent the material surfaces, we used reduced models of single monolayers that can capture the essential nature of the real NiO and Ni surfaces. Figure 1 shows the top views of the (100) single-monolayer rectangular surfaces of fcc NiO (8.89 × 8.89 Å2) and Ni (7.48 × 4.98 Å2). In each case, the total net charge was assumed to be zero. In the geometry optimization of an hfacH adsorption process performed in this study, the positions of all surface atoms were fixed, while the geometric structure of an hfacH molecule was optimized. In the reality, as an adsorbate molecule approaches the surface, the positions of surface atoms interacting with the adsorbate can change. However, we assumed that such changes in the positions of surface atoms were sufficiently small. The geometry optimization was performed with gradient minimization techniques.
As mentioned earlier, the goal of this study is to present qualitative insight into the mechanisms of the adsorption of hfacH on Ni and NiO surfaces. As the first step, we looked at adsorption processes of an hfacH molecule with Ni and NiO simple single-layer surface models mentioned above. To simplify the simulation further, we have employed several approximations for DFT calculations, as discussed earlier. The validity of the surface models as well as the use of the approximations for DFT calculations must be validated eventually. However, our preliminary DFT calculations of the same adsorption processes with larger and more realistic surface models based on different DFT calculation methods including dispersion interactions have also shown that the results are qualitatively consistent with what is presented in this article. Therefore, we are confident that our simulation results presented in this article help us understand, at least qualitatively, the mechanisms on how an hfacH molecule interacts with an Ni or NiO surface. A systematic study of the accuracy of simulation results presented in this article is deferred to a future work.
III. RESULTS AND DISCUSSION
First, we optimized the structures of the keto and enol types of hfacH using our simulation, which are shown in Fig. 2. The enol-type hfacH is known to be more stable than the keto type. The energy difference obtained from our simulation is −27.48 kJ/mol (−0.285 eV), which is consistent with the results of earlier studies.21–25 This result implies that the majority of hfacH in the gas phase exists in the enol form and, therefore, in the simulations herein, we used only the enol-type hfacH as an adsorbate.
To obtain partial charge distributions, we used a natural bond orbital (NBO) population analysis. Figure 3 shows the NBO charge distributions of both types of hfacH. Figure 3(a) shows that the oxygen atoms are the most negatively charged with a charge of −0.448 e (with e being the elementary charge) in the keto type. Furthermore, Fig. 3(b) shows that the oxygen bonded with the hydrogen atom has a charge of −0.615 e, while the other oxygen atom has a charge of −0.570 e for the enol type. It should be noted that the fluorine atoms are negatively charged (−0.335 to −0.356 e), whereas the carbon atoms of trifluoromethyl (−CF3) groups are highly positively charged (+1.050 to +1.075 e), which makes hfacH a highly polarized molecule.
We also examined the charge distributions of the model surfaces, as shown in Fig. 4. As expected, the simulation shows that NiO is an ionic crystal that consists of positively charged Ni and negatively charged O atoms. In contrast, the Ni atoms of the metallic Ni surface are nearly charge neutral.
Now, we examine how an hfacH molecule interacts with the model NiO surface of Fig. 4(a) when it approaches the surface. There are many possibilities for the orientation of an hfacH molecule when it approaches a surface. In this study, we studied three different orientations of an hfacH molecule: (1) an hfacH molecule vertically aligned toward the NiO surface with the oxygen atoms closer to the surface, as shown in Fig. 5 (this atomic configuration was called “NiO + hfacH”); (2) an hfacH molecule vertically aligned toward the surface but in the opposite direction, with the F atoms being closer to the surface, as shown in Fig. 6 (“NiO + inv. hfacH”); and (3) an hfacH molecule tilted at an angle of 56° from the configuration, as shown in Fig. 7 (“NiO ∼ 56° hfacH”). The tilting angle was measured from the surface normal, so that the vertical hfacH was defined as tilted at an angle of 0°. In all three cases, each O atom of the hfacH molecule was placed nearly above one of the two central Ni atoms of the surface with the same distance from the NiO surface plane.
In the simulation, we first evaluated the “interaction energy” between an hfacH molecule of a fixed geometric configuration and the fixed model surface as a function of distance between the molecule and the model surface. The distance between the hfacH molecule and the surface was defined as the position of the O atoms of the hfacH molecule measured from the NiO surface plane. In this configuration, the two O atoms were placed in parallel to the NiO surface plane. The “interaction energy” was defined as
where Et is the total energy of the system, E1 is the total energy of the hfacH molecule, and E2 is the total energy of the NiO model surface, all at zero temperature. Figure 8 gives the interaction energy from Eq. (1) for an hfacH molecule with an atomic configuration of Fig. 2(b) and an orientation of Figs. 5, 6, or 7, and the model (100) NiO surface of Fig. 4(a), as functions of the distance between the hfacH molecule and the surface. No geometry optimization was performed for this evaluation. In this sense, these calculation results only present a rough idea of how hfacH molecules in the three different orientations given in Figs. 5–7 interact with the NiO surface. Figure 8 suggests that the tilted hfacH molecule interacts with the NiO surface most stably with the lowest interaction energy among the three cases examined here.
Based on this observation, we placed a 56° tilted hfacH molecule at a distance of 2.6 Å (i.e., the lowest-energy position in Fig. 8) above the NiO surface as the initial configuration and performed a geometry optimization. As mentioned earlier, all surface Ni and O atoms were assumed to be immobile but the positions of all atoms of the hfacH molecule were optimized. The atomic configuration of the minimum energy state of this system is shown in Fig. 9. It is seen that the hfacH molecule is deprotonated, with the hydrogen (H) atom between its two oxygen atoms transferred to one of the oxygen atoms of the NiO surface. As a result, the hexafluoroacetylacetonate anion (hfac−) is bonded with the NiO surface with an interaction of the two O atoms of hfac− with two adjacent Ni atoms of the NiO surface. At least qualitatively, the simulation result here showing that an hfacH molecule is adsorbed stably (in the form of an hfac− anion) on an NiO surface is consistent with the experimental observations of Refs. 1 and 7.
The interaction energy, as defined in Eq. (1) with the optimization of Et with fixed positions of all NiO model surface atoms [whose atomic configuration is given in Fig. 9], is −2.02 eV, with the zero-point corrected total energy.12 This is the adsorption energy, or the “reaction energy,” of an hfacH molecule’s adsorption on the NiO surface for the model system. The bond length between an O atom of hfac− and its closest Ni atom of the NiO surface is about 2.05 Å and the hfac− is tilted with an angle of about 46° against the surface normal. The adsorption energy and the bond length obtained in this study are similar to those reported on the adsorption of hfacH on a ZnO surface by Kung and Teplyakov.25 It should be noted, however, the adsorption energy, bond length, and tilting angle given above are obtained from our DFT calculations of simplified models and should not be construed as the exact values of the real system, as we cautioned in Sec. II.
The electrical charge distribution (NBO population) of the system of Fig. 9 is shown in Fig. 10. The NBO charge distribution of the H atom on the NiO surface is 0.844 e and the total charge of hfac− shown in Fig. 10 is −0.844 e. The simulation result suggests that the deprotonation of an hfacH molecule occurs because of the presence of negatively charged O atoms on the surface. In addition, it is clearly seen that the negatively charged O atoms of hfac− are bonded with positively charged Ni atoms. In other words, the presence of ionic bonds; i.e., the fact that Ni atoms are positively charged and O atoms are negatively charged in the NiO surface; is the main reason for the deprotonation and stable adsorption of an hfacH molecule. A similar mechanism should also function for the adsorption of polarized molecules such as diketones on a surface with ionic bonds, such as an oxidized or halogenated metal surface.
We now discuss hfacH adsorption on a metallic Ni surface. It is know that hfacH molecules are dissociated when they approach a metallic Ni surface at room temperature or higher, which has been observed experimentally.1,7 Figure 11 shows the optimized structure of an hfacH molecule placed at the center of the 7.48 × 4.98 Å2 Ni monolayer surface given in Fig. 1(b). The molecule was initially placed at a distance of 2.2 Å above the surface plane with a tilt angle of 56°, in a fashion similar to Fig. 9. As before, the distance was defined as the position of the O atoms of the hfacH molecule measured from the surface plane. The two O atoms were originally placed above the two Ni atoms near the center of the rectangular Ni monolayer surface. It is seen in Fig. 11 that, after a geometry optimization (with the positions of all Ni atoms being fixed), the hfacH molecule moved closer to the Ni surface and was decomposed without deprotonation. Dissociated F atoms went to the opposite side of the monolayer surface, which of course may not occur in an actual Ni surface. The reaction energy for the atomic configurations given in Fig. 11 is −8.70 eV, with the zero-point corrected total energy.12 This large absolute value of the reaction energy arises from the formation of Ni—F and Ni—C bonds after some F atoms are removed from the hfacH molecule. On a metallic Ni surface, such reaction energy varies depending on how the adsorbed hfacH molecule is decomposed on the surface.
The NBO charge distribution of the system of Fig. 11 is shown in Fig. 12. It is seen that Ni atoms of the surface are nearly charge neutral except for those interacting closely with the F and C atoms. The simulation result suggests that, in the case of a metal surface, the negatively charged O atoms and the positively charged H atom between them in an hfacH molecule have no particular surface atoms with which to bond because all surface metal atoms are charge neutral. Therefore, deprotonation of hfacH is unlikely to occur and, if the hfacH molecule is tilted down further, some of its F and C atoms can form bonds with Ni atoms on the surface, as seen in Fig. 11.
IV. SUMMARY AND CONCLUSIONS
In this study, DFT-based first-principles QM simulations were used to clarify why hfacH can be adsorbed stably on an NiO surface while it can be dissociated on a metallic Ni surface. The simulation results showed that, when a stable enol hfacH molecule approaches an NiO surface, its negatively charged O atoms are attracted to Ni atoms of the NiO surface, where all Ni atoms are positively charged. If the distance between the hfacH molecule and the NiO surface is sufficiently short, the system becomes more stable once the H atom located between the two O atoms of hfacH is transferred to a nearby surface O atom. In this process, the dissociated H atom is positively charged and the remaining part, i.e., hfac−, is negatively charged (i.e., deprotonation of hfacH). Because of this charge transfer, the two negatively charged O atoms of hfac− are bonded more strongly with positively charged surface Ni atoms. Once this structure is established, the adsorbed hfac− anion may be slightly tilted but will not fall down on the NiO surface. This prevents the strong interaction of its C and F atoms with surface atoms and the breakdown of the hfac− anion on the surface, as our simulation indicates. In this sense, the adsorption of hfac− anion on an NiO surface is stable.
Conversely, on an Ni metallic surface, every Ni surface is equally charged or charge neutral if the surface is charge neutral. Therefore, when an hfacH molecule approaches a metallic Ni surface, there are no particular preferential sites to which it can stably attach and it is less likely to transfer its H atom or H+ cation to the surface (i.e., to deprotonate). At room temperature or higher, accidental tilting of an hfacH molecule allows its C and F atoms to be bonded with surface Ni atoms, resulting in the dissociation of hfacH on a metallic surface.
These results are qualitatively consistent with experimental observations obtained in earlier studies.1,5,7 The results of this work indicate that the adsorption of hfacH or other diketones (or even other polarized molecules) on a surface is promoted by the presence of ionic bonds (or polarization) on the surface. On an oxidized or halogenated metal surface, surface atoms are polarized, and approaching diketones may similarly deprotonate and have their O atoms bonded with positively charged metal atoms on the surface. On a metal surface, however, such polarized molecules cannot find any preferred adsorption sites and may disintegrate when they approach the surface closely.
As discussed earlier, small surface models were employed herein to represent Ni and NiO surface for analysis via first-principles QM simulations. Furthermore, we have not yet tested various DFT functionals with dispersion interactions or basis sets with the diffuse functions. Therefore, the values of adsorption energies and other physical quantities obtained in this study should not be construed to represent the true values of a real system, where the surface size is far greater than the hfacH molecule. More accurate evaluation of the reaction energies via DFT-based first-principles QM simulations for larger systems is a subject of a future study.
Stable adsorption of hfacH is only the first step of the thermal ALE process of a metal surface. At an elevated surface temperature, adsorbed hfac− anions on a metal oxide surface form volatile organometallic complexes such as Ni(hfac)2. Analysis of the formation process of such organometallic complexes is also beyond the scope of the present study and is deferred to a future work.
ACKNOWLEDGMENTS
A.H.B. appreciates support from Professional Development Consortium for Computational Materials Scientists (PCoMS), which allowed him to visit M.K. and W.W. at KIT, Germany, during this study. A.H.B and S.H. are also grateful to Yoshitada Morikawa and Ikutaro Hamada of Osaka University for fruitful discussion. This work is partially supported by the Japan Society of the Promotion of Science (JSPS) Grants-in-Aid for Scientific Research (S) (No. 15H05736), JSPS Core-to-Core Program, and by Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) under Germany`s Excellence Strategy—2082/1-390761711 and GRK 2450.