We determine, by means of density functional theory, the stability and the structure of graphene nanoribbon (GNR) edges in presence of molecules such as oxygen, water, ammonia, and carbon dioxide. As in the case of hydrogen-terminated nanoribbons, we
find that the most stable armchair and zigzag configurations are characterized by a non-metallic/non-magnetic nature, and are compatible with Clars sextet rules, well known in organic chemistry. In particular, we predict that, at thermodynamic equilibrium, neutral GNRs in oxygen-rich atmosphere should preferentially be along the armchair direction, while water-saturated GNRs should present zigzag edges. Our results promise to be particularly useful to GNRs synthesis, since the most recent and advanced experimental routes are most effective in water and/or ammonia-containing solutions.
