Structural and chemical embrittlement of grain boundaries by impurities: A general theory and first-principles calculations for copper

First-principles calculations of the Sigma 5(310)[001] symmetric tilt grain boundary in Cu with Bi, Na, and Ag substitutional impurities provide evidence that in the phenomenon of Bi embrittlement of Cu grain boundaries electronic effects do not play a major role; on the contrary, the embrittlement is mostly a structural or "size" effect. Na is predicted to be nearly as good an embrittler as Bi, whereas Ag does not embrittle the boundary in agreement with experiment. While we reject the prevailing view that "electronic" effects (i.e., charge transfer) are responsible for embrittlement, we do not exclude the role of chemistry. However, numerical results show a striking equivalence between the alkali metal Na and the semimetal Bi, small differences being accounted for by their contrasting "size" and "softness" (defined here). In order to separate structural and chemical effects unambiguously if not uniquely, we model the embrittlement process by taking the system of grain boundary and free surfaces through a sequence of precisely defined gedanken processes; each of these representing a putative mechanism. We thereby identify three mechanisms of embrittlement by substitutional impurities, two of which survive in the case of embrittlement or cohesion enhancement by interstitials. Two of the three are purely structural and the third contains both structural and chemical elements that by their very nature cannot be further unraveled. We are able to take the systems we study through each of these stages by explicit computer simulations and assess the contribution of each to the net reduction in intergranular cohesion. The conclusion we reach is that embrittlement by both Bi and Na is almost exclusively structural in origin; that is, the embrittlement is a size effect.

On thin ice: surface order and disorder during pre-melting

The effect of temperature on the structure of the ice Ih (0001) surface is considered through a series of molecular dynamics simulations on an ice slab. At relatively low temperatures (200K) a small fraction of surface self-interstitials (i.e. admolecules) appear that are formed exclusively from molecules leaving the outermost bilayer. At higher temperatures (ca. 250 K), vacancies start to appear in the inner part of the outermost bilayer exposing the underlying bilayer and providing sites with a high concentration of the dangling hydrogen bonds. Around 250-260 K aggregates of molecules formed on top of the outermost bilayer from self-interstitials become more mobile and have diffusivities approaching that of liquid water. At similar to 270-280 K the inner bilayer of one surface noticeably destructures and it appears that at above 285 K both surfaces are melting. The observed disparity in the onset of melting between the two sides of the slab is rationalised by considering the relationship between surface energy and the spatial distribution of protons at the surface; thermodynamic stability is conferred on the surface by maximising separations between dangling protons at the crystal exterior. Local hotspots associated with a high dangling proton density are suggested to be susceptible to pre-melting and may be more efficient at trapping species at the external surface than regions with low concentrations of protons thus potentially helping ice particles to catalyse reactions. A preliminary conclusion of this work is that only about 10-20 K below the melting temperature of the particular water potential employed is major disruption of the crystalline lattice noted which could be interpreted as being "liquid", the thickness of this film being about a nanometre.