▎ 摘　　要

It has been shown that the broken bonds of an unreconstructed graphene edge generate compressive edge stresses leading to edge warping. Here, we investigate edge energies and edge stresses of graphene nanoribbons with arbitrary orientations from armchair to zigzag, considering both flat and warped edge shapes in the presence and absence of hydrogen. We use the second generation reactive empirical bond order potential to calculate the edge energies and stresses for clean and hydrogenated edges. Using these energies, we perform a Wulff construction to determine the equilibrium shapes of flat graphene flakes as a function of hydrogen chemical potential. While edge stresses for clean, flat edges are compressive, they become tensile if allowed to warp. Conversely, we find that edge energies change little (similar to 1%) with edge warping. Hydrogenation of the edges virtually eliminates both the edge energy and edge stresses. For warped edges an approximately linear relationship is found between amplitudes and wavelengths. The equilibrium shape of a graphene flake is determined by the value of the hydrogen chemical potential. For very small (and large) values of it the flakes have a nearly hexagonal (dodecagon) shape with zigzag oriented edges, while for intermediate values graphene flakes are found with complex shapes.