A primer on selecting grain boundary sets for comparison of interfacial fracture properties in molecular dynamics simulations

All grain boundaries are not equal in their predisposition for fracture due to the complex coupling between lattice geometry, interfacial structure, and mechanical properties. The ability to understand these relationships is crucial to engineer materials resilient to grain boundary fracture. Here, a methodology is presented to isolate the role of grain boundary structure on interfacial fracture properties, such as the tensile strength and work of separation, using atomistic simulations. Instead of constructing sets of grain boundary models within the misorientation/structure space by simply varying the misorientation angle around a fixed misorientation axis, the proposed method creates sets of grain boundary models by means of isocurves associated with important fracture-related properties of the adjoining lattices. Such properties may include anisotropic elastic moduli, the Schmid factor for primary slip, and the propensity for simultaneous slip on multiple slip systems. This approach eliminates the effect of lattice properties from the comparative analysis of interfacial fracture properties and thus enables the identification of structure-property relationships for grain boundaries. As an example, this methodology is implemented to study crack propagation along Ni grain boundaries. Segregated H is used as a means to emphasize differences in the selected grain boundary structures while keeping lattice properties fixed.