Dislocation-pairing transitions in hot grain boundaries.

We report the finding of a novel grain-boundary structural phase transition in both molecular-dynamics and phase-field-crystal simulations of classical models of bcc Fe. This transition is characterized by pairing of individual dislocations with mixed screw and edge components. We demonstrate that this type of transition is driven by a combination of factors including elastic softening, core interaction, and core disordering. At high homologous temperatures the occurrence of this transition is shown to prevent premelting at misorientation angles where it would otherwise be expected.

Atomistic simulations of nonequilibrium crystal-growth kinetics from alloy melts.

Nonequilibrium kinetic properties of alloy crystal-melt interfaces are calculated by molecular-dynamics simulations. The relationships between the interface velocity, thermodynamic driving force, and solute partition coefficient are computed and analyzed within the framework of kinetic theories accounting for solute trapping and solute drag. The results show a transition to complete solute trapping at high growth velocities, establish appreciable solute drag at low growth velocities, and provide insights into the nature of crystalline anisotropies and solute effects on interface mobilities.

Atomistic simulation study of the structure and dynamics of a faceted crystal-melt interface.

A detailed analysis of the structure and dynamics of the crystal-melt interface region in silicon, modeled with the Stillinger-Weber potential, is performed via molecular dynamics simulations. The focus is on the faceted (111) crystal-melt interface, but properties of the rough (100) interface are also determined. We find an intrinsic 10-90 interface width of 0.681+/-0.001 nm for the coarse-grained density profile at the (111) interface and a 0.570+/-0.005 nm width at the (100) interface. Coarse-grained profiles of a suitably defined local order parameter are found to show a smaller width anisotropy between (111) and (100) interfaces while the order profiles exhibit a 0.20-0.25 nm shift in position toward the crystal phase relative to the corresponding density profiles. The structural analysis of the layer of melt adjacent to the (111) facet of the crystal finds ordered clusters with average lifetimes of 16 ps , as determined from autocorrelations of time-dependent layer structure factors, and cluster radii of gyration from 0.2 nm for the smallest cells to as large as 1.5 nm .

Structural disjoining potential for grain-boundary premelting and grain coalescence from molecular-dynamics simulations.

We describe a molecular-dynamics framework for the direct calculation of the short-ranged structural forces underlying grain-boundary premelting and grain coalescence in solidification. The method is applied in a comparative study of (i) a Sigma9115120 degrees twist and (ii) a Sigma9110{411} symmetric tilt boundary in a classical embedded-atom model of elemental Ni. Although both boundaries feature highly disordered structures near the melting point, the nature of the temperature dependence of the width of the disordered regions in these boundaries is qualitatively different. The former boundary displays behavior consistent with a logarithmically diverging premelted layer thickness as the melting temperature is approached from below, while the latter displays behavior featuring a finite grain-boundary width at the melting point. It is demonstrated that both types of behavior can be quantitatively described within a sharp-interface thermodynamic formalism involving a width-dependent interfacial free energy, referred to as the disjoining potential. The disjoining potential for boundary (i) is calculated to display a monotonic exponential dependence on width, while that of boundary (ii) features a weak attractive minimum. The results of this work are discussed in relation to recent simulation and theoretical studies of the thermodynamic forces underlying grain-boundary premelting.

Kinetic coefficient of steps at the Si(111) crystal-melt interface from molecular dynamics simulations.

Nonequilibrium molecular dynamics simulations are applied to the investigation of step-flow kinetics at crystal-melt interfaces of silicon, modeled with the Stillinger-Weber potential [Phys. Rev. B 31, 5262 (1985)]. Step kinetic coefficients are calculated from crystallization rates of interfaces that are vicinals of the faceted (111) orientation. These vicinal interfaces contain periodic arrays of bilayer steps, and they are observed to crystallize in a step-flow growth mode at undercoolings lower than 40 K. Kinetic coefficients for both [110] and [121] oriented steps are determined for several values of the average step separation, in the range of 7.7-62.4 A. The values of the step kinetic coefficients are shown to be highly isotropic, and are found to increase with increasing step separation until they saturate at step separations larger than approximately 50 A. The largest step kinetic coefficients are found to be in the range of 0.7-0.8 m(sK), values that are more than five times larger than the kinetic coefficient for the rough (100) crystal-melt interface in the same system. The dependence of step mobility on step separation and the relatively large value of the step kinetic coefficient are discussed in terms of available theoretical models for crystal growth kinetics from the melt.

Mixtures of lattice polymers with structured monomers.

The influence of monomer structure on the thermodynamic properties of lattice model polymer blends is investigated through Monte Carlo computations. The model of lattice polymers with monomer structure has been used extensively in the context of the lattice cluster theory (LCT), a thermodynamic theory for polymer mixtures in the liquid state. The Monte Carlo computations provide the first unequivocal test of the accuracy of the LCT predictions for binary mixtures of polymers with structured monomers. Four types of monomer structures are analyzed, corresponding to to the monomers of polyethylene, polypropylene, polyethylethylene, and polyisobutylene (PIB). Most computations use chains with M=12 and 24 beads and the total volume fraction of the beads is phi=0.6. Both structurally symmetric and asymmetric blends are investigated. For the symmetric case, the predictions of the LCT for the energies of mixing and the liquid-liquid coexistence curves are in qualitative agreement with the Monte Carlo computations, except for the PIB/PIB symmetric blend. For structurally asymmetric blends, the LCT does not capture contributions to the energy of mixing arising solely from structural differences between the components. Computational estimates of the nonideal entropy of mixing indicate that the LCT also underestimates the entropic cost of mixing chains with different structures, thus explaining some discrepancies between the theoretical and the Monte Carlo liquid--liquid coexistence curves.