Phase-field-crystal study of grain boundary premelting and shearing in bcc iron
Abstract
We use the phase-field-crystal (PFC) method to investigate the equilibrium premelting and nonequilibrium shearing behaviors of [001] symmetric tilt grain boundaries (GBs) at high homologous temperature over the complete range of misorientation 0<θ<90∘ in classical models of bcc Fe. We characterize the dependence of the premelted layer width W as a function of temperature and misorientation. In addition, we compute the thermodynamic disjoining potential whose derivative with respect to W represents the structural force between crystal-melt interfaces due to the spatial overlap of density waves. The disjoining potential is also computed by molecular dynamics (MD) simulations, for quantitative comparison with PFC simulations, and coarse-grained amplitude equations (AE) derived from PFC that provide additional analytical insights. We find that, for GBs over an intermediate range of misorientation (θmin<θ<θmax), W diverges as the melting temperature is approached from below, corresponding to a purely repulsive disjoining potential, while for GBs outside this range (θ<θmin or θmax<θ<90∘) W remains finite at the melting point. In the latter case, W corresponds to a shallow attractive minimum of the disjoining potential. The misorientation range where W diverges predicted by PFC simulations is much smaller than the range predicted by MD simulations when the small dimensionless parameter ɛ of the PFC model is matched to liquid structure factor properties. However, it agrees well with MD simulations with a lower ɛ value chosen to match the ratio of bulk modulus and solid-liquid interfacial free energy, consistent with the amplitude-equation prediction that θmin and 90∘-θmax scale as ∼ɛ1/2. The incorporation of thermal fluctuations in PFC is found to have a negligible effect on this range. In response to a shear stress parallel to the GB plane, GBs in PFC simulations exhibit coupled motion normal to this plane or sliding. Furthermore, the coupling factor exhibits a discontinuous change as a function of θ that reflects a transition between two coupling modes. Sliding is only observed over a range of misorientation that is a strongly increasing function of temperature for T/TM≥0.8 and matches roughly the range where W diverges at the melting point. The coupling factor for the two coupling modes is in excellent quantitative agreement with previous theoretical predictions [Cahn, Mishin, and Suzuki, Acta Mater.ACMAFD1359-645410.1016/j.actamat.2006.08.004 54, 4953 (2006)].
- Publication:
-
Physical Review B
- Pub Date:
- January 2013
- DOI:
- 10.1103/PhysRevB.87.024110
- arXiv:
- arXiv:1211.0908
- Bibcode:
- 2013PhRvB..87b4110A
- Keywords:
-
- 61.72.Mm;
- 64.70.D-;
- 81.30.Fb;
- 05.40.Ca;
- Grain and twin boundaries;
- Solid-liquid transitions;
- Solidification;
- Noise;
- Condensed Matter - Materials Science;
- Condensed Matter - Soft Condensed Matter
- E-Print:
- doi:10.1103/PhysRevB.87.024110