Efficient sampling of atomic configurational spaces.

We describe a method to explore the configurational phase space of chemical systems. It is based on the nested sampling algorithm recently proposed by Skilling (AIP Conf. Proc. 2004, 395; J. Bayesian Anal. 2006, 1, 833) and allows us to explore the entire potential energy surface (PES) efficiently in an unbiased way. The algorithm has two parameters which directly control the trade-off between the resolution with which the space is explored and the computational cost. We demonstrate the use of nested sampling on Lennard-Jones (LJ) clusters. Nested sampling provides a straightforward approximation for the partition function; thus, evaluating expectation values of arbitrary smooth operators at arbitrary temperatures becomes a simple postprocessing step. Access to absolute free energies allows us to determine the temperature-density phase diagram for LJ cluster stability. Even for relatively small clusters, the efficiency gain over parallel tempering in calculating the heat capacity is an order of magnitude or more. Furthermore, by analyzing the topology of the resulting samples, we are able to visualize the PES in a new and illuminating way. We identify a discretely valued order parameter with basins and suprabasins of the PES, allowing a straightforward and unambiguous definition of macroscopic states of an atomistic system and the evaluation of the associated free energies.

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

The environmental and economic sustainability of potential bioethanol from willow in the UK.

Life cycle assessment has been used to investigate the environmental and economic sustainability of a potential operation in the UK in which bioethanol is produced from the hydrolysis and subsequent fermentation of coppice willow. If the willow were grown on idle arable land in the UK, or, indeed, in Eastern Europe and imported as wood chips into the UK, it was found that savings of greenhouse gas emissions of 70-90%, when compared to fossil-derived gasoline on an energy basis, would be possible. The process would be energetically self-sufficient, as the co-products, e.g. lignin and unfermented sugars, could be used to produce the process heat and electricity, with surplus electricity being exported to the National Grid. Despite the environmental benefits, the economic viability is doubtful at present. However, the cost of production could be reduced significantly if the willow were altered by breeding to improve its suitability for hydrolysis and fermentation.

Absence of branches from xylan in Arabidopsis gux mutants reveals potential for simplification of lignocellulosic biomass.

As one of the most abundant polysaccharides on Earth, xylan will provide more than a third of the sugars for lignocellulosic biofuel production when using grass or hardwood feedstocks. Xylan is characterized by a linear β(1,4)-linked backbone of xylosyl residues substituted by glucuronic acid, 4-O-methylglucuronic acid or arabinose, depending on plant species and cell types. The biological role of these decorations is unclear, but they have a major influence on the properties of the polysaccharide. Despite the recent isolation of several mutants with reduced backbone, the mechanisms of xylan synthesis and substitution are unclear. We identified two Golgi-localized putative glycosyltransferases, GlucUronic acid substitution of Xylan (GUX)-1 and GUX2 that are required for the addition of both glucuronic acid and 4-O-methylglucuronic acid branches to xylan in Arabidopsis stem cell walls. The gux1 gux2 double mutants show loss of xylan glucuronyltransferase activity and lack almost all detectable xylan substitution. Unexpectedly, they show no change in xylan backbone quantity, indicating that backbone synthesis and substitution can be uncoupled. Although the stems are weakened, the xylem vessels are not collapsed, and the plants grow to normal size. The xylan in these plants shows improved extractability from the cell wall, is composed of a single monosaccharide, and requires fewer enzymes for complete hydrolysis. These findings have implications for our understanding of the synthesis and function of xylan in plants. The results also demonstrate the potential for manipulating and simplifying the structure of xylan to improve the properties of lignocellulose for bioenergy and other uses.

Transposon mutagenesis in a hyper-invasive clinical isolate of Campylobacter jejuni reveals a number of genes with potential roles in invasion.

Transposon mutagenesis has been applied to a hyper-invasive clinical isolate of Campylobacter jejuni, 01/51. A random transposon mutant library was screened in an in vitro assay of invasion and 26 mutants with a significant reduction in invasion were identified. Given that the invasion potential of C. jejuni is relatively poor compared to other enteric pathogens, the use of a hyper-invasive strain was advantageous as it greatly facilitated the identification of mutants with reduced invasion. The location of the transposon insertion in 23 of these mutants has been determined; all but three of the insertions are in genes also present in the genome-sequenced strain NCTC 11168. Eight of the mutants contain transposon insertions in one region of the genome (approximately 14 kb), which when compared with the genome of NCTC 11168 overlaps with one of the previously reported plasticity regions and is likely to be involved in genomic variation between strains. Further characterization of one of the mutants within this region has identified a gene that might be involved in adhesion to host cells.

Surface atomic properties of tetrahedral amorphous carbon

The surface energy and surface atomic structure of tetrahedral amorphous carbon has been calculated by an ab-initio method. The surface atoms are found to reconstruct into sp2 sites often bonded in graphitic rings. Placing the dangling bonds on adjacent surface atoms lower their energy by π-bonding and this is the source of the low surface energy. The even lower surface energy of hydrogenated amorphous carbon (a-C:H) is due to the hydrogenation of all broken surface bonds. © 2005 Elsevier B.V. All rights reserved.

Assessing the road-damaging potential of heavy vehicles

The relative influence of various heavy vehicle design features on road-damaging potential is discussed. Testing procedures that could be used to measure the road-damaging potential of heavy vehicles are examined. A validated vehicle simulation is used to examine some of the characteristics of dynamic tyre forces generated by typical leaf sprung and air sprung articulated heavy vehicles for typical highway conditions. The proposed EC suspension test is simulated and the results compared with dynamic tyre forces generated under highway conditions. It is concluded that the road-damaging potential of a vehicle cannot be assessed by the simplistic parametric measurement of the proposed EC test. It is questionable whether a vehicle that passes the test will be any more 'road friendly' than one that fails.