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. Our analytic approximation to the superbasin is transferable with respect to grain size and elliptical shapes and can be applied in simulations with constant temperature as well as constant heating rate. The resulting thermal desorption spectra are in close agreement with direct kinetic Monte Carlo simulations, but the calculations are computationally much more efficient. Our approach is thus applicable to much larger system sizes and provides a first step towards an atomistic understanding of the influence of structural features on the position and shape of peaks in thermal desorption spectra.This article is part of the themed issue 'The challenges of hydrogen and metals'.

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. The parametrizations are optimized for the equilibrium bcc structure and tested for atomic environments far from equilibrium that had not been included in the fitting procedure. These tests include structural energy differences for competing crystal structures; tetragonal, trigonal, hexagonal and orthorhombic deformation paths; formation energies of point defects as well as phonon dispersion relations. Our tests show good agreement with available experimental and theoretical data. In practice, we obtain the energetic ordering of vacancy, [1 1 1], [1 1 0], and [1 0 0] self-interstitial atom in agreement with density functional theory calculations.

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. We use established parametrizations for the bcc refractory metals W and Mo and test structural energy differences; tetragonal, trigonal, hexagonal and orthorhombic deformation paths; formation energies of point defects as well as phonon dispersion relations. We find that the numerical and analytic BOPs generally are in very good agreement for the calculation of energies. Different from the numerical BOPs, the forces in the analytic BOPs correspond exactly to the negative gradients of the energy. This makes it possible to use the analytic BOPs in dynamical simulations and leads to improved predictions of defect energies and phonons as compared to the numerical BOPs.

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.

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.

Spontaneous L1(2) order at Ni(90)Al(10)(110) surfaces: an x-ray and first-principles-calculation study.

We have combined x-ray diffraction studies with first-principles calculations to study the interplay between segregation and ordering at the (110) surface of Ni(90)Al(10). We find a L1(2)-ordered monolayer at the surface. The observed ordering as well as recently reported Al segregation at the surface are explained in a consistent picture. A delicate competition between the tendency for Al segregation and ordering in the Ni-Al system induced by the symmetry break at the surface stabilizes a long-range ordered surface in the entire concentration range c(Ni)>0.75.