Source author record

Graeme J. Ackland

Graeme 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

13works
8topics
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

13 published item(s)

preprint2022arXiv

Detectability of core level crossing and electronic topological transformations: the case of Osmium

Osmium, the least compressible metal, has recently been observed to undergo abrupt changes in the c/a ratio at extreme pressures. These are claimed to provide evidence for two unusual electronic behaviors: a crossing of the semicore 4f and 5p levels, and an electronic topological transition. We demonstrate that these two electronic phenomena are readily reproduced and understood in density functional theory, but that neither perturbs the trend in c/a ratio against pressure. Hence the observed anomalies in c/a must have another cause. Osmium is also notable for its high yield stress: the c/a anomalies lie well within the differential strains which osmium can support. We propose that observed c/a changes can arise from mechanical yield of crystallites with strong preferred orientation under high deviatoric stress in the experimental data. We discuss what evidence remains for the more general hypothesis that core-level overlap under pressure can have measurable effects on the crystal structure in any material.

preprint2021arXiv

Equation of state and strength of diamond in high pressure ramp loading

Diamond is used extensively as a component in high energy density experiments, but existing equation of state (EOS) models do not capture its observed response to dynamic loading. In particular, in contrast with first principles theoretical EOS models, no solid-solid phase changes have been detected, and no general-purpose EOS models match the measured ambient isotherm. We have performed density functional theory (DFT) calculations of the diamond phase to ~10TPa, well beyond its predicted range of thermodynamic stability, and used these results as the basis of a Mie-Greuneisen EOS. We also performed DFT calculations of the elastic moduli, and calibrated an algebraic elasticity model for use in simulations. We then estimated the flow stress of diamond by comparison with the stress-density relation measured experimentally in ramp-loading experiments. The resulting constitutive model allows us to place a constraint on the Taylor-Quinney factor (the fraction of plastic work converted to heat) from the observation that diamond does not melt on ramp compression.

preprint2020arXiv

Phase behaviour of the quantum Lennard-Jones solid

The Lennard-Jones potential is perhaps one of the most widely-used models for the interaction of uncharged particles, such as noble gas solids. The phase diagram of the classical LJ solid is known to exhibit transitions between hcp and fcc phases. However, the phase behaviour of the quantum Lennard-Jones solid remains unknown. Thermodynamic integration based on path integral molecular dynamics and lattice dynamics calculations are used to study the phase stability of the hcp and fcc Lennard-Jones solids. The hcp phase is shown to be stabilized by quantum effects in PIMD while fcc is shown to be favoured by lattice dynamics, which suggests a possible re-entrant low pressure hcp phase for highly quantum systems. Implications for the phase stability of noble gas solids are discussed. For parameters equating to Helium, the expansion due to zero-point vibrations is associated with quantum melting: neither crystal structure is stable at zero pressure.

preprint2020arXiv

The raman signal of a hindered rotor

We present a method for calculation of Raman modes of the quantum solid phase I solid hydrogen and deuterium. We use the mean-field assumption that the quantised excitations are localized on one molecule. This is done by explicit solution of the time-dependent Schroedinger equation in an angle-dependent potential, and direct calculation of the polarisation. We show that in the free-rotor limit, the H$_2$ and D$_2$ frequencies differ by a factor of 2, which evolves toward $\sqrt{2}$ as the modes acquire librational character due to stronger interactions. The ratio overshoots $\sqrt{2}$ if anharmonic terms weaken the harmonic potential. We also use density functional theory and molecular dynamics to calculate the E$_{2_g}$ optical phonon frequency and the Raman linewidths. The molecular dynamics shows that the molecules are not free rotors except at very low pressure and high temperature, and become like oscillators as phase II is approached. We fit the interaction strengths to experimental frequencies, but good agreement for intensities requires us to also include strong preferred-orientation and stimulated Raman effects between S$_0$(1) and S$_0$(0) contributions. The experimental Raman spectrum for phase II cannot be reproduced, showing that the mean-field assumption is invalid in that case.

preprint2020arXiv

Thermodynamic anomalies and three distinct liquid-liquid transitions in warm dense liquid hydrogen

The properties of hydrogen at high pressure have wide implications in astrophysics and high-pressure physics. Its phase change in the liquid is variously described as a metallization, H2-dissociation, density discontinuity or plasma phase transition. It has been tacitly assumed that these phenomena coincide at a first-order liquid-liquid transition (LLT). In this work, the relevant pressure-temperature conditions are thoroughly explored with first-principles molecular dynamics. We show there is a large dependency on exchange-correlation functional and significant finite size effects. We use hysteresis in a number of measurable quantities to demonstrate a first-order transition up to a critical point, above which molecular and atomic liquids are indistinguishable. At higher temperature beyond the critical point, H2-dissociation becomes a smooth cross-over in the supercritical region that can be modelled by a pseudo-transition, where the H2-2H transformation is localized and does not cause a density discontinuity at metallization. Thermodynamic anomalies and counter-intuitive transport behavior of protons are also discovered even far beyond the critical point, making this dissociative transition highly relevant to the interior dynamics of Jovian planets. Below the critical point, simulation also reveals a dynamic H2-2H chemical equilibrium with rapid interconversion, showing that H2 and H are miscible. The predicted critical temperature lies well below the ionization temperature. Our calculations unequivocally demonstrate that there are three distinct regimes in the liquid-liquid transition of warm dense hydrogen.

preprint2019arXiv

Quadrupole Arrangements and the Ground State of Solid Hydrogen

The electric quadrupole-quadrupole ($\mathcal{E}_{qq}$) interaction is believed to play an important role in the broken symmetry transition from Phase I to II in solid hydrogen. To evaluate this, we study structures adopted by purely classical quadrupoles using Markov Chain Monte Carlo simulations of fcc and hcp quadrupolar lattices. Both undergo first-order phase transitions from rotationally ordered to disordered structures, as indicated by a discontinuity in both quadrupole interaction energy ($\mathcal{E}_{qq}$) and its heat capacity. Cooling fcc reliably induced a transition to the P$a3$ structure, whereas cooling hcp gave inconsistent, frustrated and $c/a$-ratio-dependent broken symmetry states. Analysing the lowest-energy hcp states using simulated annealing, we found P$6_3/m$ and P$ca2_1$ structures found previously as minimum-energy structures in full electronic structure calculations. The candidate structures for hydrogen Phases III-V were not observed. This demonstrates that $\mathcal{E}_{qq}$ is the dominant interaction determining the symmetry breaking in Phase II. The disorder transition occurs at significantly lower temperature in hcp than fcc, showing that the $\mathcal{E}_{qq}$ cannot be responsible for hydrogen Phase II being based on hcp.

preprint2019arXiv

Structures of solid hydrogen at 300K

We present results predicting experimentally measurable structural quantities from molecular dynamics studies of hydrogen. In doing this, we propose a paradigm shift for experimentalists -- that the predictions from such calculations should be seen as the most likely hypotheses. Specifically, the experimental results should be aiming to distinguish between the candidate low-energy structures, rather than aiming to solve the simplest structure consistent with the data. We show that the room temperature X-ray diffraction patterns for hydrogen phases I, III, IV and V are very similar, with only small peaks denoting symmetry-breaking from the hcp Phase I. Because they incorporate atomic displacements the XRD patterns implied by molecular dynamics calculations are very different from those arising from the static minimum enthalpy structures found by structure searching. Simulations also show that within Phase I the molecular becomes increasingly confined to the basal plane and suggest the possibility of an unusual critical point terminating the Phase I-III boundary line.

preprint2016arXiv

Relating melting trends and elasticity in simple metals: an empirical potential approach

We demonstrate that the melting points and other thermodynamic quantities of the alkali metals can be calculated based on static crystalline properties. To do this we derive analytic interatomic potentials for the alkali metals fitted precisely to cohesive and vacancy energies, elastic moduli, lattice parameter and crystal stability. These potentials are then used to calculate melting points by simulating the equilibration of solid and liquid samples in thermal contact at ambient pressure. With the exception of lithium, remarkably good agreement is found with experimental values. The instability of the bcc structure in Li and Na at low temperatures is also reproduced, and, unusually, is not due to a soft T1N phonon mode. No forces or finite temperature properties are included in the fit, so this demonstrates a surprisingly high level of intrinsic transferrability in the simple potentials. Currently, there are few potentials available for the alkali metals, so in, addition to demonstrating trends in behaviour, we expect that the potentials will be of broad general use.

preprint2015arXiv

Charge density wave in hydrogen at high pressure

We present extensive molecular dynamics (MD) simulations investigating numerous candidate crystal structures for hydrogen in conditions around the present experimental frontier (400GPa). Spontaneous phase transitions in the simulations reveal a new structure candidate comprising twofold coordinated chains of hydrogen atoms. We explain the electronic structure of this phase in terms of a charge density wave and calculate its experimental signature. In detailed tests of the accuracy of our calculation, we find that k-point sampling is far more important in MD than in static calculations, because of the freedom it give the atoms to rearrange themselves optimally for the given sampling.

preprint2013arXiv

Simulating Raman Spectra using molecular dynamics, and identification of high-pressure phases III and IV in hydrogen

We present a technique for extracting Raman intensities from ab initio molecular dynamics (MD) simulations at high temperature. The method is applied to the highly anharmonic case of dense hydrogen up to 500 K for pressures ranging from 180 GPa to 300 GPa. On heating or pressurizing we find first-order phase transitions at the experimental conditions of the phase III - IV boundary. Direct comparison of Raman vibrons with experiment provides excellent discrimination between subtly different structures, found in MD. We find candidate structures whose Raman spectra are in good agreement with experiment. The new phase obtained in high temperature simulations adopts a dynamic, simple hexagonal structure with three layer types: freely rotating hydrogen molecules, static hexagonal trimers and rotating hexagonal trimers. We show that previously calculated structures for phase IV are inconsistent with experiment, and their appearance in simulation is due to finite size effects.

preprint2012arXiv

Binding of an oxide layer to a metal: the case of Ti(10-10)/TiO2(100)

We study the chemical nature of the bonding of an oxide layer to the parent metal. In order to disentangle chemical effects from strain/misfit, Ti(10$\bar{1}$0)/TiO$_{2}$(100) interface has been chosen. We use the density functional pseudopotential method which gives good agreement with experiment for known properties of bulk and surface Ti and TiO$_2$. Two geometries, a film-like model (with free surface in the structure) and a bulk-like model (with no free surface in the structure) are used to simulate the interface, in each case with different terminations of Ti and TiO$_2$. For the single-oxygen interfaces, the interface energy obtained using these two models agree with each other; however for the double-oxygen ones, the relative stability is quite different. The disturbance to the electronic structure is confined within a few atomic layers of the interface. The interfacial bonding is mainly ionic, and surprisingly there is more charge transfer from Ti to O in the interface than in the bulk. In consequence the Ti/TiO$_2$ interface has stronger binding than the bulk of either material. This helps to explain why the oxide forms a stable, protective layer on Ti and Ti alloys.

preprint2005arXiv

The effect of pack formation at the 2005 world orienteering championships

At the 2005 world championships there was considerable discussion and a formal protest in the long distance race arising from a perceived advantage obtained by some athletes running together. It is shown that a statistical model presented in previous work [1] is applicable to this event, giving predictions of the final times to with 2-3 minutes. Using the model, we show that pack formation was inevitable in this format. The statistical benefit gained at the elite level from running with other competitors appears to derive both from increase speed through the terrain, and the elimination of large navigational errors.

preprint2000arXiv

Theoretical and computational study of high pressure structures in barium

Recent high pressure work has suggested that elemental barium forms a high pressure self-hosting structure (Ba IV) involving two `types' of barium atom. Uniquely among reported elemental structures it cannot be described by a single crystalline lattice, instead involving two interpenetrating incommensurate lattices. In this letter we report pseudopotential calculations demonstrating the stability and the potentially disordered nature of the `guest' structure. Using band structures and nearly-free electron theory we relate the appearance of Ba IV to an instability in the close-packed structure, demonstrate that it has a zero energy vibrational mode, and speculate about the structure's stability in other divalent elements.