Researcher profile

Chris J. Pickard

Chris J. Pickard contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
9works
0followers
8topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

9 published item(s)

preprint2022arXiv

Ephemeral data derived potentials for random structure search

Structure prediction has become a key task of the modern atomistic sciences, and depends on the rapid and reliable computation of the energy landscape. First principles density functional based calculations are highly reliable, faithfully describing the entire energy landscape. They are, however, computationally intensive and slow compared to interatomic potentials. Great progress has been made in the development of machine learning, or data derived, potentials, which promise to describe the entire energy landscape at first principles quality. However, compared to first principles approaches, their preparation can be time consuming and delay searching. Ab initio random structure searching (AIRSS) is a straightforward and powerful approach to structure prediction, based on the stochastic generation of sensible initial structures, and their repeated local optimisation. Here, a scheme, compatible with AIRSS, for the rapid construction of disposable, or ephemeral, data derived potentials (EDDPs) is described. These potentials are constructed using a homogeneous, separable manybody environment vector, and iterative neural network fits, sparsely combined through non-negative least squares. The approach is first tested on methane, boron nitride, elemental boron and urea. In the case of boron, an EDDP generated using data from small unit cells is used to rediscover the complex $γ$-boron structure without recourse to symmetry or fragments. Finally, an EDDP generated for silane (SiH$_4$) at 500 GPa enables the discovery of an extremely complex, dense, structure which significantly modifies silane's high pressure phase diagram. This has implications for the theoretical exploration for high temperature superconductivity in the dense hydrides, which have so far largely depended on searches in smaller unit cells.

preprint2022arXiv

Unified chemical theory of structure and bonding in elemental metals

Most elemental metals under ambient conditions adopt simple structures such as BCC, FCC and HCP in specific groupings across the Periodic Table, and on compression, many of these elements undergo transitions to surprisingly complex structures, including open and low-symmetry phases not expected from conventional free-electron based theories of metals. First-principles calculations have been able to reproduce many observed structures and transitions, but a unified, predictive theory of bonding that underlies this behavior is not yet in hand. We propose a remarkably simple theory based on large-scale high-throughput calculations of the elements over a broad range of thermodynamic conditions. The results broaden the conventional concept of metallic bonding with a new perspective that both explains the stability of different simple structures as well as lower symmetry phases arising from electron localization at interstitial sites in these structures. The success of this simple framework reveals the important role of chemical interactions in governing the structures and electronic properties of simple metals.

preprint2021arXiv

Design Principles for High Temperature Superconductors with Hydrogen-based Alloy Backbone at Moderate Pressure

Hydrogen-based superconductors provide a route to the long-sought goal of room-temperature superconductivity, but the high pressures required to metallize these materials limit their immediate application. For example, carbonaceous sulfur hydride, the first room-temperature superconductor, can reach a critical temperature (Tc) of 288 K only at the extreme pressure of 267 GPa. The next recognized challenge is the realization of room-temperature superconductivity at significantly lower pressures. Here, we propose a strategy for the rational design of high-temperature superconductors at low pressures by alloying small-radius elements and hydrogen to form ternary hydride superconductors with alloy backbones. We identify a hitherto unknown fluorite-type backbone in compositions of the form AXH8, which exhibit high temperature superconductivity at moderate pressures. The Fm-3m phase of LaBeH8, with a fluorite-type H-Be alloy backbone, is predicted to be metastable and superconducting with a Tc ~ 191 K at 50 GPa; a substantially lower pressure than that required by the geometrically similar clathrate hydride LaH10 (170 GPa). Our approach paves the way for finding high-Tc ternary hydride superconductors at conditions close to ambient pressures.

preprint2021arXiv

High-pressure phase behaviors of titanium dioxide revealed by a $Δ$-learning potential

Titanium dioxide has been extensively studied in the rutile or anatase phases, while its high-pressure phases are less well understood, despite that many are thought to have interesting optical, mechanical and electrochemical properties. First-principles methods such as density functional theory (DFT) are often used to compute the enthalpies of TiO$_2$ phases at 0~K, but they are expensive and thus impractical for long time-scale and large system-size simulations at finite temperatures. On the other hand, cheap empirical potentials fail to capture the relative stablities of the various polymorphs. To model the thermodynamic behaviors of ambient and high-pressure phases of TiO$_2$, we design an empirical model as a baseline, and then train a machine learning potential based on the difference between the DFT data and the empirical model. This so-called $Δ$-learning potential contains long-range electrostatic interactions, and predicts the 0~K enthalpies of stable TiO$_2$ phases that are in good agreement with DFT. We construct a pressure-temperature phase diagram of TiO$_2$ in the range $0<P<70$~GPa and $100<T<1500$~K. We then simulate dynamic phase transition processes, by compressing anatase at different temperatures. At 300~K, we observe predominantly anatase-to-baddeleyite transformation at about 20~GPa, via a martensitic two-step mechanism with highly ordered and collective atomic motion. At 2000~K, anatase can transform into cotunnite around 45-55~GPa in a thermally-activated and probabilistic manner, accompanied by diffusive movement of oxygen atoms. The pressures computed for these transitions show good agreement with experiments. Our results shed light on how to synthesize and stabilize high-pressure TiO$_2$ phases, and our method is generally applicable to other functional materials with multiple polymorphs.

preprint2021arXiv

Partially Diffusive Helium-Silica Compound in the Deep Interiors of Giant Planets

Helium is the second most abundant element in the universe, and together with silica, they are major components of giant planets. Exploring the reactivity and state of helium and silica under high pressure is of fundamental importance for developing and understanding of the evolution and internal structure of giant planets. Here, using first-principles calculations and crystal structure predictions, we identify four stable phases of a helium-silica compound with seven/eight-coordinated silicon atoms at pressure range of 600-4000 GPa, corresponding to the interior condition of the outer planets in the solar system. The density of HeSiO2 agrees with current structure models of the planets. This helium-silica compound exhibits a superionic-like helium diffusive state at the high pressure and high temperature conditions along the isentropes of Saturn, a metallic fluid state in Jupiter, and a solid state in the deep interiors of Uranus and Neptune. The reaction of helium and silica may lead to the erosion of the rocky core of giant planets and form a diluted core region. These results highlight the reactivity of helium under high pressure to form new compounds, and also provides evidence to help build more sophisticated interior models of giant planets.

preprint2020arXiv

Efficient prediction of Nucleus Independent Chemical Shifts for polycyclic aromatic hydrocarbons

Nuclear Magnetic Resonance (NMR) is one of the most powerful experimental techniques to characterize the structure of molecules and confined liquids. Nevertheless, the complexity of the systems under investigation usually requires complementary computational studies to interpret the NMR results. In this work we focus on polycyclic aromatic hydrocarbons (PAHs), an important class of organic molecules which have been commonly used as simple analogues for the spectroscopic properties of more complex systems, such as porous disordered carbons. We use Density Functional Theory (DFT) to calculate 13C chemical shifts and Nucleus Independent Chemical Shifts (NICS) for 34 PAHs. The results show a clear molecular size dependence of the two quantities, as well as the convergence of the 13C NMR shifts towards the values observed for graphene. We then present two computationally cheap models for the prediction of NICS in simple PAHs. We show that while a simple dipolar model fails to produce accurate values, a perturbative tight-binding approach can be successfully applied for the prediction of NICS in this series of molecules, including some non-planar ones containing 5- and 7-membered rings. This model, one to two orders of magnitudes faster than DFT calculations, is very promising and can be further refined in order to study more complex systems.

preprint2020arXiv

Predicting the phase diagram of titanium dioxide with random search and pattern recognition

Predicting phase stabilities of crystal polymorphs is central to computational materials science and chemistry. Such predictions are challenging because they first require searching for potential energy minima and then performing arduous free-energy calculations to account for entropic effects at finite temperatures. Here, we develop a framework that facilitates such predictions by exploiting all the information obtained from random searches of crystal structures. This framework combines automated clustering, classification and visualisation of crystal structures with machine-learning estimation of their enthalpy and entropy. We demonstrate the framework on the technologically important system of TiO2, which has many polymorphs, without relying on prior knowledge of known phases. We find a number of new phases and predict the phase diagram and metastabilities of crystal polymorphs at 1600 K, benchmarking the results against full free-energy calculations.

preprint2020arXiv

Stability and superconductivity of lanthanum and yttrium decahydrides

Rare-earth hydrides can exhibit high-temperature superconductivity under high pressure. Here, we apply a crystal structure prediction method to the current record-holding $T_c$ material, LaH$_{10}$, and a candidate for even higher $T_c$, YH$_{10}$. We find a pressure-induced phase transition from the experimentally observed cubic LaH$_{10}$ phase to a new hexagonal phase at around 420 GPa. This hexagonal phase could explain experimental observations of hcp impurities in fcc samples. We find that YH$_{10}$ forms similar structures to LaH$_{10}$ and discuss the sensitivity of superconductivity calculations to the computational parameters used.

preprint2020arXiv

Structure-specific, mode-resolved phonon coherence and specularity at graphene grain boundaries

In spite of their importance for understanding phonon transport phenomena in thin films and polycrystalline solids, the effects of boundary roughness scattering on phonon specularity and coherence are poorly understood because there is no general method for predicting their dependence on phonon momentum, frequency, branch and boundary morphology. Using the recently formulated atomistic S-matrix method, we develop a theory of boundary roughness scattering to determine the mode-resolved phonon coherence and specularity parameters from the scattering amplitudes. To illustrate the theory, we apply it to phonon scattering in realistic nonsymmetric graphene grain boundary (GB) models derived from atomic structure predictions. The method is validated by comparing its predictions with frequency-resolved results from lattice dynamics-based calculations. We prove that incoherent scattering is almost perfectly diffusive. We show that phonon scattering at the graphene GB is not diffuse although coherence and specularity are significantly reduced for long-wavelength flexural acoustic phonons. Our approach can be generalized to other atomistic boundary models.