Abstract
Metal oxide nanoparticles (NPs) are regarded as good candidates for many technological applications, where their functional environment is often an aqueous solution. The correct description of metal oxide electronic structure is still a challenge for local and semilocal density functionals, whereas hybrid functional methods provide an improved description, and local atomic function-based codes such as CRYSTAL17 outperform plane wave codes when it comes to hybrid functional calculations. However, the computational cost of hybrids are still prohibitive for systems of real sizes, in a real environment. Therefore, we here present and critically assess the accuracy of our electrostatic embedding quantum mechanical/molecular mechanical (QM/MM) coupling between CRYSTAL17 and AMBER16, and demonstrate some of its capabilities via the case study of TiO2 NPs in water. First, we produced new Lennard–Jones (LJ) parameters that improve the accuracy of water–water interactions in the B3LYP/TIP3P coupling. We found that optimizing LJ parameters based on water tri- to deca-mer clusters provides a less overstructured QM/MM liquid water description than when fitting LJ parameters only based on the water dimer. Then, we applied our QM/MM coupling methodology to describe the interaction of a 1 nm wide multilayer of water surrounding a spherical TiO2 nanoparticle (NP). Optimizing the QM/MM water–water parameters was found to have little to no effect on the local NP properties, which provide insights into the range of influence that can be attributed to the LJ term in the QM/MM coupling. The effect of adding additional water in an MM fashion on the geometry optimized nanoparticle structure is small, but more evident effects are seen in its electronic properties. We also show that there is good transferability of existing QM/MM LJ parameters for organic molecules–water interactions to our QM/MM implementation, even though these parameters were obtained with a different QM code and QM/MM implementation, but with the same functional.
Highlights
Some of the authors involved in this work have previously demonstrated how electrostatic embedding quantum mechanical/molecular mechanical (QM/molecular mechanic (MM)) can be used to predict solvation dynamics of bond formation in a model catalyst system, which was later confirmed by experiment [9,10]
Since the method is focused on describing molecules and larger systems such as nanoparticles in aqueous solution, we start by focusing on neat water
We have employed two strategies for optimizing QM/MM LJ parameters: The Dimer Fit (DF) is based on optimizing the LJ parameters to minimize the difference between the multiscale interaction energy curves of the water dimer and the average of the pure QM and pure MM interaction energy curves, a methodology akin to previous work on other systems [33]
Summary
Some of the authors involved in this work have previously demonstrated how electrostatic embedding QM/MM can be used to predict solvation dynamics of bond formation in a model catalyst system, which was later confirmed by experiment [9,10]. The MPP-CRYSTAL17 [12,13,14] code has an efficient implementation of the exact exchange term [12]. This capability is crucial for successfully modeling semiconducting oxide systems, such as TiO2 nanoparticles (NPs) in aqueous solution [15]
Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have