Source author record

G. J. Ackland

G. J. Ackland appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

15works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

15 published item(s)

preprint2016arXiv

Rapid Equilibration by algorithmic quenching the ringing mode in molecular dynamics

Long wavelength acoustic phonons are normally weakly coupled to other vibrational modes in a crystalline system. This is particularly problematic in molecular dynamics calculations where vibrations at the system size scale are typically excited at initiation. The equilibration time for these vibrations depends on the strength of coupling to other modes, so is typically very long. A very simple deterministic method which avoids the use of is presented which removes this problem.

preprint2014arXiv

Efficacious calculation of Raman spectra in high pressure hydrogen

We present and evaluate an efficient method for simulating Raman spectra from molecular dynamics (MD) calculations {\it without} defining normal modes. We apply the method to high pressure hydrogen in the high-temperature "Phase IV": a plastic crystal in which the conventional picture of fixed phonon eigenmodes breaks down. Projecting trajectories onto in-phase molecular stretches is shown to be many orders of magnitude faster than polarisability calculations, allowing statistical averaging at high-temperature. The simulations are extended into metastable regimes and identify several regimes associated with symmetry-breaking on different timescales, which are shown to exhibit features in the Raman spectra at the current experimental limit of resolvability. In this paper we have concentrated on the methodology, a fuller description of the structure of Phase IV hydrogen is given in a previous paper

preprint2014arXiv

Lattice-switch Monte Carlo: the fcc-bcc problem

Lattice-switch Monte Carlo is an efficient method for calculating the free energy difference between two solid phases, or a solid and a fluid phase. Here, we provide a brief introduction to the method, and list its applications since its inception. We then describe a lattice switch for the fcc and bcc phases based on the Bain orientation relationship. Finally, we present preliminary results regarding our application of the method to the fcc and bcc phases in the Lennard-Jones system. Our initial calculations reveal that the bcc phase is unstable, quickly degenerating into some as yet undetermined metastable solid phase. This renders conventional lattice-switch Monte Carlo intractable for this phase. Possible solutions to this problem are discussed.

preprint2014arXiv

Probing atomic environments in alloys by electron spectroscopy

In alloys exhibiting substitutional disorder, the variety of atomic environments manifests itself as a `disorder broadening' in their core level binding energy spectra. Disorder broadening can be measured experimentally, and in principle can be used to deduce information about specific atomic environments within a sample. However, progress in this endeavor is hampered by the lack of a model for this phenomenon which can treat complex systems. In this work we describe such a model. The model is used to elucidate the relationship between charge transfer, atomic environment, and disorder broadening in complex systems, with a focus on the problem of characterizing the interface quality of CuNi multilayers. We also validate the model against the results of ab initio electronic structure calculations. Several counterintuitive aspects of the disorder broadening phenomenon are uncovered, an understanding of which is essential for the correct interpretation of experimental results. For instance, it is shown that systems with inhomogeneous concentration profiles can exhibit disorder broadenings significantly larger than random alloys. Furthermore in some systems a `disorder narrowing' is even possible.

preprint2014arXiv

Transition metal solute interactions with point defects in austenitic iron from first principles

We present a comprehensive set of first principles electronic structure calculations to study transition metal solutes and their interactions with point defects in austenite. Clear trends were observed across the series. Solute-defect interactions were strongly correlated to the solute size factors, consistent with local strain field effects. Strong correlations with results in ferrite show insensitivity to the underlying crystal structure in Fe. Oversized solutes act as strong traps for vacancy and self-interstitial defects and as nucleation sites for the development of proto-voids and small self-interstitial loops. The reduction in defect mobility and net defect concentrations explains the observed radiation-damage resistance in austenitic steels doped with oversized solutes. Oversized solutes remaining dissolved in oxide dispersion-strengthened (ODS) steels could contribute to their radiation-damage resistance. Ni and Co diffuse more slowly than Fe, along with any vacancy flux produced under irradiation below a critical temperature, which is 400 K for Co and their concentrations should be enhanced at defect sinks. Cr and Cu diffuse more quickly than Fe, against a vacancy flux and will be depleted at defect sinks. Oversized solutes early in the transition metal series form highly-stable solute-centred divacancy (SCD) defects with a nearest-neighbour vacancy. The vacancy-mediated diffusion of these solutes is dominated by the dissociation and reassociation of the SCDs, with a lower activation energy than for self-diffusion, which has important implications for the nucleation and growth of complex oxide nanoparticles containing these solutes in ODS steels. Interstitial-mediated solute diffusion is energetically disfavoured for all except Cr, Mn, Co and Ni. The central role that solute size plays in the results presented here means they should apply to other solvent metals and alloys.

preprint2013arXiv

Mechanism for Radiation Damage Resistance in Yttrium Oxide Dispersion Strengthened Steels

ODS steels based on yttrium oxide have been suggested as potential fusion reactor wall materials due to their observed radiation resistance properties. Presumably this radiation resistance can be related to the interaction of the particle with vacancies,self-interstitial atoms (SIAs) and other radiation damage debris. Density functional theory has been used to investigate this at the atomic scale. Four distinct interfaces, some based on HRTEM observations, between iron and yttrium oxide were investigated. It is been shown that the Y$_2$O$_3$-Fe interface acts as a strong trap with long-range attraction for both interstitial and vacancy defects, allowing recombination without altering the interface structure. The catalytic elimination of defects without change to the microstructure explains the improved behaviour of ODS steels with respect to radiation creep and swelling.

preprint2013arXiv

Molecular Dynamics and OKMC Study of Radiation Induced Motion of Voids and He Bubbles in BCC Iron

We show that voids adjacent to radiation damage cascades can be moved in their entirety by several lattice spacings. This is done using molecular dynamics cascade simulations in iron at energies of 1-5 keV. The effect of this mechanism is studied further using an OKMC code and shows enhancement of void diffusion by 2 orders of magnitude from 1x10^-22 cm^2/s to 3x10^-20 cm^2/s. Repeating the study on He bubbles shows that the movement is damped by the presence of helium in the void.

preprint2012arXiv

Temperature dependence in interatomic potentials and an improved potential for Ti

The process of deriving an interatomic potentials represents an attempt to integrate out the electronic degrees of freedom from the full quantum description of a condensed matter system. In practice it is the derivatives of the interatomic potentials which are used in molecular dynamics, as a model for the forces on a system. These forces should be the derivative of the free energy of the electronic system, which includes contributions from the entropy of the electronic states. This free energy is weakly temperature dependent, and although this can be safely neglected in many cases there are some systems where the electronic entropy plays a significant role. Here a method is proposed to incorporate electronic entropy in the Sommerfeld approximation into empirical potentials. The method is applied as a correction to an existing potential for titanium. Thermal properties of the new model are calculated, and a simple method for fixing the melting point and solid-solid phase transition temperature for existing models fitted to zero temperature data is presented.

preprint2012arXiv

The Structure of Eu-III

Previous x-ray diffraction studies have reported Eu to transform from the hcp structure to a new phase, Eu-III, at 18 GPa. Using x-ray powder diffraction we have determined that Eu remains hcp up to 33 GPa, and that the extra peaks that appear at 18 GPa are from an impurity phase with space group R-3c . Above 33 GPa the diffraction pattern becomes very much more complex, signalling a transition to a phase with a distorted hcp structure.

preprint2011arXiv

First Principles Calculations of Defects in Unstable Crystals: Austenitic Iron

First principles calculations have given a new insight into the energies of point defects in many different materials, information which cannot be readily obtained from experiment. Most such calculation are done at zero Kelvin, with the assumption that finite temperature effects on defect energies and barriers are small. In some materials, however, the stable crystal structre of interest is mechanically unstable at 0K. In such cases, alternate approaches are needed. Here we present results of first principles calculations of austenitic iron using the VASP code. We determine an appropriate reference state for collinear magnetism to be the antiferromagnetic double-layer (AFM-d) which is both stable and lower in energy than other possible models for the low temperature limit of paramagnetic fcc iron. We then consider the energetics of dissolving typical alloying impurities (Ni, Cr) in the materials, and their interaction with point defects typical of the irradiated environment. We show that using standard methods there is a very strong dependence of calculated defect formation energies on the reference state chosen. Furthermore, there is a correlation between local free volume magnetism and energetics. The effect of substitutional Ni and Cr on defect properties is weak, rarely more than tenths of eV, so it is unlikely that small amounts of Ni and Cr will have a significant effect on the radiation damage in austenitic iron at high temperatures.

preprint2010arXiv

Rescaled potential for transition metal solutes in a-iron

We present empirical potentials for dilute transition metal solutes in a-iron. It is in the Finnis-Sinclair form and is therefore suitable for billion atom molecular dynamics simulations. First principles calculation shows that there are clear trends across the transition metal series which enable us to relate the rescaling parameters to principal quantum number and number of d electrons. The potential has been developed using a rescaling technique to provide solute-iron and solute-solute interactions from an existing iron potential.

preprint2004arXiv

Development of an interatomic potential for phosphorus impurities in iron

We present the derivation of an interatomic potential for the iron phosphorus system based primarily on {\it ab initio} data. Transferrability in this system is extremely problematic, and the potential is intended specifically to address the problem of radiation damage and point defects in iron containing low concentrations of phosphorus atoms. Some preliminary molecular dynamics calculations show that P strongly affects point defect migration.

preprint1995arXiv

Ab initio studies of structural instabilities in magnesium silicate perovskite

Density-functional simulations are used to calculate structural properties and high-symmetry phonons of the hypothetical cubic phase, the stable orthorhombic phase and an intermediate tetragonal phase of magnesium silicate perovskite. We show that the structure of the stable phase is well described by freezing in a small number of unstable phonons into the cubic phase. We use the frequencies of these unstable modes to estimate transition temperatures for cubic--tetragonal and tetragonal--orthorhombic phase transitions. These are investigated further to find that the coupling with the strain suggests that phonons give a better representation than rigid unit modes. The phonons of an intermediate tetragonal phase were found to be stable except for two rotational modes. The eigenvectors of the most unstable mode of each of the cubic and tetragonal phases account for all the positional parameters of the orthorhombic phase. The phase boundary for the orthorhombic--tetragonal transition intersects possible mantle geotherms, suggesting that the tetragonal phase may be present in the lower mantle.