The Reactive empirical bond order (REBO) potential developed by Brenner et al. 1,2 for molecular dynamics (MD) simulations of hydrocarbons, and recently extended to include interactions with oxygen atoms by Ni et al. 3 is modified for graphene-oxide (GO). Based on DFT calculations, we optimized the REBO-CHO potential to improve its ability to calculate the binding energy of an oxygen atom to graphene and the equilibrium C-O bond distances. In this work, the approach towards the optimization is based on modifying the bond order term. The modified REBO-CHO potential is applied to investigate the properties of some GO samples.