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.

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 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.

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.

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.

On the relationship between hidden Markov Models and convex functional transforms for simulating scanning probe Microscopy

Models for simulating Scanning Probe Microscopy (SPM) may serve as a reference point for validating experimental data and practice. Generally, simulations use a microscopic model of the sample-probe interaction based on a first-principles approach, or a geometric model of macroscopic distortions due to the probe geometry. Examples of the latter include use of neural networks, the Legendre Transform, and dilation/erosion transforms from mathematical morphology. Dilation and the Legendre Transform fall within a general family of functional transforms, which distort a function by imposing a convex solution.In earlier work, the authors proposed a generalized approach to modeling SPM using a hidden Markov model, wherein both the sample-probe interaction and probe geometry may be taken into account. We present a discussion of the hidden Markov model and its relationship to these convex functional transforms for simulating and restoring SPM images.©2009 SPIE.

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.

Shock Propagation and MPT Noise From a Transonic Rotor in Nonuniform Flow

One of the major challenges in high-speed fan stages used in compact, embedded propulsion systems is inlet distortion noise. A body-force-based approach for the prediction of multiple-pure-tone (MPT) noise was previously introduced and validated. In this paper, it is employed with the objective of quantifying the effects of nonuniform flow on the generation and propagation of MPT noise. First-of-their-kind back-to-back coupled aero-acoustic computations were carried out using the new approach for conventional and serpentine inlets. Both inlets delivered flow to the same NASA/GE R4 fan rotor at equal corrected mass flow rates. Although the source strength at the fan is increased by 38 dB in sound power level due to the nonuniform inflow, far-field noise for the serpentine inlet duct is increased on average by only 3.1 dBA overall sound pressure level in the forward arc. This is due to the redistribution of acoustic energy to frequencies below 11 times the shaft frequency and the apparent cut-off of tones at higher frequencies including blade-passing tones. The circumferential extent of the inlet swirl distortion at the fan was found to be two blade pitches, or 1/11th of the circumference, suggesting a relationship between the circumferential extent of the inlet distortion and the apparent cut-off frequency perceived in the far field. A first-principles-based model of the generation of shock waves from a transonic rotor in nonuniform flow showed that the effects of nonuniform flow on acoustic wave propagation, which cannot be captured by the simplified model, are more dominant than those of inlet flow distortion on source noise. It demonstrated that nonlinear, coupled aerodynamic and aero-acoustic computations, such as those presented in this paper, are necessary to assess the propagation through nonuniform mean flow. A parametric study of serpentine inlet designs is underway to quantify these propagation effects. © 2013 American Society of Mechanical Engineers.

Shock propagation and MPT noise from a transonic rotor in non-uniform flow

One of the major challenges in hig4h-speed fan stages used in compact, embedded propulsion systems is inlet distortion noise. A body-force-based approach for the prediction of multiple-pure-tone (MPT) noise was previously introduced and validated. In this paper, it is employed with the objective of quantifying the effects of non-uniform flow on the generation and propagation of MPT noise. First-of-their-kind back-to-back coupled aero-acoustic computations were carried out using the new approach for conventional and serpentine inlets. Both inlets delivered flow to the same NASA/GE R4 fan rotor at equal corrected mass flow rates. Although the source strength at the fan is increased by 45 dB in sound power level due to the non-uniform inflow, farfield noise for the serpentine inlet duct is increased on average by only 3.1 dBA overall sound pressure level in the forward arc. This is due to the redistribution of acoustic energy to frequencies below 11 times the shaft frequency and the apparent cut-off of tones at higher frequencies including blade-passing tones. The circumferential extent of the inlet swirl distortion at the fan was found to be 2 blade pitches, or 1/11th of the circumference, suggesting a relationship between the circumferential extent of the inlet distortion and the apparent cut-off frequency perceived in the far field. A first-principles-based model of the generation of shock waves from a transonic rotor in non-uniform flow showed that the effects of non-uniform flow on acoustic wave propagation, which cannot be captured by the simplified model, are more dominant than those of inlet flow distortion on source noise. It demonstrated that non-linear, coupled aerodynamic and aeroacoustic computations, such as those presented in this paper, are necessary to assess the propagation through non-uniform mean flow. A parametric study of serpentine inlet designs is underway to quantify these propagation effects. Copyright © 2011 by ASME.

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.