Molecular-dynamics method for the simulation of bulk-solid interfaces at high temperatures

Method for computing short-range forces between solid-liquid interfaces driving grain boundary premelting

We present a molecular dynamics based method for accurately computing short-range structural forces resulting from the overlap of spatially diffuse solid-liquid interfaces at wetted grain boundaries close to the melting point. The method is based on monitoring the fluctuations of the liquid layer width at different temperatures to extract the excess interfacial free energy as a function of this width. The method is illustrated for a high-energy Sigma9 twist boundary in pure Ni. The short-range repulsion driving premelting is found to be dominant in comparison to long-range dispersion and entropic forces and consistent with previous experimental findings that nanometer-scale layer widths may be observed only very close to the melting point.

Molecular dynamics simulation of a glissile dislocation interface propagating a martensitic transformation

The method of Parrinello and Rahman is generalized to include slip in addition to deformation of the simulation cell. Equations of motion are derived, and a microscopic expression for traction is introduced. Lagrangian constraints are imposed so that the combination of deformation and slip conform to the invariant plane shear characteristic of martensites. Simulation of a model transformation demonstrates the nucleation and propagation of a glissile dislocation interface.