Search

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. 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.

Computational modeling of titanium dioxide nanoparticles of realistic size is extremely relevant for the direct comparison with experiments but it is also a rather demanding task. We have recently worked on a multistep/scale procedure to obtain global optimized minimum structures for chemically stable spherical titania nanoparticles of increasing size, with diameter from 1.5 nm (~300 atoms) to 4.4 nm (~4000 atoms). We use first self-consistent-charge density functional tight-binding (SCC-DFTB) methodology to perform thermal annealing simulations to obtain globally optimized structures and then hybrid density functional theory (DFT) to refine them and to achieve high accuracy in the description of structural and electronic properties. This allows also to assess SCC-DFTB performance in comparison with DFT(B3LYP) results. As a further step, we investigate photoexcitation and photoemission processes involving electron/hole pair formation, separation, trapping and recombination in the nanosphere of medium size by hybrid DFT. Finally, we show how a recently defined new set of parameters for SCC-DFTB allows for a proper description of titania/water multilayers interface, which paves the way for modeling large realistic nanoparticles in aqueous environment.

TiO2 nanoparticles (NPs) are nowadays considered fundamental building blocks for many technological applications. Morphology is found to play a key role with spherical NPs presenting higher binding properties and chemical activity. From the experimental point of view, the characterization of these nano-objects is extremely complex, opening a large room for computational investigations. In this work, TiO2 spherical NPs of different size (from 300 to 4000 atoms) have been studied with a two-scale computational approach. Global optimization to obtain stable and equilibrated NSs was performed with a self-consistent charge density functional tight-binding (SCC-DFTB) simulated annealing process, causing a considerable atomic rearrangement within the nanospheres. Those SCC-DFTB relaxed structures have been then optimized at DFT(B3LYP) level of theory. We present a systematic and comparative SCC-DFTB vs DFT(B3LYP) study of the structural properties, with particular emphasis on the surface-to-bulk sites ratio, coordination distribution of surface sites and surface energy. From the electronic point of view, we compare HOMO-LUMO and Kohn-Sham gaps, total and projected density of states. Overall, the comparisons between DFTB and hybrid DFT show that DFTB provides a rather accurate geometrical and electronic description of these nanospheres of realistic size (up to a diameter of 4.4 nm) at an extremely reduced computational cost. This opens for new challenges in simulations of very large systems and more extended molecular dynamics.

Nanoparticles of very small size (below 10 nm) of TiO2 material are nowadays the functional building blocks of many developing technological applications. Nano is clearly different from bulk or extended systems as regards surface area, molecular binding properties, charge separation efficiency, electron/hole transport, photochemical conversion properties, etc. In this work, we investigate the life path of energy (excitons) and charge (electrons and holes) carriers in anatase TiO2 nanoparticles of different size (2-3 nm) and shape (faceted vs spherical), by means of a wide set of hybrid density functional theory calculations. The attention is focused on the exciton/charge carriers formation, separation, recombination, self-trapping processes, which are analyzed in terms of structural deformations, energy gain or cost, charge localization/delocalization and electronic transitions involved. The computational models are corroborated by an extensive comparison with available experimental data based on photoluminescence measurements, electron paramagnetic resonance and transient absorption spectroscopies. Peculiar differences are observed for spherical nanoparticles with respect to faceted ones because of the higher disorder and larger diversity of coordination sites present on the surface. For example, charge delocalization on several lattice sites is more competitive with self-trapping processes in faceted than in spherical nanoparticles. This relates to the fact that selective compression or elongation of Ti-O bonds play a key role in determining the effectiveness of trapping sites, with spherical nanoparticles being more flexible. Moreover, hydroxyl groups on surface five-fold coordinated Ti sites are also found to be good hole trapping sites.

Boron-doped graphene was reported to be the best non-metal doped graphene electrocatalyst for the oxygen reduction reaction (ORR) working at an onset potential of 0.035 V [JACS 136 (2014) 4394]. In the present DFT study, intermediates and transition structures along the possible reaction pathways are determined. Both Langmuir-Hinschelwood and Eley-Rideal mechanisms are discussed. Molecular oxygen binds the positively charged B atom and forms an open shell end-on dioxygen intermediate. The associative path is favoured with respect to the dissociative one. The free energy diagrams along the four-reduction steps are investigated with the methodology by N{\o}rskov and co. [JPC B 108 (2004) 17886] in both acidic and alkaline conditions. The pH effect on the stability of the intermediates of reduction is analyzed in terms of the Pourbaix diagram. At pH = 14 we compute an onset potential value for the electrochemical ORR of U = 0.05 V, which compares very well with the experimental value in alkaline conditions.