Abstract
We determined the ground state of clusters adsorbed on one side of graphene for selected cluster sizes in the range from to . For all investigated clusters variational and diffusion Monte Carlo simulations were performed at K, and in addition for a selected subset finite temperature path integral Monte Carlo. At K the liquid or solid character of each cluster was investigated by restricting the phase using corresponding importance sampling trial-wave functions. The -graphene interaction was modeled as a sum of individual -C interactions, where both isotropic and anisotropic models were tested; also the effect of the substrate-mediated McLachlan interaction was investigated. We have found homogeneous crystallization in models of anisotropic interactions, starting from clusters with atoms in simulations without the McLachlan interaction, and between and 61 when it is included. The atoms become increasingly delocalized as one moves from the center of the cluster to the perimeter, evidenced by the Lindemann parameter. On the other hand, in the case of the isotropic interaction model, a liquidlike structure is more favorable for all considered cluster sizes. We use a liquid-drop model to extrapolate the energy per particle to the limit, and the results are compared with the values obtained in studies of bulk on graphene. Low-temperature path integral Monte Carlo simulations are in agreement with ground-state results.
2 More- Received 22 January 2016
- Revised 14 June 2016
DOI:https://doi.org/10.1103/PhysRevB.94.045428
©2016 American Physical Society