Hydride-Induced Reconstruction of Pd Electrode Surfaces: A Combined Computational and Experimental Study.

Designing electrocatalysts with optimal activity and selectivity relies on a thorough understanding of the surface structure under reaction conditions. In this study, experimental and computational approaches are combined to elucidate reconstruction processes on low-index Pd surfaces during H-insertion following proton electroreduction. While electrochemical scanning tunneling microscopy clearly reveals pronounced surface roughening and morphological changes on Pd(111), Pd(110), and Pd(100) surfaces during cyclic voltammetry, a complementary analysis using inductively coupled plasma mass spectrometry excludes Pd dissolution as the primary cause of the observed restructuring.

Efficient Parametrization of Transferable Atomic Cluster Expansion for Water.

We present a highly accurate and transferable parametrization of water using the atomic cluster expansion (ACE). To efficiently sample liquid water, we propose a novel approach that involves sampling static calculations of various ice phases and utilizing the active learning (AL) feature of the ACE-based D-optimality algorithm to select relevant liquid water configurations, bypassing computationally intensive ab initio molecular dynamics simulations. Our results demonstrate that the ACE descriptors enable a potential initially fitted solely on ice structures, which is later upfitted with few configurations of liquid, identified with AL to provide an excellent description of liquid water.

Atomic Cluster Expansion for Quantum-Accurate Large-Scale Simulations of Carbon.

We present an atomic cluster expansion (ACE) for carbon that improves over available classical and machine learning potentials. The ACE is parametrized from an exhaustive set of important carbon structures over extended volume and energy ranges, computed using density functional theory (DFT). Rigorous validation reveals that ACE accurately predicts a broad range of properties of both crystalline and amorphous carbon phases while being several orders of magnitude more computationally efficient than available machine learning models.

Role of pre-ordered liquid in the selection mechanism of crystal polymorphs during nucleation.

We investigate the atomistic mechanism of homogeneous nucleation during solidification in molybdenum employing transition path sampling. The mechanism is characterized by the formation of a pre-structured region of high bond-orientational order in the supercooled liquid followed by the emergence of the crystalline bulk phase within the center of the growing solid cluster. This precursor plays a crucial role in the process as it provides a diffusive interface between the liquid and crystalline core, which lowers the interfacial free energy and facilitates the formation of the bulk phase.

Identification of a multi-dimensional reaction coordinate for crystal nucleation in NiAl.

Nucleation during solidification in multi-component alloys is a complex process that comprises competition between different crystalline phases as well as chemical composition and ordering. Here, we combine transition interface sampling with an extensive committor analysis to investigate the atomistic mechanisms during the initial stages of nucleation in NiAl. The formation and growth of crystalline clusters from the melt are strongly influenced by the interplay between three descriptors: the size, crystallinity, and chemical short-range order of the emerging nuclei.

Experimental and Theoretical Investigation on Phase Formation and Mechanical Properties in Cr-Co-Ni Alloys Processed Using a Novel Thin-Film Quenching Technique.

The Cr-Co-Ni system was studied by combining experimental and computational methods to investigate phase stability and mechanical properties. Thin-film materials libraries were prepared and quenched from high temperatures up to 700 °C using a novel quenching technique. It could be shown that a wide A1 solid solution region exists in the Cr-Co-Ni system.

Fast diffusion mechanism in LiPS a concerted process of interstitial Li ions.

The synthesis of Li superionic conductor LiPS may be accompanied by the formation of a detrimental LiPS phase due to a high mixing sensitivity of precursor materials. This phase exhibits a poor ionic conductivity whose origins are not fully understood. Recently Dietrich investigated the energetics of Li ion migration in LiPS with nudged elastic band (NEB) calculations.

Atomistically Informed Extended Gibbs Energy Description for Phase-Field Simulation of Tempering of Martensitic Steel.

In this study we propose a unified multi-scale chemo-mechanical description of the BCT (Body-Centered Tetragonal) to BCC (Body-Centered Cubic) order-disorder transition in martensitic steel by adding the mechanical degrees of freedom to the standard CALPHAD (CALculation of PHAse Diagrams) type Gibbs energy description. The model takes into account external strain, the effect of carbon composition on the lattice parameter and elastic moduli. The carbon composition effect on the lattice parameters and elastic constants is described by a sublattice model with properties obtained from DFT (Density Functional Theory) calculations; the temperature dependence of the elasticity parameters is estimated from available experimental data.

A kinetic Monte Carlo approach to diffusion-controlled thermal desorption spectroscopy.

Atomistic simulations of thermal desorption spectra for effusion from bulk materials to characterize binding or trapping sites are a challenging task as large system sizes as well as extended time scales are required. Here, we introduce an approach where we combine kinetic Monte Carlo with an analytic approximation of the superbasins within the framework of absorbing Markov chains. We apply our approach to the effusion of hydrogen from BCC iron, where the diffusion within bulk grains is coarse grained using absorbing Markov chains, which provide an exact solution of the dynamics within a superbasin.

Atomistic insight into the non-classical nucleation mechanism during solidification in Ni.

Nucleation is a key step during crystallization, but a complete understanding of the fundamental atomistic processes remains elusive. We investigate the mechanism of nucleation during solidification in nickel for various undercoolings using transition path sampling simulations. The temperature dependence of the free energy barriers and rate constants that we obtain is consistent with the predictions of classical nucleation theory and experiments.

Thermodynamic and kinetic solid-liquid interface properties from transition path sampling.

We perform transition path sampling simulations to determine two of the key quantities in solidification, the solid-liquid interface energy and velocity, in a Lennard-Jones system. Our approach is applicable to a wide range of temperature and pressure conditions, at the melting temperature and out-of-equilibrium. We show that small system sizes are sufficient for good values of interface energies and velocities.

Collective Atomic Displacements during Complex Phase Boundary Migration in Solid-Solid Phase Transformations.

The A15 to bcc phase transition is simulated at the atomic scale based on an interatomic potential for molybdenum. The migration of the phase boundary proceeds via long-range collective displacements of entire groups of atoms across the interface. To capture the kinetics of these complex atomic rearrangements over extended time scales we use the adaptive kinetic Monte Carlo approach.

Martensitic transformation between competing phases in Ti-Ta alloys: a solid-state nudged elastic band study.

A combined density functional theory and solid-state nudged elastic band study is presented to investigate the martensitic transformation between β → (α″, ω) phases in the Ti-Ta system. The minimum energy paths along the transformation are calculated and the transformation mechanisms as well as relative stabilities of the different phases are discussed for various compositions. The analysis of the transformation paths is complemented by calculations of phonon spectra to determine the dynamical stability of the β, α″, and ω phase.

Non-collinear magnetism with analytic Bond-Order Potentials.

The theory of analytic Bond-Order Potentials as applied to non-collinear magnetic structures of transition metals is extended to take into account explicit rotations of Hamiltonian and local moment matrix elements between locally and globally defined spin-coordinate systems. Expressions for the gradients of the energy with respect to the Hamiltonian matrix elements, the interatomic forces and the magnetic torques are derived. The method is applied to simulations of the rotation of magnetic moments in α iron, as well as α and β manganese, based on d-valent orthogonal tight-binding parametrizations of the electronic structure.

Vacancy mobility and interaction with transition metal solutes in Ni.

Interaction of Re, Ta, W and Mo solutes with vacancies and their diffusion in fcc Ni is investigated by density-functional theory in combination with kinetic Monte Carlo simulations. Interaction energies are calculated for the first six neighbor shells around the solutes and a complete set of diffusion barriers for these shells is provided. Further, diffusion coefficients for the four elements in Ni as well as for vacancies in the presence of these elements are calculated.

Analytic bond-order potentials for the bcc refractory metals Nb, Ta, Mo and W.

Bond-order potentials (BOPs) are based on the tight-binding approximation for determining the energy of a system of interacting atoms. The bond energy and forces are computed analytically within the formalism of the analytic BOPs. Here we present parametrizations of the analytic BOPs for the bcc refractory metals Nb, Ta, Mo and W.

Parameterized electronic description of carbon cohesion in iron grain boundaries.

We employ a recently developed iron-carbon orthogonal tight-binding model in calculations of carbon in iron grain boundaries. We use the model to evaluate the properties of carbon near and on the Σ5 (3 1 0)[0 0 1] symmetric tilt grain boundary (GB) in iron, and calculations show that a carbon atom lowers the GB energy by 0.29 eV/atom in accordance with DFT.

Modelling the lattice dynamics in Si(x)Ge(1-x) alloys.

The development of simplified models for the simulation of thermodynamic and thermal transport properties in random alloys is of great importance. In this paper we show how a simple second nearest neighbour model can reliably capture the lattice dynamics of Si(x)Ge(1-x) alloys. The model parameters are extracted from DFT-calculated force constant matrices for pure Si, pure Ge and the Si0.

Comparison of analytic and numerical bond-order potentials for W and Mo.

Bond-order potentials (BOPs) are derived from the tight-binding approximation and provide a linearly-scaling computation of the energy and forces for a system of interacting atoms. While the numerical BOPs involve the numerical integration of the response (Green's) function, the expressions for the energy and interatomic forces are analytical within the formalism of the analytic BOPs. In this paper we present a detailed comparison of numerical and analytic BOPs.

Environmental tight-binding modeling of nickel and cobalt clusters.

Tight-binding models derived from density functional theory potentially provide a systematic approach to the development of accurate and transferable models of multicomponent systems. We introduce a systematic methodology for environmental tight binding in which both the overlap and environmental contributions to the electronic structure are included. The parameters of the model are determined directly from ab initio considerations, thus providing a formal conceptual link to density functional approaches.

High throughput density functional investigations of the stability, electronic structure and thermoelectric properties of binary silicides.

The structural stabilities of binary Mg-X (X = Si, Ge, Sn) and 4d transition metal silicides Mo-Si and Ru-Si are investigated. The convex hulls of stable alloys are in overall good agreement with the known experimental phase diagrams. It is shown how the Si-rich Ru-Si structures have band gaps at the Fermi-level and how the Ru(2)Si(3) structure is stabilized compared to the corresponding Fe(2)Si(3) structure.

Tight-binding simulation of transition-metal alloys.

In order to perform atomistic simulations of steel, it is necessary to have a detailed understanding of the complex interatomic interactions in transition metals and their alloys. The tight-binding approximation provides a computationally efficient, yet accurate, method to investigate such interactions. In the present work, an orthogonal tight-binding model for Fe, Mn and Cr, with the explicit inclusion of magnetism, has been parameterized from ab initio density-functional calculations.

New superconducting and semiconducting Fe-B compounds predicted with an ab initio evolutionary search.

New candidate ground states at 1:4, 1:2, and 1:1 compositions are identified in the well-known Fe-B system via a combination of ab initio high-throughput and evolutionary searches. We show that the proposed oP12-FeB2 stabilizes by a break up of 2D boron layers into 1D chains while oP10-FeB4 stabilizes by a distortion of a 3D boron network. The uniqueness of these configurations gives rise to a set of remarkable properties: oP12-FeB2 is expected to be the first semiconducting metal diboride and oP10-FeB4 is shown to have the potential for phonon-mediated superconductivity with a T(c) of 15-20 K.

Cluster expansions in multicomponent systems: precise expansions from noisy databases.

We have performed a systematic analysis of the numerical errors contained in the databases used in cluster expansions of multicomponent alloys. Our results underscore the importance of numerical noise in the determination of the effective cluster interactions and in the expansion determination. The relevance of the size of and information contained in the input database is highlighted.

Competition between order and phase separation in Au-Ni.

We have measured and theoretically analyzed the diffuse scattering in the binary alloy system Au-Ni, which has been proposed as a testing ground for theories of alloy phase stability. We found strong evidence that in the alloys Au3Ni and Au3Ni2, fluctuations of both ordering- and clustering-type are competing with each other. Our results resolve a long-standing controversy on the balance of relaxation and mixing energies in this alloy system and explain recent findings of ordering in thin Au-Ni films.