Self-diffusivity and interdiffusivity of molten aluminum-copper alloys under pressure, derived from molecular dynamics.

We use molecular dynamics (MD) to simulate diffusion in molten aluminum-copper (AlCu) alloys. The self-diffusivities and Maxwell-Stefan diffusivities are calculated for AlCu mixtures using the Green-Kubo formulas at temperatures from 1000 to 4000 K and pressures from 0 to 25 GPa, along with additional points at higher temperatures and pressures. The diffusivities are corrected for finite-size effects. The Maxwell-Stefan diffusivity is compared to the diffusivity calculated from the self-diffusivities using a generalization of the Darken equation. We find that the effects of cross-correlation are small. Using the calculated self-diffusivities, we have assessed whether dilute hard-sphere and dilute Lennard-Jones models apply to the molten mixture. Neither of the two dilute gas diffusivities describes the diffusivity in molten Al and Cu. We report generalized analytic models for the self-diffusivities and interdiffusivity (mutual diffusivity) that fit the MD results well. The MD-derived transport coefficients are in good agreement with the available experimental data. We also report MD calculations of the viscosity and an analytic fit to those results. The ionic thermal conductivity is discussed briefly.

Energetics and kinetics of vacancy diffusion and aggregation in shocked aluminium via orbital-free density functional theory.

A possible mechanism for shock-induced failure in aluminium involves atomic vacancies diffusing through the crystal lattice and agglomerating to form voids, which continue to grow, ultimately resulting in ductile fracture. We employ orbital-free density functional theory, a linear-scaling first-principles quantum mechanics method, to study vacancy formation, diffusion, and aggregation in aluminium under shock loading conditions of compression and tension. We calculate vacancy formation and migration energies, and find that while nearest-neighbor vacancy pairs are unstable, next-nearest-neighbor vacancy pairs are stable. As the number of nearby vacancies increases, we predict that vacancy clusters preferentially grow through next-nearest-neighbor vacancies. The energetics are found to be greatly affected by expansion and compression, leading to insight as to how vacancies behave under shock conditions.

Finding transition states for crystalline solid-solid phase transformations.

We present a method to identify transition states and minimum energy paths for martensitic solid-solid phase transformations, thereby allowing quantification of the activation energies of such transformations. Our approach is a generalization of a previous method for identifying transition states for chemical reactions, namely the climbing image-nudged elastic band algorithm, where here the global deformation of the crystalline lattice (volume and shape fluctuations) becomes the reaction coordinate instead of atomic motion. We also introduce an analogue to the Born-Oppenheimer approximation that allows a decoupling of nuclear motion and lattice deformation, where the nuclear positions along the path are determined variationally according to current deformation state. We then apply this technique to characterize the energetics of elemental lithium phase transformations as a function of applied pressure, where we see a validation of the Born-Oppenheimer-like approximation, small energy barriers (expected for martensitic transformations), and a pronounced pressure dependence of various properties characterizing the phase transitions.

Importance of shear in the bcc-to-hcp transformation in iron.

Iron shows a pressure-induced martensitic phase transformation from the ground state ferromagnetic bcc phase to a nonmagnetic hcp phase at approximately 13 GPa. The exact transformation pressure (TP) and pathway are not known. Here we present a multiscale model containing a quantum-mechanics-based multiwell energy function accounting for the bcc and hcp phases of Fe and a construction of kinematically compatible and equilibrated mixed phases. This model suggests that shear stresses have a significant influence on the bcc<-->hcp transformation. In particular, the presence of modest shear accounts for the scatter in measured TPs. The formation of mixed phases also provides an explanation for the observed hysteresis in TP.