Next Article in Journal
Trophic Transfer of Single-Walled Carbon Nanotubes at the Base of the Food Chain and Toxicological Response
Next Article in Special Issue
Distributed Polarizability Model for Covalently Bonded Fullerene Nanoaggregates: Origins of Polarizability Exaltation
Previous Article in Journal
Engineering the Exchange Spin Waves in Graded Thin Ferromagnetic Films
Previous Article in Special Issue
Influence of Group-IVA Doping on Electronic and Optical Properties of ZnS Monolayer: A First-Principles Study
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

First-Principles Molecular Dynamics Simulations on Water–Solid Interface Behavior of H2O-Based Atomic Layer Deposition of Zirconium Dioxide

1
Key Laboratory of Carbon Materials of Zhejiang Province, College of Chemistry and Materials Engineering, Wenzhou University, Wenzhou 325035, China
2
National Laboratory of Solid State Microstructures, College of Engineering and Applied Sciences, Nanjing University, Nanjing 210093, China
*
Authors to whom correspondence should be addressed.
Nanomaterials 2022, 12(24), 4362; https://doi.org/10.3390/nano12244362
Submission received: 14 November 2022 / Revised: 24 November 2022 / Accepted: 28 November 2022 / Published: 7 December 2022
(This article belongs to the Special Issue First-Principles Investigations of Low-Dimensional Nanomaterials)

Abstract

:
As an important inorganic material, zirconium dioxide (ZrO2) has a wide range of applications in the fields of microelectronics, coating, catalysis and energy. Due to its high dielectric constant and thermodynamic stability, ZrO2 can be used as dielectric material to replace traditional silicon dioxide. Currently, ZrO2 dielectric films can be prepared by atomic layer deposition (ALD) using water and zirconium precursors, namely H2O-based ALD. Through density functional theory (DFT) calculations and first-principles molecular dynamics (FPMD) simulations, the adsorption and dissociation of water molecule on the ZrO2 surface and the water–solid interface reaction were investigated. The results showed that the ZrO2 (111) surface has four Lewis acid active sites with different coordination environments for the adsorption and dissociation of water. The Zr atom on the surface can interacted with the O atom of the water molecule via the p orbital of the O atom and the d orbital of the Zr atom. The water molecules could be dissociated via the water–solid interface reaction of the first or second layer of water molecules with the ZrO2 (111) surface. These insights into the adsorption and dissociation of water and the water–solid interface reaction on the ZrO2 surface could also provide a reference for the water–solid interface behavior of metal oxides, such as H2O-based ALD.

Graphical Abstract

1. Introduction

Zirconium dioxide (ZrO2) is an important inorganic material and has wide applications in the fields of microelectronics, coating, catalysis, energy storage and conversion [1,2,3,4,5,6]. With the further scaling-down of metal-oxide semiconductor field-effect transistor (MOSFET) devices, the traditional dielectric material, silicon dioxide (SiO2), fails to overcome the quantum mechanical tunneling effect in the devices [7,8,9,10,11,12,13]. Due to its high dielectric constant and thermodynamic stability, ZrO2 has been used as the dielectric material to replace SiO2. Currently, ZrO2 films can be prepared by physical vapor deposition (PVD), chemical vapor deposition (CVD) and atomic layer deposition (ALD) methods. Among these techniques, ALD has an obvious advantage for the gate dielectric ultrathin film for the nanodevices. ALD can accurately control the composition and thickness of the deposited thin films at the atomic level by alternately pumping different precursors into the reactor. It can utilize the self-limiting and self-saturation features of surface reactions, which result in the advantages of excellent three-dimensional conformity, large-area uniformity and precise thickness control of thin films [14,15,16,17,18]. Usually, water can be used as an oxygen source for the ALD preparation of metal oxides, such as Al2O3, TiO2, ZrO2 and HfO2. Thus, the water–solid interface is an important scientific issue for H2O-based ALD of metal oxides. Meanwhile, the water–solid interface also performs very important roles in various fields of nature including microelectronics, coating and heterogeneous catalysis [19,20,21,22].
ZrO2 has three crystalline phases, namely, common monoclinic, tetragonal and cubic crystals. The adsorption behavior of water molecules on the surface of ZrO2, such as cubic ZrO2 (111) and (110) and monoclinic ZrO2 (001), ( 1 ¯ 01) and ( 1 ¯ 11), have been studied [23,24,25,26,27,28,29,30]. For example, Korhonen et al. studied the adsorption of water on the most stable ( 1 ¯ 01) and ( 1 ¯ 11) surfaces of monoclinic ZrO2 using density functional theory (DFT) and infrared spectroscopy [26,27]. Xia et al. examined the hydration of the m-ZrO2 ( 1 ¯ 11) surface using ab initio molecular dynamics simulations [29]. By obtaining the phase diagrams of surface dehydration, the representative partially hydrated m-ZrO2 ( 1 ¯ 11) surfaces at various temperatures are illustrated. The Fourier transform infrared (FTIR) spectroscopy has demonstrated that the m-ZrO2 surfaces have amounts of acidic and basic hydroxyl groups and undissociated coordinated water molecules [30]. Herein, we focused on the adsorption and dissociation of water molecules on the relatively stable surface of monoclinic ZrO2 (111) and its water–solid interface reactions using DFT calculations and first-principles molecular dynamics (FPMD) simulations. The structure and density of states (DOS) before and after adsorption and dissociation of water are also investigated. The results show that the Zr atom of the ZrO2 (111) surface is the active site for the adsorption and dissociation of water. The water molecules close to the surface will be strongly adsorbed on the substrate and form Zr–O bonding sites. Some water molecules will be directly dissociated and form O–H bonds with the O atoms on the substrate surface. The dissociation of water can also occur via a proton exchange reaction between two water molecules on the surface of the monoclinic ZrO2 (111) surface. These modeling and insights into the water–solid interface reaction of ZrO2 via FPMD simulations at the atomic level can provide theoretical references for water–solid interface studies of metal oxides, such as H2O-based ALD.

2. Computational Details

First, the monoclinic ZrO2 crystal (space group P21/c) was optimized. The optimized lattice parameters are 5.16, 5.23 and 5.35 Å, respectively, which are in accordance with the experimental values of 5.15, 5.21 and 5.31 Å [31]. The average length of the Zr–O bond is 2.295 Å. To model the surface, a monoclinic ZrO2 (111) (1 × 2) surface with 13.37 × 7.40 Å2 was adopted. The model includes two, three, four and five layers of ZrO2. For the structural optimization, the atoms of the lowest layer of the surface were fixed. To ensure that there is enough space to adsorb water molecules and avoid the influence between two interfaces, the vacuum between the two surfaces was set to 20 Å. For the FPMD simulations, three layers of O–Zr–O including 72 atoms were used. At the same time, 55 water molecules were placed in the 20 Å vacuum.
All geometric optimizations were performed using DFT calculations based on the generalized gradient approximation with the Perdew–Burke–Ernzerhof functional and implemented in the Vienna Ab initio Simulation Package (VASP 6.2) [32]. The valence electronic states were expanded based on plane waves with the core–valence interaction represented using the projector augmented plane wave approach [33,34]. An energy cutoff of 520 eV was used. A Γ-centered k-mesh of 7 × 6 × 6 was used for the calculations of the crystal. A Γ-centered k-mesh of 2 × 4 × 1 was used for the calculations of the surface. The convergence criteria for total energies and the forces for ionic relaxation were 1 × 10−5 eV and 0.02 eV/Å, respectively. For the FPMD simulations, the NVT ensemble was used. The temperature was controlled at 373 K using a Nose-Hoover thermostat. The time step was set to 1.0 fs. The total time of the FPMD simulation was 15 ps.

3. Results and Discussion

3.1. Surface Model

Surface model tests were performed using two, three and four layers of ZrO2 surface, shown in Figure 1. The water molecule was adsorbed on the same site of the ZrO2 surface models. The structural parameters and the adsorption energy are shown in Table 1. Here, the adsorption energy can be defined as:
Eads = E(substrate + adsorbate)EsubstrateEadsorbate
where E(substrate + adsorbate), Esubstrate and Eadsorbate represent the energies of the system, surface and adsorbate, respectively. The results show that there is almost no difference in the adsorption configuration of the water molecule on the surface models of two, three and four layers. The three adsorption energies are about −18.5 kcal/mol. Due to Zr atoms on the surface are strong Lewis acid sites, water can be strongly adsorbed on the ZrO2 surface, indicating that this reaction is driven by strong adsorption of adsorbates at strong Lewis acid sites. Meanwhile, all of this indicates that after the surface model reaches a certain thickness, the increase in the substrate thickness has little effect on the configuration and adsorption energy of water molecules on the surface. Therefore, to balance the computational accuracy and cost, the three-layer model of the ZrO2 surface was adopted for the following calculations.

3.2. Adsorption and Dissociation of a Water Molecule on the ZrO2 (111) Surface

As shown in Figure 2, Zr atoms of monoclinic ZrO2 bulk have seven-fold coordination with O atoms and O atoms have three- or four-fold coordination with Zr atoms. Due to the difference in the relative spatial positions of O atoms coordinated to Zr atoms, there are four nonequivalent Zr atoms with six-fold coordination on the surface of monoclinic ZrO2 (111). In other words, there are four Lewis acid active sites of surface Zr atoms with different coordination environments on the monoclinic ZrO2 (111) surface, labeled A, B, C and D, respectively. Considering the four orientations of a water molecule at a site, 16 initial structures of a water molecule adsorbed on the surface were investigated. The most stable adsorption structures at the four sites are shown in Figure 3. The corresponding adsorption energies were −20.8, −18.6, −19.5 and −34.0 kcal/mol, respectively.
Table 2 lists the bond lengths and bond angles of the four stable adsorption structures. At site A, the distances between O and H in the water molecule after adsorption are 0.979 and 1.004 Å, and ∠HOH is about 104.4°, whereas the bond length of H–O in the isolated water molecule is 0.97 Å and ∠HOH is about 103.6°. The distance between the Zr atom at site A and the O atom in the water molecule is 2.364 Å. There is an interaction between the Zr atom on the surface and the O atom in the water molecule. Similar trends were observed in the variation of structural parameters at the B, C and D. In these adsorption structures, the structure of the water molecule is distorted in value by about 5%. The distances between the Zr atom on the surface and the water molecules are all about 2.4 Å. All of these indicate that there are strong interactions between the Zr atoms on the surface and the O atoms in the water molecules.
Based on the four most stable adsorption structures, the dissociation states of the water molecule on the ZrO2 surface were also investigated. As shown in Figure 3, the water molecule can be dissociated into a H atom and the hydroxyl group, where the hydroxyl group is adsorbed on the Zr atom on the surface and the H atom interacts with the O atom on the surface to form a new hydroxyl group. Table 3 lists the bond lengths and bond angles of the dissociated structures. At the A site, after the dissociation of the water molecule, the distances of O and H in the water molecule are 0.972 and 1.642 Å, and the ∠HOH is about 109.9°. The distance between the O atom of the water molecule and the Zr atom on the surface decreases from 2.364 Å to 2.080 Å. Similarly, at the B (0.969/1.691 Å and 136.5°, 2.085 Å), C (0.973/1.774 Å and 106.7°, 2.073 Å) and D (0.972/1.465 Å and 110.8°, 2.091 Å) sites, the water molecules are dissociated. The O atom in the water molecule and the Zr atom on the surface are closer together. The distances between Zr and O atoms are below 2.1 Å. All these indicate that the water molecule is dissociated on the ZrO2 surface and a new Zr–O bond is formed.
For comparison, it can be seen that the adsorption energies of the water molecule at the four sites are −20.8, −18.6, −19.5 and −34.0 kcal/mol, respectively. The energies of the dissociated structures of the water molecules are −19.3, −16.5, −12.3 and −33.9 kcal/mol, respectively. These indicate that the adsorption energy and stability of the adsorbed structure are almost the same as that of the dissociated structure at the four sites. Among them, the energy of the dissociation structure at the C site is higher by 6.8 kcal/mol than that of the adsorption structure. This is because the different coordination environments of Zr atoms on the surface affect the adsorption and dissociation of small molecules. At the D site, the adsorption energy of water molecules is the largest, indicating that the adsorption structure at the D site is the most stable. The bonding characteristics between the water molecule and Zr atom at the D site on the substrate were further analyzed. The changes in the projected DOS of adsorbed water and dissociation of a water molecule on the D site are shown in Figure 4. The adsorption of a water molecule on the surface of ZrO2 (111) is mainly through the interaction between the O atom in the water molecule and the Zr atom on the surface. The highest occupied molecular orbital of the water molecule is mainly contributed by the 2p orbital of the O atom and the lowest unoccupied molecular orbital of the O atom is mainly contributed by the 2s orbital of the O atom. The electrons of the d orbitals of Zr atoms on the surface fill the nonbonding orbital of the O atom. From Figure 4, it can be seen that the DOS of the water molecule is shifted from lower to higher energy levels when the adsorbed state changes to the dissociated state. These are also in accordance with the results of the adsorption and dissociation of water on the HfO2 surface [35]. This shows that the water molecule gets more charge in the adsorbed state than in the dissociated adsorbed state. Moreover, the peak of the Zr atom at the Fermi energy level becomes lower in the dissociated state, indicating that the electrons are transferred from the surface to the water molecule.

3.3. The Water–Solid Interface Reaction of ZrO2

To further investigate the adsorption and dissociation behaviors of a large number of water molecules on the ZrO2 surface, FPMD simulations of the water–solid interface of ZrO2 were conducted at 373 K. Figure 5 shows the three-layer model structure of the water–solid interface of ZrO2, in which 55 water molecules are added into the vacuum slab as the explicit solvent corresponding to the density (1 g/mL) of the water.
To understand the distribution and movement rules of the particles on the surface, the radial distribution functions (RDFs) and their integrated RDFs for H and O in water molecules and for the Zr atom on the ZrO2 surface were calculated. The RDF g(r) is the distribution probability of a given particle α over a particle β when the distance is r. As shown in Figure 6, the first peak in g(r)H–H appears around 1.5 Å, showing H–H interactions between the same water molecules. The second peak in g(r)H–H appears around 2.3 Å, showing the interaction between H atoms of different water molecules. The first peak in g(r)O*–H is 1.0 Å and corresponds to the O–H bond. The second peak in g(r)O*–H appears at 1.6–2.1 Å, which shows strong hydrogen bonding between the H and O atoms in different water molecules. The first peak of g(r)Zr*–O* is 2.15 Å and corresponds to the interaction between the Zr atom on the surface and the O atom in a water molecule, which also indicates that the water molecules are strongly adsorbed and dissociated on the ZrO2 surface. The value of the integral RDF corresponding to the first peak of g(r)Zr*–O* is 0.9, indicating that the first coordination layer of the Zr atom on the surface has 0.9 O atoms and the Zr atom has approximately one water molecule adsorbed on it. The peak of g(r)O**–H between O atoms on the surface of ZrO2 and H atoms in the water molecule is about 1.0 Å and further indicates hydroxyl (–OH) formation on the ZrO2 surface. The value of the integrated RDF corresponding to the first peak of g(r)O**–H indicates that the first coordination layer of the O atom has 0.4 H atoms on the surface. All these results above show that the water molecules mainly interact with each other via hydrogen bonding. The water molecules close to the surface will be strongly adsorbed at the Zr atoms on the surface to form Zr–O bonding sites. Meanwhile, the H atoms of some water molecules are dissociated to form O–H bonds with the O atoms on the surface.
The adsorption and dissociation behaviors of H2O molecules on the ZrO2 surface can be obtained by performing FPMD simulations of the water–solid interface of ZrO2. Figure 7 and Figure 8 show the evolutions of O–H and Zr–O bonds of the interface reaction between the first layer of water and the ZrO2 surface, which is based on the four-membered ring (4MR) pathway. Initially, at the surface A site, the distance between the O(1) atom in the water molecule and the Zr(1) atom on the surface is around 6.1 Å. With the evolution of the FPMD simulation, the water molecules can go closer to the surface. At about 1000 fs, the distance between the Zr(1) and O(1) atoms decreases to around 2.141 Å, which implies the formation of Zr–O bonds. Similarly, the distance between the O(1) atom and the H(1) atom in the water molecule increases from 1.0 Å to 1.7 Å and the distance between the H(1) atom and the O(2) atom on the ZrO2 substrate decreases from 5.0 Å to about 1.0 Å. All these indicate the breakage of the O–H bond in the water molecule and the formation of the H–O bond between the H atom in the water molecule and the O atom on the substrate surface. The distance between the O(1) atom and the H(1) atom remains around 1.7 Å, indicating that strong hydrogen bonds still exist between them. To verify the validity of the three-layer ZrO2 model, we also performed FPMD simulations with a five-layer ZrO2 surface model, shown in Figure 9 and Figure S1. The simulation time of the FPMD was 1 ps. From Figure 9 and Figure S1, we obtained similar behavior of the water–solid interface between three-layer and five-layer ZrO2 surface models.
In addition to the dissociation of the first layer of water, the second layer of water on the surface also participates in the interface reaction. This water–solid interface reaction is based on the six-membered ring (6MR) pathway. Figure 10 and Figure 11 show the evolutions of the bond lengths of O–H and Zr–O bonds at the reaction center. The whole reaction process involves the ZrO2 surface at site C and two water molecules, H2O(a)(H(1)O(1)H) and H2O(b)(H(2)O(2)H), via a 6MR pathway. In the beginning, the distances of the O(1) atom and H(1) atom, the H(2) atom and the O(2) atom are both 1.0 Å. The distances between the H(1) atom and O(2) atom, and H(2) and O(3) atom are 3.3 and 3.7 Å, respectively. The distance between Zr(1) and O(1) atoms is about 2.5 Å. At 500 fs, the H(1)–O(2) and H(2)–O(3) distances gradually become shorter and their distances both decrease to about 1.6 Å. The Zr(1)–O(1) distance also decreases to 2.1 Å. All these indicate that H(1)–O(2), H(2)–O(3) and Zr(1)–O(1) bonds have the tendency to form. At 1400 fs, the H(1)–O(2) and H(2)–O(3) distances decrease to 1.0 Å. In contrast, the O(1)–H(1) and O(2)–H(2) distances increase to 1.7 Å. Afterward, the lengths of O(1)–H(1) and O(2)–H(2) bonds oscillate between 1.5 and 2.5 Å. At 2200 fs, the O(1)–H(1) and O(2)–H(2) distances oscillate around 3.0 Å. The distance between O(2) and H(1) is reduced to about 1.0 Å.
All these show that the H2O (a) molecule can be dissociated to form the hydroxyl group, which can be strongly adsorbed to the Zr(1) atom on the substrate surface. This whole process of distance change shows that the dissociated H(2) atom in H2O (b) dissociated from the H2O (b) molecule can form a H(2)–O(3) bond with the O(3) atom on the surface. The remaining hydroxyl portion of the H2O (b) molecule can combine with the dissociated H(1) atom in H2O (a) to form a new water molecule. DFT calculations and FPMD simulation results show that water molecules can not only be dissociated directly on the ZrO2 surface but also exchange the proton with the help of adjacent water and finally be dissociated.

4. Conclusions

In summary, the adsorption and dissociation behaviors of water molecules on the surface of monoclinic ZrO2 (111) were investigated using DFT static calculations and FPMD simulations. Four Lewis acid active sites of surface Zr atoms with different coordination environments have different adsorption and dissociation states and energies of water molecule on the ZrO2 (111) surface. DOS analysis shows that the O atom of the water molecule prefers to bond with the Zr atom on the ZrO2 (111) surface. There are strong interactions between the p orbitals of the O atom and the d orbitals of the Zr atom. The electrons of the d orbitals of the Zr atom on the surface fill the nonbonding orbitals of O atoms. H atoms of the water molecule prefer to bond with O atoms on the surface. Through FPMD simulations, the water–solid interface reaction of the ZrO2 (111) surface was analyzed. The water molecules mainly interact with each other by hydrogen bonding. Through the 4MR or 6MR pathways of water–solid interface reactions, the first and second water molecules on the surface can be adsorbed and dissociated at the Zr sites on the surface to form Zr–O and O–H bonds. These insights can provide the theoretical significance for the adsorption and dissociation of water and the water–solid interface behavior of metal oxides, such as H2O-based ALD.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/nano12244362/s1, Figure S1: Radial distribution functions (RDFs) and their integrated RDFs of the interface reaction between the first layer water and ZrO2 surface on the five-layer ZrO2 model.; Table S1: Atomic coordinates of 2-layer ZrO2 substrate; Table S2: Atomic coordinates of 3-layer ZrO2 substrate; Table S3: Atomic coordinates of 4-layer ZrO2 substrate; Table S4: Atomic coordinates of 5-layer ZrO2 substrate; Table S5: Atomic coordinates of 3-layer ZrO2 substrate 55 waters.

Author Contributions

Conceptualization, G.F.; methodology, G.F.; software, G.F. and L.X.; validation, G.F. and R.X.; investigation, R.X., Z.Z. and Y.W.; data curation, R.X.; writing—original draft preparation, R.X.; writing—review and editing G.F.; visualization R.X.; supervision, H.X., L.X., Y.D., X.L. and A.L.; funding acquisition, G.F. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by National Natural Science Foundation of China (21873073 and 52073142), Zhejiang Provincial Natural Science Foundation (LY22B030007), Open Project (M30031 and M35047) of National Laboratory of Solid State Microstructures at Nanjing University.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Zaera, F. The surface chemistry of thin film atomic layer deposition (ALD) processes for electronic device manufacturing. J. Mater. Chem. 2008, 18, 3521–3526. [Google Scholar] [CrossRef]
  2. Rolison, D.R.; Long, J.W.; Lytle, J.C.; Fischer, A.E.; Rhodes, C.P.; McEvoy, T.M.; Bourg, M.E.; Lubers, A.M. Multifunctional 3D nanoarchitectures for energy storage and conversion. Chem. Soc. Rev. 2008, 38, 226–252. [Google Scholar] [CrossRef]
  3. Marichy, C.; Bechelany, M.; Pinna, N. Atomic Layer Deposition of Nanostructured Materials for Energy and Environmental Applications. Adv. Mater. 2012, 24, 1017–1032. [Google Scholar] [CrossRef]
  4. O’Neill, B.J.; Jackson, D.H.K.; Lee, J.; Canlas, C.; Stair, P.C.; Marshall, C.L.; Elam, J.W.; Kuech, T.F.; Dumesic, J.A.; Huber, G.W. Catalyst Design with Atomic Layer Deposition. ACS Catal. 2015, 5, 1804–1825. [Google Scholar] [CrossRef] [Green Version]
  5. Palmstrom, A.F.; Santra, P.K.; Bent, S.F. Atomic layer deposition in nanostructured photovoltaics: Tuning optical, electronic and surface properties. Nanoscale 2015, 7, 12266–12283. [Google Scholar] [CrossRef]
  6. Asundi, A.S.; Raiford, J.A.; Bent, S.F. Opportunities for Atomic Layer Deposition in Emerging Energy Technologies. ACS Energy Lett. 2019, 4, 908–925. [Google Scholar] [CrossRef]
  7. Gaskell, J.M.; Jones, A.C.; Chalker, P.R.; Werner, M.; Aspinall, H.C.; Taylor, S.; Taechakumput, P.; Heys, P.N. Deposition of lanthanum zirconium dioxide high-κ films by liquid injection ALD and MOCVD. Chem. Vap. Depos. 2007, 13, 684–690. [Google Scholar] [CrossRef]
  8. Dezelah, C.L., IV; Niinistö, J.; Kukli, K.; Munnik, F.; Lu, J.; Ritala, M.; Leskelä, M.; Niinistö, L. The atomic layer deposition of HfO2 and ZrO2 using advanced metallocene precursors and H2O as the oxygen source. Chem. Vap. Depos. 2008, 14, 358–365. [Google Scholar] [CrossRef]
  9. Kaipio, M.; Blanquart, T.; Banerjee, M.; Xu, K.; Niinistö, J.; Longo, V.; Mizohata, K.; Devi, A.; Ritala, M.; Leskelä, M. Atomic layer deposition of TiO2 and ZrO2 thin films using heteroleptic guanidinate precursors. Chem. Vap. Depos. 2014, 20, 209–216. [Google Scholar] [CrossRef]
  10. Jung, J.-S.; Lee, S.-K.; Hong, C.-S.; Shin, J.-H.; Kim, J.-M.; Kang, J.-G. Atomic layer deposition of ZrO2 thin film on Si(100) using {η51-Cp(CH2)3NMe}Zr(NMe2)2/O3 as precursors. Thin Solid Film. 2015, 589, 831–837. [Google Scholar] [CrossRef]
  11. Kanomata, K.; Tokoro, K.; Imai, T.; Pansila, P.; Miura, M.; Ahmmad, B.; Kubota, S.; Hirahara, K.; Hirose, F. Room-temperature atomic layer deposition of ZrO2 using tetrakis(ethylmethylamino)zirconium and plasma-excited humidified argon. Appl. Surf. Sci. 2016, 387, 497–502. [Google Scholar] [CrossRef] [Green Version]
  12. Mahuli, N.; Cavanagh, A.S.; George, S.M. Atomic layer deposition of hafnium and zirconium oxyfluoride thin films. J. Vac. Sci. Technol. A Vac. Surf. Film. 2021, 39, 022403. [Google Scholar] [CrossRef]
  13. Xu, W.; Lemaire, P.C.; Sharma, K.; Gasvoda, R.J.; Hausmann, D.M.; Agarwal, S. Mechanism for growth initiation on aminosilane-functionalized SiO2 during area-selective atomic layer deposition of ZrO2. J. Vac. Sci. Technol. A Vac. Surf. Film. 2021, 39, 032402. [Google Scholar] [CrossRef]
  14. Walter, E.J.; Lewis, S.P.; Rappe, A.M. First principles study of carbon monoxide adsorption on zirconia-supported copper. Surf. Sci. 2001, 495, 44–50. [Google Scholar] [CrossRef] [Green Version]
  15. Grau-Crespo, R.; Hernández, N.C.; Sanz, J.F.; de Leeuw, N.H. Theoretical Investigation of the Deposition of Cu, Ag, and Au Atoms on the ZrO2(111) Surface. J. Phys. Chem. C 2007, 111, 10448–10454. [Google Scholar] [CrossRef]
  16. Gennard, S.; Corà, F.; Catlow, C.R.A. Comparison of the Bulk and Surface Properties of Ceria and Zirconia by ab Initio Investigations. J. Phys. Chem. B 1999, 103, 10158–10170. [Google Scholar] [CrossRef]
  17. Jung, C.; Ishimoto, R.; Tsuboi, H.; Koyama, M.; Endou, A.; Kubo, M.; Del Carpio, C.A.; Miyamoto, A. Interfacial properties of ZrO2 supported precious metal catalysts: A density functional study. Appl. Catal. A Gen. 2006, 305, 102–109. [Google Scholar] [CrossRef]
  18. Jung, K.T.; Bell, A.T. The effects of synthesis and pretreatment conditions on the bulk structure and surface properties of zirconia. J. Mol. Catal. A Chem. 2000, 163, 27–42. [Google Scholar] [CrossRef]
  19. Camellone, M.F.; Ribeiro, F.N.; Szabová, L.; Tateyama, Y.; Fabris, S. Catalytic Proton Dynamics at the Water/Solid Interface of Ceria-Supported Pt Clusters. J. Am. Chem. Soc. 2016, 138, 11560–11567. [Google Scholar] [CrossRef]
  20. Bu, Y.; Cui, T.; Zhao, M.; Zheng, W.; Gao, W.; Jiang, Q. Evolution of Water Structures on Stepped Platinum Surfaces. J. Phys. Chem. C 2018, 122, 604–611. [Google Scholar] [CrossRef]
  21. Zhang, X.; Hu, X.; Lei, H.; Hu, J.; Zhang, Y. Mechanical force-induced polymerization and depolymerization of F-actin at water/solid interfaces. Nanoscale 2016, 8, 6008–6013. [Google Scholar] [CrossRef] [PubMed]
  22. Huang, J.; Cheng, F.; Binks, B.P.; Yang, H. pH-responsive gas–water–solid interface for multiphase catalysis. J. Am. Chem. Soc. 2015, 137, 15015–15025. [Google Scholar] [CrossRef] [PubMed]
  23. Shaat, M. Viscosity of Water Interfaces with Hydrophobic Nanopores: Application to Water Flow in Carbon Nanotubes. Langmuir 2017, 33, 12814–12819. [Google Scholar] [CrossRef] [PubMed]
  24. Yali, Y.; Chunhai, L.; Juan, H.; Yi, L.; Wenkai, C. Density functional theory study of H2O adsorption and decomposition. Chin. J. Catal. 2009, 30, 328. [Google Scholar]
  25. Iskandarova, I.; Knizhnik, A.; Rykova, E.; Bagatur’Yants, A.; Potapkin, B.; Korkin, A. First-principle investigation of the hydroxylation of zirconia and hafnia surfaces. Microelectron. Eng. 2003, 69, 587–593. [Google Scholar] [CrossRef]
  26. Korhonen, S.T.; Calatayud, M.; Krause, A.O.I. Stability of Hydroxylated (1‾11) and (1‾01) Surfaces of Monoclinic Zirconia: A Combined Study by DFT and Infrared Spectroscopy. J. Phys. Chem. C 2008, 112, 6469–6476. [Google Scholar] [CrossRef]
  27. Korhonen, S.T.; Calatayud, M.; Krause, A.O.I. Structure and Stability of Formates and Carbonates on Monoclinic Zirconia: A Combined Study by Density Functional Theory and Infrared Spectroscopy. J. Phys. Chem. C 2008, 112, 16096–16102. [Google Scholar] [CrossRef]
  28. Piskorz, W.; Gryboś, J.; Zasada, F.; Cristol, S.; Paul, J.-F.; Adamski, A.; Sojka, Z. Periodic DFT and Atomistic Thermodynamic Modeling of the Surface Hydration Equilibria and Morphology of Monoclinic ZrO2 Nanocrystals. J. Phys. Chem. C 2011, 115, 24274–24286. [Google Scholar] [CrossRef]
  29. Xia, G.-J.; Wang, Y.-G. Dynamic Simulation on Surface Hydration and Dehydration of Monoclinic Zirconia. Chin. J. Chem. Phys. 2022, 35, 629–638. [Google Scholar] [CrossRef]
  30. Cerrato, G.; Bordiga, S.; Barbera, S.; Morterra, C. Surface characterization of monoclinic ZrO2: I. Morphol-ogy, FTIR spectral features, and computer modelling. Appl. Surf. Sci. 1997, 115, 53–65. [Google Scholar] [CrossRef]
  31. Whittle, K.; Lumpkin, G.; Ashbrook, S. Neutron diffraction and MAS NMR of Cesium Tungstate defect pyrochlores. J. Solid State Chem. 2006, 179, 512–521. [Google Scholar] [CrossRef]
  32. Vienna Ab initio Simulation Package (VASP) 6.2. Available online: https://www.vasp.at/ (accessed on 12 October 2021).
  33. Kresse, G.; Furthmüller, J. Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. Phys. Rev. B 1996, 54, 11169–11186. [Google Scholar] [CrossRef] [PubMed]
  34. Kresse, G.; Furthmüller, J. Efficiency of ab-initio total energy calculations for metals and semiconductors using a plane-wave basis set. Comput. Mater. Sci. 1996, 6, 15–50. [Google Scholar] [CrossRef]
  35. Li, L.; Li, Y.; Guo, X.; Zhang, Y.-F.; Chen, W.-K. Adsorption and dissociation of water on HfO2 (111) and (110) surfaces. Acta Phys.-Chim. Sin. 2013, 29, 937–945. [Google Scholar]
Figure 1. Structure of a water molecule adsorbed on two-layer (a), three-layer (b) and four-layer (c) models of the monoclinic ZrO2 (111) surface.
Figure 1. Structure of a water molecule adsorbed on two-layer (a), three-layer (b) and four-layer (c) models of the monoclinic ZrO2 (111) surface.
Nanomaterials 12 04362 g001
Figure 2. Four surface Zr atoms (A, B, C and D) with different coordination environments on the monoclinic ZrO2 (111) surface.
Figure 2. Four surface Zr atoms (A, B, C and D) with different coordination environments on the monoclinic ZrO2 (111) surface.
Nanomaterials 12 04362 g002
Figure 3. Stable adsorption structures of a water molecule at A, B, C and D sites on the monoclinic ZrO2 (111) surface and the corresponding dissociation structures.
Figure 3. Stable adsorption structures of a water molecule at A, B, C and D sites on the monoclinic ZrO2 (111) surface and the corresponding dissociation structures.
Nanomaterials 12 04362 g003
Figure 4. Density of states (DOS) of H2O adsorbed on ZrO2 (111) surfaces and the corresponding dissociative state.
Figure 4. Density of states (DOS) of H2O adsorbed on ZrO2 (111) surfaces and the corresponding dissociative state.
Nanomaterials 12 04362 g004
Figure 5. The three-layer model structure of the water–solid interface reactions of ZrO2.
Figure 5. The three-layer model structure of the water–solid interface reactions of ZrO2.
Nanomaterials 12 04362 g005
Figure 6. Radial distribution functions (RDFs) and their integrated RDFs. g(r)H–H between H atoms in the water molecule (a), g(r)O*–H between O atoms in the water molecule and H atoms in the water molecule (b), g(r)Zr*–O* between Zr atoms on the surface of ZrO2 and O atoms in the water molecule (c) and g(r)O**–H between O atoms on the surface of ZrO2 and H atoms in the water molecule (d).
Figure 6. Radial distribution functions (RDFs) and their integrated RDFs. g(r)H–H between H atoms in the water molecule (a), g(r)O*–H between O atoms in the water molecule and H atoms in the water molecule (b), g(r)Zr*–O* between Zr atoms on the surface of ZrO2 and O atoms in the water molecule (c) and g(r)O**–H between O atoms on the surface of ZrO2 and H atoms in the water molecule (d).
Nanomaterials 12 04362 g006
Figure 7. Evolutions of O–H and Zr–O bonds of the interface reaction between the first layer of water and the ZrO2 surface at site A.
Figure 7. Evolutions of O–H and Zr–O bonds of the interface reaction between the first layer of water and the ZrO2 surface at site A.
Nanomaterials 12 04362 g007
Figure 8. FPMD simulation snapshots at 0 (a), 1000 (b) and 3000 fs (c) of the water–solid interface reaction of ZrO2 at site A.
Figure 8. FPMD simulation snapshots at 0 (a), 1000 (b) and 3000 fs (c) of the water–solid interface reaction of ZrO2 at site A.
Nanomaterials 12 04362 g008
Figure 9. Evolutions of O–H and Zr–O bonds in the interface reaction between the first layer of water and the ZrO2 surface with the five-layer ZrO2 model.
Figure 9. Evolutions of O–H and Zr–O bonds in the interface reaction between the first layer of water and the ZrO2 surface with the five-layer ZrO2 model.
Nanomaterials 12 04362 g009
Figure 10. Evolutions of the bond lengths of O–H and Zr–O of two water molecules and the ZrO2 interface reaction at site C.
Figure 10. Evolutions of the bond lengths of O–H and Zr–O of two water molecules and the ZrO2 interface reaction at site C.
Nanomaterials 12 04362 g010
Figure 11. FPMD simulation snapshots at 0 (a), 500 (b), 1400 (c) and 2200 fs (d) of the water–solid interface reaction of ZrO2 at site C.
Figure 11. FPMD simulation snapshots at 0 (a), 500 (b), 1400 (c) and 2200 fs (d) of the water–solid interface reaction of ZrO2 at site C.
Nanomaterials 12 04362 g011
Table 1. Structural parameters and adsorption energies of a water molecule on different models of the ZrO2 surface.
Table 1. Structural parameters and adsorption energies of a water molecule on different models of the ZrO2 surface.
LayerdO–Zr/(Å)dO–H/(Å)∠HOH/(°)Eads/(kcal/mol)
22.4180.976/0.994103.6−18.3
32.4040.976/0.995103.6−18.5
42.4020.977/0.993103.8−18.7
Table 2. Structural parameters of water adsorption on ZrO2 surfaces.
Table 2. Structural parameters of water adsorption on ZrO2 surfaces.
SitedO–Zr/(Å)dO–H/(Å)∠HOH/(°)Eads/(kcal/mol)
A2.3640.979/1.004104.4−20.8
B2.3940.974/0.998107.8−18.6
C2.3540.977/1.000106.4−19.5
D2.2340.974/1.067107.9−34.0
Table 3. Structural parameters and related energies of water dissociation on ZrO2 surfaces.
Table 3. Structural parameters and related energies of water dissociation on ZrO2 surfaces.
SitedO–Zr/(Å)dO–H/(Å)∠HOH/(°)Eads/(kcal/mol)
A2.0800.972/1.642109.9−19.3
B2.0850.969/1.691136.5−16.5
C2.0730.973/1.774106.7−12.3
D2.0910.972/1.465110.8−33.9
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Xu, R.; Zhou, Z.; Wang, Y.; Xiao, H.; Xu, L.; Ding, Y.; Li, X.; Li, A.; Fang, G. First-Principles Molecular Dynamics Simulations on Water–Solid Interface Behavior of H2O-Based Atomic Layer Deposition of Zirconium Dioxide. Nanomaterials 2022, 12, 4362. https://doi.org/10.3390/nano12244362

AMA Style

Xu R, Zhou Z, Wang Y, Xiao H, Xu L, Ding Y, Li X, Li A, Fang G. First-Principles Molecular Dynamics Simulations on Water–Solid Interface Behavior of H2O-Based Atomic Layer Deposition of Zirconium Dioxide. Nanomaterials. 2022; 12(24):4362. https://doi.org/10.3390/nano12244362

Chicago/Turabian Style

Xu, Rui, Zhongchao Zhou, Yingying Wang, Hongping Xiao, Lina Xu, Yihong Ding, Xinhua Li, Aidong Li, and Guoyong Fang. 2022. "First-Principles Molecular Dynamics Simulations on Water–Solid Interface Behavior of H2O-Based Atomic Layer Deposition of Zirconium Dioxide" Nanomaterials 12, no. 24: 4362. https://doi.org/10.3390/nano12244362

APA Style

Xu, R., Zhou, Z., Wang, Y., Xiao, H., Xu, L., Ding, Y., Li, X., Li, A., & Fang, G. (2022). First-Principles Molecular Dynamics Simulations on Water–Solid Interface Behavior of H2O-Based Atomic Layer Deposition of Zirconium Dioxide. Nanomaterials, 12(24), 4362. https://doi.org/10.3390/nano12244362

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop