Search

Monolayer WSe2 has aroused strong interest as a semiconducting 2D material. The tuning of its electronic properties is important for expanding the range of applications. In this work, by hybrid density functional theory calculations, we investigate the effect of Na adsorption on the electronic and magnetic properties of stoichiometric and defective WSe2 monolayers and their WSe2xS2(1-x) alloys. Our results suggest that Na adsorption on defective WSe2xS2(1-x) monolayers would split and partially fill the empty defect states in the gap region inducing a local magnetic moment of 1 μ B per alkaline atom. The efficient tuning of the electronic and magnetic properties of WSe2 monolayer by Na adsorption combined with the Se vacancy and S alloying would benefit its wide applications.

High-quality nitrogen-doped graphene on nickel is prepared by exploiting both the catalytic properties of nickel and the solubility of nitrogen atoms into its bulk. Following the standard chemical vapor deposition procedure, a previously nitrogen-doped nickel substrate is exposed to carbon-containing precursors so that nitrogen atoms, segregating to the surface, remain trapped in the growing graphene network. Morphological and chemical characterization by scanning tunneling microscopy and X-ray photoelectron spectroscopy demonstrates that the process yields a flat, wide, continuous nitrogen-doped graphene layer. Experimental results are combined with a thorough density functional theory investigation of possible structural models, to obtain a clear description at the atomic scale of the various configurations of the nitrogen atoms observed in the graphene mesh. This growth method is potentially scalable and suitable for the production of high-performance nano-devices with well-defined nitrogen centers, to be exploited as metal-free carbon-based catalysts in several applicative fields such as electrochemistry, energy storage, gas storage/sensing or wastewater treatment., 21 pages, including 4 figures, plus SI

Under near-ambient pressure conditions, carbon monoxide molecules intercalate underneath an epitaxial graphene monolayer grown on Ni(111), getting trapped into the confined region at the interface. On the basis of ab-initio density functional theory calculations, we provide here a full investigation of the intercalated CO pattern, highlighting the modifications induced on the graphene electronic structure. For a CO coverage as low as 0.14 monolayer (ML), the graphene layer is spatially decoupled from the metallic substrate, with a significant C 1s core level shift towards lower binding energies. The most relevant signature of the CO intercalation is a clear switching of the graphene doping state, which changes from n-type, when strongly interacting with the metal surface, to p-type. The shift of the Dirac cone linearly depends on the CO coverage, reaching about 0.9 eV for the saturation value of 0.57 ML. Theoretical predictions are compared with the results of scanning tunnelling microscopy, low-energy electron diffraction and photoemission spectroscopy experiments, which confirm the proposed scenario for the nearly saturated intercalated CO system. This result opens the way to the application of the graphene/Ni(111) interface as gas sensor to easily detect and quantify the presence of carbon monoxide.

Single Atom Catalysts (SAC) in graphene have been recently gaining more and more attention. They are usually non-noble transition metal (TM) adatoms getting trapped at the carbon vacancies during the fabrication of the graphene layer, which then act as active centers for catalysis and adsorption. In this work we present a systematic and comparative investigation, by means of dispersion–corrected density functional theory (DFT) calculations, of Fe, Co, Ni, and Cu as possible SACs when they become trapped at graphene C vacancies. The stability of these TM atoms is further increased by introducing pyridinic nitrogen (N) atoms and transforming graphene into a giant porphyrin-like macrocyclic ligand. The structural, electronic and energetics properties of these systems, even under the effect of a metal substrate (weakly interacting Cu (111) or strongly interacting Ni (111)), are comparatively examined in great detail by means of crystal/ligand field theories and through ad-hoc energy decomposition analysis to highlight trends and peculiar behaviors. The position of the TM d-orbitals with respect to the Fermi level of the whole system is of considerable importance for designing prospective device applications in catalysis, electrocatalysis and sensors. To this purpose, we also examine how the reactivity of the SACs in graphene towards the hydrogen evolution reaction (HER) can be tuned with N-doping and with different substrates.

In many potential applications, nanoparticles are typically in an aqueous medium. This has strong influence on the stability, optical properties and reactivity, in particular for their functionalization. Therefore, the understanding of the chemistry at the interface between the solvent and the nanoparticle is of utmost importance. In this work, we present a comparative ReaxFF reactive molecular dynamics investigation on spherical TiO2 nanoparticles (NSs) of realistic size, with diameters from 2.2 to 4.4 nm, immersed in a large drop of bulk water. After force field validation for its use for a curved anatase TiO2 surface/water interface, we performed several simulations of the TiO2 nanoparticles of increasing size in a water drop. We found that water can be adsorbed jointly in a molecular and dissociative way on the surface. A Langmuir isotherm indicating an adsorption/desorption mechanism of water on the NS is observed. Regarding the dissociative adsorption, atomistic details reveal two different mechanisms, depending on the water concentration around the NS. At low coverage, the first mechanism involves direct dissociation of a single water molecule, whereas, at higher water coverage, the second mechanism is a proton transfer reaction involving two water molecules, also known as Grotthuss-like mechanism. Thermal annealing simulations show that several water molecules remain on the surface in agreement with the experimental reports. The capacity of adsorption is higher for the 2.2 and 3.0 nm NSs than for the 4.4 nm NS. Finally, a comparative investigation with flat surfaces indicates that NSs present a higher water adsorption capacity (undissociated and dissociated) than flat surfaces, which can be rationalized considering that NSs present many more low-coordinated Ti atoms available for water adsorption. This journal is.

Metal-supported hexagonal boron nitride monolayers (h-BN/M) are emerging as new potential electrocatalysts for various energy-related oxidation or reduction process. So far, several preparation methods have been developed to introduce, in a controlled way, defects such as vacancies or substitutional heteroatoms. Herein, we investigate by dispersion-corrected density functional theory (DFT) calculations, defective and metal-doped h-BN/Cu(111) systems as electrocatalysts for the hydrogen evolution reaction (HER). By calculating the hydrogen binding energy (ΔG*H) at different coverage conditions, we observe how the interaction between the defective/metal-doped h-BN layer and the Cu(111) substrate plays a key role in tuning the reactivity, leading to a thermoneutral hydrogen adsorption step (i.e., ΔG*H ≈ 0). These results could be generalized to other h-BN/M interfaces and may help their rational design for an improved H2-evolving electrocatalysis.

Polymeric graphitic carbon nitride (gCN) compounds are promising materials in photoactivated electrocatalysis thanks to their peculiar structure of periodically spaced voids exposing reactive pyridinic N atoms. These are excellent sites for the adsorption of isolated transition metal atoms or small clusters that can highly enhance the catalytic properties. However, several polymorphs of gCN can be obtained during synthesis, differing for their structural and electronic properties that ultimately drive their potential as catalysts. The accurate characterization of the obtained material is critical for the correct rationalization of the catalytic results, however, an unambiguous experimental identification of the actual polymer is challenging, especially without any reference spectroscopic features for the assignment. In this work, we optimized several models of melem-based gCN, taking into account different degrees of polymerization and arrangement of the monomers, and we present a thorough computational characterization of their simulated XRD, XPS, and NEXAFS spectroscopic properties, based on state-of-the-art density functional theory calculations. Through this detailed study, we could identify the peculiar fingerprints of each model and correlate them with its structural and/or electronic properties. Theoretical predictions were compared with the experimental data whenever they were available.

We report the study of two-dimensional graphitic carbon nitride (GCN) functionalized with copper single atoms as a catalyst for the reduction of CO2 (CO2RR). The correct GCN structure, as well as the adsorption sites and the coordination of the Cu atoms, was carefully determined by combining experimental techniques, such as X-ray diffraction, transmission electron microscopy, X-ray absorption, and X-ray photoemission spectroscopy, with DFT theoretical calculations. The CO2RR products in KHCO3 and phosphate buffer solutions were determined by rotating ring disk electrode measurements and confirmed by 1H-NMR and gas chromatography. Formate was the only liquid product obtained in bicarbonate solution, whereas only hydrogen was obtained in phosphate solution. Finally, we demonstrated that GCN is a promising substrate able to stabilize metal atoms, since the characterization of the Cu-GCN system after the electrochemical work did not show the aggregation of the copper atoms.

Inorganic nanoparticles are gaining increasing attention as drug carriers because they respond to external physical stimuli, allowing therapy to be combined with diagnosis. Their drawback is low drug loading capacity, which can be improved by proper and efficacious functionalization. In this computational study, we take TiO2 spherical nanoparticles as prototype photoresponsive inorganic nanoparticles and we fully decorate them with two different types of bifunctional ligands: TETTs and DOPACs, which present different surface anchoring groups (silanol or catechol) but the same drug tethering COOH group, although in different concentrations (3 vs. 1), thus causing different steric hindrances. Then, we put these two types of nanocarriers in bulk water and in the presence of several DOX molecules and let the systems evolve through molecular dynamics (MD) simulations, clearly observing drug loading on the nanocarriers. This comparative MD study allows the investigation of the loading mechanism, performance of a conformational analysis and establishment of the guiding interactions through an energy decomposition analysis. We learn that DOX mostly interacts with the functionalized NPs through electrostatics, as a consequence of the protonated amino group, although several H-bonds are also established both with the ligands and with the oxide surface. Different ligands induce a different electrostatic potential around the NP; therefore, those which lead to the formation of more negative hotspots (here TETTs) are found to favour DOX binding. The leading role of electrostatics can provide a rational explanation for a pH-dependent drug release mechanism that is often invoked for DOX when reaching diseased cells because under anomalous acidic conditions both the NP surface and the carboxylate groups of the ligands are expected to get protonated, which of course would weaken, if not totally quench, the interaction of the nanocarrier with protonated DOX., A comparative MD study unravels the loading mechanism of the doxorubicin drug to TiO2 nanoparticles functionalized with two different organic bifunctional ligands: TETT and DOPAC.

Donor-modified TiO2 nanoparticles are interesting hybrid systems shifting the absorption edge of this semiconductor from the ultra-violet to the visible or infrared light spectrum, which is a benefit for several applications ranging from photochemistry, photocatalysis, photovoltaics, or photodynamic therapy. Here, we investigate the absorption properties of two catechol-like molecules, that is, dopamine and DOPAC ligands, when anchored to a spherical anatase TiO2 nanoparticle of realistic size (2.2 nm), by means of time-dependent density functional theory calculations. By the differential absorbance spectra with the bare nanoparticle, we show how it is possible to determine the injection mechanism. Since new low-energy absorption peaks are observed, we infer a direct charge transfer injection, which, unexpectedly, does not involve the lowest energy conduction band states. We also find that the more perpendicular the molecular benzene ring is to the surface, the more intense is the absorption, which suggests aiming at high molecular packing in the synthesis. Through a comparative investigation with a flat TiO2 surface model, we unravel both the curvature and coverage effects.

Iron oxide magnetic nanoparticles (NPs) are excellent systems in catalysis and in nanomedicine, where they are mostly immersed in aqueous media. Even though the NP solvation by water is expected to play an active role, the detailed structural insight at the nanostructure oxide/water interface is still missing. Here, based on our previous efforts to obtain accurate models of dehydrated Fe3O4 NPs and of their magnetic properties and through multiscale molecular dynamics simulations combining the density functional tight binding method and force field, we unravel the atomistic details of the short range (chemical) and long range (physical) interfacial effects when magnetite nanoparticles are immersed in water. The influence of the first hydration shell on the structural, electronic and magnetic properties of Fe3O4 NPs is revealed by high-level hybrid density functional calculations. Hydrated Fe3O4 NPs possess larger magnetic moment than dehydrated ones. This work bridges the large gap between experimental studies on solvated Fe3O4 NPs and theoretical investigations on flat Fe3O4 surfaces covered with water and paves the way for further study of Fe3O4 NPs in biological environments. This journal is

The accurate description of iron oxides/water interfaces requires reliable force field parameters that can be developed through the 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 tri-layer 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 the water adsorption on the curved surface of a spherical Fe3O4 nanoparticle of realistic size (2.5 nm).

Molecules intercalate at the graphene/metal interface even though defect-free graphene is impermeable to any atomic and molecular species in the gas and liquid phase, except hydrogen. The mechanism of molecular intercalation is still a big open question. In this Letter, by means of a combined experimental (STM, XPS, and LEED) and theoretical (DFT) study, we present a proof of how CO molecules succeed in permeating the graphene layer and get into the confined zone between graphene and the Ni(111) surface. The presence of N-dopants in the graphene layer is found to highly facilitate the permeation process, reducing the CO threshold pressure by more than one order of magnitude, through the stabilization of multiatomic vacancy defects that are the open doors to the bidimensional nanospace, with crucial implications for the catalysis under cover and for the graphene-based electrochemistry.

Nanoparticle functionalization is a modern strategy in nanotechnology to build up devices for several applications. Modeling fully decorated metal oxide nanoparticles of realistic size (few nanometers) in an aqueous environment is a challenging task. In this work, we present a case study relevant for solar-light exploitation and for biomedical applications, i.e., a dopamine-functionalized TiO2 nanoparticle (1700 atoms) in bulk water, for which we have performed an extensive comparative investigation with both MM and QM/MM approaches of the structural properties and of the conformational dynamics. We have used a combined multiscale protocol for a more efficient exploration of the complex conformational space. On the basis of the results of this study and of some QM and experimental data, we have defined strengths and limitations of the existing force field parameters. Our findings will be useful for an improved modeling and simulation of many other similar hybrid bioinorganic nanosystems in an aqueous environment that are pivotal in a broad range of nanotechnological applications.

In this work, we present a quantum mechanical investigation, based on the self-consistent charge density functional tight-binding (SCC-DFTB) method, of the functionalization with silane-type ligands (TETT) of a spherical TiO2 nanoparticle of realistic size (2.2 nm containing 700 atoms) to create an efficient nanosystem for simultaneous photodynamic therapy and drug transport. We determine the mechanism of the TETT ligand anchoring and its stability under thermal treatment, through molecular dynamics simulations at 300 K. Then, we build a medium and a full coverage model (22 and 40 TETTs, respectively) and analyze the interaction among TETT ligands and between the ligands and the surface. Finally, on the fully covered nanoparticle, we succeed in localizing two minimum energy structures for an attached doxorubicin anticancer molecule (DOX) and provide the atomistic details for both the covalent and the non-covalent (electrostatic) types of interaction. A future development of this work will be the investigation of the loading capacity of this drug delivery system and of the pH effect of the surrounding aqueous environment.

In the confined zone between a bidimensional material and a metal surface, unexpected effects can take place. In this study, we show that when a nonregular two-dimensional h-BN layer is grown on a Cu(111) surface, metal adatoms spontaneously pop up from the bulk to fill the holes in the structure. We provide ample theoretical support to our findings based on a large set of dispersion-corrected density functional theory calculations and on a detailed analysis of the electronic properties and of the chemical processes at this peculiar interface. The observation can be rationalized in terms of a high affinity of Cu adatoms toward N-donor species. Defective h-BN, exposing N-terminated edges, behaves like a giant multi-N-donor macrocyclic ligand that can encapsulate metal atoms as a consequence of a huge stabilization deriving from the Cu–N bond formation. Our conclusions could apply to other metal surfaces and could even stimulate the idea of trapping different metal atoms from those of the underlying surface (e.g., more precious but more active metals) for catalytic purposes

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 TiO2 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/TiO2 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 TiO2 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 TiO2 surface, showing a prominent cover effect. Reduced TiO2 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/TiO2 interface and thus its application in photocatalysis, electronic devices, and nanomedicine.

Formaldehyde (CH2O) adsorption on the anatase TiO2(101) surface was studied with a combination of experimental and theoretical methods. Scanning tunneling microscopy, noncontact atomic force microscopy, temperature-programmed desorption, and X-ray photoelectron spectroscopy were employed on the experimental side. Density functional theory was used to calculate formaldehyde adsorption configurations and energy barriers for transitions between them. At low coverages (

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 O2 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 O2 release from the catalyst, which is often observed to be a major limit for other potential electrocatalysts.

Functionalizing nanoparticles (NPs) with biological molecules is apromising modern strategy in bionanotechnology to build up smart bioinorganicdevices for medical applications. Bifunctional linkers provide an interesting andductile bioconjugation approach especially because they behave not only asanchoring and tethering agents but also as spacers between the NP and thebiomolecules, which helps in maintaining their 3D structural and functionalproperties. In this work, by means of a wide set of density functional theory(DFT) electronic structure calculations and density functional tight binding(DFTB) molecular dynamics simulations, we provide an all-round investigation ofthe functionalization of realistic curved TiO2 NPs (2−3 nm size with∼800atoms) with a catechol derivative, such as dopamine and DOPAC. We span fromsingle-molecule adsorption to the full coverage regime. For the low coverage, weachieve a detailed description of the mechanisms of molecular adsorption, of theinterfacial electronic charge-transfer effects, and of the processes following visiblelight irradiation (exciton formation, trapping, charge carrier diffusion, or recombination). We then consider a growing molecularlayer on the NP and analyze the self-assembling mechanism and the effects on the electronic properties of the complex. Finally,for the maximum coverage (46 molecules per NP) we perform molecular dynamics runs at 300 K to compare the molecularconfiguration and electronic properties of the NP/linker complex interface before and after thermal treatment to better accountfor the competition between molecule/surface and molecule/molecule interactions. The use of curved NP surfaces combinedwith dopamine, with respect to aflat one and DOPAC, respectively, is found to be more effective for bioconjugation

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.

Sodium tantalate, NaTaO3, has attracted interest as one of the few photocatalysts capable to perform overall water splitting, that is, simultaneously produce oxygen and hydrogen from water. In particular, an interesting and not fully understood observation is that the efficiency of NaTaO3 increases dramatically in the presence of cation doping. To obtain better insight into the origin of this effect, we use first-principles calculations to investigate the fundamental structural, electronical, and chemical properties of pristine and Sr-doped NaTaO3, a system for which several experimental studies have recently become available. Our results show that Sr donor-acceptor codoping at Na and Ta sites significantly reduces the formation energy of the Sr dopants. Further study of the energetics of the oxygen evolution reaction (OER) shows a substantial reduction of the OER overpotential for the codoped material, consistent with recent suggestions that codoping is crucial for increasing NaTaO3's efficiency. The detailed insights provided by our work could benefit the design and preparation of new efficient catalysts based on NaTaO3

TiO2 nanoparticles (NPs) are attracting materials for biomedical applications, provided that they are coated with polymers to improve solubility, dispersion and biocompatibility. Conformation, coverage density and solvent effects largely influence their functionality and stability. In this work, we use atomistic molecular dynamics simulations to study polyethylene glycol (PEG) grafting to highly curved TiO2 NPs (2-3 nm) in different solvents. We compare the coating polymer conformations on NPs with those on (101) flat surfaces. In water, the transition from mushroom to brush conformation starts only at high density ({\sigma} = 2.25 chains/nm2). In dichloromethane (DCM), at low-medium coverage ({\sigma} < 1.35 chains/nm2), several interactions between the PEG chains backbone and undercoordinated Ti atoms are established, whereas at {\sigma} = 2.25 chains/nm2 the conformation clearly becomes brush-like. Finally, we demonstrate that these spherical brushes, when immersed in water, but not in DCM, follow the Daoud-Cotton (DC) classical scaling model for the polymer volume fraction dependence with the distance from the center of star-shaped systems.

Graphene-h-BN hybrid nanostructures are grown in one step on the Pt(111) surface by ultra-high vacuum chemical vapor deposition using a single precursor, the dimethylamino borane complex. By varying the deposition conditions, different nanostructures ranging from a fully continuous hybrid monolayer to well-separated Janus nanodots can be obtained. The growth starts with heterogeneous nucleation on morphological defects such as Pt step edges and proceeds by the addition of small clusters formed by the decomposition of the dimethylamino borane complex. Scanning tunneling microscopy measurements indicate that a sharp zigzag in-plane boundary is formed when graphene grows aligned with the Pt substrate and consequently with the h-BN layer as well. When graphene is rotated by 30°, the graphene armchair edges are seamlessly connected to h-BN zigzag edges. This is confirmed by a thorough density functional theory (DFT) study. Angle resolved photoemission spectroscopy (ARPES) data suggests that both h-BN and graphene present the typical electronic structure of self-standing non-interacting materials.[Figure not available: see fulltext.].

In this work, we present an overview on how density functional theory calculations can be used to design novel electrocatalytic materials for fuel cells. In particular, we focus the attention on non-metal doped graphene systems, which were reported to present excellent performances as electrocatalysts for the oxygen reduction reaction (ORR) at the cathodic electrode of fuel cells and are, thus, considered promising substitutes of platinum or platinum alloys electrodes. The methodology, originally proposed by Nørskov et al. (J. Phys. Chem. B 2004, 108, 17886) for electrochemical processes at metal electrodes, is revisited and applied specifically to doped graphene. Finite molecular models of graphene are found to perform as well as periodic models for localized properties or reactions. Therefore, the sophisticated molecular quantum mechanics methodologies can be safely used to compute reliable Gibbs free energies of reaction in an aqueous environment for the various steps of reduction (at the cathode) or of oxidation (at the anode). Details of the reaction mechanisms and accurate cell onset- or over-potentials can be derived from the Gibbs free energy diagrams. The latter are computational quantities which can be directly compared to experimentally obtained cell overpotentials. Modeling electrocatalysis at fuel cells is, thus, an extremely powerful and convenient tool to improve our understanding of how fuel cells work and to design novel potentially active electrocatalytic materials. In this work, we present two specific applications of B-doped graphene, as electrocatalysts for the ORR at a half-cell cathode and for the methanol oxidation reaction (MOR) at a half-cell anode.

2D tungsten diselenide (2D-WSe2) is one of the most successful bidimensional materials for optoelectronic and photonic applications, thanks to its strong photoluminescence properties and to a characteristic large exciton binding energy. Although these optical properties are widely recognized by the scientific community, there is no general understanding of the atomistic details of the excitonic species giving rise to them. In this work, we present a density functional theory investigation of excitons in 2D-WSe2, where we compare results obtained by standard generalized gradient approximation (GGA) methods (including spin-orbit coupling) with those by hybrid density functionals. Our study provides information on the size of the self-trapped exciton, the number and type of atoms involved, the structural reorganization, the self-trapping energy, and the photoluminescence energy, whose computed value is in good agreement with experimental measurements in the literature. Moreover, based on the comparative analysis of the self-trapping energy for the exciton with that for isolated charge carriers (unbound electrons and holes), we also suggest a simplified approach for the theoretical estimation of the excitonic binding energy, which can be compared with previous estimates from different approaches or from experimental data.

Despite the many successful syntheses and applications of dopamine-functionalized TiO 2 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 2 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 2 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 3+ 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.

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.

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.

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.

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.

Understanding magnetism in defective graphene is paramount to improve and broaden its technological applications. A single vacancy in graphene is expected to lead to a magnetic moment with both a σ (1 μB) and a π (1 μB) component. Theoretical calculations based on standard LDA or GGA functional on periodic systems report a partial quenching of the π magnetization (0.5 μB) due to the crossing of two spin split bands at the Fermi level. In contrast, STS experiments ( Phys. Rev. Lett. 2016, 117, 166801) have recently proved the existence of two defect spin states that are separated in energy by 20–60 meV. In this work, we show that self-interaction corrected hybrid functional methods (B3LYP-D*) are capable of correctly reproducing this finite energy gap and, consequently, provide a π magnetization of 1 μB. The crucial role played by the exact exchange is highlighted by comparison with PBE-D2 results and by the magnetic moment dependence with the exact exchange portion in the functional used. The ground state ferromagnetic planar solution is compared to the antiferromagnetic and to the diamagnetic ones, which present an out-of-plane distortion. Periodic models are then compared to graphene nanoflakes of increasing size (up to C383H48). For large models, the triplet spin configuration (total magnetization 2 μB) is the most stable, independently of the functional used, which further corroborates the conclusions of this work and puts an end to the long-debated issue of the magnetic properties of an isolated C monovacancy in graphene

The nature of the interaction of water with the WO3 surface is of crucial importance for the use of this semiconductor oxide in photocatalysis. In this work, we investigate water adsorption and dissociation on both clean and O-deficient (001) WO3 surfaces by means of an accurate DFT approach. The O vacancy formation energy (computed with respect to O2) has been evaluated for all possible surface configurations, and the removal of the terminal O atom along the c axis is found to be preferred, costing about half the corresponding energy in the bulk. The presence of oxygen vacancies leads to a semiconductor to metal transition, confirming the experimental evidence of n-type conductivity in defective WO3 films. H2O preferably adsorbs on WO3 in a molecular undissociated form, due to the presence of W ions at the surface that act as Lewis acid sites. This interaction, about -1 eV per H2O molecule, is not very strong. Contrary to what is usually expected, the presence of oxygen vacancies does not significantly affect H2O adsorption. Finally, we investigated the H2O desorption from a hydroxylated surface. This suggests that the exposure of WO3 to H2 directly results in a hydroxylated surface and the corresponding H2O desorption turns out to be a very efficient mechanism to generate a reduced oxide surface, with important consequences on the electronic structure of this oxide.

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.

Nanoparticles of very small size (below 10 nm) of TiO 2 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 TiO 2 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.

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.

The mechanism of water photooxidation on TiO2 surfaces is still controversial. Here we report a first-principles density functional study based on a hybrid functional method in which an adsorbed water molecule is found to directly interact with a self-trapped hole at a bridging oxygen site and to transform into an OH(•) radical species through a concerted proton/hole transfer. This study analyzes both the thermodynamics and kinetics of this step of the reaction, which is generally considered to be the rate determining one. The fate of the OH(•) radical is then investigated in terms of its reactivity with different surface species, with a second OH(•) radical, or with a second water molecule coming from the environment. We find that OH(•) radicals can either acquire a hydrogen from surrounding water molecules or, if they meet, couple to form hydrogen peroxide with highly associated energy gain.

On-surface synthesis of thin organic and organometallic films in a bottom-up fashion has become a promising approach for the development of new nanotechnologies. In this work we studied 5,11-dibromotetracene (C18H10Br2) as a prototypical case of rodlike polyaromatic molecules functionalized with two bromine atoms on the sides. The adsorption and temperature-stimulated transformations of dibromotetracene assemblies on Ag(110) have been investigated by a combination of synchrotron radiation X-ray photoemission spectroscopy (XPS), near-edge X-ray absorption spectroscopy (NEXAFS), scanning tunneling microscopy (STM), and density functional theory (DFT) calculations. Upon the contact with the Ag substrate, the Br-C bonds are promptly cleaved at room temperature, and Ag-coordinated protopolymers are formed along the [001] substrate direction. The organometallic diners and trimers remain on the surface up to 523 K. The stabilization of the protopolymers is driven by the substrate anisotropy and weak interactions with nearby Br atoms. The short oligomers formed at elevated temperatures are weakly bounded to the substrate and desorb before covalent structures can be formed.

In biomedical applications, TiO2 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 TiO2 (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 [-(OCH2CH2)n], on one side, and of a short PEG dimer molecule [H(OCH2CH2)2OH], 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. (Chemical Equation Presented).

Tungsten trioxide (WO$_3$) is a paradigmatic electrochromic material, whose peculiar optical properties in the presence of oxygen vacancies or intercalated alkali atoms have been observed and investigated for a long time. In this paper we propose a rationalization of experiments based on first-principles calculations of optical and electrical properties of oxygen deficient (reduced) WO$_3$. Our approach is based on a parameter-free dielectric-dependent hybrid density functional methodology, used in combination with the charge transition levels formalism, for studying excitation mechanisms in the presence of defects. Our results indicate that oxygen vacancies lead to a different physics in $\gamma$-monoclinic WO$_3$, depending on the orientation of the W-O-W chain where the vacancy is created, thus evidencing strong anisotropic effects rooted in the peculiar structural properties of the original nondefective monoclinic cell. Different types of oxygen vacancies can hence be classified on the basis of the calculated ground state properties, electronic structure, and excitation/emission energies, giving a satisfactory explanation to a range of experimental observations made on oxygen deficient WO$_3$., Comment: Accepted for publication in J. Phys. Chem. C

Graphene/TiO2 nanocomposites are successfully applied in both photocatalysis and photovoltaics. The enhanced performances are attributed to their improved interfacial charged transfer and charge separation, which reduce the recombination rate of the photoexcited charge carriers. Here, it is shown that only density functional methods which provide corrections for the spurious self-interaction and for the van der Waals forces can correctly describe the electronic structure, the adhesion energy, and the atomic equilibrium distances. It is also proven that residual O atoms at the interface largely enhances the binding energy and causes further electronic states hybridization between G and TiO2, which is expected to favor interfacial electron transfers. Finally, evidence that electrons are preferentially trapped at subsurface layers of TiO2, while holes are preferentially delocalized on the G sheet, is provided. This opposite tendency is proposed to be at the basis of the reduced recombination leading to the observed improved outcomes in photocatalytic and photovoltaic applications. Separation of photoexcited e-/h+ pairs at the interface between graphene and (101) anatase TiO2 surface is investigated through dispersion corrected hybrid density functional calculations. Electrons are preferentially trapped at second layer Ti3+ sites while holes tend to fully delocalize in the carbon layer with relevant consequences for photocatalysis and photovoltaics

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 (C16H8Br2). 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.

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.

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. Boron found to matter: Advanced density functional studies on the enhancement of oxygen reactivity through non-metal doping and/or interfacing graphene with metal (copper) and semiconductor (TiO2) surfaces are reviewed. Oxidized boron species are formed in the graphene sheets. These are bound to the underlying support through direct covalent bonds, which enhances the intimate connection between materials with different chemical and physical properties, thus creating an interactive hybrid interface.

No abstract available]

International audience; In order to assign the defect responsible for the experimental electron paramagnetic resonance (EPR) signal with trigonal symmetry (T center), we have studied the properties of different paramagnetic centers in yttria-stabilized cubic zirconia by computing the EPR g-tensor from density functional perturbation theory. We have considered reduced vacancy-zirconium complexes and reduced Ti impurities. These first-principles calculations allow us to discard the experimental assignment of the T center to an extrinsic Ti3+ ion nearest neighbor to a single vacancy. Instead, the calculated EPR g tensors of both a Zr3+ or a Ti3+ ion at the center of a divacancy aligned along the directions are compatible with the experimental EPR signal. However, since the EPR signal of the T center is correlated experimentally with an optical absorption band at 370 nm, calculated optical excitations allow us to decide in favor of the Ti3+ divacancy complex.

The formation of Pd dimers on the surface of MgO has been studied by means of density functional theory (DFT) cluster model calculations. The following surface sites have been considered: regular five-coordinated anions at the (100) terraces, monoatomic steps, OH groups, and neutral vacancies (F centers). We discuss the energy balance of forming a dimer at a given site with respect to two isolated Pd atoms, one adsorbed at the defect and one on a regular terrace site. We found that all the defects considered lead to an energy gain when the dimer is formed, suggesting that they can be involved in nucleation and growth processes of metal clusters on the MgO surface. The dimerization energy is moderate for steps (approximate to0.8 eV), large for OH groups (approximate to1.3 eV) and rather small (

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

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.

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 ${\mathrm{Fe}}_{3}{\mathrm{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 ${\mathrm{Fe}}_{3}{\mathrm{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.