Electronic and magnetic properties of Ti(2)O(3), Cr(2)O(3), and Fe(2)O(3) calculated by the screened exchange hybrid density functional.

The electronic and magnetic properties of the transition metal sesqui-oxides Cr(2)O(3), Ti(2)O(3), and Fe(2)O(3) have been calculated using the screened exchange (sX) hybrid density functional. This functional is found to give a band structure, bandgap, and magnetic moment in better agreement with experiment than the local density approximation (LDA) or the LDA+U methods. Ti(2)O(3) is found to be a spin-paired insulator with a bandgap of 0.22 eV in the Ti d orbitals. Cr(2)O(3) in its anti-ferromagnetic phase is an intermediate charge transfer Mott-Hubbard insulator with an indirect bandgap of 3.31 eV. Fe(2)O(3), with anti-ferromagnetic order, is found to be a wide bandgap charge transfer semiconductor with a 2.41 eV gap. Interestingly sX outperforms the HSE functional for the bandgaps of these oxides.

A hybrid density functional view of native vacancies in gallium nitride.

We investigated the transition energy levels of the vacancy defects in gallium nitride by means of a hybrid density functional theory approach (DFT). We show that, in contrast to predictions from a recent study on the level of purely local DFT, the inclusion of screened exchange stabilizes the triply positive charge state of the nitrogen vacancy for Fermi energies close to the valence band. On the other hand, the defect levels associated with the negative charge states of the nitrogen vacancy hybridize with the conduction band and turn out to be energetically unfavorable, except for high n-doping. For the gallium vacancy, the increased magnetic splitting between up-spin and down-spin bands due to stronger exchange interactions in sX-LDA pushes the defect levels deeper into the band gap and significantly increases the associated charge transition levels. Based on these results, we propose the ϵ(0| - 1) transition level as an alternative candidate for the yellow luminescence in GaN.

Between scylla and charybdis: Hydrophobic graphene-guided water diffusion on hydrophilic substrates

The structure of water confined in nanometer-sized cavities is important because, at this scale, a large fraction of hydrogen bonds can be perturbed by interaction with the confining walls. Unusual fluidity properties can thus be expected in the narrow pores, leading to new phenomena like the enhanced fluidity reported in carbon nanotubes. Crystalline mica and amorphous silicon dioxide are hydrophilic substrates that strongly adsorb water. Graphene, on the other hand, interacts weakly with water. This presents the question as to what determines the structure and diffusivity of water when intercalated between hydrophilic substrates and hydrophobic graphene. Using atomic force microscopy, we have found that while the hydrophilic substrates determine the structure of water near its surface, graphene guides its diffusion, favouring growth of intercalated water domains along the C-C bond zigzag direction. Molecular dynamics and density functional calculations are provided to help understand the highly anisotropic water stripe patterns observed.

QM/MM simulation of liquid water with an adaptive quantum region.

The simulation of complex chemical systems often requires a multi-level description, in which a region of special interest is treated using a computationally expensive quantum mechanical (QM) model while its environment is described by a faster, simpler molecular mechanical (MM) model. Furthermore, studying dynamic effects in solvated systems or bio-molecules requires a variable definition of the two regions, so that atoms or molecules can be dynamically re-assigned between the QM and MM descriptions during the course of the simulation. Such reassignments pose a problem for traditional QM/MM schemes by exacerbating the errors that stem from switching the model at the boundary. Here we show that stable, long adaptive simulations can be carried out using density functional theory with the BLYP exchange-correlation functional for the QM model and a flexible TIP3P force field for the MM model without requiring adjustments of either. Using a primary benchmark system of pure water, we investigate the convergence of the liquid structure with the size of the QM region, and demonstrate that by using a sufficiently large QM region (with radius 6 Å) it is possible to obtain radial and angular distributions that, in the QM region, match the results of fully quantum mechanical calculations with periodic boundary conditions, and, after a smooth transition, also agree with fully MM calculations in the MM region. The key ingredient is the accurate evaluation of forces in the QM subsystem which we achieve by including an extended buffer region in the QM calculations. We also show that our buffered-force QM/MM scheme is transferable by simulating the solvated Cl(-) ion.

Channel selective tunnelling through a nanographene assembly.

We report selective tunnelling through a nanographene intermolecular tunnel junction achieved via scanning tunnelling microscope tip functionalization with hexa-peri-hexabenzocoronene (HBC) molecules. This leads to an offset in the alignment between the energy levels of the tip and the molecular assembly, resulting in the imaging of a variety of distinct charge density patterns in the HBC assembly, not attainable using a bare metallic tip. Different tunnelling channels can be selected by the application of an electric field in the tunnelling junction, which changes the condition of the HBC on the tip. Density functional theory-based calculations relate the imaged HBC patterns to the calculated molecular orbitals at certain energy levels. These patterns bear a close resemblance to the π-orbital states of the HBC molecule calculated at the relevant energy levels, mainly below the Fermi energy of HBC. This correlation demonstrates the ability of an HBC functionalized tip as regards accessing an energy range that is restricted to the usual operating bias range around the Fermi energy with a normal metallic tip at room temperature. Apart from relating to molecular orbitals, some patterns could also be described in association with the Clar aromatic sextet formula. Our observations may help pave the way towards the possibility of controlling charge transport between organic interfaces.

Calculation of point defects in rutile TiO 2 by the screened-exchange hybrid functional

The formation energies of the oxygen vacancy and titanium interstitial in rutile TiO 2 were calculated by the screened-exchange (sX) hybrid density functional method, which gives a band gap of 3.1 eV, close to the experimental value. The oxygen vacancy gives rise to a gap state lying 0.7 eV below the conduction band edge, whose charge density is localized around the two of three Ti atoms next to the vacancy. The Ti interstitial (Ti int) generates four defect states in the gap, whose unpaired electrons lie on the interstitial and the adjacent Ti 3d orbitals. The formation energy for the neutral oxygen vacancy is 1.9 eV for the O-poor chemical potential. The neutral Ti interstitial has a lower formation energy than the O vacancy under O-poor conditions. This indicates that both the O vacancy and Ti int are relevant for oxygen deficiency in rutile TiO 2 but the O vacancy will dominate under O-rich conditions. This resolves questions about defect localization and defect predominance in the literature. © 2012 American Physical Society.

Accurate screened exchange band structures for the transition metal monoxides MnO, FeO, CoO and NiO.

We report calculations of the band structures and density of states of the four transition metal monoxides MnO, FeO, CoO and NiO using the hybrid density functional sX-LDA ('screened exchange local density approximation'). Late transition metal oxides are prototypical examples of strongly correlated materials, which pose challenges for electronic structure methods. We compare our results with available experimental data and show that our calculations generally yield accurate predictions for the fundamental band gaps and valence bands, in favourable agreement with previously reported theoretical studies. For MnO, the band gaps are still underestimated, suggesting additional many-body effects that are not captured by our screened hybrid functional approach.

Communication: energy benchmarking with quantum Monte Carlo for water nano-droplets and bulk liquid water.

We show the feasibility of using quantum Monte Carlo (QMC) to compute benchmark energies for configuration samples of thermal-equilibrium water clusters and the bulk liquid containing up to 64 molecules. Evidence that the accuracy of these benchmarks approaches that of basis-set converged coupled-cluster calculations is noted. We illustrate the usefulness of the benchmarks by using them to analyze the errors of the popular BLYP approximation of density functional theory (DFT). The results indicate the possibility of using QMC as a routine tool for analyzing DFT errors for non-covalent bonding in many types of condensed-phase molecular system.