Damien J. Carter

Benchmarking Calculated Lattice Parameters and Energies of Molecular Crystals Using van der Waals Density Functionals.

The development of new functionals and methods to accurately describe van der Waals forces in density functional theory (DFT) has become popular in recent years, with the vast majority of studies assessing the accuracy of the energetics of collections of molecules, and to a lesser extent molecular crystalline systems. As the energies are a function of the atom positions, we assess the accuracy of DFT calculations from both a geometric and energetics point of view for the C21 reference data set of Otero-de-la-Roza and Johnson for molecular crystals, and a set of monosaccharide molecular crystals. In particular, we examine the performance of exchange-correlation functionals designed to handle van der Waals forces, including the vdW-DF, vdW-DF2, and XDM methods. We also assess the effect of using small and large basis sets, the choice of basis functions (local atomic orbitals using the SIESTA code versus planewaves using the Quantum ESPRESSO code), and the effect of corrections for basis set superposition errors. Finally, we examine the geometries and energies of the S22 reference set of molecular complexes. Overall, the most accurate geometries for both choices of basis functions are obtained with the vdW-DF2 functional, while the most accurate lattice energies are obtained using vdW-DF2 with local atomic orbitals and XDM with planewaves with mean absolute errors of less than 4 kJ/mol.

Noncovalent Interactions in SIESTA Using the vdW-DF Functional: S22 Benchmark and Macrocyclic Structures.

We investigate the performance of the vdW-DF functional of Dion et al. implemented in the SIESTA code. In particular, the S22 data set and several calixarene-based host-guest structures are examined to assess the performance of the functional. The binding energy error statistics for the S22 data set reveal that the vdW-DF functional performs very well when compared to a range of other methods of treating dispersion in density functional theory, and to vdW-DF implementations in other codes. For the calixarene host-guest structures, the structural properties and binding energies are compared to previous experimental and computational studies, and in most cases we find that vdW-DF provides superior results to other computational studies.

Phosphorus δ-doped silicon: mixed-atom pseudopotentials and dopant disorder effects

Within a full density functional theory framework we calculate the band structure and doping potential for phosphorus δ-doped silicon. We compare two different representations of the dopant plane; pseudo-atoms in which the nuclear charge is fractional between silicon and phosphorus, and explicit arrangements employing distinct silicon and phosphorus atoms. While the pseudo-atom approach offers several computational advantages, the explicit model calculations differ in a number of key points, including the valley splitting, the Fermi level and the width of the doping potential. These findings have implications for parameters used in device modelling.

Bottom-up assembly of metallic germanium

Extending chip performance beyond current limits of miniaturisation requires new materials and functionalities that integrate well with the silicon platform. Germanium fits these requirements and has been proposed as a high-mobility channel material, a light emitting medium in silicon-integrated lasers, and a plasmonic conductor for bio-sensing. Common to these diverse applications is the need for homogeneous, high electron densities in three-dimensions (3D). Here we use a bottom-up approach to demonstrate the 3D assembly of atomically sharp doping profiles in germanium by a repeated stacking of two-dimensional (2D) high-density phosphorus layers. This produces high-density (1019 to 1020 cm−3) low-resistivity (10−4Ω · cm) metallic germanium of precisely defined thickness, beyond the capabilities of diffusion-based doping technologies. We demonstrate that free electrons from distinct 2D dopant layers coalesce into a homogeneous 3D conductor using anisotropic quantum interference measurements, atom probe tomography, and density functional theory.

Quantum confinement effects in gallium nitride nanostructures: ab initio investigations

We present ab initio density functional investigations of the atomic and electronic structure of gallium nitride nanodots and nanowires. With increasing diameter, the average Ga-N bond length in the nanostructures increases, as does the relative stability (heat of formation), approaching the values for bulk GaN. As the diameter decreases, the band gap increases, with the variation for the nanodots greater than that for the nanowires, in qualitative accordance with expectations based on simple geometrical quantum confinement considerations. Interestingly, in contrast to nanowires, the lowest unoccupied states of the nanodots exhibit an extended delocalized (Ga-derived) character, weighted in the centre of the nanodot.

Incorporation of cyano transition metal complexes in KCl crystals: Experimental and computational studies

Experimental and computational studies of the incorporation of hexacyanoferrate(II), hexacyanocobaltate(III), and hexacyanoferrate(III) into potassium chloride crystals are described. The experimental results showed that the extent of incorporation follows the trend hexacyanoferrate(II) » hexacyanoferrate(III) > hexacyanocobaltate(III). Computational modelling produced replacement energies that match the experimental trend. The calculated geometry of the incorporated complexes was also found to match well with previous experimental results.

Valley splitting in a silicon quantum device platform.

By suppressing an undesirable surface Umklapp process, it is possible to resolve the two most occupied states (1Γ and 2Γ) in a buried two-dimensional electron gas (2DEG) in silicon. The 2DEG exists because of an atomically sharp profile of phosphorus dopants which have been formed beneath the Si(001) surface (a δ-layer). The energy separation, or valley splitting, of the two most occupied bands has critical implications for the properties of δ-layer derived devices, yet until now, has not been directly measurable. Density functional theory (DFT) allows the 2DEG band structure to be calculated, but without experimental verification the size of the valley splitting has been unclear. Using a combination of direct spectroscopic measurements and DFT we show that the measured band structure is in good qualitative agreement with calculations and reveal a valley splitting of 132 ± 5 meV. We also report the effective mass and occupation of the 2DEG states and compare the dispersions and Fermi surface with DFT.

van der Waals corrected density functional calculations of the adsorption of benzene on the Cu (111) surface

We investigate the performance of several van der Waals (vdW) functionals at calculating the interactions between benzene and the copper (111) surface, using the local orbital approach in the SIESTA code. We demonstrate the importance of using surface optimized basis sets to calculate properties of pure surfaces, including surface energies and the work function. We quantify the errors created using (3 × 3) supercells to study adsorbate interactions using much larger supercells, and show non‐negligible errors in the binding energies and separation distances. We examine the eight high‐symmetry orientations of benzene on the Cu (111) surface, reporting the binding energies, separation distance, and change in work function. The optimized vdW‐DF(optB88‐vdW) functional provides superior results to the vdW‐DF(revPBE) and vdW‐DF2(rPW86) functionals, and closely matches the experimental and experimentally deduced values. This work demonstrates that local orbital methods using appropriate basis sets combined with a vdW functional can model adsorption between metal surfaces and organic molecules.

Chemical evolution via beta decay: a case study in strontium-90

Using (90)Sr as a representative isotope, we present a framework for understanding beta decay within the solid state. We quantify three key physical and chemical principles, namely momentum-induced recoil during the decay event, defect creation due to physical displacement, and chemical evolution over time. A fourth effect, that of electronic excitation, is also discussed, but this is difficult to quantify and is strongly material dependent. The analysis is presented for the specific cases of SrTiO(3) and SrH(2). By comparing the recoil energy with available threshold displacement data we show that in many beta-decay situations defects such as Frenkel pairs will not be created during decay as the energy transfer is too low. This observation leads to the concept of chemical evolution over time, which we quantify using density functional theory. Using a combination of Bader analysis, phonon calculations and cohesive energy calculations, we show that beta decay leads to counter-intuitive behavior that has implications for nuclear waste storage and novel materials design.

Silver(I), gold(I) and palladium(II) complexes of a NHC-pincer ligand with an aminotriazine core: a comparison with pyridyl analogues

Dinuclear silver, di- and tetra-nuclear gold, and mononuclear palladium complexes with chelating C,N,C diethylaminotriazinyl-bridged bis(NHC) pincer ligands were prepared and characterised. The silver and gold complexes exist in a twisted, helical conformation in both the solution- and the solid state. In contrast, an analogous dinuclear gold complex with pyridyl-bridged NHCs exists in a linear conformation. Computational studies have been performed to rationalise the formation of twisted/helical vs. linear forms.