First-principles prediction of doped graphane as a high-temperature electron-phonon superconductor.

We predict by first-principles calculations that p-doped graphane is an electron-phonon superconductor with a critical temperature above the boiling point of liquid nitrogen. The unique strength of the chemical bonds between carbon atoms and the large density of electronic states at the Fermi energy arising from the reduced dimensionality give rise to a giant Kohn anomaly in the optical phonon dispersions and push the superconducting critical temperature above 90 K. As evidence of graphane was recently reported, and doping of related materials such as graphene, diamond, and carbon nanostructures is well established, superconducting graphane may be feasible.

First-principles energetics of water clusters and ice: a many-body analysis.

Standard forms of density-functional theory (DFT) have good predictive power for many materials, but are not yet fully satisfactory for cluster, solid, and liquid forms of water. Recent work has stressed the importance of DFT errors in describing dispersion, but we note that errors in other parts of the energy may also contribute. We obtain information about the nature of DFT errors by using a many-body separation of the total energy into its 1-body, 2-body, and beyond-2-body components to analyze the deficiencies of the popular PBE and BLYP approximations for the energetics of water clusters and ice structures. The errors of these approximations are computed by using accurate benchmark energies from the coupled-cluster technique of molecular quantum chemistry and from quantum Monte Carlo calculations. The systems studied are isomers of the water hexamer cluster, the crystal structures Ih, II, XV, and VIII of ice, and two clusters extracted from ice VIII. For the binding energies of these systems, we use the machine-learning technique of Gaussian Approximation Potentials to correct successively for 1-body and 2-body errors of the DFT approximations. We find that even after correction for these errors, substantial beyond-2-body errors remain. The characteristics of the 2-body and beyond-2-body errors of PBE are completely different from those of BLYP, but the errors of both approximations disfavor the close approach of non-hydrogen-bonded monomers. We note the possible relevance of our findings to the understanding of liquid water.

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.

First-principles energetics of water: a many-body analysis

Standard forms of density-functional theory (DFT) have good predictive power for many materials, but are not yet fully satisfactory for solid, liquid and cluster forms of water. We use a many-body separation of the total energy into its 1-body, 2-body (2B) and beyond-2-body (B2B) components to analyze the deficiencies of two popular DFT approximations. We show how machine-learning methods make this analysis possible for ice structures as well as for water clusters. We find that the crucial energy balance between compact and extended geometries can be distorted by 2B and B2B errors, and that both types of first-principles error are important.

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.

Shifting atomic patterns: on the origin of the different atomic-scale patterns of graphite as observed using scanning tunnelling microscopy.

We present an in-depth study of the myriad atomically resolved patterns observed on graphite using the scanning tunnelling microscope (STM) over the past three decades. Through the use of highly resolved atomic resolution images, we demonstrate how the interactions between the different graphene layers comprising graphite affect the local surface atomic charge density and its resulting symmetry orientation, with particular emphasis on interactions that are thermodynamically unstable. Moreover, the interlayer graphene coupling is controlled experimentally by varying the tip-surface interaction, leading to associated changes in the atomic patterns. The images are corroborated by first-principles calculations, further validating our claim that surface graphene displacement, coming both from lateral and vertical displacement of the top graphene layer, forms the basis of the rich variety of atomic patterns observed in STM experiments on graphite.

Growth of straight InAs-on-GaAs nanowire heterostructures.

One of the main motivations for the great interest in semiconductor nanowires is the possibility of easily growing advanced heterostructures that might be difficult or even impossible to achieve in thin films. For III-V semiconductor nanowires, axial heterostructures with an interchange of the group III element typically grow straight in only one interface direction. In the case of InAs-GaAs heterostructures, straight nanowire growth has been demonstrated for growth of GaAs on top of InAs, but so far never in the other direction. In this article, we demonstrate the growth of straight axial heterostructures of InAs on top of GaAs. The heterostructure interface is sharp and we observe a dependence on growth parameters closely related to crystal structure as well as a diameter dependence on straight nanowire growth. The results are discussed by means of accurate first principles calculations of the interfacial energies. In addition, the role of the gold seed particle, the effect of its composition at different stages during growth, and its size are discussed in relation to the results observed.

Tunable Rashba effect in two-dimensional LaOBiS2 films: Ultrathin candidates for spin field effect transistors

Rashba spin splitting is a two-dimensional (2D) relativistic effect closely related to spintronics. However, so far there is no pristine 2D material to exhibit enough Rashba splitting for the fabrication of ultrathin spintronic devices, such as spin field effect transistors (SFET). On the basis of first-principles calculations, we predict that the stable 2D LaOBiS2 with only 1 nm of thickness can produce remarkable Rashba spin splitting with a magnitude of 100 meV. Because the medium La2O2 layer produces a strong polar field and acts as a blocking barrier, two counter-helical Rashba spin polarizations are localized at different BiS 2 layers. The Rashba parameter can be effectively tuned by the intrinsic strain, while the bandgap and the helical direction of spin states sensitively depends on the external electric field. We propose an advanced Datta-Das SFET model that consists of dual gates and 2D LaOBiS2 channels by selecting different Rashba states to achieve the on-off switch via electric fields. © 2013 American Chemical Society.

Electrical control of nanoscale functionalization in graphene by the scanning probe technique

Functionalized graphene is a versatile material that has well-known physical and chemical properties depending on functional groups and their coverage. However, selective control of functional groups on the nanoscale is hardly achievable by conventional methods utilizing chemical modifications. We demonstrate electrical control of nanoscale functionalization of graphene with the desired chemical coverage of a selective functional group by atomic force microscopy (AFM) lithography and their full recovery through moderate thermal treatments. Surprisingly, our controlled coverage of functional groups can reach 94.9% for oxygen and 49.0% for hydrogen, respectively, well beyond those achieved by conventional methods. This coverage is almost at the theoretical maximum, which is verified through scanning photoelectron microscope measurements as well as first-principles calculations. We believe that the present method is now ready to realize 'chemical pencil drawing' of atomically defined circuit devices on top of a monolayer of graphene. © 2014 Nature Publishing Group All rights reserved.

Accuracy and transferability of GAP models for tungsten

We introduce interatomic potentials for tungsten in the bcc crystal phase and its defects within the Gaussian Approximation Potential (GAP) framework, fitted to a database of first principles density functional theory (DFT) calculations. We investigate the performance of a sequence of models based on databases of increasing coverage in configuration space and showcase our strategy of choosing representative small unit cells to train models that predict properties only observable using thousands of atoms. The most comprehensive model is then used to calculate properties of the screw dislocation, including its structure, the Peierls barrier and the energetics of the vacancy-dislocation interaction. All software and raw data are available at www.libatoms.org.

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.

Nucleated polymerization with secondary pathways. I. Time evolution of the principal moments.

Self-assembly processes resulting in linear structures are often observed in molecular biology, and include the formation of functional filaments such as actin and tubulin, as well as generally dysfunctional ones such as amyloid aggregates. Although the basic kinetic equations describing these phenomena are well-established, it has proved to be challenging, due to their non-linear nature, to derive solutions to these equations except for special cases. The availability of general analytical solutions provides a route for determining the rates of molecular level processes from the analysis of macroscopic experimental measurements of the growth kinetics, in addition to the phenomenological parameters, such as lag times and maximal growth rates that are already obtainable from standard fitting procedures. We describe here an analytical approach based on fixed-point analysis, which provides self-consistent solutions for the growth of filamentous structures that can, in addition to elongation, undergo internal fracturing and monomer-dependent nucleation as mechanisms for generating new free ends acting as growth sites. Our results generalise the analytical expression for sigmoidal growth kinetics from the Oosawa theory for nucleated polymerisation to the case of fragmenting filaments. We determine the corresponding growth laws in closed form and derive from first principles a number of relationships which have been empirically established for the kinetics of the self-assembly of amyloid fibrils.