Direct calculation of lattice Green function with arbitrary interactions for general crystals

Efficient computation of lattice defect geometries such as point defects, dislocations, disconnections, grain boundaries, interfaces and free surfaces requires accurate coupling of displacements near the defect to the long-range elastic strain. Flexible boundary condition methods embedded a defect in infinite harmonic bulk through the lattice Green function. We demonstrate an efficient and accurate calculation of the lattice Green function from the force-constant matrix for general crystals with an arbitrary basis by extending a method for Bravais lattices. New terms appear due to the presence of optical modes and the possible loss of inversion symmetry. By separately treating poles and discontinuities in reciprocal space, numerical accuracy is controlled at all distances. We compute the lattice Green function for a two-dimensional model with broken symmetry to elucidate the role of different coupling terms. The algorithm is generally applicable in two and three dimensions, to crystals with arbitrary number of atoms in the unit cell, symmetry, and interactions.

Prediction of thermal cross-slip stress in magnesium alloys from direct first principles data

We develop a first-principles model of thermally-activated cross-slip in magnesium in the presence of a random solute distribution. Electronic structure methods provide data for the interaction of solutes with prismatic dislocation cores and basal dislocation cores. Direct calculations of interaction energies are possible for solutes---K, Na, and Sc---that lower the Mg prismatic stacking fault energy to improve formability. To connect to thermally activated cross-slip, we build a statistical model for the distribution of activation energies for double kink nucleation, barriers for kink migration, and roughness of the energy landscape to be overcome by an athermal stress. These distributions are calculated numerically for a range of concentrations, as well as alternate approximate analytic expressions for the dilute limit. The analytic distributions provide a simplified model for the maximum cross-slip softening for a solute as a function of temperature. The direct interaction calculations predict lowered forming temperatures for Mg-0.7at.%Sc, Mg-0.4at.%K, and Mg-0.6at.%Na of approximately 250C.

First-principles data for solid-solution strengthening of magnesium: From geometry and chemistry to properties

Solid-solution strengthening results from solutes impeding the glide of dislocations. Existing theories of strength rely on solute-dislocation interactions, but do not consider dislocation core structures, which need an accurate treatment of chemical bonding. Here, we focus on strengthening of Mg, the lightest of all structural metals and a promising replacement for heavier steel and aluminum alloys. Elasticity theory, which is commonly used to predict the requisite solute-dislocation interaction energetics, is replaced with quantum-mechanical first-principles calculations to construct a predictive mesoscale model for solute strengthening of Mg. Results for 29 different solutes are displayed in a strengthening design map as a function of solute misfits that quantify volumetric strain and slip effects. Our strengthening model is validated with available experimental data for several solutes, including Al and Zn, the two most common solutes in Mg. These new results highlight the ability of quantum-mechanical first-principles calculations to predict complex material properties such as strength.