Fredrik
Grote
and
Alexander P.
Lyubartsev
*
Department of Materials and Environmental Chemistry, Stockholm University, Svante Arrhenius väg 16 C, 106 91 Stockholm, Sweden. E-mail: alexander.lyubartsev@mmk.su.se
First published on 3rd October 2022
Water structure, dynamics and reactivity at the surface of a small TiO2-nanoparticle fully immersed in water was investigated by an ab initio molecular dynamics simulation. Several modes of water binding were identified by assigning each atom to an atom type, representing a distinct chemical environment in the ab initio ensemble, and then computing radial distribution functions between the atom types. Surface reactivity was investigated by monitoring how populations of atom types change during the simulation. In order to acquire further insight, electron densities for a set of representative system snapshots were analyzed using an atoms-in-molecules approach. Our results reveal that water dissociation, where a water molecule splits at a bridging oxygen site to form a hydroxyl group and a protonated oxygen bridge, can occur by a mechanism involving transfer of a proton over several water molecules. The hydroxyl group and protonated oxygen bridge formed in the process persist (on a 10 ps time scale) and the hydroxyl group undergoes exchange using a mechanism similar to the one responsible for water dissociation. Rotational and translational dynamics of water molecules around the nanoparticle were analyzed in terms of reorientational time correlation functions and mean square displacement. While reorientation of water O–H vectors decreases quickly in the proximity of the nanoparticle surface, translational diffusion slows down more gradually. Our results give new insight into water structure, dynamics and reactivity on TiO2-nanoparticle surfaces and suggest that water dissociation on curved TiO2-nanoparticle surfaces can occur via more complex mechanisms than those previously identified for flat defect-free surfaces.
Recently, properties of nanoclusters, that is small nanoparticles under 2 nm in size, have become a hot topic of research due to the possibility of controlling a nanocluster structure at up to atomistic resolutions.10 The photocatalytic efficiency of TiO2 can be improved by nanostructuring where the size and shape of nanoclusters are controlled in order to optimize their catalytic properties.11–13 Here, it is important to better understand the reactivity of small TiO2-clusters with well-defined structures. Despite the significant progress in nanotechnology during the last decades, developing effective and green catalysts for energy applications and reducing the hazards of nanoparticles to humans and the environment remain central challenges for nanoscience. Since many of these processes take place in aqueous media, a detailed understanding of nanoparticle–water interactions is essential. Here, water dissociation and proton transfers modifying the nanoparticle surface and its interaction with the environment play an important role.
Experimentally, water reactivity on TiO2-surfaces has been studied using X-ray photoelectron spectroscopy14 (XPS) and X-ray absorption spectroscopy15 (XAS), which have provided valuable insights into electronic structures and adsorption thermodynamics. The TiO2–water interface has also been studied using nuclear magnetic resonance (NMR) spectroscopy16 providing information about the structure and mobility of water around TiO2-nanoparticles. These results have suggested that water molecules are organized in layers where the mobility of water molecules decreases with their proximity to the nanoparticle surface. Furthermore, Wang et al.17 have conducted in situ molecular beam scattering STM studies, which provided a detailed picture of the energetics governing water dissociation on the rutile (110) surface. Their results provide evidence of water dissociation at bridging oxygen sites (2-coordinated oxygen atoms) leading to the formation of protonated oxygen bridges and hydroxyl groups. They reported that the energy difference between molecularly adsorbed water and dissociated water is on the order of kBT at room temperature and that the energy barrier for dissociation is roughly 10kBT. This suggests that water dissociation on this surface occurs even at ambient conditions.
TiO2–water interactions have also been studied using computational techniques, including density functional theory (DFT),18 which emphasize the role of the local surface structure, in particular structural defects including step edges and oxygen vacancies. Classical molecular dynamics (MD) simulations, using standard molecular-mechanical force fields, have provided additional insights into water structures and dynamics at TiO2-surfaces.19–22 These simulations are, however, limited by the quality of force fields describing atomic interactions, which do not allow for the breaking or formation of chemical bonds as well as the absence of polarizability due to the use of a fixed point charge representation of the charge distribution. Huang et al. has partially resolved the first of these problems by carrying out MD simulations of TiO2 and water using a reactive force field23 allowing chemical bonds to break and form. Although this approach is very attractive from a computational point of view, it is also highly dependent on the parametrization of the interaction potential. Also, proton transfer reactions at the anatase (101) aqueous interface were probed by MD simulations with an ab initio trained neural network model.24 With the continual increase of computer power, ab initio molecular dynamics simulations25 (AIMD) have gained popularity. Here, the force acting on each atom is obtained by solving the electronic structure problem at each MD step making it possible for chemical bonds to break and form and allowing the electronic cloud to deform in response to electric fields. In practice, however, this approach is currently limited to system sizes of up to 1000 atoms and a 10 ps time scale. The interaction between water and a number of flat defect-free TiO2-surfaces has been studied using AIMD simulations.26 These simulations gave important mechanistic insights into water structure and water dissociation on the investigated surfaces showing that water molecules split via a number of mechanisms involving one or two water molecules.
It can be suggested that a real nanoparticle represents some mixture of different surfaces with various defects, i.e. at edges between the surfaces. Defects can be associated with certain effects that are not captured by the planar defect-free geometry of slab systems. On a curved surface, water molecules can arrange in ways that are not possible on a planar surface due to the absence of curvature. Surface atoms at the place of defects are often more under-coordinated compared to atoms at an ideally planar surface. Furthermore, reactive events often take place on edges, which are not present in a periodically repeated defect-free slab. Attempts to consider these effects were made by ab initio simulations of smaller TiO2-clusters in the presence of water molecules, focusing on wetting thermodynamics.27 However, the size of these systems was rather limited (nanoparticle diameter of ∼1 nm), and the number of water molecules nH2O ≤ 30 corresponded to no more than one layer of water at the nanoparticle surface. To access effects of an irregular surface on water hydration, in this work we have carried out an AIMD simulation of a larger nanoparticle (∼1.5 nm) fully immersed in water in order to gain detailed atomistic insights into water structure, dynamics and reactivity on TiO2-surfaces that are different from the planar defect-free systems that have been studied previously. Recently the TiO2–water interface has also attracted interest in the context of the decoupling of translational and rotational diffusion similar to the anomalies that have been observed in supercooled water.28 This also makes these systems interesting from a theoretical point of view. We have therefore carried out an extensive analysis of translational and reorientational dynamics of water molecules as a function of the distance to the nanoparticle surface. The next section, “Methods and models”, describes details of the simulated system and the analysis of the obtained trajectory and electron density. In the “Results and discussion”, we present and discuss the results from this analysis. The last section, “Conclusion”, provides our concluding remarks.
Equilibration was monitored by following how temperature, potential energy, root mean square displacement (RMSD) of the nanoparticle atoms, and radial distribution functions (RDF) were changing as functions of time (Section S2, Fig. S3 and S6†). As shown in the ESI (Fig. S3†), temperature and potential energy reached stable levels within 1 ps. The RMSD of the nanoparticle atoms (Fig. S4†), obtained with respect to the original crystalline structure for both surface and interior atoms, shows a slightly increasing trend but remains low, indicating that no structural transformation occurred, and the structure remains close to the anatase phase. Similar conclusions were made from previous ab initio simulations of weakly hydrated TiO2-nanoparticles.27 RDF computed between Ti and O atoms showed a stable shape after 1 ps of the simulations (Fig. S5†). We also computed RDF between water hydrogen and surface oxygen atoms (O2c) during different time windows (Fig. S6†), which confirmed that the structure of the interface reached equilibrium. This short equilibration time can be explained by the fact that the starting structure of the AIMD simulation was obtained as a result of a 50 ns simulation with a classical force field, which provided an almost equilibrium structure needing only a short time to adjust to ab initio interactions. In the analysis of the simulation, we excluded the first 3 ps of the trajectory.
Each atom in the system can be mapped onto one of these atom types at every point during the simulation by determining the coordination number of atoms by counting other atoms that fall within a certain cut-off radius. However, the atom types determined by such geometrical criteria will frequently change due to fast bond oscillations and will not always represent actual chemical transformations from one type to another. In order for the change in atom types to reflect the change of local environment due to reactive events, we must average out these fast bond vibrations. For this reason, we assign atom types in 100 fs windows to determine coordination numbers of atoms by averaging the number of neighboring atoms found within a sphere of radius R, which represents the maximal distance where atoms can still be considered as bonded to each other, here taken to be the first minimum of the respective RDF. For Ti–O we use R = 3.0 Å and for O–H we take R = 1.3 Å. Populations of atom types obtained using this algorithm for the start configuration (generated by a simulation with the classical force field) and averaged over the whole trajectory of the AIMD simulation are presented in Table 1. Note that the original (non-hydrated) nanoparticle contained 21, 12 and 20 of Ti6c, Ti5c and Ti4c atoms, respectively, but upon hydration and adsorption of water molecules during a classical MD simulation the number of 4- and 5-coordinated Ti atoms was greatly reduced, and in particular no 4-coordinated Ti atoms remained. We observed that during almost 50 ps of the AIMD simulation only small changes in the populations take place: a small increase of Ti5c (5-coordinated) Ti atoms, and the appearance of small amounts of OH (hydroxyl) and O2cH (protonated bridge) types of oxygen atoms due to the splitting of water molecules (discussed in more detail below). This situation is different from that of the previous AIMD simulation of a TiO2 anatase nanoparticle at low hydration where about 30% of water molecules were split to form surface hydroxyl groups,27 and also from simulations of plain anatase TiO2 slabs, which showed no water splitting or reactivity at some surfaces ((100) and (101)) while there was a high level of reactivity at the (001) surface where hydroxyl groups were formed at about one-third of all surface Ti sites.26 Experimentally, high numbers of OH– groups were found on the anatase (101) surface at low hydration39 while sum frequency generation spectrum measurements and AIMD simulations indicated a stable bilayer of intact water molecules on the anatase (101) surface at full hydration.40
Atom type | q NP (e) | q slab (e) | N start | N mean |
---|---|---|---|---|
Ti6c | 2.16 ± 0.05 | — | 50 | 48.69 |
Ti5c | 2.08 ± 0.04 | 2.2 | 3 | 4.31 |
O3c | −1.15 ± 0.03 | — | 58 | 58.00 |
O2c | −1.06 ± 0.03 | −1.0 | 46 | 45.53 |
Ow | −0.75 ± 0.03 | −0.753 | 377 | 378.30 |
Ow1c | −0.81 ± 0.03 | — | 49 | 47.23 |
OH | −0.91 ± 0.03 | — | 0 | 0.46 |
O2cH | −1.00 ± 0.02 | — | 0 | 0.47 |
H | 0.38 ± 0.02 | 0.43 | 852 | 852.00 |
In order to feel further confidence that this choice of atom types is physically motivated, we carried out an atoms-in-molecules analysis where a set of snapshots from the trajectory were extracted and from their electron densities net atomic charges were computed using the DDEC6 method.41 The result presented in Table 1 shows that the difference of net atomic charges between atom types is larger or about the same as the spread within the types, which is also seen in the distributions of net atomic charges presented in Fig. S7 in Section S3 of the ESI.† Furthermore, the charges computed for the nanoparticle are very similar to previously computed net charges for a planar surface of the same atom types. This suggests that the atom types to which we have assigned atoms represent distinct chemical environments, and they are transferable over different types of surfaces. These atom types can be useful for discussions of surface chemistry and other processes taking place on TiO2-nanoparticle surfaces including biomolecular adsorption and catalysis. They also form a natural choice of atom types for use in classical MD simulations, which is needed in order to simulate longer lengths and time scales. In addition, the distribution of atom types presented here can provide important guidance for setting up starting configurations for such simulations.
![]() | ||
Fig. 2 Water density distribution as a function of distance from the nanoparticle surface with water layers L1, L2, L3a, L3b and L4. |
While the water density distribution gives an overall picture of the structure of the TiO2–water interface, a more detailed description can be obtained by analyzing correlations between individual atom types. RDFs are important functions for characterization of structure, in particular for systems where there is no long-range order and where only local correlations occur. In order to obtain detailed structural information about the simulated system, we calculated RDFs between the atom types defined in the previous section. Of particular interest are RDFs between atom types belonging to the nanoparticle and water since these contain valuable information about the structure of the TiO2–water interface. Fig. 3 shows RDFs between different oxygen and hydrogen atom types. The RDF between oxygen in bulk water (Ow type) and hydrogen, as shown in Fig. 3(a), has three peaks typical for O–H RDF in pure water. A high and narrow peak appears just above 1 Å and corresponds to the covalent O–H bond, with the broader one, centered at just below 2 Å, being assigned to the hydrogen of a hydrogen-bonded water molecule, while the third peak at roughly 3.5 Å corresponds to another hydrogen bonding of that molecule to a second solvation shell. The RDF in Fig. 3(b) between oxygen in molecularly adsorbed water (Owlc type) and hydrogen atoms shows 3 peaks in similar positions but of different magnitudes. In particular, the peak at roughly 2 Å, corresponding to hydrogen bonding, is significantly reduced due to coordination of oxygen to the nanoparticle surface, which hinders hydrogen bonding to another water molecule. Fig. 3(c) shows the RDF between bridging (2-coordinated) oxygen atoms and hydrogen. This RDF has a well-expressed hydrogen bond peak just below 2 Å and a second broader peak at about 2.5 Å. It can be noted that correlations between bridging oxygens and hydrogen persist for longer when compared with those between oxygen atoms in bulk water and hydrogen. This can be attributed to an ordering effect of the rigid nanoparticle structure. Fig. 3(d) shows the RDF between 3-coordinated oxygen atoms and hydrogen. In this graph the characteristic peak due to hydrogen bonding is completely absent and the RDF has complex structure beyond 3 Å since O3c atoms are found throughout the interior of the nanoparticle at different distances from the surface. This shows that only bridging oxygen atoms are involved in hydrogen bonding while 3-coordinated oxygen atoms are either found in the interior of the nanoparticle or in buried surface sites where formation of hydrogen bonds is not possible. Consistently, RDFs between oxygen atom types O3c–OH and O2c–OH presented in Fig. S8 in Section S4 of the ESI† show a pronounced hydrogen bond peak at about 3 Å for O2c, which is absent for O3c. In Fig. 3(e), the RDF between oxygen as a protonated bridging oxygen (O2cH) and hydrogen is shown. This RDF is essentially identical to that for the Ow–H bond but is noisier since there is only one atom belonging to atom type O2cH. Fig. 3(f) shows the RDF between oxygen in a hydroxyl group (OH) and hydrogen. This RDF shows that the highest hydrogen bond peak is attributed to strong hydrogen bonding between the hydroxyl oxygen and a hydrogen in water.
![]() | ||
Fig. 3 Oxygen–hydrogen RDFs between atom types (a) Ow–H, (b) Owlc–H, (c) O2c–H, (d) O3c–H, (e) O2cH–H, and (f) OH–H. |
Fig. S9 in Section S4 of the ESI† show RDFs between titanium atom types and the oxygen atom of water. The RDFs for Ti6c–Owlc and Ti5c–Owlc bonds are presented in Fig. S9(a) and S9(b)†. These RDFs have a high peak just above 2 Å corresponding to the bond between titanium atoms and oxygen atoms in molecularly adsorbed water molecules. The peaks occurring at longer distances correspond to titanium atoms and oxygen in molecularly adsorbed water separated by several bonds. Fig. S9(c) and S9(d)† show the corresponding RDFs between titanium atom types and oxygen atoms in water molecules in the bulk nanoparticle. Here the first maximum occurs at 4 Å since water molecules located closer than this distance are those which are adsorbed to the surface. A third shell of water molecules gives rise to a maximum at roughly 7 Å showing that structural correlations in the water surrounding the nanoparticle persist out to at least this distance.
RDFs between titanium atom and oxygen atom types in the nanoparticle are shown in Fig. S10 in Section S4 of the ESI.† These RDFs all have a sharp peak at about 2 Å corresponding to the bond between titanium and oxygen. At roughly 4 Å a broader peak occurs, which is due to titanium–oxygen neighbors that are separated by more than one bond.
The RDFs presented here give a detailed insight into the structure of the TiO2-nanoparticle and the TiO2–water interface and can be important target data in parametrization of classical force fields for these systems.
As illustrated in Fig. 4a these events take place on an edge/corner of the nanoparticle, which is shown in close up in Fig. 4b together with three coordinated water molecules that play key roles in the set of reactions that lead to the splitting of a water molecule. We now turn to a more detailed analysis of how water dissociates on the surface. The atoms involved in the reaction are shown in Fig. 4c where we have also marked two distances that will be used in our discussion. Fig. 5 shows the O75–H948 and O430–H431 distances as functions of time. In order to illustrate the mechanism by which water dissociation occurs, a number of snapshots have been extracted along the reaction pathway. In snapshot 1, three water molecules are arranged in a hydrogen-bonded constellation around a bridging oxygen atom. The water molecule closest to the oxygen bridge then shares one of its protons with the bridging oxygen, as shown in snapshot 2, forming the configuration in snapshot 3. This structure is not stable, and the proton is quickly captured back by the water molecule (snapshot 4). However, the configuration in snapshot 3 is then formed again, as shown in snapshot 5. Then a synchronous proton transfer over three water molecules occurs leading to the dissociation of a water molecule and the formation of a hydroxyl group and a protonated oxygen bridge. As illustrated in snapshot 6, a water molecule adsorbed to a Ti-atom donates a proton to a neighboring water molecule, which donates one of its proton to the water molecule at the oxygen bridge, which in turn donates one of its protons to the oxygen bridge and then diffuses out into the solution. This mechanism involves three water molecules, where one is bridge protonating (BP-water), the second is hydroxyl group forming (HF-water), and the third is proton mediating (PM-water) and transfers a proton between the BP- and the HF-water. This result suggests that water dissociation on curved TiO2 surfaces can occur by mechanisms involving complex proton transfers over several molecules.
![]() | ||
Fig. 5 Distances between atoms O75–H948 and O430–H431 as a function of time during water dissociation. |
Although interatomic distances are easy to interpret and provide important information about breaking and formation of chemical bonds, it should be pointed out that the spatial proximity of two atoms does not unambiguously quantify the amount of covalent interaction between them. For this reason, we have repeated the analysis presented here using bond orders, computed from the electron density of the snapshots in Fig. 5, instead of interatomic distances. Bond orders provide information about the amount of shared electron density between atoms due to ‘dressed exchange’ and are considered good indicators of the extent of covalent interactions between atoms. In the ESI (ection S5, Fig. S11†), we show that the analysis using bond orders leads to the same results as the approach presented here.
In previous works,42 proton transfer was described using two coordinates: (1) the donor–acceptor distance and (2) δ = RA–H − RD–H, where RA–H is the distance between the acceptor atom and the proton and RD–H is the distance between the donor atom and the proton. In order to obtain further insight, we computed the distribution of these coordinates during water dissociation. The free energy profiles along δ for short and intermediate donor–acceptor distances have been reported to be symmetric double-well separated by a barrier. As the donor–acceptor distance decreases the barrier becomes smaller facilitating proton transfer between the donor and acceptor. This picture is consistent with the distribution shown in Fig. S12 in Section S6 of the ESI† being nearly symmetric about δ = 0 with two maxima corresponding to the proton bonded to a donor and an acceptor, respectively. From this graph, it is also evident that the proton is transferred when the donor–acceptor distance is short, i.e., 2.5–2.6 Å.
Another observed event involving the rearranging of chemical bonds was hydroxyl group exchange. We carried out a similar analysis of atomic distances, defined in Fig. 4d, as was done for the water dissociation. Fig. 6 shows the O430–Ti11 and O439–H720 distances as functions of time and a set of snapshots extracted along the reaction pathway. Snapshot 1 shows a configuration with a water molecule adsorbed to a Ti-atom with a hydroxyl group and between them there is a second water hydrogen-bonded to the hydroxyl group. Because of thermal fluctuations, the second water molecule can donate its proton to the hydroxyl group (see snapshot 2), and simultaneously accept a proton (as shown in snapshot 3) from the adsorbed water, which then forms a hydroxyl group. This results in the structure shown in snapshot 4 where the hydroxyl group has been exchanged. This suggests that hydroxyl group exchange can occur by a mechanism similar to the one presented for water dissociation. One water protonates the hydroxyl group, while another water mediates a proton from an adsorbed water molecule. Just as in the mechanism for water dissociation this requires several water molecules to come together in a specific arrangement where proton transfer occurs. This suggests that the local surface structure (i.e. curvature, edges and defects) is of considerable importance for the surface chemistry of these materials. Our results show that water dissociation on curved TiO2-surfaces with edges can be more complex than water dissociation previously discovered for flat surfaces.
![]() | ||
Fig. 6 Distances between atoms O75–H948 and O430–H431 as a function of time during hydroxyl group exchange. |
C(τ) = 〈P2[û(t)·û(t + τ)]〉 | (1) |
In order to analyze water diffusion near the TiO2 surface we computed mean square displacement (MSD) of water in layers L1–L4. MSD was computed for water oxygen atoms present in the region both at initial time t and later time t + τ. The resulting MSD curves (see Fig. 7(b)) show that the mobility of water molecules is strongly reduced as they come into the proximity of the surface. In contrast to the reorientational dynamics, which decreases abruptly at R = 3 Å, it is observed that a decrease of translational diffusion starts at longer distances from the surface. Fig. 7(b) also shows the total MSD (computed over all water molecules), which increases faster than MSD-curves for R = 2–5 Å due to interlayer diffusion and a contribution from water molecules located at R > 5 Å. Diffusion coefficients were computed by fitting straight lines to the MSD plots presented in Fig. 7(b) in the range where the MSD is linear (4 ps ≤ τ ≤ 20 ps) according to eqn (2):
![]() | (2) |
We address further the question of a possible decoupling of translational and rotational diffusion. Fig. 7(c) shows the TCF of the water O–H vector as a function of the MSD taken over the same time. From this graph, one observes an indication of how far water molecules diffuse per rotation for each of the shells L1–L4, i.e., a faster decay means that the molecule diffuses over a shorter distance per rotation. While the TCFs for L3 and L4 show a similar decay, it is observed that for L2 TCF decreases more slowly and for L1 it decreases more quickly. This means that the displacement per rotation in layer L2 is larger than it is in L3 and L4, which in turn is larger than it is for L1, thus supporting the decoupling of rotational and translational diffusion in the proximity of the nanoparticle surface.
The motion of water molecules changing hydration shell can also be characterized by survival probability, which is the probability that a water molecule initially residing inside a selected region (shell) remains inside that region after a given time, and is defined according to eqn (3):
![]() | (3) |
Taken together, our results indicate that both translational and rotational dynamics of water molecules substantially decrease with shorter distances to the nanoparticle surface. However, water molecules in the layer closest to the nanoparticle (L1) undergo faster reorientation compared to those in the second hydration shell. It is possible that the extent to which these water molecules participate in hydrogen bonding is related to this effect. The decoupling of rotational and translational diffusion, where translation slows down more than rotation, is observed where water molecules in the second hydration layer (L2) translate further per rotation compared to water molecules in layers L3 and L4. Our analysis also indicates that water reorientation slows down abruptly with decreasing distance to the nanoparticle surface between shells L3 and L2, while a decrease of translational diffusion occurs more gradually. It can be noted that in the inner region where R < 3 Å (i.e., L1 and L2), both rotational and reorientational dynamics are substantially slower compared to those in the bulk. For 3 Å < R < 4 Å (i.e., L3), reorientation is as fast as in the bulk, but diffusion still occurs at a reduced rate. Analysis of survival probability in the various hydration shells showed that water molecules located within 3 Å of the surface undergo slow exchange while water molecules beyond 3 Å are replaced much more quickly.
We have also investigated the dynamics of water around the TiO2-nanoparticle by calculating orientational time correlation functions and the mean square displacement of water molecules at different distances from the nanoparticle surface. We found that water reorientation rates quickly decrease at a distance of around 3 Å from the surface, while diffusion slows down more gradually when starting from longer distances. Our results also indicated that the solvation shells located 2–5 Å from the surface consist of two regions: (1) an inner region with tightly held water molecules having slow reorientational and translational dynamics; and (2) an outer region where reorientation of water molecules is as fast as in the bulk nanoparticle while diffusion occurs at a reduced rate. Overall, local atom environment analysis, and RDFs and MSDs calculated in this work give a detailed picture of the structure, dynamics and reactivity at a curved TiO2–water interface. Furthermore, the computed water properties are important target data for parametrization and validation of classical force fields allowing for extended simulations of the aqueous interface with TiO2-nanoparticles over much longer time and length scales.
Footnote |
† Electronic supplementary information (ESI) available: Supplementary method details and supplementary figures (.pdf). See DOI: https://doi.org/10.1039/d2nr02354g |
This journal is © The Royal Society of Chemistry 2022 |