Duane D. Johnson

A generalized solid-state nudged elastic band method

A generalized solid-state nudged elastic band (G-SSNEB) method is presented for determining reaction pathways of solid-solid transformations involving both atomic and unit-cell degrees of freedom. We combine atomic and cell degrees of freedom into a unified description of the crystal structure so that calculated reaction paths are insensitive to the choice of periodic cell. For the rock-salt to wurtzite transition in CdSe, we demonstrate that the method is robust for mechanisms dominated either by atomic motion or by unit-cell deformation; notably, the lowest-energy transition mechanism found by our G-SSNEB changes with cell size from a concerted transformation of the cell coordinates in small cells to a nucleation event in large cells. The method is efficient and can be applied to systems in which the force and stress tensor are calculated using density functional theory.

Recycling Krylov Subspaces for Sequences of Linear Systems

Many problems in science and engineering require the solution of a long sequence of slowly changing linear systems. We propose and analyze two methods that significantly reduce the total number of matrix-vector products required to solve all systems. We consider the general case where both the matrix and right-hand side change, and we make no assumptions regarding the change in the right-hand sides. Furthermore, we consider general nonsingular matrices, and we do not assume that all matrices are pairwise close or that the sequence of matrices converges to a particular matrix. Our methods work well under these general assumptions, and hence form a significant advancement with respect to related work in this area. We can reduce the cost of solving subsequent systems in the sequence by recycling selected subspaces generated for previous systems. We consider two approaches that allow for the continuous improvement of the recycled subspace at low cost. We consider both Hermitian and non-Hermitian problems, and we analyze our algorithms both theoretically and numerically to illustrate the effects of subspace recycling. We also demonstrate the effectiveness of our algorithms for a range of applications from computational mechanics, materials science, and computational physics.

Generalized planar fault energies and twinning in Cu–Al alloys

We report ab initio density functional theory calculations of generalized planar fault energies of fcc Cu–xAl (x=0, 5.0, and 8.3at.%) alloys. We investigate the effects of substitutional solute Al on the unstable intrinsic γus and twin γut stacking fault energies (SFEs). Our results reveal an increased tendency of Cu–Al to deform preferentially by twinning with increasing Al content, consistent with experiment. We attribute this mechanical behavior to appreciable lowering of the twinning barrier γut, along with the stable intrinsic and twin SFEs.

Recent advances in the study of structural materials compatibility with hydrogen.

Hydrogen is a ubiquitous element that enters materials from many different sources. It almost always has a deleterious effect on mechanical properties. In non-hydride-forming systems, research to date has identified hydrogen-enhanced localized plasticity and hydrogen-induced decohesion as two viable mechanisms for embrittlement. However, a fracture prediction methodology that associates macroscopic parameters with the degradation mechanisms at the microscale has not been established, as of yet. In this article, we report recent work on modeling and simulation of hydrogen-induced crack initiation and growth. Our goal is to develop methodologies to relate characteristics of the degradation mechanisms from microscopic observations and first-principles calculations with macroscopic indices of embrittlement. The approach we use involves finite element analysis of the coupled hydrogen transport problem with hydrogen-assisted elastoplastic deformation, thermodynamic theories of decohesion, and ab initio density functional theory calculations of the hydrogen effect on grain boundaries.

A Comprehensive Search for Stable Pt-Pd Nanoalloy Configurations and Their Use as Tunable Catalysts

Using density-functional theory, we predict stable alloy configurations (ground states) for a 1 nm Pt-Pd cuboctahedral nanoparticle across the entire composition range and demonstrate their use as tunable alloy catalysts via hydrogen-adsorption studies. Unlike previous works, we use simulated annealing with a cluster expansion Hamiltonian to perform a rapid and comprehensive search that encompasses both high and low-symmetry configurations. The ground states show Pt(core)-Pd(shell) type configurations across all compositions but with specific Pd patterns. For catalysis studies at room temperatures, the ground states are more realistic structural models than the commonly assumed random alloy configurations. Using the ground states, we reveal that the hydrogen adsorption energy increases (decreases) monotonically with at. % Pt for the {111} hollow ({100} bridge) adsorption site. Such trends are useful for designing tunable Pd-Pt nanocatalysts for the hydrogen evolution reaction.

Site-preference and valency for rare-earth sites in (R-Ce)2Fe14B magnets

Rare-earth (R) permanent magnets of R2Fe14B have technological importance due to their high energy products, and they have two R-sites (Wyckoff 4f and 4g, with four-fold multiplicity) that affect chemistry and valence. Designing magnetic behavior and stability via alloying is technologically relevant to reduce critical (expensive) R-content while retaining key properties; cerium, an abundant (cheap) R-element, offers this potential. We calculate magnetic properties and Ce site preference in (R1−xCex) 2Fe14B [R = La,Nd] using density functional theory (DFT) methods—including a DFT + U scheme to treat localized 4f-electrons. Fe moments compare well with neutron data—almost unaffected by Hubbard U, and weakly affected by spin-orbit coupling. In La2Fe14B, Ce alloys for 0≤x≤1 and prefers smaller R(4f) sites, as observed, a trend we find unaffected by valence. Whereas, in Nd2Fe14B, Ce is predicted to have limited alloying (x≤0.3) with a preference for larger R(4g) sites, resulting in weak partial ordering and s...

Dirac node arcs in PtSn4

A combination of detailed photoelectron spectroscopy measurements and numerical simulations reveal the presence of so-called Dirac node arcs in the electronic structure of PtSn4.

Reliable First-Principles Alloy Thermodynamics via Truncated Cluster Expansions

In alloys cluster expansions (CE) are increasingly used to combine first-principles electronic-structure calculations and Monte Carlo methods to predict thermodynamic properties. As a basis-set expansion in terms of lattice geometrical clusters and effective cluster interactions, the CE is exact if infinite, but is tractable only if truncated. Yet until now a truncation procedure was not well defined and did not guarantee a reliable truncated CE. We present an optimal truncation procedure for CE basis sets that provides reliable thermodynamics. We then exemplify its importance in Ni3V, where the CE has failed unpredictably, and now show agreement to a range of measured values, predict new low-energy structures, and explain the cause of previous failures.

Thermal stability of MnBi magnetic materials

MnBi has attracted much attention in recent years due to its potential as a rare-earth-free permanent magnet material. It is unique because its coercivity increases with increasing temperature, which makes it a good hard phase material for exchange coupling nanocomposite magnets. MnBi phase is difficult to obtain, partly because the reaction between Mn and Bi is peritectic, and partly because Mn reacts readily with oxygen. MnO formation is irreversible and harmful to magnet performance. In this paper, we report our efforts toward developing MnBi permanent magnets. To date, high purity MnBi (>90%) can be routinely produced in large quantities. The produced powder exhibits 74.6 emu g(-1) saturation magnetization at room temperature with 9 T applied field. After proper alignment, the maximum energy product (BH)max of the powder reached 11.9 MGOe, and that of the sintered bulk magnet reached 7.8 MGOe at room temperature. A comprehensive study of thermal stability shows that MnBi powder is stable up to 473 K in air.

Energy pathways and directionality in deformation twinning

We present ab initio density functional theory calculations of twinning energy pathways for two opposite twinning modes, (111)[112¯] and (111)[1¯1¯2], in fcc materials to examine the directional nature of twinning which cannot be explained by classical twin nucleation models or the “twinnability” criterion. By accounting for these energy pathways in a multiscale model, we quantitatively predict the critical twinning stress for the (111)[1¯1¯2] mode to be substantially higher compared to the favorable (111)[112¯] mode (whose predicted stresses are in agreement with experiment), thus, ruling out twinning in the (111)[1¯1¯2] mode.