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.

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.

Calculation of metallic and insulating phases of V2O3 by hybrid density functionals.

The electronic structure of vanadium sesquioxide V2O3 in its different phases has been calculated using the screened exchange hybrid density functional. The hybrid functional accurately reproduces the experimental electronic properties of all three phases, the paramagnetic metal (PM) phase, the anti-ferromagnetic insulating phase, and the Cr-doped paramagnetic insulating (PI) phase. We find that a fully relaxed supercell model of the Cr-doped PI phase based on the corundum structure has a monoclinic-like local strain around the substitutional Cr atoms. This is found to drive the PI-PM transition, consistent with a Peierls-Mott transition. The PI phase has a calculated band gap of 0.15 eV, in good agreement with experiment.

Development of a classical force field for the oxidized Si surface: application to hydrophilic wafer bonding.

We have developed a classical two- and three-body interaction potential to simulate the hydroxylated, natively oxidized Si surface in contact with water solutions, based on the combination and extension of the Stillinger-Weber potential and of a potential originally developed to simulate SiO(2) polymorphs. The potential parameters are chosen to reproduce the structure, charge distribution, tensile surface stress, and interactions with single water molecules of a natively oxidized Si surface model previously obtained by means of accurate density functional theory simulations. We have applied the potential to the case of hydrophilic silicon wafer bonding at room temperature, revealing maximum room temperature work of adhesion values for natively oxidized and amorphous silica surfaces of 97 and 90 mJm(2), respectively, at a water adsorption coverage of approximately 1 ML. The difference arises from the stronger interaction of the natively oxidized surface with liquid water, resulting in a higher heat of immersion (203 vs 166 mJm(2)), and may be explained in terms of the more pronounced water structuring close to the surface in alternating layers of larger and smaller densities with respect to the liquid bulk. The computed force-displacement bonding curves may be a useful input for cohesive zone models where both the topographic details of the surfaces and the dependence of the attractive force on the initial surface separation and wetting can be taken into account.

Graphite intercalation compounds under pressure

Motivated by recent experimental work, we use first-principles density functional theory methods to conduct an extensive search for low enthalpy structures of C$_6$Ca under pressure. As well as a range of buckled structures, which are energetically competitive over an intermediate range of pressures, we show that the high pressure system ($\gtrsim 18$ GPa) is unstable towards the formation of a novel class of layered structures, with the most stable compound involving carbon sheets containing five- and eight-membered rings. As well as discussing the energetics of the different classes of low enthalpy structures, we comment on the electronic structure of the high pressure compound and its implications for superconductivity.

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.

Gaussian Approximation Potential: an interatomic potential derived from first principles Quantum Mechanics

Simulation of materials at the atomistic level is an important tool in studying microscopic structure and processes. The atomic interactions necessary for the simulation are correctly described by Quantum Mechanics. However, the computational resources required to solve the quantum mechanical equations limits the use of Quantum Mechanics at most to a few hundreds of atoms and only to a small fraction of the available configurational space. This thesis presents the results of my research on the development of a new interatomic potential generation scheme, which we refer to as Gaussian Approximation Potentials. In our framework, the quantum mechanical potential energy surface is interpolated between a set of predetermined values at different points in atomic configurational space by a non-linear, non-parametric regression method, the Gaussian Process. To perform the fitting, we represent the atomic environments by the bispectrum, which is invariant to permutations of the atoms in the neighbourhood and to global rotations. The result is a general scheme, that allows one to generate interatomic potentials based on arbitrary quantum mechanical data. We built a series of Gaussian Approximation Potentials using data obtained from Density Functional Theory and tested the capabilities of the method. We showed that our models reproduce the quantum mechanical potential energy surface remarkably well for the group IV semiconductors, iron and gallium nitride. Our potentials, while maintaining quantum mechanical accuracy, are several orders of magnitude faster than Quantum Mechanical methods.