Insight into the interface between Fe3O4 (001) surface and water overlayers through multiscale molecular dynamics simulations

In this work, we investigate the Fe₃O₄ (001) surface/water interface by combining several theoretical approaches, ranging from a hybrid functional method (HSE06) to density-functional tight-binding (DFTB) to molecular mechanics (MM). First, we assess the accuracy of the DFTB method to correctly reproduce HSE06 results on structural details and energetics and available experimental data for adsorption of isolated water, dimers, and trimers up to a water monolayer. Second, we build two possible configurations of a second and a third overlayer and perform molecular dynamics simulations with DFTB, monitoring the water orientation, the H-bond network, and the ordered water structure formation. To make our models more realistic, we then build a 12 nm-thick water multilayer on top of the Fe₃O₄ (001) surface slab model, which we investigate through MM-molecular dynamics (MD). The water layer structuring, revealed by the analysis of the atomic positions from a long MM-MD run for this large MM model, extends up to about 6-7 Å and nicely compares with that observed for a water trilayer model. However, MM and DFTB MD simulations show some discrepancy due to the poor description of the Fe⋯OH₂ distance in MM that calls for further work in the parameterization of the model.

Shaping Magnetite Nanoparticles from First Principles

Iron oxide magnetic nanoparticles (NPs) are stimuli-responsive materials at the forefront of nanomedicine. Their realistic finite temperature simulations are a formidable challenge for first-principles methods. Here, we use density functional tight binding to open up the required time and length scales and obtain global minimum structures of Fe_{3}O_{4} NPs of realistic size (1400 atoms, 2.5 nm) and of different shapes, which we then refine with hybrid density functional theory methods to accomplish proper electronic and magnetic properties, which have never been accurately described in simulations. On this basis, we develop a general empirical formula and prove its predictive power for the evaluation of the total magnetic moment of Fe_{3}O_{4} NPs. By converting the total magnetic moment into the macroscopic saturation magnetization, we rationalize the experimentally observed dependence with shape. We also reveal interesting reconstruction mechanisms and unexpected patterns of charge ordering.

Optimizing PEGylation of TiO2 Nanocrystals through a Combined Experimental and Computational Study

PEGylation of metal oxide nanoparticles is the common approach to improve their biocompatibility and in vivo circulation time. In this work, we present a combined experimental and theoretical study to determine the operating condition that guarantee very high grafting densities, which are desirable in any biomedical application. Moreover, we present an insightful conformational analysis spanning different coverage regimes and increasing polymer chain lengths. Based on 13C NMR measurements and molecular dynamics simulations, we show that classical and popular models of polymer conformation on surfaces fail in determining the mushroom-to-brush transition point and prove that it actually takes place only at rather high grafting density values.

Rational design of nanosystems for simultaneous drug delivery and photodynamic therapy by quantum mechanical modeling

Drug delivery systems are based on reversible interactions between carriers and drugs. Spacers are often introduced to tailor the type of interaction and to keep drugs intact. Here, we model a drug delivery system based on a functionalized curved TiO₂ nanoparticle of realistic size (700 atoms - 2.2 nm) by the neurotransmitter dopamine to carry the anticancer chemotherapeutic agent doxorubicin (DOX). The multiscale quantum chemical study aims at unraveling the nature and mechanism of the interactions between the components and the electronic properties of the composite system. We simulate the temperature effect through molecular dynamics runs of thermal annealing. Dopamine binds preferentially to low coordinated Ti sites on the nanoparticle through dissociated bidentate and chelate modes involving the diol groups. DOX is tethered by H-bonds, π-π stacking, dipole-dipole interactions and dispersion forces. Comparing different coverage densities of the spacer on the nanoparticle surface, we assess the best conditions for an effective drug transport and release: only at full coverage, DOX does not slip among the dopamine molecules to reach the nanoparticle surface, which is crucial to avoid the formation of stable coordinative bonds with under-coordinated Ti atoms. Finally, given the strong absorption properties and fluorescence of DOX and of the TiO₂ photocatalyst, we model the effect of light irradiation through excited state calculations to localize excitons and to follow the charge carrier's life path. This fundamental study on the nature and mechanism of drug/carrier interaction provides a solid ground for the rational design of new experimental protocols for a more efficient drug transport and release and its combination with photodynamic therapy.

Computational Electrochemistry of Water Oxidation on Metal-Doped and Metal-Supported Defective h-BN

Metal-doped and metal-supported two-dimensional materials are attracting a lot of interest as potentially active electrocatalysts for reduction and oxidation processes. Previously, when a non-regular 2 D h-BN layer was grown on a Cu(111) surface, metal adatoms were found to spontaneously emerge from the bulk to fill the atomic holes in the structure and become available for surface catalysis. Herein, computational electrochemistry is used to investigate and compare the performance of Cu-doped and Cu-supported pristine and defective h-BN systems for the electrocatalytic water oxidation reaction. For the various model systems, the intermediate species of this multistep oxidation process are identified and the free-energy variations for each step of reaction are computed, even for those steps that do not involve an electron or a proton transfer. Both associative and O₂ direct evolution mechanisms are considered. On this thermodynamic basis, the potential-determining step, the thermodynamic-determining step, and consequently the theoretical overpotential are determined for comparison with experiments. Small Cu clusters (tetramers) trapped in the h-BN defective lattice on a Cu(111) support are found to be very active for the water oxidation reaction since such systems are characterized by a low overpotential and by a small energy cost for O₂ release from the catalyst, which is often observed to be a major limit for other potential electrocatalysts.

Proton Transfers at a Dopamine-Functionalized TiO2 Interface

Despite the many successful syntheses and applications of dopamine-functionalized TiO₂ nanohybrids, there has not yet been an atomistic understanding of the interaction of this 1,2-dihydroxybenzene derivative ligand with the titanium dioxide surfaces. In this work, on the basis of a wide set of dispersion-corrected hybrid density functional theory (DFT) calculations and density functional tight binding (DFTB) molecular dynamics simulations, we present a detailed study of the adsorption modes, patterns of growth, and configurations of dopamine on the anatase (101) TiO₂ surface, with reference to the archetype of 1,2-dihydroxybenzene ligands, i.e., catechol. At low coverage, the isolated dopamine molecule prefers to bend toward the surface, coordinating the NH₂ group to a Ti_5c ion. At high coverage, the packed molecules succeed in bending toward the surface only in some monolayer configurations. When they do, we observe a proton transfer from the surface to the ethyl-amino group, forming terminal NH₃ ⁺ species, which highly interact with the O atoms of a neighboring dopamine molecule. This strong Coulombic interaction largely stabilizes the self-assembled monolayer. On the basis of these results, we predict that improving the probability of dopamine molecules being free to bend toward the surface through thermodynamic versus kinetic growth conditions will lead to a monolayer of fully protonated dopamine molecules.

An efficient way to model complex magnetite: Assessment of SCC-DFTB against DFT

Magnetite has attracted increasing attention in recent years due to its promising and diverse applications in biomedicine. Theoretical modelling can play an important role in understanding magnetite-based nanomaterials at the atomic scale for a deeper insight into the experimental observations. However, calculations based on density functional theory (DFT) are too costly for realistically large models of magnetite nanoparticles. Classical force field methods are very fast but lack of precision and of the description of electronic effects. Therefore, a cheap and efficient quantum mechanical simulation method with comparable accuracy to DFT is highly desired. Here, a less computationally demanding DFT-based method, i.e., self-consistent charge density functional tight-binding (SCC-DFTB), is adopted to investigate magnetite bulk and low-index (001) surfaces with newly proposed parameters for Fe-O interactions. We report that SCC-DFTB with on-site Coulomb correction provides results in quantitatively comparable agreement with those obtained by DFT + U and hybrid functional methods. Therefore, SCC-DFTB is valued as an efficient and reliable method for the description magnetite. This assessment will promote SCC-DFTB computational studies on magnetite-based nanostructures that attract increasing attention for medical applications.

Interfacing CRYSTAL/AMBER to Optimize QM/MM Lennard⁻Jones Parameters for Water and to Study Solvation of TiO₂ Nanoparticles

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

Curved TiO2 Nanoparticles in Water: Short (Chemical) and Long (Physical) Range Interfacial Effects

In most technological applications, nanoparticles are immersed in a liquid environment. Understanding nanoparticles/liquid interfacial effects is extremely relevant. This work provides a clear and detailed picture of the type of chemistry and physics taking place at the prototypical TiO₂ nanoparticles/water interface, which is crucial in photocatalysis and photoelectrochemistry. We present a multistep and multiscale investigation based on hybrid density functional theory (DFT), density functional tight-binding, and quantum mechanics/molecular mechanics calculations. We consider increasing water partial pressure conditions from ultra-high vacuum up to the bulk water environment. We first investigate single water molecule adsorption modes on various types of undercoordinated sites present on a realistic curved nanoparticle (2-3 nm) and then, by decorating all the adsorption sites, we study a full water monolayer to identify the degree of water dissociation, the Brønsted-Lowry basicity/acidity of the nanoparticle in water, the interface effect on crystallinity, surface energy, and electronic properties, such as the band gap and work function. Furthermore, we increase the water coverage by adding water multilayers up to a thickness of 1 nm and perform molecular dynamics simulations, which evidence layer structuring and molecular orientation around the curved nanoparticle. Finally, we clarify whether these effects arise as a consequence of the tension at the water drop surface around the nanosphere by simulating a bulk water up to a distance of 3 nm from the oxide surface. We prove that the nanoparticle/water interfacial effects go rather long range since the dipole orientation of water molecules is observed up to a distance of 5 Å, whereas water structuring extends at least up to a distance of 8 Å from the surface.

Bulk-terminated or reconstructed Fe3O4(001) surface: water makes a difference

Surfaces and their interaction with water play an important role in most of materials' applications. Magnetite has attracted continued interest in the fields of catalysis, spintronic devices, magnetic resonance imaging (MRI) and drug delivery. In this work, water adsorption and its effect on the stability diagram and on the electronic structure of the Fe3O4(001) surface are investigated by hybrid density functional theory calculations combined with an ab initio atomistic thermodynamic approach. We span a wide range of gaseous O2 and vapor H2O partial pressures. At low water pressure, a reconstructed SCV surface model is confirmed to be the most stable model at common working O2 partial pressures. However, at high water coverage, an unexpected stability inversion is observed that makes the hydrated bulk-terminated DBT surface the most favored. These results open up new horizons in Fe3O4 surface chemistry when working in an aqueous environment and are of key importance to develop rational strategies to surface engineering for high performance Fe3O4 nanomaterials.

Water on Graphene-Coated TiO2: Role of Atomic Vacancies

Beyond two-dimensional (2D) materials, interfaces between 2D materials and underlying supports or 2D-coated metal or metal oxide nanoparticles exhibit excellent properties and promising applications. The hybrid interface between graphene and anatase TiO₂ shows great importance in photocatalytic, catalytic, and nanomedical applications due to the excellent and complementary properties of the two materials. Water, as a ubiquitous and essential element in practical conditions and in the human body, plays a significant role in the applications of graphene/TiO₂ composites for both electronic devices and nanomedicine. Carbon vacancies, as common defects in chemically prepared graphene, also need to be considered for the application of graphene-based materials. Therefore, the behavior of water on top and at the interface of defective graphene on anatase TiO₂ surface was systematically investigated by dispersion-corrected hybrid density functional calculations. The presence of the substrate only slightly enhances the on-top adsorption and reduces the on-top dissociation of water on defective graphene. However, at the interface, dissociated water is largely preferred compared with undissociated water on bare TiO₂ surface, showing a prominent cover effect. Reduced TiO₂ may further induce oxygen diffusion into the bulk. Our results are helpful to understand how the presence of water in the surrounding environment affects structural and electronic properties of the graphene/TiO₂ interface and thus its application in photocatalysis, electronic devices, and nanomedicine.

Water-Assisted Hole Trapping at the Highly Curved Surface of Nano-TiO2 Photocatalyst

Heterogeneous photocatalysis is vital in solving energy and environmental issues that this society is confronted with. Although photocatalysts are often operated in the presence of water, it has not been yet clarified how the interaction with water itself affects charge dynamics in photocatalysts. Using water-coverage-controlled steady and transient infrared absorption spectroscopy and large-model (∼800 atoms) ab initio calculations, we clarify that water enhances hole trapping at the surface of TiO2 nanospheres but not of well-faceted nanoparticles. This water-assisted effect unique to the nanospheres originates from water adsorption as a ligand at a low-coordinated Ti-OH site or through robust hydrogen bonding directly to the terminal OH at the highly curved nanosphere surface. Thus, the interaction with water at the surface of nanospheres can promote photocatalytic reactions of both oxidation and reduction by elongating photogenerated carrier lifetimes. This morphology-dependent water-assisted effect provides a novel and rational basis for designing and engineering nanophotocatalyst morphology to improve photocatalytic performances.

Synthesis of corrugated C-based nanostructures by Br-corannulene oligomerization

The structure and electronic properties of carbon-based nanostructures obtained by metal surface assisted synthesis is highly dependent on the nature of the precursor molecule.

Band Gap in Magnetite above Verwey Temperature Induced by Symmetry Breaking

Magnetite exhibits a famous phase transition, called Verwey transition, at the critical temperature T_V of about 120 K. Although numerous efforts have been devoted to the understanding of this interesting transition, up to now, it is still under debate whether a charge ordering and a band gap exist in magnetite above T_V. Here, we systematically investigate the charge ordering and the electronic properties of magnetite in its cubic phase using different methods based on density functional theory: DFT+U and hybrid functionals. Our results show that, upon releasing the symmetry constraint on the density but not on the geometry, charge disproportionation (Fe²⁺/Fe³⁺) is observed, resulting in a band gap of around 0.2 eV at the Fermi level. This implies that the Verwey transition is probably a semiconductor-to-semiconductor transition and that the conductivity mechanism above T_V is small polaron hopping.

Methanol on Anatase TiO2 (101): Mechanistic Insights into Photocatalysis

The photoactivity of methanol adsorbed on the anatase TiO2 (101) surface was studied by a combination of scanning tunneling microscopy (STM), temperature-programmed desorption (TPD), X-ray photoemission spectroscopy (XPS), and density functional theory (DFT) calculations. Isolated methanol molecules adsorbed at the anatase (101) surface show a negligible photoactivity. Two ways of methanol activation were found. First, methoxy groups formed by reaction of methanol with coadsorbed O2 molecules or terminal OH groups are photoactive, and they turn into formaldehyde upon UV illumination. The methoxy species show an unusual C 1s core-level shift of 1.4 eV compared to methanol; their chemical assignment was verified by DFT calculations with inclusion of final-state effects. The second way of methanol activation opens at methanol coverages above 0.5 monolayer (ML), and methyl formate is produced in this reaction pathway. The adsorption of methanol in the coverage regime from 0 to 2 ML is described in detail; it is key for understanding the photocatalytic behavior at high coverages. There, a hydrogen-bonding network is established in the adsorbed methanol layer, and consequently, methanol dissociation becomes energetically more favorable. DFT calculations show that dissociation of the methanol molecule is always the key requirement for hole transfer from the substrate to the adsorbed methanol. We show that the hydrogen-bonding network established in the methanol layer dramatically changes the kinetics of proton transfer during the photoreaction.

Water at the Interface Between Defective Graphene and Cu or Pt (111) Surfaces

The presence of defects in the graphenic layers deposited on metal surfaces modifies the nature of the interaction. Unsaturated carbon atoms, due to vacancies in the lattice, form strong organometallic bonds with surface metal atoms that highly enhance the binding energy between the two materials. We investigate by means of a wide set of dispersion-corrected density functional theory calculations how such strong chemical bonds affect both the electronic properties of these hybrid interfaces and the chemical reactivity with water, which is commonly present in the working conditions. We compare different metal substrates (Cu vs Pt) that present a different type of interaction with graphene and with defective graphene. This comparative analysis allows us to unravel the controlling factors of water reactivity, the role played by the carbon vacancies and by the confinement or "graphene cover effect". Water is capable of breaking the C-Cu bond by dissociating at the undercoordinated carbon atom of the vacancy, restoring the weak van der Waals type of interaction between the two materials that allows for an easy detachment of graphene from the metal, but the same is not true in the case of Pt, where C-Pt bonds are much stronger. These conclusions can be used to rationalize water reactivity at other defective graphene/metal interfaces.

Water Multilayers on TiO2 (101) Anatase Surface: Assessment of a DFTB-Based Method

A water/(101) anatase TiO₂ interface has been investigated with the DFT-based self-consistent-charge density functional tight-binding theory (SCC-DFTB). By comparison of the computed structural, energetic, and dynamical properties with standard DFT-GGA and experimental data, we assess the accuracy of SCC-DFTB for this prototypical solid–liquid interface. We tested different available SCC-DFTB parameters for Ti-containing compounds and, accordingly, combined them to improve the reliability of the method. To better describe water energetics, we have also introduced a modified hydrogen-bond-damping function (HBD). With this correction, equilibrium structures and adsorption energies of water on (101) anatase both for low (0.25 ML) and full (1 ML) coverages are in excellent agreement with those obtained with a higher level of theory (DFT-GGA). Furthermore, Born–Oppenheimer molecular dynamics (MD) simulations for mono-, bi-, and trilayers of water on the surface, as computed with SCC-DFTB, evidence similar ordering and energetics as DFT-GGA Car–Parrinello MD results. Finally, we have evaluated the energy barrier for the dissociation of a water molecule on the anatase (101) surface. Overall, the combined set of parameters with the HBD correction (SCC-DFTB+HBD) is shown to provide a description of the water/water/titania interface, which is very close to that obtained by standard DFT-GGA, with a remarkably reduced computational cost. Hence, this study opens the way to the future investigations on much more extended and realistic TiO₂/liquid water systems, which are extremely relevant for many modern technological applications.

Ab Initio Investigation of Polyethylene Glycol Coating of TiO2 Surfaces

In biomedical applications, TiO₂ nanoparticles are generally coated with polymers to prevent agglomeration, improve biocompatibility, and reduce cytotoxicity. Although the synthesis processes of such composite compounds are well established, there is still a substantial lack of information on the nature of the interaction between the titania surface and the organic macromolecules. In this work, the adsorption of polyethylene glycol (PEG) on the TiO₂ (101) anatase surface is modeled by means of dispersion-corrected density functional theory (DFT-D2) calculations. The two extreme limits of an infinite PEG polymer [-(OCH₂CH₂) _n ], on one side, and of a short PEG dimer molecule [H(OCH₂CH₂)₂OH], on the other, are analyzed. Many different molecular configurations and modes of adsorption are compared at increasing surface coverage densities. At low and medium coverage, PEG prefers to lay down on the surface, while at full coverage, the adsorption is maximized when PEG molecules bind perpendicularly to the surface and interact with each other through lateral dispersions, following a mushroom to brush transition. Finally, we also consider the adsorption of competing water molecules at different coverage densities, assessing whether PEG would remain bonded to the surface or desorb in the presence of the aqueous solvent.

Synthesis of graphene nanoribbons with a defined mixed edge-site sequence by surface assisted polymerization of (1,6)-dibromopyrene on Ag(110)

By a combination of scanning tunneling microscopy, X-ray spectroscopic techniques and density functional theory calculations, we prove the formation of extended patterns of parallel, graphene nanoribbons with alternate zig-zag and armchair edges and selected width by surface-assisted Ullmann coupling polymerization and dehydrogenation of 1,6-dibromopyrene (C₁₆H₈Br₂). Besides the relevance of these nanostructures for their possible application in nanodevices, we demonstrate the peculiarity of halogenated pyrene derivatives for the formation of nanoribbons, in particular on Ag(110). These results open the possibility of tuning the shape and dimension of nanoribbons (and hence the correlated electronic properties) by choosing suitably tailored or on-purpose designed molecular precursors.

Magnetic properties of nitrogen-doped ZrO2: Theoretical evidence of absence of room temperature ferromagnetism

N-dopants in bulk monoclinic ZrO2 and their magnetic interactions have been investigated by DFT calculations, using the B3LYP hybrid functional. The electronic and magnetic properties of the paramagnetic N species, substitutionals and interstitials, are discussed. Their thermodynamic stability has been estimated as a function of the oxygen partial pressure. At 300 K, N prefers interstitial sites at any range of oxygen pressure, while at higher temperatures (700-1000 K), oxygen poor-conditions facilitate substitutional dopants. We have considered the interaction of two N defects in various positions in order to investigate the possible occurrence of ferromagnetic ordering. A very small magnetic coupling constant has been calculated for several 2N-ZrO2 configurations, thus demonstrating that magnetic ordering can be achieved only at very low temperatures, well below liquid nitrogen. Furthermore, when N atoms replace O at different sites, resulting in slightly different positions of the corresponding N 2p levels, a direct charge transfer can occur between the two dopants with consequent quenching of the magnetic moment. Another mechanism that contributes to the quenching of the N magnetic moments is the interplay with oxygen vacancies. These effects contribute to reduce the concentration of magnetic impurities, thus limiting the possibility to establish magnetic ordering.

Catalysis under Cover: Enhanced Reactivity at the Interface between (Doped) Graphene and Anatase TiO2

The "catalysis under cover" involves chemical processes which take place in the confined zone between a 2D material, such as graphene, h-BN, or MoS2, and the surface of an underlying support, such as a metal or a semiconducting oxide. The hybrid interface between graphene and anatase TiO2 is extremely important for photocatalytic and catalytic applications because of the excellent and complementary properties of the two materials. We investigate and discuss the reactivity of O2 and H2O on top and at the interface of this hybrid system by means of a wide set of dispersion-corrected hybrid density functional calculations. Both pure and boron- or nitrogen-doped graphene are interfaced with the most stable (101) anatase surface of TiO2 in order to improve the chemical activity of the C-layer. Especially in the case of boron, an enhanced reactivity toward O2 dissociation is observed as a result of both the contribution of the dopant and of the confinement effect in the bidimensional area between the two surfaces. Extremely stable dissociation products are observed where the boron atom bridges the two systems by forming very stable B-O covalent bonds. Interestingly, the B defect in graphene could also act as the transfer channel of oxygen atoms from the top side across the C atomic layer into the G/TiO2 interface. On the contrary, the same conditions are not found to favor water dissociation, proving that the "catalysis under cover" is not a general effect, but rather highly depends on the interfacing material properties, on the presence of defects and impurities and on the specific reaction involved.

Theoretical Studies of Oxygen Reactivity of Free-Standing and Supported Boron-Doped Graphene

Graphene inertness towards chemical reactivity can be considered as an accepted postulate by the research community. This limit has been recently overcome by chemically and physically modifying graphene through non-metal doping or interfacing with acceptor/donor materials (metals or semiconductors). As a result, outstanding performances as catalytic, electrocatalytic, and photocatalytic material have been observed. In this critical Review we report computational work performed, by our group, on the reactivity of free-standing, metal- and semiconductor-supported B-doped graphene towards oxygen, which is at the basis of extremely important energy-related chemical processes, such as the oxygen reduction reaction. It appears that a combination of doping and interfacing approaches for the activation of graphene can open unconventional and unprecedented reaction paths, thus boosting the potential of modified graphene in many chemical applications.