Parametrization of the Fe-Owater cross-interaction for a more accurate Fe3O4/water interface model and its application to a spherical Fe3O4 nanoparticle of realistic size.

The accurate description of iron oxides/water interfaces requires reliable force field parameters that can be developed through comparison with sophisticated quantum mechanical calculations. Here, a set of CLASS2 force field parameters is optimized to describe the Fe-Owater cross-interaction through comparison with hybrid density functional theory (HSE06) calculations of the potential energy function for a single water molecule adsorbed on the Fe3O4 (001) surface and with density functional tight binding (DFTB+U) molecular dynamics simulations for a water trilayer on the same surface. The performance of the new parameters is assessed through the analysis of the number density profile of a water bulk (12 nm) sandwiched between two magnetite slabs of large surface area. Their transferability is tested for water adsorption on the curved surface of a spherical Fe3O4 nanoparticle of realistic size (2.5 nm).

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

In this work, we investigate the Fe3O4 (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 Fe3O4 (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⋯OH2 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.

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.

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. Here, we report on a combined scanning tunneling microscopy, soft X-ray spectroscopy and density functional theory investigation on the surface assisted polymerization of Br-corannulene at Ag(110) and on the possibility of building a mesh of π-conjugated polymers starting from buckyball shaped molecules. Indeed, the corannulene units form one-molecule-wide ribbons in which the natural concavity of the precursor molecule is maintained. These C-based nanostructures are corrugated and merge into a covalent network on the surface.

Interfacing CRYSTAL/AMBER to Optimize QM/MM Lennard⁻Jones Parameters for Water and to Study Solvation of TiO2 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 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.

Modelling realistic TiO2 nanospheres: A benchmark study of SCC-DFTB against hybrid DFT.

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 sizes (from 300 to 4000 atoms) have been studied with a two-scale computational approach. Global optimization to obtain stable and equilibrated nanospheres 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 the 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 density functional theory 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.

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.

Spectroscopic properties of doped and defective semiconducting oxides from hybrid density functional calculations.

CONSPECTUS: Very rarely do researchers use metal oxides in their pure and fully stoichiometric form. In most of the countless applications of these compounds, ranging from catalysis to electronic devices, metal oxides are either doped or defective because the most interesting chemical, electronic, optical, and magnetic properties arise when foreign components or defects are introduced in the lattice. Similarly, many metal oxides are diamagnetic materials and do not show a response to specific spectroscopies such as electron paramagnetic resonance (EPR) spectroscopy. However, doped or defective oxides may exhibit an interesting and informative paramagnetic behavior. Doped and defective metal oxides offer an expanding range of applications in contemporary condensed matter science; therefore researchers have devoted enormous effort to the understanding their physical and chemical properties. The interplay between experiment and computation is particularly useful in this field, and contemporary simulation techniques have achieved high accuracies with these materials. In this Account, we show how the direct comparison between spectroscopic experimental and computational data for some selected and relevant materials provides ways to understand and control these complex systems. We focus on the EPR properties and electronic transitions that arise from the presence of dopants and defects in bulk metal oxide materials. We analyze and compare the effect of nitrogen doping in TiO2 and ZnO (two semiconducting oxides) and MgO (a wide gap insulator) and examine the effect of oxygen deficiency in the semiconducting properties of TiO2-x, ZnO1-x, and WO3-x materials. We chose these systems because of their relevance in applications including photocatalysis, touch screens, electrodes in magnetic random access memories, and smart glasses. Density functional theory (DFT) provides the general computational framework used to illustrate the electronic structure of these systems. However, for a more accurate description of the oxide band gap and of the electron localization of the impurity states associated with dopants or defects, we resorted to the use of hybrid functionals (B3LYP), where a portion of exact exchange in the exchange-correlation functional partly corrects for the self-interaction error inherent in DFT. In many cases, the self-interaction correction is very important, and these results can lead to a completely different physical picture than that obtained using local or semilocal functionals. We analyzed the electronic transitions in terms of their transition energy levels, which provided a more accurate comparison with experimental spectroscopic data than Kohn-Sham eigenvalues. The effects of N-doping were similar among the three oxides that we considered. The nature of the impurity state is always localized at the dopant site, which may limit their application in photocatalytic processes. Photocatalytic systems require highly delocalized photoexcited carriers within the material to effectively trigger redox processes at the surface. The nature of the electronic states associated with the oxygen deficiency differed widely in the three investigated oxides. In ZnO1-x and WO3-x the electronic states resemble the typical F-centers in insulating oxides or halides, with the excess electron density localized at the vacancy site. However, TiO2 acts as a reducible oxide, and the removal of neutral oxygen atoms reduced Ti(4+) to Ti(3+).

Control of the intermolecular coupling of dibromotetracene on Cu(110) by the sequential activation of C-Br and C-H bonds.

Dibromotetracene molecules are deposited on the Cu(110) surface at room temperature. The complex evolution of this system has been monitored at different temperatures (i.e., 298, 523, 673, and 723 K) by means of a variety of complementary techniques that range from STM and temperature-programmed desorption (TPD) to high-resolution X-ray spectroscopy (XPS) and near-edge X-ray absorption fine structure spectroscopy (NEXAFS). State-of-the-art density-functional calculations were used to determine the chemical processes that take place on the surface. After deposition at room temperature, the organic molecules are transformed into organometallic monomers through debromination and carbon-radical binding to copper adatoms. Organometallic dimers, trimers, or small oligomers, which present copper-bridged molecules, are formed by increasing the temperature. Surprisingly, further heating to 673 K causes the formation of elongated chains along the Cu(110) close-packed rows as a consequence of radical-site migration to the thermodynamically more stable molecule heads. Finally, massive dehydrogenation occurs at the highest temperature followed by ring condensation to nanographenic patches. This study is a paradigmatic example of how intermolecular coupling can be modulated by the stepwise control of a simple parameter, such as temperature, through a sequence of domino reactions.

Communication: Hole localization in Al-doped quartz SiO2 within ab initio hybrid-functional DFT.

We investigate the long-standing problem of hole localization at the Al impurity in quartz SiO2, using a relatively recent DFT hybrid-functional method in which the exchange fraction is obtained ab initio, based on an analogy with the static many-body COHSEX approximation to the electron self-energy. As the amount of the admixed exact exchange in hybrid functionals has been shown to be determinant for properly capturing the hole localization, this problem constitutes a prototypical benchmark for the accuracy of the method, allowing one to assess to what extent self-interaction effects are avoided. We obtain good results in terms of description of the charge localization and structural distortion around the Al center, improving with respect to the more popular B3LYP hybrid-functional approach. We also discuss the accuracy of computed hyperfine parameters, by comparison with previous calculations based on other self-interaction-free methods, as well as experimental values. We discuss and rationalize the limitations of our approach in computing defect-related excitation energies in low-dielectric-constant insulators.

Defect calculations in semiconductors through a dielectric-dependent hybrid DFT functional: The case of oxygen vacancies in metal oxides.

We investigate the behavior of oxygen vacancies in three different metal-oxide semiconductors (rutile and anatase TiO2, monoclinic WO3, and tetragonal ZrO2) using a recently proposed hybrid density-functional method in which the fraction of exact exchange is material-dependent but obtained ab initio in a self-consistent scheme. In particular, we calculate charge-transition levels relative to the oxygen-vacancy defect and compare computed optical and thermal excitation/emission energies with the available experimental results, shedding light on the underlying excitation mechanisms and related materials properties. We find that this novel approach is able to reproduce not only ground-state properties and band structures of perfect bulk oxide materials but also provides results consistent with the optical and electrical behavior observed in the corresponding substoichiometric defective systems.

Vitamin C Affinity to TiO2 Nanotubes: A Computational Study by Hybrid Density Functional Theory Calculations.

Titanium dioxide nanotubes (TNT) have been extensively studied because of their unique properties, which make such systems ideal candidates for biomedical application, especially for the targeted release of drugs. However, knowledge about the properties of TiO2 nanotubes with typical dimensions of the order of the nanometer is limited, especially concerning the adsorption of molecules that can be potentially loaded in actual devices. In this work, we investigate, by means of simulations based on hybrid density functional theory, the adsorption of Vitamin C molecules on different nanotubes through a comparative analysis of the properties of different structures. We consider two different anatase TiO2 surfaces, the most stable (101) and the more reactive (001)A; we evaluate the role of the curvature, the thickness and of the diameter as well as of the rolling direction of the nanotube. Different orientations of the molecule with respect to the surface are studied in order to identify any trends in the adsorption mechanism. Our results show that there is no preferential functional group of the molecule interacting with the substrate, nor any definite spatial dependency, like a rolling orientation or the concavity of the nanotube. Instead, the adsorption is driven by geometrical factors only, i.e., the favorable matching of the position and the alignment of any functional groups with undercoordinated Ti atoms of the surface, through the interplay between chemical and hydrogen bonds. Differently from flat slabs, thicker nanotubes do not improve the stability of the adsorption, but rather develop weaker interactions, due to the enhanced curvature of the substrate layers.

The effect of polymer coating on nanoparticles' interaction with lipid membranes studied by coarse-grained molecular dynamics simulations.

Nanoparticles' (NPs) permeation through cell membranes, whether it happens via passive or active transport, is an essential initial step for their cellular internalization. The NPs' surface coating impacts the way they translocate through the lipid bilayer and the spontaneity of the process. Understanding the molecular details of NPs' interaction with cell membranes allows the design of nanosystems with optimal characteristics for crossing the lipid bilayer: computer simulations are a powerful tool for this purpose. In this work, we have performed coarse-grained molecular dynamics simulations and free energy calculations on spherical titanium dioxide NPs conjugated with polymer chains of different chemical compositions. We have demonstrated that the hydrophobic/hydrophilic character of the chains, more than the nature of their terminal group, plays a crucial role in determining the NPs' interaction with the lipid bilayer and the thermodynamic spontaneity of NPs' translocation from water to the membrane. We envision that this computational work will be helpful to the experimental community in terms of the rational design of NPs for efficient cell membrane permeation.

Improving the Oxygen Evolution Reaction on Fe3O4(001) with Single-Atom Catalysts.

Doping magnetite surfaces with transition-metal atoms is a promising strategy to improve the catalytic performance toward the oxygen evolution reaction (OER), which governs the overall efficiency of water electrolysis and hydrogen production. In this work, we investigated the Fe3O4(001) surface as a support material for single-atom catalysts of the OER. First, we prepared and optimized models of inexpensive and abundant transition-metal atoms, such as Ti, Co, Ni, and Cu, trapped in various configurations on the Fe3O4(001) surface. Then, we studied their structural, electronic, and magnetic properties through HSE06 hybrid functional calculations. As a further step, we investigated the performance of these model electrocatalysts toward the OER, considering different possible mechanisms, in comparison with the pristine magnetite surface, on the basis of the computational hydrogen electrode model developed by Nørskov and co-workers. Cobalt-doped systems were found to be the most promising electrocatalytic systems among those considered in this work. Overpotential values (∼0.35 V) were in the range of those experimentally reported for mixed Co/Fe oxide (0.2-0.5 V).

Modeling Zeta Potential for Nanoparticles in Solution: Water Flexibility Matters.

Nonequilibrium molecular dynamics simulations were performed to study the electrokinetic properties of five mainstream TIPxP water models (namely, TIP3P-FB, TIP3Pm, TIP4P-FB, TIP4P-Ew, and TIP4P/2005) in NaCl aqueous solutions in the presence of a negatively charged TiO2 surface. The impact of solvent flexibility and system geometry on the electro-osmotic (EO) mobility and flow direction was systematically assessed and compared. We found that lack of water flexibility decelerates the forward EO flow of aqueous solutions at moderate (0.15 M) or high (0.30 M) NaCl concentrations, in some special cases to such an extent that EO flow reversal occurs. Zeta potential (ZP) values were then determined from the bulk EO mobilities using the Helmholtz-Smoluchowski formula. The straight comparison against available experimental data strongly suggests that water flexibility improves the ZP determination of NaCl solutions adjacent to a realistic TiO2 surface under neutral pH conditions.

Metadynamics simulations for the investigation of drug loading on functionalized inorganic nanoparticles.

Inorganic nanoparticles show promising properties that allow them to be efficiently used as drug carriers. The main limitation in this type of application is currently the drug loading capacity, which can be overcome with a proper functionalization of the nanoparticle surface. In this study, we present, for the first time, a computational approach based on metadynamics to estimate the binding free energy of the doxorubicin drug (DOX) to a functionalized TiO2 nanoparticle under different pH conditions. On a thermodynamic basis, we demonstrate the robustness of our approach to capture the overall mechanism behind the pH-triggered release of DOX due to environmental pH changes. Notably, binding free energy estimations align well with what is expected for a pH-sensitive drug delivery system. Based on our results, we envision the use of metadynamics as a promising computational tool for the rational design and in silico optimization of organic ligands with improved drug carrier properties.

Chemistry of the Interaction and Retention of TcVII and TcIV Species at the Fe3O4(001) Surface.

The pertechnetate ion TcVIIO4 - is a nuclear fission product whose major issue is the high mobility in the environment. Experimentally, it is well known that Fe3O4 can reduce TcVIIO4 - to TcIV species and retain such products quickly and completely, but the exact nature of the redox process and products is not completely understood. Therefore, we investigated the chemistry of TcVIIO4 - and TcIV species at the Fe3O4(001) surface through a hybrid DFT functional (HSE06) method. We studied a possible initiation step of the TcVII reduction process. The interaction of the TcVIIO4 - ion with the magnetite surface leads to the formation of a reduced TcVI species without any change in the Tc coordination sphere through an electron transfer that is favored by the magnetite surfaces with a higher FeII content. Furthermore, we explored various model structures for the immobilized TcIV final products. TcIV can be incorporated into a subsurface octahedral site or adsorbed on the surface in the form of TcIVO2·xH2O chains. We propose and discuss three model structures for the adsorbed TcIVO2·2H2O chains in terms of relative energies and simulated EXAFS spectra. Our results suggest that the periodicity of the Fe3O4(001) surface matches that of the TcO2·2H2O chains. The EXAFS analysis suggests that, in experiments, TcO2·xH2O chains were probably not formed as an inner-shell adsorption complex with the Fe3O4(001) surface.

Molecular Dynamics for the Optimal Design of Functionalized Nanodevices to Target Folate Receptors on Tumor Cells.

Atomistic details on the mechanism of targeting activity by biomedical nanodevices of specific receptors are still scarce in the literature, where mostly ligand/receptor pairs are modeled. Here, we use atomistic molecular dynamics (MD) simulations, free energy calculations, and machine learning approaches on the case study of spherical TiO2 nanoparticles (NPs) functionalized with folic acid (FA) as the targeting ligand of the folate receptor (FR). We consider different FA densities on the surface and different anchoring approaches, i.e., direct covalent bonding of FA γ-carboxylate or through polyethylene glycol spacers. By molecular docking, we first identify the lowest energy conformation of one FA inside the FR binding pocket from the X-ray crystal structure, which becomes the starting point of classical MD simulations in a realistic physiological environment. We estimate the binding free energy to be compared with the existing experimental data. Then, we increase complexity and go from the isolated FA to a nanosystem decorated with several FAs. Within the simulation time framework, we confirm the stability of the ligand-receptor interaction, even in the presence of the NP (with or without a spacer), and no significant modification of the protein secondary structure is observed. Our study highlights the crucial role played by the spacer, FA protonation state, and density, which are parameters that can be controlled during the nanodevice preparation step.