On the consistency of the exfoliation free energy of graphenes by molecular simulations

Monolayer graphene is now produced at significant yields, by liquid phase exfoliation of graphites in solvents. This has increased the interest in molecular simulation studies to give new insights in the field. We use decoupling simulations to compute the exfoliation free energy of graphenes in a liquid environment. Starting from a bilayer graphene configuration, we decouple the Van der Waals interactions of a graphene monolayer in the presence of saline water. Then, we introduce the monolayer back into water by coupling its interactions with water molecules and ions. A different approach to compute the graphene exfoliation free energy is to use umbrella sampling. We apply umbrella sampling after pulling the graphene monolayer on the shear direction up to a distance from a bilayer. We show that the decoupling and umbrella methods give highly consistent free energy results for three bilayer graphene samples with different size. This strongly suggests that the systems in both methods remain closely in equilibrium as we move between the states before and after the exfoliation. Therefore, the amount of nonequilibrium work needed to peel the two layers apart is minimized efficiently.