A combination of fundamental measure density functional theory and Monte~Carlo computer simulation is used to determine the orientation-resolved interfacial tension and stiffness for the equilibrium hard-sphere crystal-fluid interface. Microscopic density functional theory is in quantitative agreement with simulations and predicts a tension of 0.66kBT/$\sigma${}2 with a small anisotropy of about 0.025kBT and stiffnesses with, e.g., 0.53kBT/$\sigma${}2 for the (001) orientation and 1.03kBT/$\sigma${}2 for the (111) orientation. Here kBT is denoting the thermal energy and $\sigma${} the hard-sphere diameter. We compare our results with existing experimental findings.
Nucleation in undercooled Ni is investigated by a combination of differential scanning calorimetry (DSC) experiments and Monte~Carlo (MC) simulation. By systematically varying the sample size in the DSC experiments, nucleation rates J over a range of 8 orders of magnitude are obtained. Evidence is given that these rates correspond to homogeneous nucleation. Free energy barriers $\Delta${}G*, as extracted from the measured J, are in very good agreement with those from the MC simulation. The MC simulation indicates a nonspherical geometry of crystalline clusters, fluctuating between prolate and oblate shape at a given size. Nevertheless, the temperature dependence of $\Delta${}G* is well described by classical nucleation theory.
Crystal-liquid interfaces in nickel are investigated by molecular-dynamics computer simulations. Inhomogeneous systems of size Lx×Ly×Lz with Lz=5Ly are prepared where the crystal fcc phase at different orientations coexists with the liquid phase, separated by planar interfaces in the xy-plane. The lateral dimensions are varied, using two different geometries with Lx=Ly and with LyLx. In the framework of capillary wave theory (CWT), anisotropic interfacial stiffnesses and tensions are determined using different predictions of CWT with respect to the spectrum, finite-size broadening and different geometries. From a parameterization in terms of cubic harmonics up to 8th order, the anisotropic interfacial free energy is obtained.
We compare experimental results from a quasi-two-dimensional colloidal hard sphere fluid to a Monte Carlo simulation of hard disks with small particle displacements. The experimental short-time self-diffusion coefficient scaled by the diffusion coefficient at infinite dilution, , strongly depends on the area fraction, pointing to significant hydrodynamic interactions at short times in the experiment, which are absent in the simulation. In contrast, the area fraction dependence of the experimental long-time self-diffusion coefficient is in quantitative agreement with obtained from the simulation. This indicates that the reduction in the particle mobility at short times due to hydrodynamic interactions does not lead to a proportional reduction in the long-time self-diffusion coefficient. Furthermore, the quantitative agreement between experiment and simulation at long times indicates that hydrodynamic interactions effectively do not affect the dependence of on the area fraction. In light of this, we discuss the link between structure and long-time self-diffusion in terms of a configurational excess entropy and do not find a simple exponential relation between these quantities for all fluid area fractions.
3 November 2017 in The Journal of Physical Chemistry C
A combination of quasielastic neutron scattering (QNS) and molecular-dynamics (MD) simulation is used to modify an embedded atom method (EAM) potential for Ti with respect to the description of melt properties and crystallization from the melt. In the MD simulation, the EAM model is optimized such that agreement with accurate data of self-diffusion coefficients from QNS is achieved. As a result the density and the melting temperature are in good agreement with experiment.