First principles interatomic potential for tungsten based on Gaussian process regression

An accurate description of atomic interactions, such as that provided by first principles quantum mechanics, is fundamental to realistic prediction of the properties that govern plasticity, fracture or crack propagation in metals. However, the computational complexity associated with modern schemes explicitly based on quantum mechanics limits their applications to systems of a few hundreds of atoms at most. This thesis investigates the application of the Gaussian Approximation Potential (GAP) scheme to atomistic modelling of tungsten - a bcc transition metal which exhibits a brittle-to-ductile transition and whose plasticity behaviour is controlled by the properties of $\frac{1}{2} \langle 111 \rangle$ screw dislocations. We apply Gaussian process regression to interpolate the quantum-mechanical (QM) potential energy surface from a set of points in atomic configuration space. Our training data is based on QM information that is computed directly using density functional theory (DFT). To perform the fitting, we represent atomic environments using a set of rotationally, permutationally and reflection invariant parameters which act as the independent variables in our equations of non-parametric, non-linear regression. We develop a protocol for generating GAP models capable of describing lattice defects in metals by building a series of interatomic potentials for tungsten. We then demonstrate that a GAP potential based on a Smooth Overlap of Atomic Positions (SOAP) covariance function provides a description of the $\frac{1}{2} \langle 111 \rangle$ screw dislocation that is in agreement with the DFT model. We use this potential to simulate the mobility of $\frac{1}{2} \langle 111 \rangle$ screw dislocations by computing the Peierls barrier and model dislocation-vacancy interactions to QM accuracy in a system containing more than 100,000 atoms.

Variance prediction in sea

A method is presented for predicting the variance of the energy levels in a built-up system to accompany the mean values predicted by SEA. Closed form expressions for the variance are obtained in terms of the standard SEA parameters and an additional set of parameters αk that describe the nature of the power input to each subsystem k, and αks that describe the nature of the coupling between subsystems k and s.

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.

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.

Nanoengineering coaxial carbon nanotube-dual-polymer heterostructures.

We describe studies of new nanostructured materials consisting of carbon nanotubes wrapped in sequential coatings of two different semiconducting polymers, namely, poly(3-hexylthiophene) (P3HT) and poly(9,9'-dioctylfluorene-co-benzothiadiazole) (F8BT). Using absorption spectroscopy and steady-state and ultrafast photoluminescence measurements, we demonstrate the role of the different layer structures in controlling energy levels and charge transfer in both solution and film samples. By varying the simple solution processing steps, we can control the ordering and proportions of the wrapping polymers in the solid state. The resulting novel coaxial structures open up a variety of new applications for nanotube blends and are particularly promising for implementation into organic photovoltaic devices. The carbon nanotube template can also be used to optimize both the electronic properties and morphology of polymer composites in a much more controlled fashion than achieved previously, offering a route to producing a new generation of polymer nanostructures.

Oxygen vacancy defects in Ta2O5 showing long-range atomic re-arrangements

The structure, formation energy, and energy levels of the various oxygen vacancies in Ta2O5 have been calculated using the λ phase model. The intra-layer vacancies give rise to unusual, long-range bonding rearrangements, which are different for each defect charge state. The 2-fold coordinated intra-layer vacancy is the lowest cost vacancy and forms a deep level 1.5 eV below the conduction band edge. The 3-fold intra-layer vacancy and the 2-fold inter-layer vacancy are higher cost defects, and form shallower levels. The unusual bonding rearrangements lead to low oxygen migration barriers, which are useful for resistive random access memory applications. © 2014 AIP Publishing LLC.